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 ($\lambda_{R}$) profiles. We observe a spread in the $\lambda_{R}$ profiles of our sample dEs, finding that the $\lambda_{R}$ profiles of half of them are as high as those of low-mass field galaxies. The remaining dEs exhibit $\lambda_{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-$\lambda_{R}$ dEs in our sample were processed by their previous host halo, prior to their infall to Virgo, and that the high-$\lambda_{R}$ dEs may be experiencing ram pressure stripping in Virgo.

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 $\lambda_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 $\lambda_R$ decreases towards lower stellar mass, with a corresponding increase in the proportion of slowly rotating galaxies in this regime. The decrease of $\lambda_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.

We present the study of the south-west group in the Fornax cluster centred on the brightest group galaxy (BGG) Fornax A, observed as part of the Fornax Deep Survey (FDS). This includes the analysis of the bright group members (mB < 16 mag) and the intra-group light (IGL). The main objective of this work is to investigate the assembly history of the Fornax A group and to compare its physical quantities as a function of the environment to that of the Fornax cluster core. For all galaxies, we extract the azimuthally averaged surface brightness profiles in three optical bands (g, r, i) by modelling the galaxy's isophotes. We derive their colour profiles and structural parameters in all respective bands. The long integration time and large covered area of the FDS allow us to also estimate the amount of IGL. The majority of galaxies in the Fornax A group are late-type galaxies (LTGs), spanning a range of stellar mass of $8 < log (M_* M_{\odot}) < 10.5$. Six out of nine LTGs show a Type III (up-bending) break in their light profiles, which is either suggestive of strangulation halting star-formation in their outskirts or their HI-richness causing enhanced star-formation in their outer-discs. The estimated luminosity of the IGL is $6 \pm 2 \times 10^{10} L_{\odot}$ in g-band, which corresponds to about 16% of the total light in the group. The Fornax A group appears to be in an early stage of assembly with respect to the cluster core. The environment of the Fornax A group is not as dense as that of the cluster core, with all galaxies except the BGG showing similar morphology, comparable colours and stellar masses, and Type III disc-breaks, without any clear trend of these properties with group-centric distances. The main contribution to the IGL is from the minor merging in the outskirts of the BGG NGC1316 and, probably, the disrupted dwarf galaxies close to the group centre.

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

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}/\sigma \sim$ 0.4, and blue GCs have weak rotation with $V_{\rm max}/\sigma \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.

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}+\sigma_{\star_e}^{2}$, which combines rotation velocity $V_{R_e}$, and velocity dispersion $\sigma_{\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}-\sigma_{\star_e}$ ones. We also explore a more general Universal Fundamental Plane described by the equation $log(\Upsilon_{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, $\Upsilon_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, $\Upsilon_{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.

Globular cluster (GC) systems of massive galaxies often show a bimodal colour distribution. This has been interpreted as a metallicity bimodality, created by a two-stage galaxy formation where the red, metal-rich GCs were formed in the parent halo and the blue metal-poor GCs were accreted. This interpretation, however, crucially depends on the assumption that GCs are exclusively old stellar systems with a linear colour-metallicity relation (CZR). The shape of the CZR and range of GC ages are currently under debate, because their study requires high quality spectra to derive reliable stellar population properties. We determined metallicities with full spectral fitting from a sample of 187 GCs with high spectral signal-to-noise ratio in 23 galaxies of the Fornax cluster that were observed as part of the Fornax 3D project. The derived CZR from this sample is non-linear and can be described by a piecewise linear function with a break point at ($g - z$) $\sim$ 1.1 mag. The less massive galaxies in our sample ($M_\ast < 10^{10} M_\odot$) appear to have slightly younger GCs, but the shape of the CZR is insensitive to the GC ages. Although the least massive galaxies lack red, metal-rich GCs, a non-linear CZR is found irrespective of the galaxy mass, even in the most massive galaxies ($M_\ast \geq 10^{11} M_\odot$). Our CZR predicts narrow unimodal GC metallicity distributions for low mass and broad unimodal distributions for very massive galaxies, dominated by a metal-poor and metal-rich peak, respectively, and bimodal distributions for galaxies with intermediate masses (10$^{10}$ $\leq$ $M_\ast < 10^{11} M_\odot$) as a consequence of the relative fraction of red and blue GCs. The diverse metallicity distributions challenge the simple differentiation of GC populations solely based on their colour.

Globular clusters (GCs) are found ubiquitously in massive galaxies and due to their old ages, they are regarded as fossil records of galaxy evolution. Spectroscopic studies of GC systems are often limited to the outskirts of galaxies, where GCs stand out against the galaxy background and serve as bright tracers of galaxy assembly. In this work, we use the capabilities of the Multi Unit Explorer Spectrograph (MUSE) to extract a spectroscopic sample of 722 GCs in the inner regions ($\lesssim 3 R_\text{eff}$) of 32 galaxies in the Fornax cluster. These galaxies were observed as part of the Fornax 3D project, a MUSE survey that targets early and late-type galaxies within the virial radius of Fornax. After accounting for the galaxy background in the GC spectra, we extracted line-of-sight velocities and determined metallicities of a sub-sample of 238 GCs. We found signatures of rotation within GC systems, and comparing the GC kinematics and that of the stellar body shows that the GCs trace the spheroid of the galaxies. While the red GCs prove to closely follow the metallicity profile of the host galaxy, the blue GCs show a large spread of metallicities but they are generally more metal-poor than the host.

We investigate the shape of the Rotation Curves (RCs) of $z\sim 1$ star-forming galaxies and compare them with the local star-forming galaxies. For this purpose, we have used $409$ galaxies from the K-band Multi-Object Spectrograph (KMOS) for Redshift One Spectroscopic Survey (KROSS). This sample covers the redshift range $0.57\leq z \leq 1.04$, effective radii $0.69 \leq R_e \left[kpc \right] \leq 7.73$, absolute H-band magnitude $-24.46 \leq M_H \leq -18.85$ with median stellar mass $log\left(M_* \ \left[M_\odot \right]\right)=9.95$ and median total star-formation rate $log\left(SFR_{tot} \ \left[M_\odot \ yr^{-1} \right]\right)=1.49$. Using $^{3D}$BAROLO (Barolo), we extract $H\alpha$ kinematic maps and corresponding Rotation Curves (RCs). The main advantage of Barolo is that it incorporates the beam smearing in the 3D observational space, which provide us with the intrinsic rotation velocity even in the low spatial resolution data. Using Asymmetric Drift Correction (ADC), we have corrected the RCs for the pressure gradient effect, which seems to be a more dominant effect than beam smearing in high-z galaxies. Nearly all objects ($0.1 <v/ \sigma< 15$) are affected by the pressure gradient, and we noticed that ADC improves the rotation velocity of these systems by $\sim 10-87\%$. Only a combination of the three techniques (3D-kinematic modelling + 3D-Beamsmearing correction + ADC ) yields the intrinsic RC of an individual galaxy. Further, we present the co-added RCs constructed out of 237 high-quality objects to obtain intrinsic RC shapes out to 6.4 $\times$ disk scale length. We do not see any change in the shape of RCs with respect to the local star-forming disk-type galaxies. In contrast, we do find a significant evolution in the stellar-disk length ($R_D$) of the galaxies. Therefore, we conclude stellar disk evolves over cosmic time while total mass distribution stays constant.

This work is based on deep multi-band (g, r, i) data from the Fornax Deep Survey with VST. We analyse the surface brightness profiles of the 19 bright ETGs inside the virial radius of the Fornax cluster. The main aim of this work is to identify signatures of accretion onto galaxies by studying the presence of outer stellar halos, and understand their nature and occurrence. Our analysis also provides a new and accurate estimate of the intra-cluster light inside the virial radius of Fornax. We performed multi-component fits to the azimuthally averaged surface brightness profiles available for all sample galaxies. This allows to quantify the relative weight of all components in the galaxy structure that contribute to the total light. In addition, we derived the average g-i colours in each component identified by the fit, as well as the azimuthally averaged g-i colour profiles, to correlate them with the stellar mass of each galaxy and the location inside the cluster. We find that in the most massive and reddest ETGs the fraction of light in, probably accreted, halos is much larger than in the other galaxies. Less-massive galaxies have an accreted mass fraction lower than 30%, bluer colours and reside in the low-density regions of the cluster. Inside the virial radius of the cluster, the total luminosity of the intra-cluster light, compared with the total luminosity of all cluster members, is about 34%. Inside the Fornax cluster there is a clear correlation between the amount of accreted material in the stellar halos of galaxies and the density of the environment in which those galaxies reside. By comparing this quantity with theoretical predictions and previous observational estimates, there is a clear indication that the driving factor for the accretion process is the total stellar mass of the galaxy, in agreement with the hierarchical accretion scenario.

Extragalactic planetary nebulae (PNe) are detectable through relatively strong nebulous [OIII] emission and act as direct probes into the local stellar population. Because they have an apparently universal invariant magnitude cut-off, PNe are also considered to be a remarkable standard candle for distance estimation. Through detecting PNe within the galaxies, we aim to connect the relative abundances of PNe to the properties of their host galaxy stellar population. By removing the stellar background components from FCC 167 and FCC 219, we aim to produce PN luminosity functions (PNLF) of these galaxies, and thereby also estimate the distance modulus to these two systems. Finally, we test the reliability and robustness of our novel detection and analysis method. It detects unresolved point sources by their [OIII] 5007{\AA} emission within regions that have previously been unexplored. We model the [OIII] emissions in the spatial and spectral dimensions together, as afforded to us by the Multi Unit Spectroscopic Explorer (MUSE), and we draw on data gathered as part of the Fornax3D survey. For each source, we inspect the properties of the nebular emission lines to remove other sources that might hinder the safe construction of the PNLF, such as supernova remnants and HII regions. As a further step, we characterise any potential limitations and draw conclusions about the reliability of our modelling approach through a set of simulations. By applying this novel detection and modelling approach to integral field unit observations, we report for the distance estimates and luminosity-specific PNe frequency values for the two galaxies. Furthermore, we include an overview of source contamination, galaxy differences, and possible effects on the PNe populations in the dense stellar environments.

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).

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$\alpha$ 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.

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 $\lambda_z$ vs. radius $r$, the intrinsic stellar population distribution in age $t$ vs. metallicity $Z$, and the correlation between orbits' circularity $\lambda_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. $\sigma_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.

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).

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.

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.

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.

We present an orbit-based method of combining stellar and cold gas kinematics to constrain the dark matter profile of early-type galaxies. We apply this method to early-type galaxy NGC 2974, using Pan-STARRS imaging and SAURON stellar kinematics to model the stellar orbits, and introducing HI kinematics from VLA observation as a tracer of the gravitational potential. The introduction of the cold gas kinematics shows a significant effect on the confidence limits of especially the dark halo properties: we exclude more than $95\%$ of models within the $1-\sigma$ confidence level of Schwarzschild modelling with only stellar kinematics, and reduce the relative uncertainty of the dark matter fraction significantly to $10\%$ within $5 R_\mathrm{e}$. Adopting a generalized-NFW dark matter profile, we measure a shallow cuspy inner slope of $0.6^{+0.2}_{-0.3}$ when including the cold gas kinematics in our model. We cannot constrain the inner slope with the stellar kinematics alone.

We perform spatially resolved stellar population analysis for a sample of 69 early-type galaxies (ETGs) from the CALIFA integral field spectroscopic survey, including 48 ellipticals and 21 S0's. We generate and quantitatively characterize profiles of light-weighted mean stellar age and metallicity within $\lesssim 2R_e$, as a function of radius and stellar-mass surface density $\mu_*$. We study in detail the dependence of profiles on galaxies' global properties, including velocity dispersion $\sigma_e$, stellar mass, morphology. ETGs are universally characterized by strong, negative metallicity gradients ($\sim -0.3\,\text{dex}$ per $R_e$) within $1\,R_e$, which flatten out moving towards larger radii. A quasi-universal local $\mu_*$-metallicity relation emerges, which displays a residual systematic dependence on $\sigma_e$, whereby higher $\sigma_e$ implies higher metallicity at fixed $\mu_*$. Age profiles are typically U-shaped, with minimum around $0.4\,R_e$, asymptotic increase to maximum ages beyond $\sim 1.5\,R_e$, and an increase towards the centre. The depth of the minimum and the central increase anti-correlate with $\sigma_e$. A possible qualitative interpretation of these observations is a two-phase scenario. In the first phase, dissipative collapse occurs in the inner $1\,R_e$, establishing a negative metallicity gradient. The competition between the outside-in quenching due to feedback-driven winds and some form of inside-out quenching, possibly caused by central AGN feedback or dynamical heating, determines the U-shaped age profiles. In the second phase, the accretion of ex-situ stars from quenched and low-metallicity satellites shapes the flatter stellar population profiles in the outer regions.

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.

[Abridged] We present the apparent stellar angular momentum of 300 galaxies across the Hubble sequence, using integral-field spectroscopic data from the CALIFA survey. Adopting the same $\lambda_\mathrm{R}$ parameter previously used to distinguish between slow and fast rotating early-type (elliptical and lenticular) galaxies, we show that spiral galaxies as expected are almost all fast rotators. Given the extent of our data, we provide relations for $\lambda_\mathrm{R}$ measured in different apertures, including conversions to long-slit 1D apertures. Our sample displays a wide range of $\lambda_\mathrm{Re}$ values, consistent with previous IFS studies. The fastest rotators are dominated by relatively massive and highly star-forming Sb galaxies, which preferentially reside in the main star-forming sequence. These galaxies reach $\lambda_\mathrm{Re}$ values of $\sim$0.85, are the largest galaxies at a given mass, and display some of the strongest stellar population gradients. Compared to the population of S0 galaxies, our findings suggest that fading may not be the dominant mechanism transforming spirals into lenticulars. Interestingly, we find that $\lambda_\mathrm{Re}$ decreases for late-type Sc and Sd spiral galaxies, with values than in occasions puts them in the slow-rotator regime. While for some of them this can be explained by their irregular morphologies and/or face-on configurations, others are edge-on systems with no signs of significant dust obscuration. The latter are typically at the low-mass end, but this does not explain their location in the classical ($V/\sigma$,$\varepsilon$) and ($\lambda_\mathrm{Re}$,$\varepsilon$) diagrams. Our initial investigations, based on dynamical models, suggest that these are dynamically hot disks, probably influenced by the observed important fraction of dark matter within R$_\mathrm{e}$.

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.

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.

Using new long-slit spectroscopy obtained with X-Shooter at ESO-VLT, we study, for the first time, radial gradients of optical and Near-Infrared IMF-sensitive features in a representative sample of galaxies at the very high-mass end of the galaxy population. The sample consists of seven early-type galaxies (ETGs) at $z\sim0.05$, with central velocity dispersion in the range $300<\sigma<350$km/s. Using state-of-art stellar population synthesis models, we fit a number of spectral indices, from different chemical species (including TiO's and Na indices), to constrain the IMF slope (i.e. the fraction of low-mass stars), as a function of galactocentric distance, over a radial range out to $\sim4$kpc. ETGs in our sample show a significant correlation of IMF slope and surface mass density. The bottom-heavy population (i.e. an excess of low-mass stars in the IMF) is confined to central galaxy regions with surface mass density above $\sim 10^{10} M_\odot kpc^{-2}$, or, alternatively, within a characteristic radius of $\sim2$~kpc. Radial distance, in physical units, and surface mass density, are the best correlators to IMF variations, with respect to other dynamical (e.g. velocity dispersion) and stellar population (e.g. metallicity) properties. Our results for the most massive galaxies suggest that there is no single parameter} that fully explains variations in the stellar IMF, but IMF radial profiles at z$\sim$0 rather result from the complex formation and mass accretion history of galaxy inner and outer regions.

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.

The discovery around the turn of the millenium of a population of very massive (M$_\star$ > 2$\times$10$^6$ M$_\odot$) compact stellar systems (CSS) with physical properties (radius, velocity dispersion, stellar mass etc.) that are intermediate between those of the classical globular cluster (GC) population and galaxies led to questions about their exact nature. Recently a consensus has emerged that these objects, usually called ultra compact dwarfs (UCDs), are a mass-dependent mixture of high mass star clusters and remnant nuclei of tidally disrupted galaxies. The existence of genuine star clusters with stellar masses >10$^7$ M$_\odot$ naturally leads to questions about the upper mass limit of the star cluster formation process. In this work we compile a comprehensive catalog of compact stellar systems, and reinforce the evidence that the true ancient star cluster population has a maximum mass of M$_\star$ ~ 5$\times$10$^7$ M$_\odot$, corresponding to a stellar mass at birth of close to 10$^8$ M$_\odot$. We then discuss several physical and statistical mechanisms potentially responsible for creating this limiting mass.

Stellar feedback plays a significant role in modulating star formation, redistributing metals, and shaping the baryonic and dark structure of galaxies -- however, the efficiency of its energy deposition to the interstellar medium is challenging to constrain observationally. Here we leverage HST and ALMA imaging of a molecular gas and dust shell ($M_{H2} \sim 2\times 10^{5} ~{\rm M}_{\odot}$) in an outflow from the nuclear star forming ring of the galaxy NGC 3351, to serve as a boundary condition for a dynamical and energetic analysis of the outflowing ionised gas seen in our MUSE TIMER survey. We use \texttt{STARBURST99} models and prescriptions for feedback from simulations to demonstrate that the observed star formation energetics can reproduce the ionised and molecular gas dynamics -- provided a dominant component of the momentum injection comes from direct photon pressure from young stars, on top of supernovae, photoionisation heating and stellar winds. The mechanical energy budget from these sources is comparable to low luminosity AGN, suggesting that stellar feedback can be a relevant driver of bulk gas motions in galaxy centres - although here $\lesssim 10^{-3}$ of the ionized gas mass is escaping the galaxy. We test several scenarios for the survival/formation of the cold gas in the outflow, including in-situ condensation and cooling. Interestingly, the geometry of the molecular gas shell, observed magnetic field strengths and emission line diagnostics are consistent with a scenario where magnetic field lines aided survival of the dusty ISM as it was initially launched (with mass loading factor $\lesssim 1$) from the ring by stellar feedback. This system's unique feedback driven morphology can hopefully serve as a useful litmus test for feedback prescriptions in magnetohydrodynamical galaxy simulations.

We construct Schwarzschild orbit-based models of NGC 7457, known as a peculiar low-mass lenticular galaxy. Our best-fitting model successfully retrieves most of the unusual kinematics behaviours of this galaxy, in which, the orbital distribution of stars is dominated by warm and hot orbits. The reconstructed surface brightness of the hot component matches fairly well the photometric bulge and the reconstructed LOSVD map of this component shows clear rotation around the major photometric axis of the galaxy. In the absence of a dominant cold component, the outer part of our model is dominated by warm orbits, representing an exponential thick disc. Our orbital analysis also confirm the existence of a counter-rotating orbital substructure in the very centre, reported in previous observational studies. By comparing our model with a variety of simulation studies, and considering the stellar kinematics and populations properties of this galaxy, we suggest that the thick disc is most likely a dynamically heated structure, formed through the interactions and accretion of satellite(s) with near-polar initial inclination. We also suggest a merger-driven process as the most plausible scenario to explain the observed and dynamically-modelled properties of the bulge of NGC 7457. We conclude that both the high level of cylindrical rotation and unusually low velocity dispersion reported for the NGC 7457 have most-likely external origins. Therefore, NGC 7457 could be considered as a candidate for merger-driven cylindrical rotation in the absence of a strong bar in disc galaxies.

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. The large NSC (r$_\rm{eff} = 66$ pc) and rich GC system of the early-type Fornax cluster galaxy FCC47 (NGC1336) render this galaxy an ideal laboratory to constrain NSC formation. Using MUSE science verification data assisted by adaptive optics, we obtained maps for the stellar kinematics and for stellar-population properties of FCC47. We extracted the spectra of the central NSC and determined line-of-sight velocities of 24 GCs and metallicities of five. FCC47 shows two decoupled components (KDCs): a rotating disk and the NSC. Our orbit-based dynamical Schwarzschild model revealed that the NSC is a distinct kinematic feature and it constitutes the peak of metallicity and old ages in the galaxy. The main body consists of two counter-rotating populations and is dominated by a more metal-poor population. The GC system is bimodal with a dominant metal-poor population and the total GC system mass is $\sim 17\%$ of the NSC mass ($\sim$ 7 $\times$ $10^8 M_\odot$). The rotation, high metallicity and high mass of the NSC cannot be uniquely explained by GC-inspiral and most likely requires additional, but quickly quenched, in-situ formation. The presence of two KDCs most probably are evidence of a major merger that has altered the structure of FCC47 significantly, indicating the important role of galaxy mergers in forming the complex kinematics in the galaxy-NSC system.

We present the study of a magnitude limited sample (mB < 16.6 mag) of 13 late type galaxies (LTGs), observed inside the virial radius, Rvir 0.7 Mpc, of the Fornax cluster within the Fornax Deep Survey (FDS). The main objective is to use surface brightness (SB) profiles and g-i colour maps to obtain information on the internal structure of these galaxies and find signatures of the mechanisms that drive their evolution in high-density environment, which is inside the virial radius of the cluster. By modelling galaxy isophotes, we extract the azimuthally averaged surface brightness profiles in four optical bands. We also derive g-i colour profiles, and structural parameters like total magnitude and effective radius. For 10 of the galaxies in this sample, we observe a clear discontinuity in their SB profiles, derive their break radius (Br), and classify their disc-breaks into Type-II (down-bending) or Type-III (up-bending). We find that Type-II galaxies have bluer average (g-i) colour in their outer discs while Type-III galaxies are redder. Br increases with stellar mass and molecular gas mass while decreases with molecular gas-fractions. The inner and outer scale-lengths increase monotonically with absolute magnitude, as found in other works. In Fornax, galaxies with morphological type 5< T< 9 (~60 % of the sample) are located beyond the high-density, ETG-dominated regions, but no correlation found between T and the disc-break type. The main results of this work suggest that the disc-breaks of LTGs inside the virial radius of the Fornax cluster seem to have arisen through a variety of mechanisms, which is evident in their outer-disc colours and the absence of molecular gas beyond their break radius in some cases. This can result in a variety of stellar populations inside and outside the break radii.

The 31 brightest galaxies (m_B < 15 mag) inside the virial radius of the Fornax cluster were observed from the centres to the outskirts with the Multi Unit Spectroscopic Explorer on the Very Large Telescope. These observations provide detailed high-resolution maps of the line-of-sight kinematics and line strengths of the stars and ionised gas reaching 2-3 Re for 21 early-type galaxies and 1-2 Re for 10 late-type galaxies. The majority of the galaxies are regular rotators, with eight hosting a kinematically distinct core. Only two galaxies are slow rotators. The mean age, total metallicity, and [Mg/Fe] abundance ratio in the bright central region inside 0.5 Re and in the galaxy outskirts are presented. Extended emission-line gas is detected in 13 galaxies, most of them are late-type objects with wide-spread star formation. The measured structural properties are analysed in relation to the galaxies' position in the projected phase space of the cluster. This shows that the Fornax cluster appears to consist of three main groups of galaxies inside the virial radius: the old core; a clump of galaxies, which is aligned with the local large-scale structure and was accreted soon after the formation of the core; and a group of galaxies that fell in more recently.

The formation and assembly process of massive galaxies is a combination of two phases: an initial in-situ-dominated one followed by an ex-situ-dominated evolution. Separating these two contributions is therefore crucial to understand the baryonic cycle within massive halos. A recently discovered population of so-called naked red nuggets, galaxies that shortcut the ex-situ stage preserving their pristine properties, presents a unique opportunity to study in detail star formation in massive galaxies without the confounding effect of later accretion. We investigate the spatially resolved star formation histories of a sample of 12 naked red nuggets. We measure how their radial light distributions, star formation rates and central densities evolved in time. We find that, while forming stars, red nuggets become gradually more concentrated, reaching a maximum concentration at quenching. After being quenched, they kept forming stars in a more disky-like configuration. Our measurements suggest that super-massive black holes and host galaxies grow their mass in a self-regulated way until a characteristic Mbh/Mhalo is reached. Once black holes are massive enough, red nuggets get quenched and depart from the star formation main sequence. While in the main sequence, red nuggets evolve at roughly constant star formation rate. This can explain up to ~0.3 dex of the scatter of the star formation main sequence, as well as its higher normalization observed in the early Universe. Hence, our results suggest that the main sequence is composed of populations of galaxies at different evolutionary stages, and that the scatter is therefore due to secular processes.

Motivated by the recently discovered kinematic "Hubble sequence" shown by the stellar orbit-circularity distribution of 260 CALIFA galaxies, we make use of a comparable galaxy sample at z = 0 with a stellar mass range from 5E9 to 5E11 solar masses, selected from the IllustrisTNG simulation and study their stellar orbit compositions in relation to a number of other fundamental galaxy properties.We find that the TNG100 simulation broadly reproduces the observed fractions of different orbital components and their stellar mass dependencies. In particular, the mean mass dependencies of the luminosity fractions for the kinematically warm and hot orbits are well reproduced within model uncertainties of the observed galaxies. The simulation also largely reproduces the observed peak and trough features at a stellar mass of 1-2E10 solar masses, in the mean distributions of the cold- and hot-orbit fractions, respectively, indicating fewer cooler orbits and more hotter orbits in both more- and less-massive galaxies beyond such a mass range. Several marginal disagreements are seen between the simulation and observations: the average cold-orbit (counter-rotating) fractions of the simulated galaxies below (above) a stellar mass of 6E10 solar masses, are systematically higher than the observational data by < 10% (absolute orbital fraction); the simulation also seems to produce more scatter for the cold-orbit fraction and less so for the non-cold orbits at any given galaxy mass. Possible causes that stem from the adopted heating mechanisms are discussed.

We present a combination of the Schwarzschild orbit-superposition dynamical modelling technique with the spatially-resolved mean stellar age and metallicity maps to uncover the formation history of galaxies. We apply this new approach to a remarkable 5-pointing mosaic of VLT/MUSE observations obtained by Gu\'erou et al. (2016) extending to a maximum galactocentric distance of $120''$ ($\sim 5.6\ {\rm kpc}$) along the major axis, corresponding to $\sim 2.5\ {\rm R}_e$. Our method first identifies `families' of orbits from the dynamical model that represent dynamically-distinct structures of the galaxy. Individual ages and metallicities of these components are then fit for using the stellar-population information. Our results highlight components of the galaxy that are distinct in the combined stellar dynamics/populations space, which implies distinct formation paths. We find evidence for a dynamically-cold, metal-rich disk, consistent with a gradual in-situ formation. This disk is embedded in a generally-old population of stars, with kinematics ranging from dispersion-dominated in the centre to an old, diffuse, metal-poor stellar halo at the extremities. We find also a direct correlation between the dominant dynamical support of these components, and their associated age, akin to the relation observed in the Milky Way. This approach not only provides a powerful model for inferring the formation history of external galaxies, but also paves the way to a complete population-dynamical model.

We evaluate the capabilities of Schwarzschild's orbit-superposition method by applying it to galaxies from the large scale, high resolution Illustris simulation. Nine early-type galaxies with a range of triaxiality are selected, and we create mock integral field unit data for five line-of-sight projections of each galaxy. Each of the 45 mock data sets is taken as an independent observed galaxy. Using van den Bosch's 2008 triaxial Schwarzschild implementation, we assess model estimates of various galaxy properties, covering mass profiles, intrinsic shapes, stellar orbit distributions and velocity anisotropies. Total mass within $\overline{R_{\rm e}}$ is recovered well with average deviations within $\pm15$ percent. Stellar mass is underestimated by $\sim24$ percent and dark matter overestimated by $\sim38$ percent (assuming an NFW dark matter profile and allowing for degeneracy between stellar mass and dark matter mass). Using a gNFW profile, these values improve to $\sim13$ percent for stellar mass and $\sim18$ percent for dark matter. Axis ratio estimates show a moderate bias of $\Delta (b/a)=0.07$ and $\Delta (c/a)=0.14$ ($a\ge b\ge c$). Distributions of the orbit circularities $\lambda_z$ and $\lambda_x$, representing rotation about the minor and major axes, are well reconstructed. Separating orbits into thermal categories, our models match the average fractions of these categories to within $10$ percent. Velocity anisotropy is well estimated with values matching in the inner regions but becoming slightly radially biased in the outer regions. Overall, the galaxy property estimates we obtained using Schwarzschild modelling are not implausible and are representative of the simulated galaxies we modelled.

The Fornax Deep Survey Dwarf galaxy Catalog (FDSDC) includes 564 dwarf galaxies in the Fornax cluster and the in-falling Fornax A subgroup. We use the FDSDC galaxies for statistical analysis of the structural and stellar population differences in the range of galactic environments within the Fornax cluster. We present the standard scaling relations for the dwarfs and analyze trends as a function of cluster-centric radius. We find a different behavior for the bright dwarfs (-18.5 mag < M$_r$ < -16 mag) as compared to the fainter ones (M$_r$ > -16 mag): While considering galaxies in the same magnitude-bins, we find that, while for fainter dwarfs the g'-r' color is redder for lower surface brightness objects (as expected from fading stellar populations), for brighter dwarfs the color is redder for the higher surface brightness and higher S\'ersic n objects. The trend of the bright dwarfs might be explained by those galaxies being affected by harassment and by slower quenching of star formation in their inner parts. As the fraction of early-type dwarfs with respect to late-types increases toward the central parts of the cluster, the color-surface brightness trends are also manifested in the cluster-centric trends, confirming that it is indeed the environment that changes the galaxies. We also estimate the strengths of the ram-pressure stripping, tidal disruption, and harassment in the Fornax cluster, and find that our observations are consistent with the theoretically expected ranges of galaxy properties where each of those mechanisms dominate. We furthermore find that the luminosity function, color-magnitude relation, and axis-ratio distribution of the dwarfs in the center of the Fornax cluster are similar to those in the center of the Virgo cluster.

We used deep MUSE observations to perform a stellar-kinematic and population analysis of FCC 153 and FCC 177, two edge-on S0 galaxies in the Fornax cluster. The geometrical definition of the different structural components of these two galaxies allows us to describe the nature of their thick disks. These are both old, relatively metal poor and [Mg/Fe]-enhanced, and their star formation history (SFH) reveals a minor younger component whose chemical properties suggest its later accretion. Moreover, the outer regions of these geometrically defined thick disks show higher values of metallicity and lower values of [Mg/Fe]. These stars probably formed in the thin-disk region and they were dynamically heated to form the flares present in these two galaxies. We propose different formation scenarios for the three populations of these thick disks: in-situ formation, accretion and disk heating. A clear distinction in age is found between the metal poor and [Mg/Fe]-enhanced thick disks (old, $\sim 12-13$ Gyr), and the metal rich and less [Mg/Fe]-enhanced thin disks (young, $\sim 4-5$ Gyr). These two galaxies show signs of relatively recent star formation in their thin disks and nuclear regions. While the thin disks show more continuous SFHs, the nuclei display a rather bursty SFH. These two galaxies are located outside of the densest region of the Fornax cluster where FCC 170 resides. This other edge-on S0 galaxy was studied by \citet{Pinna2019}. We compare and discuss our results with this previous study. The differences between these three galaxies, at different distances from the cluster center, suggest that the environment can have a strong effect on the galaxy evolutionary path.

Photometric surveys of galaxy clusters have revealed a large number of ultra compact dwarfs (UCDs) around predominantly massive elliptical galaxies. Their origin is still debated as some UCDs are considered to be the remnant nuclei of stripped dwarf galaxies while others seem to mark the high-mass end of the star cluster population. We aim to characterise the properties of a UCD found at very close projected distance (1.1 kpc) from the centre of the low-mass (M~10^10 M_sun) early-type galaxy FCC47. This is a serendipitous discovery from MUSE adaptive optics science verification data. We explore the potential origin of this UCD as either a massive cluster or the remnant nucleus of a dissolved galaxy. We use archival Hubble Space Telescope data to study the photometric and structural properties of FCC47-UCD1. In the MUSE data, the UCD is unresolved, but we use its spectrum to determine the radial velocity and metallicity. FCC47-UCD1's surface brightness is best described by a single King profile with low concentration C = R_t/R_c ~ 10 and large effective radius (r_e = 24pc). Its integrated magnitude and a blue colour (G = -10.6 mag, g-z = 1.46 mag) combined with with a metallicity of [M/H] = -1.12+-0.10 dex and an age > 8 Gyr obtained from the full fitting of the MUSE spectrum suggests a stellar population mass of M_star = 4.87x10^6 M_sun. The low S/N of the MUSE spectrum prevents detailed stellar population analysis. Due to the limited spectral resolution of MUSE, we can only give an upper limit on the velocity dispersion (sig < 17km/s), and consequently on its dynamical mass (M_dyn < 1.3x10^7 M_sun). The origin of the UCD cannot be constrained with certainty. The low metallicity, old age and magnitude are consistent with a star cluster origin, whereas the extended size is consistent with an origin as the stripped nucleus of a dwarf galaxy with a initial stellar mass of a few 10^8 M_sun.

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.

The stellar initial mass function (IMF) regulates the baryonic cycle within galaxies, and is a key ingredient to translate observations into physical quantities. Although for decades it was assumed to be universal, there is now growing observational evidence showing that the center of massive early-type galaxies host an enhanced population of low-mass stars compared to the expectations from the Milky Way. Moreover, these variations in the IMF have been found to be related to the radial metallicity variations in massive galaxies. We present here a two-dimensional stellar population analysis of the massive lenticular galaxy FCC 167 (NGC 1380) as part of the Fornax3D project. Using a newly developed stellar population fitting scheme, we derive a full two-dimensional IMF map of an early-type galaxy. This two-dimensional analysis allows us go further than a radial analysis, showing how the metallicity changes along a disc-like structure while the IMF follows a distinct, less disky distribution. Thus, our findings indicate that metallicity cannot be the sole driver of the observed radial IMF variations. In addition, a comparison with the orbital decomposition shows suggestive evidence of a coupling between stellar population properties and the internal dynamical structure of FCC 167, where metallicity and IMF maps seem to track the distribution of cold and warm orbits, respectively.

Maintenance-mode feedback from low-accretion rate AGN, manifesting itself observationally through radio-loudness, is invoked in all cosmological galaxy formation models as a mechanism that prevents excessive star-formation in massive galaxies (M$_*$ $\gtrsim$ 3$\times$10$^{10}$ M$_{\odot}$). We demonstrate that at a fixed mass the incidence of radio-loud AGN (L $>$ 10$^{23}$ WHz$^{- 1}$) identified in the FIRST and NVSS radio surveys among a large sample of quiescent (non-star forming) galaxies selected from the SDSS is much higher in geometrically round galaxies than in geometrically flat, disk-like galaxies. As found previously, the RL AGN fraction increases steeply with stellar velocity dispersion $\sigma_*$ and stellar mass, but even at a fixed velocity dispersion of 200-250 kms$^{-1}$ this fraction increases from 0.3% for flat galaxies (projected axis ratio of q $<$ 0.4) to 5% for round galaxies (q $>$ 0.8). We rule out that this strong trend is due to projection effects in the measured velocity dispersion. The large fraction of radio-loud AGN in massive, round galaxies is consistent with the hypothesis that such AGN deposit energy into their hot gaseous halos, preventing cooling and star-formation. However, the absence of such AGN in disk-like quiescent galaxies -- most of which are not satellites in massive clusters, raises important questions: is maintenance-mode feedback a generally valid explanation for quiescence; and, if so, how does that feedback avoid manifesting at least occasionally as a radio-loud galaxy?

The halo masses $M_{halo}$ of low surface brightness (LSB) galaxies are critical measurements for understanding their formation processes. One promising method to estimate a galaxy's $M_{halo}$ is to exploit the empirical scaling relation between $M_{halo}$ and the number of associated globular clusters ($N_{\mathrm{GC}}$). We use a Bayesian mixture model approach to measure $N_{\mathrm{GC}}$ for 175 LSB ($23\leq\left\langle \mu_{e,r} \right\rangle [\mathrm{mag\ arcsec}^{-2}]\leq 28$) galaxies in the Fornax cluster using the Fornax Deep Survey (FDS) data; this is the largest sample of low mass galaxies so-far analysed for this kind of study. The proximity of the Fornax cluster means that we can measure galaxies with much smaller physical sizes ($0.3\leq r_{e,r}\ [\mathrm{kpc}]\leq 9.5$) compared to previous studies of the GC systems of LSB galaxies, probing stellar masses down to $M_{*}\sim10^{5}\mathrm{M_{\odot}}$. The sample also includes \nudg\ ultra-diffuse galaxies (UDGs), with projected $r$-band half-light radii greater than 1.5 kpc. Our results are consistent with an extrapolation of the $M_{*}-M_{halo}$ relation predicted from abundance matching. In particular, our UDG measurements are consistent with dwarf sized halos, having typical masses between $10^{10}$ and $10^{11}\mathrm{M_{\odot}}$. Overall, our UDG sample is statistically indistinguishable from smaller LSB galaxies in the same magnitude range. We do not find any candidates likely to be as rich as some of those found in the Coma cluster. We suggest that environment might play a role in producing GC-rich LSB galaxies.

The formation of two stellar bars within a galaxy has proved challenging for numerical studies. It is yet not clear whether the inner bar is born via a star formation process promoted by gas inflow along the outer bar, or whether it is dynamically assembled from instabilities in a small-scale stellar disc. Observational constraints to these scenarios are scarce. We present a thorough study of the stellar content of two double-barred galaxies observed by the MUSE TIMER project, NGC 1291 and NGC 5850, combined with a two-dimensional multi-component photometric decomposition performed on the 3.6{\mu}m images from S4G. Our analysis confirms the presence of {\sigma}-hollows appearing in the stellar velocity dispersion distribution at the ends of the inner bars. Both galaxies host inner discs matching in size with the inner bars, suggestive of a dynamical formation for the inner bars from small-scale discs. The analysis of the star formation histories for the structural components shaping the galaxies provides constraints on the epoch of dynamical assembly of the inner bars, which took place >6.5 Gyr ago for NGC 1291 and >4.5 Gyr ago for NGC 5850. This implies that inner bars are long-lived structures.

We present and discuss the stellar kinematics and populations of the S0 galaxy FCC 170 (NGC 1381) in the Fornax cluster, using deep MUSE data from the Fornax 3D survey. We show the maps of the first four moments of the stellar line-of-sight velocity distribution and of the mass-weighted mean stellar age, metallicity and [Mg/Fe] abundance ratio. The high-quality MUSE stellar kinematic measurements unveil the structure of this massive galaxy: a nuclear disk, a bar seen as a boxy bulge with a clear higher-velocity-dispersion X shape, a fast-rotating and flaring thin disk and a slower rotating thick disk. Whereas their overall old age makes it difficult to discuss differences in the formation epoch between these components, we find a clear-cut distinction between metal-rich and less [Mg/Fe]-enhanced populations in the thin-disk, boxy-bulge and nuclear disk, and more metal-poor and [Mg/Fe]-enhanced stars in the thick disk. Located in the densest region of the Fornax cluster, where signs of tidal stripping have been recently found, the evolution of FCC 170 might have been seriously affected by its environment. We discuss the possibility of its "pre-processing" in a subgroup before falling into the present-day cluster, which would have shaped this galaxy a long time ago. The thick disk displays a composite star formation history, as a significant fraction of younger stars co-exist with the main older thick-disk population. The former sub-population is characterized by even lower-metallicity and higher-[Mg/Fe] values, suggesting that these stars formed later and faster in a less chemically evolved satellite, which was subsequently accreted. Finally, we discuss evidence that metal-rich and less [Mg/Fe]-enhanced stars were brought in the outer parts of the thick disk by the flaring of the thin disk.

Galaxies continuously reprocess their interstellar material. One can therefore expect changing dust grain properties in galaxies which have followed different evolutionary pathways. Determining the intrinsic dust grain mix of a galaxy helps in reconstructing its evolutionary history. Early-type galaxies occasionally display regular dust lanes in their central regions. Due to the relatively simple geometry and composition of their stellar bodies, these galaxies are ideal to disentangle dust mix variations from geometric effects. We therefore model the various components of such a galaxy (FCC 167). We reconstruct its recent history, and investigate the possible fate of the dust lane. MUSE and ALMA observations reveal a nested ISM structure. An ionised-gas disk pervades the central regions of FCC 167, including those occupied by the main dust lane. Inward of the dust lane, we also find a disk/ring of cold molecular gas where stars are forming and HII regions contribute to the ionised-gas emission. Further in, the gas ionisation points towards an active galactic nucleus and the fuelling of a central supermassive black hole from its surrounding ionised and molecular reservoir. Observational constraints and radiative transfer models suggest the dust and gas are distributed in a ring-like geometry and the dust mix lacks small grains. The derived dust destruction timescales from sputtering in hot gas are short and we conclude that the dust must be strongly self-shielding and clumpy, or will quickly be eroded and disappear. Our findings show how detailed analysis of individual systems can complement statistical studies of dust-lane ETGs.

This paper is based on the multi-band (ugri) Fornax Deep Survey (FDS) with the VLT Survey Telescope (VST). We study bright early-type galaxies (m_B<15 mag) in the 9 square degrees around the core of the Fornax cluster, which covers the virial radius. The main goal of the present work is to provide the analysis of the light distribution for all galaxies out to unprecedented limits (in radius and surface brightness) and to release the main products resulting from this analysis in all FDS bands. From the isophote fit, we derive the azimuthally averaged surface brightness profiles, the position angle and ellipticity profiles as a function of the semi-major axis. In each band, we derive the total magnitudes, effective radii, integrated colors and stellar mass-to-light ratios. The long integration times, the arcsec-level angular resolution of OmegaCam@VST and the large covered area of FDS allow us to map the light and color distributions out to large galactocentric distances (up to about 10-15 R_e) and surface brightness levels beyond mu_r = 27 mag/arcsec^2 (mu_B > 28 mag/arcsec^2). Therefore, the new FDS data allow us to explore in great detail the morphology and structure of cluster galaxies out to the region of the stellar halo. The observations suggest that the Fornax cluster is not completely relaxed inside the virial radius. The bulk of the gravitational interactions between galaxies happens in the W-NW core region of the cluster, where most of the bright early-type galaxies are located and where the intra-cluster baryons (diffuse light and GCs) are found. We suggest that the W-NW sub-clump of galaxies results from an infalling group onto the cluster, which has modified the structure of the galaxy outskirts (making asymmetric stellar halos) and has produced the intra-cluster baryons (ICL and GCs), concentrated in this region of the cluster.

Double bars are thought to be important features for secular evolution in the central regions of galaxies. However, observational evidence about their origin and evolution is still scarce. We report on the discovery of the first Box-Peanut (B/P) structure in an inner bar detected in the face-on galaxy NGC 1291. We use the integral field data obtained from the MUSE spectrograph within the TIMER project. The B/P structure is detected as bi-symmetric minima of the $h_4$ moment of the line-of-sight velocity distribution along the major axis of the inner bar, as expected from numerical simulations. Our observations demonstrate that inner bars can follow a similar evolutionary path as outer bars, undergoing buckling instabilities. They also suggest that inner bars are long-lived structures, thus imposing tight constraints to their possible formation mechanisms

We analyze radial stellar metallicity and kinematic profiles out to 1Re in 244 CALIFA galaxies ranging from morphological type E to Sd, to study the evolutionary mechanisms of stellar population gradients. We find that linear metallicity gradients exhibit a clear correlation with galaxy morphological type - with early type galaxies showing the steepest gradients. We show that the metallicity gradients simply reflect the local mass-metallicity relation within a galaxy. This suggests that the radial stellar population distribution within a galaxys effective radius is primarily a result of the \emph{in-situ} local star formation history. In this simple picture, the dynamically derived stellar surface mass density gradient directly predicts the metallicity gradient of a galaxy. We show that this correlation and its scatter can be reproduced entirely by using independent empirical galaxy structural and chemical scaling relations. Using Schwarzschild dynamical models, we also explore the link between a galaxys local stellar populations and their orbital structures. We find that galaxies angular momentum and metallicity gradients show no obvious causal link. This suggests that metallicity gradients in the inner disk are not strongly shaped by radial migration, which is confirmed by the lack of correlation between the metallicity gradients and observable probes of radial migration in the galaxies, such as bars and spiral arms. Finally, we find that galaxies with positive metallicity gradients become increasingly common towards low mass and late morphological types - consistent with stellar feedback more efficiently modifying the baryon cycle in the central regions of these galaxies.

The Fornax Deep Survey (FDS), an imaging survey in the u', g', r', and i'-bands, has a supreme resolution and image depth compared to the previous spatially complete Fornax Cluster Catalog (FCC). Our new data allows us to study the galaxies down to r'-band magnitude m$_{r'}\approx$21 mag (M$_{r'}\approx$-10.5 mag). These data provide an important legacy dataset to study the Fornax cluster. We aim to present the Fornax Deep Survey (FDS) dwarf galaxy catalog, focusing on explaining the data reduction and calibrations, assessing the quality of the data, and describing the methods used for defining the cluster memberships for the catalog objects. As a first step we used the SExtractor fine-tuned for dwarf galaxy detection, to find galaxies from the FDS data, covering a 26 deg$^2$ area of the main cluster, and the area around the Fornax A substructure. We made 2D-decompositions of the identified galaxies using GALFIT. We used color-magnitude, luminosity-radius and luminosity-concentration relations to separate the cluster galaxies from the background galaxies. We then divided the cluster galaxies into early- and late-type galaxies according to their morphology and gave first order morphological classifications. Our final catalog includes 14,095 galaxies. We classify 590 galaxies as being likely Fornax cluster galaxies, of which 564 are dwarfs (M$_{r'}$ > -18.5 mag) consisting our Fornax dwarf catalog. Of the cluster dwarfs we classify 470 as early-types, and 94 as late-type galaxies. Our final catalog reaches its 50% completeness limit at magnitude M$_{r'}$ = -10.5 mag and surface brightness $\bar{\mu}_{e,r'}$ = 26 mag arcsec-2, which is approximately three magnitudes deeper than the FCC. Based on previous works and comparison with a spectroscopically confirmed subsample, we estimate that our final Fornax dwarf galaxy catalog has < 10% contamination from the background objects.

The Time Inference with MUSE in Extragalactic Rings (TIMER) project is a survey with the VLT-MUSE integral-field spectrograph of 24 nearby barred galaxies with prominent central structures (e.g., nuclear rings or inner discs). The main goals of the project are: (i) estimating the cosmic epoch when discs of galaxies settle, leading to the formation of bars; (ii) testing the hypothesis whereby discs in more massive galaxies are assembled first; and (iii) characterising the history of external gas accretion in disc galaxies. We present details on the sample selection, observations, data reduction, and derivation of high-level data products, including stellar kinematics, ages and metallicities. We also derive star formation histories and physical properties and kinematics of ionised gas. We illustrate how this dataset can be used for a plethora of scientific applications, e.g., stellar feedback, outflows, nuclear and primary bars, stellar migration and chemical enrichment, and the gaseous and stellar dynamics of nuclear spiral arms, barlenses, box/peanuts and bulges. Amongst our first results - based on a few selected galaxies -, we show that the dynamics of nuclear rings and inner discs is consistent with the picture in which they are formed by bars, that the central few hundred parsecs in massive disc galaxies tend to show a pronounced peak in stellar metallicity, and that nuclear rings can efficiently prevent star formation in this region. Finally, we present evidence that star-bursting nuclear rings can be fed with low-metallicity gas from low-mass companions.

State-of-the-art integral field surveys like $\mathrm{ATLAS^{3D}}$, SLUGGS, CALIFA, SAMI, and MaNGA provide large data sets of kinematical observations of early-type galaxies (ETGs), yielding constraints on the formation of ETGs. Using the cosmological hydrodynamical \textit{Magneticum Pathfinder} simulations, we investigate the paradigm of fast and slow rotating ETGs in a fully cosmological context. We show that the ETGs within the \textit{Magneticum} simulation are in remarkable agreement with the observations, revealing fast and slow rotators quantified by the angular momentum proxy $\lambda_{\mathrm{R}}$ and the flattening $\epsilon$ with the observed prevalence. Taking full advantage of the three-dimensional data, we demonstrate that the dichotomy between fast and slow rotating galaxies gets enhanced, showing an upper and lower population separated by an underpopulated region in the edge-on $\lambda_{\mathrm{R}}$-$\epsilon$ plane. Following the evolution of the $\lambda_{\mathrm{R}}$-$\epsilon$ plane through cosmic time, we find that, while the upper population is already in place at $z=2$, the lower population gets statistically significant below $z=1$ with a gradual increase. At least $50\%$ of the galaxies transition from fast to slow rotators on a short timescale, in most cases associated to a significant merger event. Furthermore, we connect the $M_{*}$-$j_{*}$ plane, quantified by the $b$-value, with the $\lambda_{\mathrm{R}}$-$\epsilon$ plane, revealing a strong correlation between the position of a galaxy in the $\lambda_{\mathrm{R}}$-$\epsilon$ plane and the $b$-value. Going one step further, we classify our sample based on features in their velocity map, finding all five common kinematic groups, also including the recently observed group of prolate rotators, populating distinct regions in the $\lambda_{\mathrm{R}}$-$b$ plane.

We used ionized gas and stellar kinematics for 667 spatially resolved galaxies publicly available from the Calar Alto Legacy Integral Field Area survey (CALIFA) 3rd Data Release with the aim of studying kinematic scaling relations as the Tully $\&$ Fisher (TF) relation using rotation velocity, $V_{rot}$, the Faber $\&$ Jackson (FJ) relation using velocity dispersion, $\sigma$, and also a combination of $V_{rot}$ and $\sigma$ through the $S_{K}$ parameter defined as $S_{K}^2 = KV_{rot}^2 + \sigma^2$ with constant $K$. Late-type and early-type galaxies reproduce the TF and FJ relations. Some early-type galaxies also follow the TF relation and some late-type galaxies the FJ relation, but always with larger scatter. On the contrary, when we use the $S_{K}$ parameter, all galaxies, regardless of the morphological type, lie on the same scaling relation, showing a tight correlation with the total stellar mass, $M_\star$. Indeed, we find that the scatter in this relation is smaller or equal to that of the TF and FJ relations. We explore different values of the $K$ parameter without significant differences (slope and scatter) in our final results with respect the case $K=0.5$ besides than a small change in the zero point. We calibrate the kinematic $S_{K}^2$ dynamical mass proxy in order to make it consistent with sophisticated published dynamical models within $0.15\ dex$. We show that the $S_{K}$ proxy is able to reproduce the relation between the dynamical mass and the stellar mass in the inner regions of galaxies. Our result may be useful in order to produce fast estimations of the central dynamical mass in galaxies and to study correlations in large galaxy surveys.

Based on the stellar orbit distribution derived from orbit-superposition Schwarzschild models, we decompose each of 250 representative present-day galaxies into four orbital components: cold with strong rotation, warm with weak rotation, hot with dominant random motion and counter-rotating (CR). We rebuild the surface brightness ($\Sigma$) of each orbital component and we present in figures and tables a quantification of their morphologies using the Sersic index \textit{n}, concentration $C = \log{(\Sigma_{0.1R_e}/\Sigma_{R_e})}$ and intrinsic flattening $q_{\mathrm{Re}}$ and $q_{\mathrm{Rmax}}$, with $R_e$ the half-light-radius and $R_{\mathrm{max}}$ the CALIFA data coverage. We find that: (1) kinematic hotter components are generally more concentrated and rounder than colder components, and (2) all components become more concentrated and thicker/rounder in more massive galaxies; they change from disk-like in low mass late-type galaxies to bulge-like in high-mass early type galaxies. Our findings suggest that Sersic \textit{n} is not a good discriminator between rotating bulges and non-rotating bulges. The luminosity fraction of cold orbits $f_{\rm cold}$ is well correlated with the photometrically-decomposed disk fraction $f_{\rm disk}$ as $f_{\mathrm{cold}} = 0.14 + 0.23f_{\mathrm{\mathrm{disk}}}$. Similarly, the hot orbit fraction $f_{\rm hot}$ is correlated with the bulge fraction $f_{\rm bulge}$ as $f_{\mathrm{hot}} = 0.19 + 0.31f_{\mathrm{\mathrm{bulge}}}$. The warm orbits mainly contribute to disks in low-mass late-type galaxies, and to bulges in high-mass early-type galaxies. The cold, warm, and hot components generally follow the same morphology ($\epsilon = 1-q_{\rm Rmax}$) versus kinematics ($\sigma_z^2/\overline{V_{\mathrm{tot}}^2}$) relation as the thin disk, thick disk/pseudo bulge, and classical bulge identified from cosmological simulations.

We present a comparative study of molecular and ionized gas kinematics in nearby galaxies. These results are based on observations from the EDGE survey, which measured spatially resolved $^{12}$CO(J=1-0) in 126 nearby galaxies. Every galaxy in EDGE has corresponding resolved ionized gas measurements from CALIFA. Using a sub-sample of 17 rotation dominated, star-forming galaxies where precise molecular gas rotation curves could be extracted, we derive CO and H$\alpha$ rotation curves using the same geometric parameters out to $\gtrsim$1 $R_e$. We find that $\sim$75% of our sample galaxies have smaller ionized gas rotation velocities than the molecular gas in the outer part of the rotation curve. In no case is the molecular gas rotation velocity measurably lower than that of the ionized gas. We suggest that the lower ionized gas rotation velocity can be attributed to a significant contribution from extraplanar diffuse ionized gas in a thick, turbulence supported disk. Using observations of the H$\gamma$ transition also available from CALIFA, we measure ionized gas velocity dispersions and find that these galaxies have sufficiently large velocity dispersions to support a thick ionized gas disk. Kinematic simulations show that a thick disk with a vertical rotation velocity gradient can reproduce the observed differences between the CO and H$\alpha$ rotation velocities. Observed line ratios tracing diffuse ionized gas are elevated compared to typical values in the midplane of the Milky Way. In galaxies affected by this phenomenon, dynamical masses measured using ionized gas rotation curves will be systematically underestimated.

The Fornax cluster provides a uniquely compact laboratory to study the detailed history of early-type galaxies and the role played by environment in driving their evolution and their transformation from late-type galaxies. Using the superb capabilities of the Multi Unit Spectroscopic Explorer on the Very Large Telescope, high-quality integral-field spectroscopic data were obtained for the inner regions of all the bright ($m_B\leq15$) galaxies within the virial radius of Fornax. The stellar haloes of early-type galaxies are also covered out to about four effective radii. State-of-the-art stellar dynamical and population modelling allows to aim in particular at better characterising the disc components of fast-rotating early-type galaxies, constraining radial variations in the stellar initial-mass functions and measuring the stellar age, metallicity, and $\alpha$-element abundance of stellar haloes in cluster galaxies. This paper describes the sample selection, observations, and overall goals of the survey, and provides initial results based on the spectroscopic data, including the detailed characterisation of stellar kinematics and populations to large radii; decomposition of galaxy components directly via their orbital structure; the ability to identify globular clusters and planetary nebulae, and derivation of high-quality emission-line diagnostics in the presence of complex ionised gas.

We present the results of a wide spectroscopic survey aimed at detecting extragalactic globular clusters (GCs) in the core of the Fornax cluster. About 4500 low resolution spectra (from 4800 to 10000 \AA ) were observed in 25 VLT/VIMOS masks covering the central 1 deg$^{2}$ around the dominant galaxy NGC 1399 corresponding to $\sim$175 kpc galactocentric radius. We describe the methodology used for data reduction and data analysis. We found a total of 387 unique physical objects (372 GCs and 15 ultra compact dwarfs) in the field covered by our observations. Most of these objects are associated with NGC 1399, with only 10% likely belonging to other giant galaxies. The new VIMOS dataset is complementary to the many GC catalogues already present in the literature and it brings the total number of tracer particles around NGC 1399 to more than 1130 objects. With this comprehensive radial velocity sample we have found that the velocity dispersion of the GC population (equally for red and blue GC populations) shows a relatively sharp increase from low velocity dispersion ($\sim250$-$350$ kms$^{-1}$) to high velocity dispersion ($\sim300$-$400$ kms$^{-1}$) at projected radius of $\approx10$ arcmin ($\sim60$ kpc) from the galaxy centre. This suggests that at a projected radius of $\approx60$ kpc both blue and red GC populations begin to be governed by the dominating Fornax cluster potential, rather than by the central NGC 1399 galaxy potential. This kinematic evidence corroborates similar results found using surface brightness analysis and planetary nebulae kinematics.

Deriving circular velocities of galaxies from stellar kinematics can provide an estimate of their total dynamical mass, provided a contribution from the velocity dispersion of the stars is taken into account. Molecular gas (e.g., CO) on the other hand, is a dynamically cold tracer and hence acts as an independent circular velocity estimate without needing such a correction. In this paper we test the underlying assumptions of three commonly used dynamical models, deriving circular velocities from stellar kinematics of 54 galaxies (S0-Sd) that have observations of both stellar kinematics from the CALIFA survey, and CO kinematics from the EDGE survey. We test the Asymmetric Drift Correction (ADC) method, as well as Jeans, and Schwarzschild models. The three methods each reproduce the CO circular velocity at 1Re to within 10%. All three methods show larger scatter (up to 20%) in the inner regions (R < 0.4Re) which may be due to an increasingly spherical mass distribution (which is not captured by the thin disk assumption in ADC), or non-constant stellar M/L ratios (for both the JAM and Schwarzschild models). This homogeneous analysis of stellar and gaseous kinematics validates that all three models can recover Mdyn at 1Re to better than 20%, but users should be mindful of scatter in the inner regions where some assumptions may break down.

In this second paper of T$h$e role of $E$nvironment in shaping $L$ow-mass $E$arly-type $N$earby g$a$laxies (hELENa) series we study [Mg/Fe] abundance distribution trends of early-type galaxies observed with the SAURON integral field unit, spanning a wide range in mass and local environment densities: 20 low-mass early-types (dEs) of Sybilska et al. (2017) and 258 massive early types (ETGs) of the $ATLAS^{3D}$ project, all homogeneously reduced and analyzed. We show that the [Mg/Fe] ratios scale with velocity dispersion ($\sigma$) at fixed [Fe/H] and that they evolve with [Fe/H] along similar paths for all early-types, grouped in bins of increasing local and global {$\sigma$}, as well as the second velocity moment $V_{rms}$, indicating a common inside-out formation pattern. We then place our dEs on the [Mg/Fe] $vs.$ [Fe/H] diagram of Local Group galaxies and show that dEs occupy the same region and show a similar trend line slope in the diagram as the high-metallicity stars of the Milky Way and the Large Magellanic Cloud. This finding extends the similar trend found for dwarf spheroidal $vs.$ dwarf irregular galaxies and supports the notion that dEs have evolved from late-type galaxies that have lost their gas at a point of their evolution, which likely coincided with them entering denser environments.

To investigate star formation and assembly processes of massive galaxies, we present here a spatially-resolved stellar populations analysis of a sample of 45 elliptical galaxies (Es) selected from the CALIFA survey. We find rather flat age and [Mg/Fe] radial gradients, weakly dependent on the effective velocity dispersion of the galaxy within half-light radius. However, our analysis shows that metallicity gradients become steeper with increasing galaxy velocity dispersion. In addition, we have homogeneously compared the stellar populations gradients of our sample of Es to a sample of nearby relic galaxies, i.e., local remnants of the high-z population of red nuggets. This comparison indicates that, first, the cores of present-day massive galaxies were likely formed in gas-rich, rapid star formation events at high redshift (z>2). This led to radial metallicity variations steeper than observed in the local Universe, and positive [Mg/Fe] gradients. Second, our analysis also suggests that a later sequence of minor dry mergers, populating the outskirts of early-type galaxies (ETGs), flattened the pristine [Mg/Fe] and metallicity gradients. Finally, we find a tight age-[Mg/Fe] relation, supporting that the duration of the star formation is the main driver of the [Mg/Fe] enhancement in massive ETGs. However, the star formation time-scale alone is not able to fully explain our [Mg/Fe] measurements. Interestingly, our results match the expected effect that a variable stellar initial mass function would have on the [Mg/Fe] ratio.

The stellar velocity ellipsoid (SVE) in galaxies can provide important information on the processes that participate in the dynamical heating of their disc components (e.g. giant molecular clouds, mergers, spiral density waves, bars). Earlier findings suggested a strong relation between the shape of the disc SVE and Hubble type, with later-type galaxies displaying more anisotropic ellipsoids and early-types being more isotropic. In this paper, we revisit the strength of this relation using an exhaustive compilation of observational results from the literature on this issue. We find no clear correlation between the shape of the disc SVE and morphological type, and show that galaxies with the same Hubble type display a wide range of vertical-to-radial velocity dispersion ratios. The points are distributed around a mean value and scatter of $\sigma_z/\sigma_R=0.7\pm 0.2$. With the aid of numerical simulations, we argue that different mechanisms might influence the shape of the SVE in the same manner and that the same process (e.g. mergers) does not have the same impact in all the galaxies. The complexity of the observational picture is confirmed by these simulations, which suggest that the vertical-to-radial axis ratio of the SVE is not a good indicator of the main source of disc heating. Our analysis of those simulations also indicates that the observed shape of the disc SVE may be affected by several processes simultaneously and that the signatures of some of them (e.g. mergers) fade over time.

Galaxy formation entails the hierarchical assembly of mass, along with the condensation of baryons and the ensuing, self-regulating star formation. The stars form a collisionless system whose orbit distribution retains dynamical memory that can constrain a galaxy's formation history. The ordered-rotation dominated orbits with near maximum circularity $\lambda_z \simeq1$ and the random-motion dominated orbits with low circularity $\lambda_z \simeq0$ are called kinematically cold and kinematically hot, respectively. The fraction of stars on `cold' orbits, compared to the fraction of stars on `hot' orbits, speaks directly to the quiescence or violence of the galaxies' formation histories. Here we present such orbit distributions, derived from stellar kinematic maps via orbit-based modelling for a well defined, large sample of 300 nearby galaxies. The sample, drawn from the CALIFA survey, includes the main morphological galaxy types and spans the total stellar mass range from $10^{8.7}$ to $10^{11.9}$ solar masses. Our analysis derives the orbit-circularity distribution as a function of galaxy mass, $p(\lambda_z~|~M_\star)$, and its volume-averaged total distribution, $p(\lambda_z)$. We find that across most of the considered mass range and across morphological types, there are more stars on `warm' orbits defined as $0.25\le \lambda_z \le 0.8$ than on either `cold' or `hot' orbits. This orbit-based "Hubble diagram" provides a benchmark for galaxy formation simulations in a cosmological context.

Super-massive black holes, with masses larger than a million times that of the Sun, appear to inhabit the centers of all massive galaxies. Cosmologically-motivated theories of galaxy formation need feedback from these super-massive black holes to regulate star formation. In the absence of such feedback, state-of-the-art numerical simulations dramatically fail to reproduce the number density and properties of massive galaxies in the local Universe. However, there is no observational evidence of this strongly coupled co-evolution between super-massive black holes and star formation, impeding our understanding of baryonic processes within galaxies. Here we show that the star formation histories (SFHs) of nearby massive galaxies, as measured from their integrated optical spectra, depend on the mass of the central super-massive black hole. Our results suggest that black hole mass growth scales with gas cooling rate in the early Universe. The subsequent quenching of star formation takes place earlier and more efficiently in galaxies hosting more massive central black holes. The observed relation between black hole mass and star formation efficiency applies to all generations of stars formed throughout a galaxy's life, revealing a continuous interplay between black hole activity and baryon cooling.

We determine abundance ratios of 37 dwarf ellipticals (dEs) in the nearby Virgo cluster. This sample is representative of the early-type population of galaxies in the absolute magnitude range -19.0 < Mr < -16.0. We analyze their absorption line-strength indices by means of index-index diagrams and scaling relations and use the stellar population models to interpret them. We present ages, metallicities and abundance ratios obtained from these dEs within an aperture size of Re/8. We calculate [Na/Fe] from NaD, [Ca/Fe] from Ca4227 and [Mg/Fe] from Mgb. We find that [Na/Fe] is under-abundant with respect to solar while [Mg/Fe] is around solar. This is exactly opposite to what is found for giant ellipticals, but follows the trend with metallicity found previously for the Fornax dwarf NGC 1396. We discuss possible formation scenarios that can result in such elemental abundance patterns and we speculate that dEs have disk-like SFH favouring them to originate from late-type dwarfs or small spirals. Na-yields appear to be very metal-dependent, in agreement with studies of giant ellipticals, probably due to the large dependence on the neutron-excess in stars. We conclude that dEs have undergone a considerable amount of chemical evolution, they are therefore not uniformly old, but have extended SFH, similar to many of the Local Group galaxies.

In the core of the Fornax cluster, on the West side of NGC1399, we have detected a previously unknown region of intra-cluster light (ICL). It is made up by several faint ($\mu_r \simeq 28 - 29$~mag/arcsec$^2$) {\it patches} of diffuse light. The bulk of the ICL is located in between the three bright galaxies in the core, NGC1387, NGC1379 and NGC1381, at $10\leq R \leq40$~arcmin ($\sim 58 - 230$~kpc) from the central galaxy NGC~1399. We show that the ICL is the counterpart in the diffuse light of the known over-density in the population of blue globular clusters (GCs). The total g-band luminosity of the ICL is $L_g\simeq 8.3 \times 10^{9}$ $L_{\odot}$, which is $\sim5\%$ of the total luminosity of NGC1399. This is consistent with the fraction of the blue GCs in the same region of the cluster. The ICL has $g-r \sim 0.7$~mag, which is similar to the colors in the halo of the bright galaxies in the cluster core. The new findings were compared with theoretical predictions for the ICL formation and they support a scenario in which the intra-cluster population detected in the core of the Fornax cluster is build up by the tidal stripping of material (stars and GCs) from galaxy outskirts in a close passage with the cD. Moreover, the diffuse form of the ICL and its location close to the core of the cluster is expected in a dynamically evolved cluster as Fornax.

We analyze the globular cluster (GC) systems in two very different galaxies, N3115 and N1399. With the papers of this series, we aim at highlighting common and different properties in the GC systems in galaxies covering a wide range of parameter space. We compare the GCs in N3115 and N1399 as derived from the analysis of u, g, and i band images taken with the VST telescope as part of the VEGAS and FDS surveys. We selected GC candidates using as reference the morpho-photometric and color properties of confirmed GCs. The surface density of GCs in N3115 reveal a morphology similar to the profile of field stars; the same is true when blue and red GCs are taken separately. The GC maps for N1399 are richer in structure and confirm the existence of an intracluster GC component. We confirm the presence of a spatial offset in the N1399 GC centroid and find that the centroid of the GCs for N3115 coincides well with the galaxy center. Both GC systems show unambiguous color bimodality in g-i and u-i. The azimuthal average of the radial density profiles in both galaxies reveals a larger spatial extent for the total GCs population with respect to the galaxy surface brightness profile. As for the specific frequency of GCs, Sn, we find it is a factor of two higher in N1399 than for N3115. By inspecting the radial behavior of the specific frequency, Sn(<r), for the total, blue, and red GCs, we find notable similarities between the trends for red GCs in the two targets. In spite of extremely different host environments, the red GCs in both cases appear closely linked to the light distribution of field stars. Blue GCs extend to larger galactocentric scales than red GCs, marking a significant difference between the two galaxies: the blue/red GCs and field stellar components of N3115 appear well thermalized with each other and the blue GCs in N1399 appear to fade into an unrelaxed intracluster GC population.

We study a sample of 28 S0 galaxies extracted from the integral-field spectroscopic (IFS) survey CALIFA. We combine an accurate two-dimensional (2D) multi-component photometric decomposition with the IFS kinematic properties of their bulges to understand their formation scenario. Our final sample is representative of S0s with high stellar masses ($M_{star}/M_{sun} > 10^{10}$). They lay mainly on the red sequence and live in relatively isolated environments similar to that of the field and loose groups. We use our 2D photometric decomposition to define the size and photometric properties of the bulges, as well as their location within the galaxies. We perform mock spectroscopic simulations mimicking our observed galaxies to quantify the impact of the underlying disc on our bulge kinematic measurements ($\lambda$ and $v/\sigma$). We compare our bulge corrected kinematic measurements with the results from Schwarzschild dynamical modelling. The good agreement confirms the robustness of our results and allows us to use bulge reprojected values of $\lambda$ and $v/\sigma$. We find that the photometric ($n$ and $B/T$) and kinematic ($v/\sigma$ and $\lambda$) properties of our field S0 bulges are not correlated. We demonstrate that this morpho-kinematic decoupling is intrinsic to the bulges and it is not due to projection effects. We conclude that photometric diagnostics to separate different types of bulges (disc-like vs classical) might not be useful for S0 galaxies. The morpho-kinematics properties of S0 bulges derived in this paper suggest that they are mainly formed by dissipation processes happening at high redshift, but dedicated high-resolution simulations are necessary to better identify their origin.

We present new proper motion (PM) measurements of the dwarf spheroidal galaxies (dSphs) Draco and Sculptor using multi-epoch images obtained with the Hubble Space Telescope ACS/WFC. Our PM results have uncertainties far lower than previous measurements, even made with the same instrument. The PM results for Draco and Sculptor are (mu_W,mu_N)_Dra = (-0.0562+/-0.0099,-0.1765+/-0.0100) mas/yr and (mu_W,mu_N)_Scl = (-0.0296+/-0.0209,-0.1358 +/-0.0214) mas/yr. The implied Galactocentric velocity vectors for Draco and Sculptor have radial and tangential components: (V_rad,V_tan)_Dra = (-88.6,161.4) +/- (4.4,5.6) km/s; and (V_rad,V_tan)_Scl = (72.6,200.2) +/- (1.3,10.8) km/s. We study the detailed orbital history of both Draco and Sculptor via numerical orbit integrations. Orbital periods of Draco and Sculptor are found to be 1-2 and 2-5 Gyrs, respectively, accounting for uncertainties in the MW mass. We also study the influence of the Large Magellanic Cloud (LMC) on the orbits of Draco and Sculptor. Overall, the inclusion of the LMC increases the scatter in the orbital results. Based on our calculations, Draco shows a rather wide range of orbital parameters depending on the MW mass and inclusion/exclusion of the LMC, but Sculptor's orbit is very well constrained with its most recent pericentric approach to the MW being 0.3-0.4 Gyr ago. Our new PMs imply that the orbital trajectories of both Draco and Sculptor are confined within the Disk of Satellites (DoS), better so than implied by earlier PM measurements, and likely rule out the possibility that these two galaxies were accreted together as part of a tightly bound group.

Schwarzschild orbit-based dynamical models are widely used to uncover the internal dynamics of early-type galaxies and globular clusters. Here we present for the first time the Schwarzschild models of late-type galaxies: an SBb galaxy NGC 4210 and an S0 galaxy NGC 6278 from the CALIFA survey. The mass profiles within $2\,R_e$ are constrained well with $1\sigma$ statistical error of $\sim 10\%$. The luminous and dark mass can be disentangled with uncertainties of $\sim 20\%$ and $\sim 50\%$ respectively. From $R_e$ to $2\,R_e$, the dark matter fraction increases from $14\pm10\%$ to $18\pm10\%$ for NGC 4210 and from $15\pm10\%$ to $30\pm20\%$ for NGC 6278. The velocity anisotropy profiles of both $\sigma_r/\sigma_t$ and $\sigma_z/\sigma_R$ are well constrained. The inferred internal orbital distributions reveal clear substructures. The orbits are naturally separated into three components: a cold component with near circular orbits, a hot component with near radial orbits, and a warm component in between. The photometrically-identified exponential disks are predominantly made up of cold orbits only beyond $\sim 1\,R_e$, while they are constructed mainly with the warm orbits inside. Our dynamical hot components are concentrated in the inner regions, similar to the photometrically-identified bulges. The reliability of the results, especially the orbit distribution, are verified by applying the model to mock data.

We present H$\alpha$ spectroscopic observations and detailed modelling of the Balmer filaments in the supernova remnant Tycho. We used Galaxy H$\alpha$ Fabry-P\'erot Spectrometer on the William Herschel Telescope with a 3.4'$\times$3.4' field-of-view, 0.2" pixel scale and $\sigma_\rm{instr}=8.1$ km/s resolution at 1" seeing for $\sim10$ hours, resulting in 82 spatial-spectral bins that resolve the narrow H$\alpha$ line in the entire Tycho's northeastern rim. For the first time, we can mitigate artificial line broadening from unresolved differential motion, and probe H$\alpha$ emission parameters in varying shock and ambient medium conditions. Broad H$\alpha$ line remains unresolved within spectral coverage of 392 km/s. We employed Bayesian inference to obtain reliable parameter confidence intervals, and quantify the evidence for models with multiple line components. The median H$\alpha$ narrow-line full-width at half-maximum of all bins and models is $W_\rm{NL}=(54.8\pm1.8)$ km/s at the $95\%$ confidence level, varying within [35, 72] km/s between bins and clearly broadened compared to the intrinsic (thermal) $\approx20$ km/s. Possible line splits are accounted for, significant in $\approx18\%$ of the filament, and presumably due to remaining projection effects. We also find wide-spread evidence for intermediate-line emission of a broad-neutral precursor, with median $W_\rm{IL}=(180\pm14)$ km/s ($95\%$ confidence). Finally, we present a measurement of the remnant's systemic velocity, $V_\rm{LSR}=-34$ km/s, and map differential line-of-sight motions. Our results confirm the existence and interplay of shock precursors in Tycho's remnant. In particular, we show that suprathermal narrow-line emission is near-universal in Tycho and that, in absence of an alternative explanation, collisionless supernova remnant shocks constitute a viable acceleration source for Galactic TeV Cosmic-Ray protons.

We present wide-field, spatially and highly resolved spectroscopic observations of Balmer filaments in the northeastern rim of Tycho's supernova remnant in order to investigate the signal of cosmic-ray (CR) acceleration. The spectra of Balmer-dominated shocks (BDSs) have characteristic narrow (FWHM $\sim$ 10 kms$^{-1}$) and broad (FWHM $\sim$ 1000 kms$^{-1}$) H$\alpha$ components. CRs affect the H$\alpha$-line parameters: heating the cold neutrals in the interstellar medium results in broadening of the narrow H$\alpha$-line width beyond 20 kms$^{-1}$, but also in reduction of the broad H$\alpha$-line width due to energy being removed from the protons in the post-shock region. For the first time we show that the width of the narrow H$\alpha$ line, much larger than 20 kms$^{-1}$, is not a resolution or geometric effect nor a spurious result of a neglected intermediate (FWHM $\sim$ 100 kms$^{-1}$) component resulting from hydrogen atoms undergoing charge exchange with warm protons in the broad-neutral precursor. Moreover, we show that a narrow line width $\gg$ 20 kms$^{-1}$ extends across the entire NE rim, implying CR acceleration is ubiquitous, and making it possible to relate its strength to locally varying shock conditions. Finally, we find several locations along the rim, where spectra are significantly better explained (based on Bayesian evidence) by inclusion of the intermediate component, with a width of 180 kms$^{-1}$ on average.

We present interferometric CO observations made with the Combined Array for Millimeter-wave Astronomy (CARMA) of galaxies from the Extragalactic Database for Galaxy Evolution survey (EDGE). These galaxies are selected from the Calar Alto Legacy Integral Field Area (CALIFA) sample, mapped with optical integral field spectroscopy. EDGE provides good quality CO data (3$\sigma$ sensitivity $\Sigma_{\rm mol}$ $\sim$ 11 M$_\odot$ pc$^{-2}$ before inclination correction, resolution $\sim1.4$ kpc) for 126 galaxies, constituting the largest interferometric CO survey of galaxies in the nearby universe. We describe the survey, the data characteristics, the data products, and present initial science results. We find that the exponential scale-lengths of the molecular, stellar, and star-forming disks are approximately equal, and galaxies that are more compact in molecular gas than in stars tend to show signs of interaction. We characterize the molecular to stellar ratio as a function of Hubble type and stellar mass, present preliminary results on the resolved relations between the molecular gas, stars, and star formation rate, and discuss the dependence of the resolved molecular depletion time on stellar surface density, nebular extinction, and gas metallicity. EDGE provides a key dataset to address outstanding topics regarding gas and its role in star formation and galaxy evolution, which will be publicly available on completion of the quality assessment.

Early-type dwarf galaxies are not simply featureless, old objects, but were found to be much more diverse, hosting substructures and a variety of stellar population properties. To explore the stellar content of faint early-type galaxies, and to investigate in particular those with recent central star formation, we study colours and colour gradients within one effective radius in optical (g-r) and near-infrared (i-H) bands for 120 Virgo cluster early types with -19 mag < $M_{r}$ < -16 mag. Twelve galaxies turn out to have blue cores, when defined as g-r colour gradients larger than 0.10 mag/$R_{\rm eff}$, which represents the positive tail of the gradient distribution. For these galaxies, we find that they have the strongest age gradients, and that even outside the blue core, their mean stellar population is younger than the mean of ordinary faint early types. The metallicity gradients of these blue-cored early-type dwarf galaxies are, however, in the range of most normal faint early types, which we find to have non-zero gradients with higher central metallicity. The blue central regions are consistent with star formation activity within the last few 100 Myr. We discuss that these galaxies could be explained by environmental quenching of star formation in the outer galaxy regions, while the inner star formation activity has continued.

We observed twelve nearby HI -detected early-type galaxies (ETGs) of stellar mass $\sim 10^{10}M\odot \leq M_* \leq \sim 10^{11}M\odot$ with the Mitchell Integral-Field Spectrograph, reaching approximately three half-light radii in most cases. We extracted line-of-sight velocity distributions for the stellar and gaseous components. We find little evidence of transitions in the stellar kinematics of the galaxies in our sample beyond the central effective radius, with centrally fast-rotating galaxies remaining fast-rotating and centrally slow-rotating galaxies likewise remaining slow-rotating. This is consistent with these galaxies having not experienced late dry major mergers; however, several of our objects have ionised gas that is misaligned with respect to their stars, suggesting some kind of past interaction. We extract Lick index measurements of the commonly-used H$\beta$, Fe5015, Mg\, b, Fe5270 and Fe5335 absorption features, and we find most galaxies to have flat H$\beta$ gradients and negative Mg, b gradients. We measure gradients of age, metallicity and abundance ratio for our galaxies using spectral fitting, and for the majority of our galaxies find negative age and metallicity gradients. We also find the stellar mass-to-light ratios to decrease with radius for most of the galaxies in our sample. Our results are consistent with a view in which intermediate-mass ETGs experience mostly quiet evolutionary histories, but in which many have experienced some kind of gaseous interaction in recent times.

We present new evidence for eight early-type galaxies (ETGs) from the CALIFA Survey that show clear rotation around their major photometric axis ("prolate rotation"). These are LSBCF560-04, NGC 0647, NGC 0810, NGC 2484, NGC 4874, NGC 5216, NGC 6173 and NGC 6338. Including NGC 5485, a known case of an ETG with stellar prolate rotation, as well as UGC 10695, a further possible candidate for prolate rotation, we report ten CALIFA galaxies in total that show evidence for such a feature in their stellar kinematics. Prolate rotators correspond to ~9% of the volume-corrected sample of CALIFA ETGs, a fraction much higher than previously reported. We find that prolate rotation is more common among the most massive ETGs. We investigate the implications of these findings by studying N-body merger simulations, and show that a prolate ETG with rotation around its major axis could be the result of a major polar merger, with the amplitude of prolate rotation depending on the initial bulge-to-total stellar mass ratio of its progenitor galaxies. Additionally, we find that prolate ETGs resulting from this formation scenario show a correlation between their stellar line-of-sight velocity and higher order moment h_3, opposite to typical oblate ETGs, as well as a double peak of their stellar velocity dispersion along their minor axis. Finally, we investigate the origin of prolate rotation in polar galaxy merger remnants. Our findings suggest that prolate rotation in massive ETGs might be more common than previously expected, and can help towards a better understanding of their dynamical structure and formation origin.

Boxy/peanut bulges are considered to be part of the same stellar structure as bars and both could be linked through the buckling instability. The Milky Way is our closest example. The goal of this letter is determining if the mass assembly of the different components leaves an imprint in their stellar populations allowing to estimate the time of bar formation and its evolution. To this aim we use integral field spectroscopy to derive the stellar age distributions, SADs, along the bar and disc of NGC 6032. The analysis shows clearly different SADs for the different bar areas. There is an underlying old (>=12 Gyr) stellar population for the whole galaxy. The bulge shows star formation happening at all times. The inner bar structure shows stars of ages older than 6 Gyrs with a deficit of younger populations. The outer bar region presents a SAD similar to that of the disc. To interpret our results, we use a generic numerical simulation of a barred galaxy. Thus, we constrain, for the first time, the epoch of bar formation, the buckling instability period and the posterior growth from disc material. We establish that the bar of NGC 6032 is old, formed around 10 Gyr ago while the buckling phase possibly happened around 8 Gyr ago. All these results point towards bars being long-lasting even in the presence of gas.

We present the first in a series of papers in T$h$e role of $E$nvironment in shaping $L$ow-mass $E$arly-type $N$earby g$a$laxies (hELENa) project. In this paper we combine our sample of 20 low-mass early types (dEs) with 258 massive early types (ETGs) from the ATLAS$^{\mathrm{3D}}$ survey - all observed with the SAURON integral field unit (IFU) - to investigate early-type galaxies' stellar population scaling relations and the dependence of the population properties on local environment, extended to the low-{\sigma} regime of dEs. The ages in our sample show more scatter at lower {\sigma} values, indicative of less massive galaxies being affected by the environment to a higher degree. The shape of the age-{\sigma} relation for cluster vs. non-cluster galaxies suggests that cluster environment speeds up the placing of galaxies on the red sequence. While the scaling relations are tighter for cluster than for the field/group objects, we find no evidence for a difference in average population characteristics of the two samples. We investigate the properties of our sample in the Virgo cluster as a function of number density (rather than simple clustrocentric distance) and find that dE ages negatively correlate with the local density, likely because galaxies in regions of lower density are later arrivals to the cluster or have experienced less pre-processing in groups, and consequently used up their gas reservoir more recently. Overall, dE properties correlate more strongly with density than those of massive ETGs, which was expected as less massive galaxies are more susceptible to external influences.

The internal dynamics of globular clusters (GCs) is strongly affected by two-body interactions that bring the systems to a state of partial energy equipartition. Using a set of Monte Carlo clusters simulations, we investigate the role of the onset of energy equipartition in shaping the mass-to-light ratio (M/L) in GCs. Our simulations show that the M/L profiles cannot be considered constant and their specific shape strongly depends on the dynamical age of the clusters. Dynamically younger clusters display a central peak up to M/L $\simeq25$ $M_\odot/L_\odot$ caused by the retention of dark remnants; this peak flattens out for dynamically older clusters. Moreover, we find that also the global values of M/L correlate with the dynamical state of a cluster quantified as either the number of relaxation times a system has experienced $n_{rel}$ or the equipartition parameter $m_{eq}$: clusters closer to full equipartition (higher $n_{rel}$ or lower $m_{eq}$) display a lower M/L. We show that the decrease of M/L is primarily driven by the dynamical ejection of dark remnants, rather than by the escape of low-mass stars. The predictions of our models are in good agreement with observations of GCs in the Milky Way and M31, indicating that differences in relaxation state alone can explain variations of M/L up to a factor of $\simeq3$. Our characterization of the M/L as a function of relaxation state is of primary relevance for the application and interpretation of dynamical models.

Context: According to numerical simulations, stars are not always kept at their birth galactocentric distances but migrate. The importance of this radial migration in shaping galactic light distributions is still unclear. However, if it is indeed important, galaxies with different surface brightness (SB) profiles must display differences in their stellar population properties. Aims: We investigate the role of radial migration on the light distribution and the radial stellar content by comparing the inner colour, age and metallicity gradients for galaxies with different SB profiles. We define these inner parts avoiding the bulge and bar regions and up to around three disc scale-lengths (type I, pure exponential) or the break radius (type II, downbending; type III, upbending). Methods: We analyse 214 spiral galaxies from the CALIFA survey covering different SB profiles. We make use of GASP2D and SDSS data to characterise their light distribution and obtain colour profiles. The stellar age and metallicity profiles are computed using a methodology based on full-spectrum fitting techniques (pPXF, GANDALF, and STECKMAP) to the IFS CALIFA data. Results: The distributions of the colour, stellar age and stellar metallicity gradients in the inner parts for galaxies displaying different SB profiles are unalike as suggested by Kolmogorov-Smirnov and Anderson-Darling tests. We find a trend in which type II galaxies show the steepest profiles of all and type III the shallowest, with type I galaxies displaying an intermediate behaviour. Conclusions: These results are consistent with a scenario in which radial migration is more efficient for type III galaxies than for type I systems with type II galaxies presenting the lowest radial migration efficiency. In such scenario, radial migration mixes the stellar content flattening the radial stellar properties and... [abriged]

As part of the Fornax Deep Survey with the ESO VLT Survey Telescope, we present new $g$ and $r$ bands mosaics of the SW group of the Fornax cluster. It covers an area of $3 \times 2$ square degrees around the central galaxy NGC1316. The deep photometry, the high spatial resolution of OmegaCam and the large covered area allow us to study the galaxy structure, to trace stellar halo formation and look at the galaxy environment. We map the surface brightness profile out to 33arcmin ($\sim 200$kpc $\sim15R_e$) from the galaxy centre, down to $\mu_g \sim 31$ mag arcsec$^{-2}$ and $\mu_r \sim 29$ mag arcsec$^{-2}$. This allow us to estimate the scales of the main components dominating the light distribution, which are the central spheroid, inside 5.5 arcmin ($\sim33$ kpc), and the outer stellar envelope. Data analysis suggests that we are catching in act the second phase of the mass assembly in this galaxy, since the accretion of smaller satellites is going on in both components. The outer envelope of NGC1316 still hosts the remnants of the accreted satellite galaxies that are forming the stellar halo. We discuss the possible formation scenarios for NGC1316, by comparing the observed properties (morphology, colors, gas content, kinematics and dynamics) with predictions from cosmological simulations of galaxy formation. We find that {\it i)} the central spheroid could result from at least one merging event, it could be a pre-existing early-type disk galaxy with a lower mass companion, and {\it ii)} the stellar envelope comes from the gradual accretion of small satellites.

Dedicated photometric and spectroscopic surveys have provided unambiguous evidence for a strong stellar mass-size evolution of galaxies within the last 10 Gyr. The likely progenitors of today's most massive galaxies are remarkably small, disky, passive and have already assembled much of their stellar mass at redshift z=2. An in-depth analysis of these objects, however, is currently not feasible due to the lack of high-quality, spatially-resolved photometric and spectroscopic data. In this paper, we present a sample of nearby compact elliptical galaxies (CEGs), which bear resemblance to the massive and quiescent galaxy population at earlier times. Hubble Space Telescope (HST) and wide-field integral field unit (IFU) data have been obtained, and are used to constrain orbit-based dynamical models and stellar population synthesis (SPS) fits, to unravel their structural and dynamical properties. We first show that our galaxies are outliers in the present-day stellar mass-size relation. They are, however, consistent with the mass-size relation of compact, massive and quiescent galaxies at redshift z=2. The compact sizes of our nearby galaxies imply high central stellar mass surface densities, which are also in agreement with the massive galaxy population at higher redshift, hinting at strong dissipational processes during their formation. Corroborating evidence for a largely passive evolution within the last 10 Gyr is provided by their orbital distribution as well as their stellar populations, which are difficult to reconcile with a very active (major) merging history. This all supports that we can use nearby CEGs as local analogues of the high-redshift, massive and quiescent galaxy population, thus providing additional constraints for models of galaxy formation and evolution.

CO observations allow estimations of the gas content of molecular clouds, which trace the reservoir of cold gas fuelling star formation, as well as to determine extinction via H$_2$ column density, N(H$_2$). Here, we studied millimetric and optical properties at 26 supernovae (SNe) locations of different types in a sample of 23 nearby galaxies by combining molecular $^{12}$C$^{16}$O (J = 1 $\rightarrow$ 0) resolved maps from the EDGE survey and optical Integral Field Spectroscopy from the CALIFA survey. We found an even clearer separation between type II and type Ibc SNe in terms of molecular gas than what we found in the optical using H$\alpha$ emission as a proxy for current SF rate, which reinforces the fact that SNe Ibc are more associated with SF-environments. While A$_V$ at SN locations is similar for SNe II and SNe Ibc, and higher compared to SNe Ia, N(H$_2$) is significantly higher for SNe Ibc than for SNe II and SNe Ia. When compared to alternative extinction estimations directly from SN photometry and spectroscopy, we find that our SNe Ibc have also redder color excess but showed standard Na I D absorption pseudo-equivalent widths ($\sim$1 \AA). In some cases we find no extinction when estimated from the environment, but high amounts of extinction when measured from SN observations, which suggests that circumstellar material or dust sublimation may be playing a role. This work serves as a benchmark for future studies combining last generation millimeter and optical IFS instruments to reveal the local environmental properties of extragalactic SNe.

The detection of intermediate mass black holes (IMBHs) in Galactic globular clusters (GCs) has so far been controversial. In order to characterize the effectiveness of integrated-light spectroscopy through integral field units, we analyze realistic mock data generated from state-of-the-art Monte Carlo simulations of GCs with a central IMBH, considering different setups and conditions varying IMBH mass, cluster distance, and accuracy in determination of the center. The mock observations are modeled with isotropic Jeans models to assess the success rate in identifying the IMBH presence, which we find to be primarily dependent on IMBH mass. However, even for a IMBH of considerable mass (3% of the total GC mass), the analysis does not yield conclusive results in 1 out of 5 cases, because of shot noise due to bright stars close to the IMBH line-of-sight. This stochastic variability in the modeling outcome grows with decreasing BH mass, with approximately 3 failures out of 4 for IMBHs with 0.1% of total GC mass. Finally, we find that our analysis is generally unable to exclude at 68% confidence an IMBH with mass of $10^3~M_\odot$ in snapshots without a central BH. Interestingly, our results are not sensitive to GC distance within 5-20 kpc, nor to mis-identification of the GC center by less than 2'' (<20% of the core radius). These findings highlight the value of ground-based integral field spectroscopy for large GC surveys, where systematic failures can be accounted for, but stress the importance of discrete kinematic measurements that are less affected by stochasticity induced by bright stars.

Galaxies in the local Universe are known to follow bimodal distributions in the global stellar populations properties. We analyze the distribution of the local average stellar-population ages of 654,053 sub-galactic regions resolved on ~1-kpc scales in a volume-corrected sample of 394 galaxies, drawn from the CALIFA-DR3 integral-field-spectroscopy survey and complemented by SDSS imaging. We find a bimodal local-age distribution, with an old and a young peak primarily due to regions in early-type galaxies and star-forming regions of spirals, respectively. Within spiral galaxies, the older ages of bulges and inter-arm regions relative to spiral arms support an internal age bimodality. Although regions of higher stellar-mass surface-density, mu*, are typically older, mu* alone does not determine the stellar population age and a bimodal distribution is found at any fixed mu*. We identify an "old ridge" of regions of age ~9 Gyr, independent of mu*, and a "young sequence" of regions with age increasing with mu* from 1-1.5 Gyr to 4-5 Gyr. We interpret the former as regions containing only old stars, and the latter as regions where the relative contamination of old stellar populations by young stars decreases as mu* increases. The reason why this bimodal age distribution is not inconsistent with the unimodal shape of the cosmic-averaged star-formation history is that i) the dominating contribution by young stars biases the age low with respect to the average epoch of star formation, and ii) the use of a single average age per region is unable to represent the full time-extent of the star-formation history of "young-sequence" regions.

We construct triaxial dynamical models for the Milky Way nuclear star cluster using Schwarzschild's orbit superposition technique. We fit the stellar kinematic maps presented in Feldmeier et al. (2014). The models are used to constrain the supermassive black hole mass M_BH, dynamical mass-to-light ratio M/L, and the intrinsic shape of the cluster. Our best-fitting model has M_BH = (3.0 +1.1 -1.3)x10^6 M_sun, M/L = (0.90 +0.76 -0.08) M_sun/L_{sun,4.5micron}, and a compression of the cluster along the line-of-sight. Our results are in agreement with the direct measurement of the supermassive black hole mass using the motion of stars on Keplerian orbits. The mass-to-light ratio is consistent with stellar population studies of other galaxies in the mid-infrared. It is possible that we underestimate M_BH and overestimate the cluster's triaxiality due to observational effects. The spatially semi-resolved kinematic data and extinction within the nuclear star cluster bias the observations to the near side of the cluster, and may appear as a compression of the nuclear star cluster along the line-of-sight. We derive a total dynamical mass for the Milky Way nuclear star cluster of M_MWNSC = (2.1 +-0.7)x10^7 M_sun within a sphere with radius r = 2 x r_eff = 8.4 pc. The best-fitting model is tangentially anisotropic in the central r = 0.5-2 pc of the nuclear star cluster, but close to isotropic at larger radii. Our triaxial models are able to recover complex kinematic substructures in the velocity map.

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.

We present a two-dimensional multi-component photometric decomposition of 404 galaxies from the CALIFA Data Release 3. They represent all possible galaxies with no clear signs of interaction and not strongly inclined in the final CALIFA data release. Galaxies are modelled in the g, r, and i SDSS images including, when appropriate, a nuclear point source, bulge, bar, and an exponential or broken disc component. We use a human-supervised approach to determine the optimal number of structures to be included in the fit. The dataset, including the photometric parameters of the CALIFA sample, is released together with statistical errors and a visual analysis of the quality of each fit. The analysis of the photometric components reveals a clear segregation of the structural composition of galaxies with stellar mass. At high masses (log(Mstar/Msun)>11), the galaxy population is dominated by galaxies modelled with a single Sersic or a bulge+disc with a bulge-to-total (B/T) luminosity ratio B/T>0.2. At intermediate masses (9.5<log(Mstar/Msun)<11), galaxies described with bulge+disc but B/T < 0.2 are preponderant, whereas, at the low mass end (log(Mstar/Msun)<9.5), the prevailing population is constituted by galaxies modelled with either pure discs or nuclear point sources+discs (i.e., no discernible bulge). We obtain that 57% of the volume corrected sample of disc galaxies in the CALIFA sample host a bar. This bar fraction shows a significant drop with increasing galaxy mass in the range 9.5<log(Mstar/Msun)<11.5. The analyses of the extended multi-component radial profile result in a volume-corrected distribution of 62%, 28%, and 10% for the so-called Type I, Type II, and Type III disc profiles, respectively. These fractions are in discordance with previous findings. We argue that the different methodologies used to detect the breaks are the main cause for these differences.

Nuclear Star Clusters (NSCs) are commonly observed in the centres of most galactic nuclei, including our own Milky Way. While their study can reveal important information about the build-up of the innermost regions of galaxies, the physical processes that regulate their formation are still poorly understood. NSCs might have been formed through gas infall and subsequent in situ star formation, and/or through the infall and merging of multiple star clusters into the centre of the galaxy. Here, we investigate the viability of the latter, by studying direct N-body simulations of inspiralling clusters to the centre of a Milky-Way-like nuclear bulge that hosts a massive black hole. We find that the NSC that forms through this process can show both morphological and kinematical properties that make it comparable with observations of the Milky Way NSC, including significant rotation- a fact that has been so far attributed mainly to gas infall. We explore its kinematic evolution, to see if and how the merger history can imprint fossil records on its dynamical structure. Moreover, we study the effect of stellar foreground contamination in the line-of-sight kinematics of the NSC. Our study shows that no fine tuning of the orientation of the infalling globular clusters is necessary to result in a rotating NSC. We suggest that cluster-inspiral is a viable mechanism for the formation of rotating NSCs.

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_{\phi}}/\sigma_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_{\phi}}/\sigma_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.

We present the stellar kinematic maps of a large sample of galaxies from the integral-field spectroscopic survey CALIFA. The sample comprises 300 galaxies displaying a wide range of morphologies across the Hubble sequence, from ellipticals to late-type spirals. This dataset allows us to homogeneously extract stellar kinematics up to several effective radii. In this paper, we describe the level of completeness of this subset of galaxies with respect to the full CALIFA sample, as well as the virtues and limitations of the kinematic extraction compared to other well-known integral-field surveys. In addition, we provide averaged integrated velocity dispersion radial profiles for different galaxy types, which are particularly useful to apply aperture corrections for single aperture measurements or poorly resolved stellar kinematics of high-redshift sources. The work presented in this paper sets the basis for the study of more general properties of galaxies that will be explored in subsequent papers of the survey.

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.

We present Fabry-P\'erot interferometric observations of the narrow H$\alpha$ component in the shock front of the historical supernova remnant Tycho (SN 1572). Using GH$\alpha$FaS (Galaxy H$\alpha$ Fabry-P\'erot 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$\alpha$-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$\alpha$ 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.

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 $\gamma = 0.5 \pm 0.3$. The metal-rich population is nearly isotropic (with $\beta_r^{red} = 0.0\pm0.1$) while the metal-poor population is tangentially anisotropic (with $\beta_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}/\sigma_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 ($\gamma = 0$) and cusp ($\gamma = 1$), and to constrain the possible internal rotation to better than $1\,\sigma$ 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.

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}/\sigma_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.

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$\mu$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 $\mu$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$\%$.

Deep observations of the dwarf elliptical (dE) galaxy NGC 1396 (M$_V = -16.60$, Mass $\sim 4\times10^8$ M$_\odot$), located in the Fornax cluster, have been performed with the VLT/ MUSE spectrograph in the wavelength region from $4750-9350$ \AA{}. In this paper we present a stellar population analysis studying chemical abundances, the star formation history (SFH) and the stellar initial mass function (IMF) as a function of galacto-centric distance. Different, independent ways to analyse the stellar populations result in a luminosity-weighted age of $\sim$ 6 Gyr and a metallicity [Fe/H]$\sim$ $-0.4$, similar to other dEs of similar mass. We find unusually overabundant values of [Ca/Fe] $\sim +0.1$, and under-abundant Sodium, with [Na/Fe] values around $-0.1$, while [Mg/Fe] is overabundant at all radii, increasing from $\sim+0.1$ in the centre to $\sim +0.2$ dex. We notice a significant metallicity and age gradient within this dwarf galaxy. To constrain the stellar IMF of NGC 1396, we find that the IMF of NGC 1396 is consistent with either a Kroupa-like or a top-heavy distribution, while a bottom-heavy IMF is firmly ruled out. An analysis of the abundance ratios, and a comparison with galaxies in the Local Group, shows that the chemical enrichment history of NGC 1396 is similar to the Galactic disc, with an extended star formation history. This would be the case if the galaxy originated from a LMC-sized dwarf galaxy progenitor, which would lose its gas while falling into the Fornax cluster.

The velocity function is a fundamental observable statistic of the galaxy population, similarly impor- tant as the luminosity function, but much more difficult to measure. In this work we present the first directly measured circular velocity function that is representative between 60 < v_circ < 320 km/s for galaxies of all morphological types at a given rotation velocity. For the low mass galaxy population (60 < v_circ < 170 km/s), we use the HIPASS velocity function. For the massive galaxy population (170 < v_circ < 320 km/s), we use stellar circular velocities from the Calar Alto Legacy Integral Field Area Survey (CALIFA). In earlier work we obtained the measurements of circular velocity at the 80% light radius for 226 galaxies and demonstrated that the CALIFA sample can produce volume- corrected galaxy distribution functions. The CALIFA velocity function includes homogeneous velocity measurements of both late and early-type rotation-supported galaxies and has the crucial advantage of not missing gas-poor massive ellipticals that HI surveys are blind to. We show that both velocity functions can be combined in a seamless manner, as their ranges of validity overlap. The resulting observed velocity function is compared to velocity functions derived from cosmological simulations of the z = 0 galaxy population. We find that dark matter-only simulations show a strong mismatch with the observed VF. Hydrodynamic simulations fare better, but still do not fully reproduce observations.

In this letter we describe how we use stellar dynamics information to constrain the shape of the stellar IMF in a sample of 27 early-type galaxies from the CALIFA survey. We obtain dynamical and stellar mass-to-light ratios, $\Upsilon_\mathrm{dyn}$ and $\Upsilon_{\ast}$, over a homogenous aperture of 0.5~$R_{e}$. We use the constraint $\Upsilon_\mathrm{dyn} \ge \Upsilon_{\ast}$ to test two IMF shapes within the framework of the extended MILES stellar population models. We rule out a single power law IMF shape for 75% of the galaxies in our sample. Conversely, we find that a double power law IMF shape with a varying high-mass end slope is compatible (within 1$\sigma$) with 95% of the galaxies. We also show that dynamical and stellar IMF mismatch factors give consistent results for the systematic variation of the IMF in these galaxies.

This paper describes the Third Public Data Release (DR3) of the Calar Alto Legacy Integral Field Area (CALIFA) survey. Science-grade quality data for 667 galaxies are made public, including the 200 galaxies of the Second Public Data Release (DR2). Data were obtained with the integral-field spectrograph PMAS/PPak mounted on the 3.5m telescope at the Calar Alto Observatory. Three different spectral setups are available, i) a low-resolution V500 setup covering the wavelength range 3749-7500 AA (4240-7140 AA unvignetted) with a spectral resolution of 6.0 AA (FWHM), for 646 galaxies, ii) a medium-resolution V1200 setup covering the wavelength range 3650-4840 AA (3650-4620 AA unvignetted) with a spectral resolution of 2.3 AA (FWHM), for 484 galaxies, and iii) the combination of the cubes from both setups (called COMBO), with a spectral resolution of 6.0 AA and a wavelength range between 3700-7500 AA (3700-7140 AA unvignetted), for 446 galaxies. The Main Sample, selected and observed according to the CALIFA survey strategy covers a redshift range between 0.005 and 0.03, spans the color-magnitude diagram and probes a wide range of stellar mass, ionization conditions, and morphological types. The Extension Sample covers several types of galaxies that are rare in the overall galaxy population and therefore not numerous or absent in the CALIFA Main Sample. All the cubes in the data release were processed using the latest pipeline, which includes improved versions of the calibration frames and an even further improved im- age reconstruction quality. In total, the third data release contains 1576 datacubes, including ~1.5 million independent spectra. It is available at http://califa.caha.es/DR3.

The space density distribution of galaxies in the absolute magnitude - rotation velocity plane: a volume-complete Tully-Fisher relation from CALIFA stellar kinematics

We observed the lenticular galaxy NGC 3998 with the Mitchell Integral-Field Spectrograph and extracted line-of-sight velocity distributions out to 3 half-light radii. We constructed collisionless orbit models in order to constrain NGC 3998's dark and visible structure, using kinematics from both the Mitchell and SAURON instruments. We find NGC 3998 to be almost axisymmetric, seen nearly face on with a flattened intrinsic shape - i.e., a face-on fast-rotator. We find an I-band mass-to-light ratio of $4.7_{-0.45}^{+0.32}$ in good agreement with previous spectral fitting results for this galaxy. Our best-fit orbit model shows a both a bulge and a disc component, with a non-negligible counter-rotating component also evident. We find that relatively little dark matter is needed to model this galaxy, with an inferred dark mass fraction of just $(7.1^{+8.1}_{-7.1})\%$ within one half-light radius.

The search of extragalactic regions with conspicuous presence of Wolf-Rayet (WR) stars outside the Local Group is challenging task due to the difficulties in detecting their faint spectral features. In this exploratory work, we develop a methodology to perform an automated search of WR signatures through a pixel-by-pixel analysis of integral field spectroscopy (IFS) data belonging to the Calar Alto Legacy Integral Field Area survey, CALIFA. This technique allowed us to build the first catalogue of Wolf-Rayet rich regions with spatially-resolved information, allowing to study the properties of these complexes in a 2D context. The detection technique is based on the identification of the blue WR bump (around He II 4686 {\AA}, mainly associated to nitrogen-rich WR stars, WN) and the red WR bump (around C IV 5808 {\AA} and associated to carbon-rich WR stars, WC) using a pixel-by-pixel analysis. We identified 44 WR-rich regions with blue bumps distributed in 25 galaxies of a total of 558. The red WR bump was identified only in 5 of those regions. We found that the majority of the galaxies hosting WR populations in our sample are involved in some kind of interaction process. Half of the host galaxies share some properties with gamma-ray burst (GRB) hosts where WR stars, as potential candidates to being the progenitors of GRBs, are found. We also compared the WR properties derived from the CALIFA data with stellar population synthesis models, and confirm that simple star models are generally not able to reproduce the observations. We conclude that other effects, such as the binary star channel (which could extend the WR phase up to 10 Myr), fast rotation or other physical processes that causes the loss of observed Lyman continuum photons, are very likely affecting the derived WR properties, and hence should be considered when modelling the evolution of massive stars.

Lenticular galaxies (S0s) represent the majority of early-type galaxies in the local Universe, but their formation channels are still poorly understood. While galaxy mergers are obvious pathways to suppress star formation and increase bulge sizes, the marked parallelism between spiral and lenticular galaxies (e.g. photometric bulge-disc coupling) seemed to rule out a potential merger origin. Here, we summarise our recent work in which we have shown, through N-body numerical simulations, that disc-dominated lenticulars can emerge from major mergers of spiral galaxies, in good agreement with observational photometric scaling relations. Moreover, we show that mergers simultaneously increase the light concentration and reduce the angular momentum relative to their spiral progenitors. This explains the mismatch in angular momentum and concentration between spirals and lenticulars recently revealed by CALIFA observations, which is hard to reconcile with simple fading mechanisms (e.g. ram-pressure stripping).

The metallicity of a supernova (SN) progenitor, together with its mass, is one of the main parameters that rules their outcome. We present a metallicity study of 115 nearby SN host galaxies (0.005<z<0.03) which hosted 142 SNe using Integral Field Spectroscopy (IFS) from the CALIFA survey. Using O3N2 we found no statistically significant differences between the gas-phase metallicities at the locations of the three main SN types (Ia, Ib/c and II) all having ~8.50$\pm$0.02 dex. The total galaxy metallicities are also very similar and we argue that this is because our sample consists only of SNe discovered in massive galaxies (log(M/Msun)>10 dex) by targeted searches. We also found no evidence that the metallicity at the SN location differs from the average metallicity at the GCD of the SNe. By extending our SN sample with published metallicities at the SN location, we studied the metallicity distributions for all SN subtypes split into SN discovered in targeted and untargeted searches. We confirm a bias toward higher host masses and metallicities in the targeted searches. Combining data from targeted and untargeted searches we found a sequence from higher to lower local metallicity: SN Ia, Ic, and II show the highest metallicity, which is significantly higher than SN Ib, IIb, and Ic-BL. Our results support the picture of SN Ib resulting from binary progenitors and, at least part of, SN Ic being the result of single massive stars stripped of their outer layers by metallicity driven winds. We studied several proxies of the local metallicity frequently used in the literature and found that the total host metallicity allows for the estimation of the metallicity at the SN location with an accuracy better than 0.08 dex and very small bias. In addition, weak AGNs not seen in total spectra may only weakly bias (by 0.04 dex) the metallicity estimate from integrated spectra. (abridged)

High-precision kinematic studies of globular clusters require an accurate knowledge of all possible sources of contamination. Amongst other sources, binary stars can introduce systematic biases in the kinematics. Using a set of Monte Carlo cluster simulations with different concentrations and binary fractions, we investigate the effect of unresolved binaries on proper-motion dispersion profiles, treating the simulations like HST proper-motion samples. Since globular clusters evolve towards a state of partial energy equipartition, more massive stars lose energy and decrease their velocity dispersion. As a consequence, on average, binaries have a lower velocity dispersion, since they are more massive kinematic tracers. We show that, in the case of clusters with high binary fraction (initial binary fraction of 50%) and high concentration (i.e., closer to energy equipartition), unresolved binaries introduce a color-dependent bias in the velocity dispersion of main-sequence stars of the order of 0.1-0.3 km s$^{-1}$ (corresponding to $1-6$% of the velocity dispersion), with the reddest stars having a lower velocity dispersion, due to the higher fraction of contaminating binaries. This bias depends on the ability to distinguish binaries from single stars, on the details of the color-magnitude diagram and the photometric errors. We apply our analysis to the HSTPROMO data set of NGC 7078 (M15) and show that no effect ascribable to binaries is observed, consistent with the low binary fraction of the cluster. Our work indicates that binaries do not significantly bias proper-motion velocity-dispersion profiles, but should be taken into account in the error budget of kinematic analyses.

Two-body interactions play a major role in shaping the structural and dynamical properties of globular clusters (GCs) over their long-term evolution. In particular, GCs evolve toward a state of partial energy equipartition that induces a mass-dependence in their kinematics. By using a set of Monte Carlo cluster simulations evolved in quasi-isolation, we show that the stellar mass dependence of the velocity dispersion $\sigma(m)$ can be described by an exponential function $\sigma^2\propto \exp(-m/m_\mathrm{eq})$, with the parameter $m_\mathrm{eq}$ quantifying the degree of partial energy equipartition of the systems. This simple parametrization successfully captures the behaviour of the velocity dispersion at lower as well as higher stellar masses, that is, the regime where the system is expected to approach full equipartition. We find a tight correlation between the degree of equipartition reached by a GC and its dynamical state, indicating that clusters that are more than about 20 core relaxation times old, have reached a maximum degree of equipartition. This equipartition$-$dynamical state relation can be used as a tool to characterize the relaxation condition of a cluster with a kinematic measure of the $m_\mathrm{eq}$ parameter. Vice versa, the mass-dependence of the kinematics can be predicted knowing the relaxation time solely on the basis of photometric measurements. Moreover, any deviations from this tight relation could be used as a probe of a peculiar dynamical history of a cluster. Finally, our novel approach is important for the interpretation of state-of-the-art Hubble Space Telescope proper motion data, for which the mass dependence of kinematics can now be measured, and for the application of modeling techniques which take into consideration multi-mass components and mass segregation.

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.

We report the discovery of a complex extended density enhancement in the Globular Clusters (GCs) in the central $\sim 0.5(^{\circ})^2$ ($\sim 0.06$ Mpc$^2$) of the Fornax cluster, corresponding to $\sim 50\%$ of the area within 1 core radius. This overdensity connects the GC system of NGC1399 to most of those of neighboring galaxies within $\sim 0.6^{\circ}$ ($\sim 210$ kpc) along the W-E direction. The asymmetric density structure suggests that the galaxies in the core of the Fornax cluster experienced a lively history of interactions that have left a clear imprint on the spatial distribution of GCs. The extended central dominant structure is more prominent in the distribution of blue GCs, while red GCs show density enhancements that are more centrally concentrated on the host galaxies. We propose that the relatively small-scale density structures in the red GCs are caused by galaxy-galaxy interactions, while the extensive spatial distribution of blue GCs is due to stripping of GCs from the halos of core massive galaxies by the Fornax gravitational potential. Our investigations is based on density maps of candidate GCs extracted from the multi-band VLT Survey Telescope (VST) survey of Fornax (FDS), identified in a three-dimensional color space and further selected based on their $g$-band magnitude and morphology.

We study the stellar properties of 44 face-on spiral galaxies from the Calar Alto Legacy Integral Field Area survey via full spectrum fitting techniques. We compare the age profiles with the surface brightness distribution in order to highlight differences between profile types (type I, exponential profile; and II, down-bending profile). We observe an upturn ("U-shape") in the age profiles for 17 out of these 44 galaxies with reliable stellar information up to their outer parts. This "U-shape" is not a unique feature for type II galaxies but can be observed in type I as well. These findings suggest that the mechanisms shaping the surface brightness and stellar population distributions are not directly coupled. This upturn in age is only observable in the light-weighted profiles while it flattens out in the mass-weighted profiles. Given recent results on the outer parts of nearby systems and the results presented in this Letter, one of the most plausible explanations for the age upturn is an early formation of the entire disc ($\sim$~10~Gyr ago) followed by an inside-out quenching of the star formation.

We investigate the compact, early-type galaxy NGC 1281 with integral field unit observations to map the stellar LOSVD out to 5 effective radii and construct orbit-based dynamical models to constrain its dark and luminous matter content. Under the assumption of mass-follows-light, the H-band stellar mass-to-light ratio (M/L) is {\Upsilon} = 2.7(+-0.1) {\Upsilon}_{sun}, higher than expected from our stellar population synthesis fits with either a canonical Kroupa ({\Upsilon} = 1.3 {\Upsilon}_{sun}) or Salpeter ({\Upsilon} = 1.7 {\Upsilon}_{sun}) stellar initial mass function. Such models also cannot reproduce the details of the LOSVD. Models with a dark halo recover the kinematics well and indicate that NGC 1281 is dark matter dominated, making up ~ 90 per cent of the total enclosed mass within the kinematic bounds. Parameterised as a spherical NFW profile, the dark halo mass is 11.5 < log(M_{DM}/M_{sun}) < 11.8 and the stellar M/L is 0.6 < {\Upsilon} < 1.1. However, this stellar M/L is lower than predicted by its old stellar population. Moreover, the halo mass within the kinematic extent is ten times larger than expected based on {\Lambda}CDM predictions, and an extrapolation yields cluster sized dark halo masses. Adopting {\Upsilon} = 1.7 {\Upsilon}_{sun} yields more moderate dark halo virial masses, but these models fit the kinematics worse. A non-NFW model might solve the discrepancy between the unphysical consequences of the best-fitting dynamical models and models based on more reasonable assumptions for the dark halo and stellar mass-to-light ratio, which are disfavoured according to our parameter estimation.

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 $\lambda_\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.

The mid-infrared is an optimal window to trace stellar mass in nearby galaxies and the 3.6$\mu m$ IRAC band has been exploited to this effect, but such mass estimates can be biased by dust emission. We present our pipeline to reveal the old stellar flux at 3.6$\mu m$ and obtain stellar mass maps for more than 1600 galaxies available from the Spitzer Survey of Stellar Structure in Galaxies (S$^{4}$G). This survey consists of images in two infrared bands (3.6 and 4.5$\mu m$), and we use the Independent Component Analysis (ICA) method presented in Meidt et al. (2012) to separate the dominant light from old stars and the dust emission that can significantly contribute to the observed 3.6$\mu m$ flux. We exclude from our ICA analysis galaxies with low signal-to-noise ratio (S/N < 10) and those with original [3.6]-[4.5] colors compatible with an old stellar population, indicative of little dust emission (mostly early Hubble types, which can directly provide good mass maps). For the remaining 1251 galaxies to which ICA was successfully applied, we find that as much as 10-30% of the total light at 3.6$\mu m$ typically originates from dust, and locally it can reach even higher values. This contamination fraction shows a correlation with specific star formation rates, confirming that the dust emission that we detect is related to star formation. Additionally, we have used our large sample of mass estimates to calibrate a relationship of effective mass-to-light ratio ($M/L$) as a function of observed [3.6]-[4.5] color: $\log(M/L)=-0.339 (\pm 0.057) \times ([3.6]-[4.5]) -0.336 (\pm 0.002)$. Our final pipeline products have been made public through IRSA, providing the astronomical community with an unprecedentedly large set of stellar mass maps ready to use for scientific applications.

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.

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\'ersic 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.

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.

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.

We present spatially resolved stellar and/or ionized gas kinematic properties for a sample of 103 interacting galaxies, tracing all merger stages: close companions, pairs with morphological signatures of interaction, and coalesced merger remnants. We compare our sample with 80 non-interacting galaxies. We measure for the stellar and the ionized gas components the major (projected) kinematic position angles (PA$_{\mathrm{kin}}$, approaching and receding) directly from the velocity fields with no assumptions on the internal motions. This method allow us to derive the deviations of the kinematic PAs from a straight line ($\delta$PA$_{\mathrm{kin}}$). Around half of the interacting objects show morpho-kinematic PA misalignments that cannot be found in the control sample. Those misalignments are present mostly in galaxies with morphological signatures of interaction. Alignment between the kinematic sides for both samples is similar, with most of the galaxies displaying small misalignments. Radial deviations of the kinematic PA from a straight line in the stellar component measured by $\delta$PA$_{\mathrm{kin}}$ are large for both samples. However, for a large fraction of interacting galaxies the ionized gas $\delta$PA$_{\mathrm{kin}}$ is larger than typical values derived from isolated galaxies (48%), making this parameter a good indicator to trace the impact of interaction and mergers in the internal motions of galaxies. By comparing the stellar and ionized gas kinematic PA, we find that 42% (28/66) of the interacting galaxies have misalignments larger than 16 degrees, compared to 10% from the control sample. Our results show the impact of interactions in the internal structure of galaxies as well as the wide variety of their velocity distributions. This study also provides a local Universe benchmark for kinematic studies in merging galaxies at high redshift.

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.

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 $\lambda_\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.

We use optical integral-field spectroscopic (IFS) data from 103 nearby galaxies at different stages of the merging event, from close pairs to merger remnants provided by the CALIFA survey, to study the impact of the interaction in the specific star formation and oxygen abundance on different galactic scales. To disentangle the effect of the interaction and merger from internal processes, we compared our results with a control sample of 80 non-interacting galaxies. We confirm the moderate enhancement (2-3 times) of specific star formation for interacting galaxies in central regions as reported by previous studies; however, the specific star formation is comparable when observed in extended regions. We find that control and interacting star-forming galaxies have similar oxygen abundances in their central regions, when normalized to their stellar masses. Oxygen abundances of these interacting galaxies seem to decrease compared to the control objects at the large aperture sizes measured in effective radius. Although the enhancement in central star formation and lower metallicities for interacting galaxies have been attributed to tidally induced inflows, our results suggest that other processes such as stellar feedback can contribute to the metal enrichment in interacting galaxies.

While studies of gas-phase metallicity gradients in disc galaxies are common, very little has been done in the acquisition of stellar abundance gradients in the same regions. We present here a comparative study of the stellar metallicity and age distributions in a sample of 62 nearly face-on, spiral galaxies with and without bars, using data from the CALIFA survey. We measure the slopes of the gradients and study their relation with other properties of the galaxies. We find that the mean stellar age and metallicity gradients in the disc are shallow and negative. Furthermore, when normalized to the effective radius of the disc, the slope of the stellar population gradients does not correlate with the mass or with the morphological type of the galaxies. Contrary to this, the values of both age and metallicity at $\sim$2.5 scale-lengths correlate with the central velocity dispersion in a similar manner to the central values of the bulges, although bulges show, on average, older ages and higher metallicities than the discs. One of the goals of the present paper is to test the theoretical prediction that non-linear coupling between the bar and the spiral arms is an efficient mechanism for producing radial migrations across significant distances within discs. The process of radial migration should flatten the stellar metallicity gradient with time and, therefore, we would expect flatter stellar metallicity gradients in barred galaxies. However, we do not find any difference in the metallicity or age gradients in galaxies with without bars. We discuss possible scenarios that can lead to this absence of difference.

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.

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\alpha\ 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.

The bar pattern speed ($\Omega_{\rm b}$) is defined as the rotational frequency of the bar, and it determines the bar dynamics. Several methods have been proposed for measuring $\Omega_{\rm b}$. The non-parametric method proposed by Tremaine \& Weinberg (1984; TW) and based on stellar kinematics is the most accurate. This method has been applied so far to 17 galaxies, most of them SB0 and SBa types. We have applied the TW method to a new sample of 15 strong and bright barred galaxies, spanning a wide range of morphological types from SB0 to SBbc. Combining our analysis with previous studies, we investigate 32 barred galaxies with their pattern speed measured by the TW method. The resulting total sample of barred galaxies allows us to study the dependence of $\Omega_{\rm b}$ on galaxy properties, such as the Hubble type. We measured $\Omega_{\rm b}$ using the TW method on the stellar velocity maps provided by the integral-field spectroscopy data from the CALIFA survey. Integral-field data solve the problems that long-slit data present when applying the TW method, resulting in the determination of more accurate $\Omega_{\rm b}$. In addition, we have also derived the ratio $\cal{R}$ of the corotation radius to the bar length of the galaxies. According to this parameter, bars can be classified as fast ($\cal{R}$ $< 1.4$) and slow ($\cal{R}$>1.4). For all the galaxies, $\cal{R}$ is compatible within the errors with fast bars. We cannot rule out (at 95$\%$ level) the fast bar solution for any galaxy. We have not observed any significant trend between $\cal{R}$ and the galaxy morphological type. Our results indicate that independent of the Hubble type, bars have been formed and then evolve as fast rotators. This observational result will constrain the scenarios of formation and evolution of bars proposed by numerical simulations.

The detection of intermediate mass black holes (IMBHs) in globular clusters has been hotly debated, with different observational methods delivering different outcomes for the same object. In order to understand these discrepancies, we construct detailed mock integral field spectroscopy (IFU) observations of globular clusters, starting from realistic Monte Carlo cluster simulations. The output is a data cube of spectra in a given field-of-view that can be analyzed in the same manner as real observations and compared to other (resolved) kinematic measurement methods. We show that the main discrepancies arise because the luminosity-weighted IFU observations can be strongly biased by the presence of a few bright stars that introduce a scatter in velocity dispersion measurements of several km/s. We show that this intrinsic scatter can prevent a sound assessment of the central kinematics, and therefore should be fully taken into account to correctly interpret the signature of an IMBH.

We present spatially resolved kinematics and global stellar populations and mass-to-light ratios for a sample of 39 dwarf early-type (dE) galaxies in the Virgo cluster studied as part of the SMAKCED stellar absorption-line spectroscopy and imaging survey. This sample is representative of the early-type population in the absolute magnitude range -19.0 < M_r < -16.0. For each dE, we measure the rotation curve and velocity dispersion profile and fit an analytic function to the rotation curve. We study the significance of the departure of the rotation curve from the best fit analytic function (poorly fit) and of the difference between the approaching and receding sides of the rotation curve (asymmetry). We find that 62 +/- 8 % (23 out of the 39) of the dEs have a significant anomaly in their rotation curve. Analysis of the images reveals photometric anomalies for most galaxies. However, there is no clear correlation between the significance of the photometric and kinematic anomalies. We measure age-sensitive and metallicity-sensitive Lick spectral indices and find a wide range of ages and metallicities. We also find that 4 dEs have emission partially filling in the Balmer absorption lines. Finally, we estimate the total masses and dark matter fractions of the dEs. They have a median total mass and dark matter fraction within the Re of log Me = 9.1 +/- 0.2 and f_DM = 46 +/- 18 %. We plot several scaling relations and show that dEs seem to be the bridge between massive early-type and dwarf spheroidal galaxies.

We analyze the stellar kinematics of 39 dwarf early-type galaxies (dEs) in the Virgo cluster. Based on the specific stellar angular momentum lambda_e and the ellipticity, we find 11 slow rotators and 28 fast rotators. The fast rotators in the outer parts of the Virgo cluster rotate significantly faster than fast rotators in the inner parts of the cluster. Moreover, 10 out of the 11 slow rotators are located in the inner 3 degrees (D < 1 Mpc) of the cluster. The fast rotators contain subtle disky structures that are visible in high-pass filtered optical images, while the slow rotators do not exhibit these structures. In addition, two of the dEs have kinematically decoupled cores and four more have emission partially filling in the Balmer absorption lines. These properties suggest that Virgo cluster dEs may have originated from late-type star-forming galaxies that were transformed by the environment after their infall into the cluster. The correlation between lambda_e and the clustercentric distance can be explained by a scenario where low luminosity star-forming galaxies fall into the cluster, their gas is rapidly removed by ram pressure stripping, although some of it can be retained in their core, their star-formation is quenched but their stellar kinematics are preserved. After a long time in the cluster and several passes through its center, the galaxies are heated up and transformed into slow rotating dEs.

H ii regions in galaxies are the sites of star formation and thus particular places to understand the build-up of stellar mass in the universe. The line ratios of this ionized gas are frequently used to characterize the ionization conditions. We use the Hii regions catalogue from the CALIFA survey (~5000 H ii regions), to explore their distribution across the classical [OIII]/Hbeta vs. [NII]/Halpha diagnostic diagram, and how it depends on the oxygen abundance, ionization parameter, electron density, and dust attenuation. We compared the line ratios with predictions from photoionization models. Finally, we explore the dependences on the properties of the host galaxies, the location within those galaxies and the properties of the underlying stellar population. We found that the location within the BPT diagrams is not totally predicted by photoionization models. Indeed, it depends on the properties of the host galaxies, their galactocentric distances and the properties of the underlying stellar population. These results indicate that although H ii regions are short lived events, they are affected by the total underlying stellar population. One may say that H ii regions keep a memory of the stellar evolution and chemical enrichment that have left an imprint on the both the ionizing stellar population and the ionized gas

This paper describes the Second Public Data Release (DR2) of the Calar Alto Legacy Integral Field Area (CALIFA) survey. The data for 200 objects are made public, including the 100 galaxies of the First Public Data Release (DR1). Data were obtained with the integral-field spectrograph PMAS/PPak mounted on the 3.5 m telescope at the Calar Alto observatory. Two different spectral setups are available for each galaxy, (i) a low-resolution V500 setup covering the wavelength range 3745-7500 \AA\ with a spectral resolution of 6.0 \AA\ (FWHM), and (ii) a medium-resolution V1200 setup covering the wavelength range 3650-4840 \AA\ with a spectral resolution of 2.3 \AA\ (FWHM). The sample covers a redshift range between 0.005 and 0.03, with a wide range of properties in the Color-Magnitude diagram, stellar mass, ionization conditions, and morphological types. All released cubes were reduced with the latest pipeline, including improved spectrophotometric calibration, spatial registration and spatial resolution. The spectrophotometric calibration is better than 6% and the median spatial resolution is 2.5". Altogether the second data release contains over 1.5 million spectra. It is available at http://califa.caha.es/DR2.

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 $\lambda_{\rm Re}-\epsilon_{\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 $\lambda_{\rm Re}$ values than any elliptical and lenticular galaxy; (b) low-mass spiral galaxies with observed $\lambda_{\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.

[Abridged] We use optical IFS of nearby SN host galaxies provided by the CALIFA Survey with the goal of finding correlations in the environmental parameters at the location of different SN types. We recover the sequence in association of different SN types to the star-forming regions by using several indicators of the ongoing and recent SF related to both the ionized gas and the stellar populations. While the total ongoing SF is on average the same for the three SN types, SNe Ibc/IIb tend to happen closer to star-forming regions and occur in higher SF density locations compared to SNe II and SNe~Ia, the latter showing the weakest correlation. SNe~Ia host galaxies have on average masses that are $\sim$0.3-0.8~dex higher than CC SNe hosts due to a larger fraction of old stellar populations in the SNe~Ia hosts. Using the recent SN~Ia delay-time distribution and the SFHs of the galaxies, we show that the SN~Ia hosts in our sample should presently produce a factor 2 more SNe~Ia than the CC~SN hosts. Since both types are in hosts with similar SF rate and hence similar CC~SN rate, this can explain the mass difference between the SN~Ia and CC~SN hosts, and reinforce the finding that at least part of SNe~Ia should originate from very old progenitors. Comparing the mean SFH of the eight least massive galaxies to that of the massive SF SN~Ia hosts we find that the low-mass galaxies formed their stars over more extended time than the massive SN~Ia hosts. We estimate that the low-mass galaxies should produce by a factor of 10 less SNe~Ia, and a factor of 3 less CC~SNe than the high-mass group. Therefore the ratio between the number of CC~SNe and SNe~Ia is expected to increase with decreasing the galaxy mass. CC~SNe tend to explode at positions with younger stellar populations than the galaxy average, but the galaxy properties at SNe~Ia locations are one average the same as the global ones.

This work provides an overall characterization of the kinematic behavior of the ionized gas of the galaxies included in the Calar Alto Legacy Integral field Area (CALIFA), offering kinematic clues to potential users of this survey for including kinematical criteria for specific studies. From the first 200 galaxies observed by CALIFA, we present the 2D kinematic view of the 177 galaxies satisfying a gas detection threshold. After removing the stellar contribution, we used the cross-correlation technique to obtain the radial velocity of the dominant gaseous component. The main kinematic parameters were directly derived from the radial velocities with no assumptions on the internal motions. Evidence of the presence of several gaseous components with different kinematics were detected by using [OIII] profiles. Most objects in the sample show regular velocity fields, although the ionized-gas kinematics are rarely consistent with simple coplanar circular motions. 35% of the objects present evidence of a displacement between the photometric and kinematic centers larger than the original spaxel radii. Only 17% of the objects in the sample exhibit kinematic lopsidedness when comparing receding and approaching sides of the velocity fields, but most of them are interacting galaxies exhibiting nuclear activity. Early-type galaxies in the sample present clear photometric-kinematic misaligments. There is evidence of asymmetries in the emission line profiles suggesting the presence of kinematically distinct gaseous components at different distances from the nucleus. This work constitutes the first determination of the ionized gas kinematics of the galaxies observed in the CALIFA survey. The derived velocity fields, the reported kinematic peculiarities and the identification of the presence of several gaseous components might be used as additional criteria for selecting galaxies for specific studies.

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.

We describe and discuss the selection procedure and statistical properties of the galaxy sample used by the Calar Alto Legacy Integral Field Area Survey (CALIFA), a public legacy survey of 600 galaxies using integral field spectroscopy. The CALIFA "mother sample" was selected from the Sloan Digital Sky Survey (SDSS) DR7 photometric catalogue to include all galaxies with an r-band isophotal major axis between 45" and 79.2" and with a redshift 0.005 < z < 0.03. The mother sample contains 939 objects, 600 of which will be observed in the course of the CALIFA survey. The selection of targets for observations is based solely on visibility and thus keeps the statistical properties of the mother sample. By comparison with a large set of SDSS galaxies, we find that the CALIFA sample is representative of galaxies over a luminosity range of -19 > Mr > -23.1 and over a stellar mass range between 10^9.7 and 10^11.4Msun. In particular, within these ranges, the diameter selection does not lead to any significant bias against - or in favour of - intrinsically large or small galaxies. Only below luminosities of Mr = -19 (or stellar masses < 10^9.7Msun) is there a prevalence of galaxies with larger isophotal sizes, especially of nearly edge-on late-type galaxies, but such galaxies form < 10% of the full sample. We estimate volume-corrected distribution functions in luminosities and sizes and show that these are statistically fully compatible with estimates from the full SDSS when accounting for large-scale structure. We also present a number of value-added quantities determined for the galaxies in the CALIFA sample. We explore different ways of characterizing the environments of CALIFA galaxies, finding that the sample covers environmental conditions from the field to genuine clusters. We finally consider the expected incidence of active galactic nuclei among CALIFA galaxies.

We use spatially and temporally resolved maps of stellar population properties of 300 galaxies from the CALIFA integral field survey to investigate how the stellar metallicity (Z*) relates to the total stellar mass (M*) and the local mass surface density ($\mu$*) in both spheroidal and disk dominated galaxies. The galaxies are shown to follow a clear stellar mass-metallicity relation (MZR) over the whole 10$^9$ to 10$^{12}$ M$_{\odot}$ range. This relation is steeper than the one derived from nebular abundances, which is similar to the flatter stellar MZR derived when we consider only young stars. We also find a strong relation between the local values of $\mu$* and Z* (the $\mu$ZR), betraying the influence of local factors in determining Z*. This shows that both local ($\mu$*-driven) and global (M*-driven) processes are important in determining the metallicity in galaxies. We find that the overall balance between local and global effects varies with the location within a galaxy. In disks, $\mu$* regulates Z*, producing a strong $\mu$ZR whose amplitude is modulated by M*. In spheroids it is M* who dominates the physics of star formation and chemical enrichment, with $\mu$* playing a minor, secondary role. These findings agree with our previous analysis of the star formation histories of CALIFA galaxies, which showed that mean stellar ages are mainly governed by surface density in galaxy disks and by total mass in spheroids.

Aims: We present the first in-depth study of the kinematic properties and derive the structural parameters of NGC 4372 based on the fit of a Plummer profile and a rotating, physical model. We explore the link between internal rotation to different cluster properties and together with similar studies of more GCs, we put these in the context of globular cluster formation and evolution. Methods: We present radial velocities for 131 cluster member stars measured from high-resolution FLAMES/GIRAFFE observations. Their membership to the GC is additionally confirmed from precise metallicity estimates. Using this kinematic data set we build a velocity dispersion profile and a systemic rotation curve. Additionally, we obtain an elliptical number density profile of NGC 4372 based on optical images using a MCMC fitting algorithm. From this we derive the cluster's half-light radius and ellipticity as r_h=3.4'+/-0.04' and e=0.08+/-0.01. Finally, we give a physical interpretation of the observed morphological and kinematic properties of this GC by fitting an axisymmetric, differentially rotating, dynamical model. Results: Our results show that NGC 4372 has an unusually high ratio of rotation amplitude to velocity dispersion (1.2 vs. 4.5 km/s) for its metallicity. This, however, puts it in line with two other exceptional, very metal-poor GCs - M 15 and NGC 4590. We also find a mild flattening of NGC 4372 in the direction of its rotation. Given its old age, this suggests that the flattening is indeed caused by the systemic rotation rather than tidal interactions with the Galaxy. Additionally, we estimate the dynamical mass of the GC M_dyn=2.0+/-0.5 x 10^5 M_Sun based on the dynamical model, which constrains the mass-to-light ratio of NGC 4372 between 1.4 and 2.3 M_Sun/L_Sun, representative of an old, purely stellar population.

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.

We present 80 stellar and ionised gas velocity maps from the Calar Alto Legacy Integral Field Area (CALIFA) survey in order to characterize the kinematic orientation of non-interacting galaxies. The study of galaxies in isolation is a key step towards understanding how fast-external processes, such as major mergers, affect kinematic properties in galaxies. We derived the global and individual (projected approaching and receding sides) kinematic position angles (PAs) for both the stellar and ionised gas line-of-sight velocity distributions. When compared to the photometric PA, we find that morpho-kinematic differences are smaller than 22 degrees in 90% of the sample for both components; internal kinematic misalignments are generally smaller than 16 degrees. We find a tight relation between the global stellar and ionised gas kinematic PA consistent with circular-flow pattern motions in both components. This relation also holds generally in barred galaxies across the bar and galaxy disk scales. Our findings suggest that even in the presence of strong bars, both the stellar and the gaseous components tend to follow the gravitational potential of the disk. As a result, kinematic orientation can be used to assess the degree of external distortions in interacting galaxies.

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)

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.

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.

Driven by gravity, massive structures like galaxies and clusters of galaxies are believed to grow continuously through hierarchical merging and accretion of smaller systems. Observational evidence of accretion events is provided by the coherent stellar streams crossing the outer haloes of massive galaxies, such as the Milky Way or Andromeda. At similar mass-scales, around $10^{11}$ solar masses in stars, further observational evidence of merging activity is also ample. Mergers of lower-mass galaxies are expected within the hierarchical process of galaxy formation, but have hitherto not been seen for galaxies with less than about $10^9$ solar masses in stars. Here, we report the kinematic detection of a stellar stream in one of the satellite galaxies of Andromeda, the dwarf spheroidal galaxy Andromeda II, which has a mass of only $10^7$ solar masses in stars. The properties of the stream show that we are observing the remnant of a merger between two dwarf galaxies. This had a dramatic influence on the dynamics of the remnant, which is now rotating around its projected major axis. The stellar stream in Andromeda II illustrates the scale-free character of the formation of galaxies, down to the lowest galactic mass scales.

We present evidence for kinematically decoupled cores (KDCs) in two dwarf early-type (dE) galaxies in the Virgo cluster, VCC 1183 and VCC 1453, studied as part of the SMAKCED stellar absorption-line spectroscopy and imaging survey. These KDCs have radii of 1.8'' (0.14 kpc) and 4.2'' (0.33 kpc), respectively. Each of these KDCs is distinct from the main body of its host galaxy in two ways: (1) inverted sense of rotation; and (2) younger (and possibly more metal-rich) stellar population. The observed stellar population differences are probably associated with the KDC, although we cannot rule out the possibility of intrinsic radial gradients in the host galaxy. We describe a statistical analysis method to detect, quantify the significance of, and characterize KDCs in long-slit rotation curve data. We apply this method to the two dE galaxies presented in this paper and to five other dEs for which KDCs have been reported in the literature. Among these seven dEs, there are four significant KDC detections, two marginal KDC detections, and one dE with an unusual central kinematic anomaly that may be an asymmetric KDC.The frequency of occurence of KDCs and their properties provide important constraints on the formation history of their host galaxies. We discuss different formation scenarios for these KDCs in cluster environments and find that dwarf-dwarf wet mergers or gas accretion can explain the properties of these KDCs. Both of these mechanisms require that the progenitor had a close companion with a low relative velocity. This suggests that KDCs were formed in galaxy pairs residing in a poor group environment or in isolation whose subsequent infall into the cluster quenched star formation.

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.

We present the largest and most homogeneous catalog of HII regions and associations compiled so far. The catalog comprises more than 7000 ionized regions, extracted from 306 galaxies observed by the CALIFA survey. We describe the procedures used to detect, select, and analyse the spectroscopic properties of these ionized regions. In the current study we focus on the characterization of the radial gradient of the oxygen abundance in the ionized gas, based on the study of the deprojected distribution of HII regions. We found that all galaxies without clear evidence of an interaction present a common gradient in the oxygen abundance, with a characteristic slope of alpha = -0.1 dex/re between 0.3 and 2 disk effective radii, and a scatter compatible with random fluctuations around this value, when the gradient is normalized to the disk effective radius. The slope is independent of morphology, incidence of bars, absolute magnitude or mass. Only those galaxies with evidence of interactions and/or clear merging systems present a significant shallower gradient, consistent with previous results. The majority of the 94 galaxies with H ii regions detected beyond 2 disk effective radii present a flattening in the oxygen abundance. The flattening is statistically significant. We cannot provide with a conclusive answer regarding the origin of this flattening. However, our results indicate that its origin is most probably related to the secular evolution of galaxies. Finally, we find a drop/truncation of the oxygen abundance in the inner regions for 26 of the galaxies. All of them are non-interacting, mostly unbarred, Sb/Sbc galaxies. This feature is associated with a central star-forming ring, which suggests that both features are produced by radial gas flows induced by resonance processes.

Over the past decade, 3D optical spectroscopy has become the preferred tool for understanding the properties of galaxies and is now increasingly used to carry out galaxy surveys. Low redshift surveys include SAURON, DiskMass, ATLAS3D, PINGS and VENGA. At redshifts above 0.7, surveys such as MASSIV, SINS, GLACE, and IMAGES have targeted the most luminous galaxies to study mainly their kinematic properties. The on-going CALIFA survey ($z\sim0.02$) is the first of a series of upcoming Integral Field Spectroscopy (IFS) surveys with large samples representative of the entire population of galaxies. Others include SAMI and MaNGA at lower redshift and the upcoming KMOS surveys at higher redshift. Given the importance of spatial scales in IFS surveys, the study of the effects of spatial resolution on the recovered parameters becomes important. We explore the capability of the CALIFA survey and a hypothetical higher redshift survey to reproduce the properties of a sample of objects observed with better spatial resolution at lower redshift. Using a sample of PINGS galaxies, we simulate observations at different redshifts. We then study the behaviour of different parameters as the spatial resolution degrades with increasing redshift.

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.

We present the first galactic-scale model of the gas dynamics of the prototype barred Seyfert 1 galaxy NGC1097. We use large scale FaNTOmM Fabry-Perot interferometric data covering the entire galactic disc and combine the distribution and kinematics maps with high resolution two-dimensional spectroscopy from the Gemini telescope. We build a dynamical model for the gravitational potential by applying the analytic solution to the equations of motion, within the epicyclic approximation. Our model reproduces all the significant kinematic and structural signatures of this galaxy. We find that the primary bar is 7.9+/-0.6 kpc long and has a pattern speed of 36 +/- 2 km s^-1 kpc^-1. This places the corotation radius at 8.6 +/-0.5 kpc, the outer Lindblad resonance at 14.9+/-0.9 kpc and two inner Lindblad resonances at 60+/-5 pc and 2.9+/-0.1 kpc. These derivations lead to a ratio of the corotation radius over bar length of 1.0--1.2, which is in agreement with the predictions of simulations for fast galaxy bars. Our model presents evidence that the circumnuclear ring in this galaxy is not located near any of the resonance radii in this galaxy. The ring might have once formed at the outer inner Lindblad resonance radius, and it has been migrating inward, toward the centre of the galactic gravitational potential.

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.

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.

We present a detailed two-dimensional stellar dynamical analysis of a sample of 44 cosmological hydrodynamical simulations of individual central galaxies and their satellites. Kinematic maps of the stellar line-of-sight velocity, velocity dispersion, and higher-order Gauss-Hermite moments $h_3$ and $h_4$ are constructed for each central galaxy and for the most massive satellites. The amount of rotation is quantified using the $\lambda_{\mathrm{R}}$-parameter. The velocity, velocity dispersion, $h_3$, and $h_4$ fields of the simulated galaxies show a diversity similar to observed kinematic maps of early-type galaxies in the ATLAS$^{\rm{3D}}$ survey. This includes fast (regular), slow, and misaligned rotation, hot spheroids with embedded cold disk components as well as galaxies with counter-rotating cores or central depressions in the velocity dispersion. We link the present day kinematic properties to the individual cosmological formation histories of the galaxies. In general, major galaxy mergers have a significant influence on the rotation properties resulting in both a spin-down as well as a spin-up of the merger remnant. Lower mass galaxies with significant in-situ formation of stars, or with additional gas-rich major mergers - resulting in a spin-up - in their formation history, form elongated fast rotators with a clear anti-correlation of $h_3$ and $v/\sigma$. An additional formation path for fast rotators includes gas poor major mergers leading to a spin-up of the remnants. This formation path does not result in anti-correlated $h_3$ and $v/\sigma$. The galaxies most consistent with the rare class of non-rotating round early-type galaxies grow by gas-poor minor mergers alone. In general, more massive galaxies have less in-situ star formation since $z \sim 2$, rotate slower and have older stellar populations. (shortened)

We study the radial structure of the stellar mass surface density ($\mu$) and stellar population age as a function of the total stellar mass and morphology for a sample of 107 galaxies from the CALIFA survey. We use the fossil record to recover the star formation history (SFH) in spheroidal and disk dominated galaxies with masses from 10$^9$ to 10$^{12}$ M$_\odot$. We derive the half mass radius, and we find that galaxies are on average 15% more compact in mass than in light. HMR/HLR decreases with increasing mass for disk galaxies, but is almost constant in spheroidal galaxies. We find that the galaxy-averaged stellar population age, stellar extinction, and $\mu$ are well represented by their values at 1 HLR. Negative radial gradients of the stellar population ages support an inside-out formation. The larger inner age gradients occur in the most massive disk galaxies that have the most prominent bulges; shallower age gradients are obtained in spheroids of similar mass. Disk and spheroidal galaxies show negative $\mu$ gradients that steepen with stellar mass. In spheroidal galaxies $\mu$ saturates at a critical value that is independent of the galaxy mass. Thus, all the massive spheroidal galaxies have similar local $\mu$ at the same radius (in HLR units). The SFH of the regions beyond 1 HLR are well correlated with their local $\mu$, and follow the same relation as the galaxy-averaged age and $\mu$; suggesting that local stellar mass surface density preserves the SFH of disks. The SFH of bulges are, however, more fundamentally related to the total stellar mass, since the radial structure of the stellar age changes with galaxy mass even though all the spheroid dominated galaxies have similar radial structure in $\mu$. Thus, galaxy mass is a more fundamental property in spheroidal systems while the local stellar mass surface density is more important in disks.

Galaxies, which often contain ionised gas, sometimes also exhibit a so-called low-ionisation nuclear emission line region (LINER). For 30 years this was attributed to a central mass-accreting supermassive black hole (AGN) of low luminosity, making LINER galaxies the largest AGN-sub-population, dominating in numbers over higher luminosity Seyfert galaxies and quasars. This, however, poses a serious problem. While the inferred energy balance is plausible, many LINERs clearly do not contain any other independent signatures of an AGN. Using integral field spectroscopic data from the CALIFA survey, we aim at comparing the observed radial surface brightness profiles with what is expected from illumination by an AGN. Essential for this analysis is a proper extraction of emission-lines, especially weak lines such as the Balmer Hb line which is superposed on an absorption trough. To accomplish this, we use the GANDALF code which simultaneously fits the underlying stellar continuum and emission lines. We show for 48 galaxies with LINER-like emission, that the radial emission-line surface brightness profiles are inconsistent with ionisation by a central point-source and hence cannot be due to an AGN alone. The most probable explanation for the excess LINER-like emission is ionisation by evolved stars during the short but very hot and energetic phase known as post-AGB. This leads us to an entirely new interpretation. Post-AGB stars are ubiquitous and their ionising effect should be potentially observable in every galaxy with gas present and stars older than ~1 Gyr, unless a stronger radiation field from young hot stars or an AGN outshines them. This means that galaxies with LINER-like emission are in fact not a class defined by a property, but rather by the absence of a property. It also explains why LINER emission is observed mostly in massive galaxies with old stars and little star formation.

Using spectroscopy from the Large Binocular Telescope and imaging from the Hubble Space Telescope we discovered the first strong galaxy lens at z(lens)>1. The lens has a secure photometric redshift of z=1.53+/-0.09 and the source is spectroscopically confirmed at z=3.417. The Einstein radius (0.35"; 3.0 kpc) encloses 7.6 x 10^10 Msol, with an upper limit on the dark matter fraction of 60%. The highly magnified (40x) source galaxy has a very small stellar mass (~10^8 Msol) and shows an extremely strong [OIII]_5007A emission line (EW_0 ~ 1000A) bolstering the evidence that intense starbursts among very low-mass galaxies are common at high redshift.

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.

We have developed spherically symmetric dynamical models of dwarf spheroidal galaxies using Schwarzschild's orbit superposition method. This type of modelling yields constraints both on the total mass distribution (e.g. enclosed mass and scale radius) as well as on the orbital structure of the system (e.g. velocity anisotropy). This method is thus less prone to biases introduced by assumptions in comparison to the more commonly used Jeans modelling, and it allows us to derive the dark matter content in a robust way. Here we present our results for the Sculptor dwarf spheroidal galaxy, after testing our methods on mock data sets. We fit both the second and fourth velocity moment profile to break the mass-anisotropy degeneracy. For an NFW dark matter halo profile, we find that the mass of Sculptor within 1 kpc is M_1kpc = 1.03 \pm 0.07 x 10^8 Msol, and that its velocity anisotropy profile is tangentially biased and nearly constant with radius. The preferred concentration (c ~ 15) is low for its dark matter mass but consistent within the scatter found in N-body cosmological simulations. When we let the value of the central logarithmic slope \alpha{} vary, we find that the best-fit model has \alpha{} = 0, although an NFW cusp or shallower is consistent at 1 \sigma{} confidence level. On the other hand, very cuspy density profiles with logarithmic central slopes \alpha{} < -1.5 are strongly disfavoured for Sculptor.

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\alpha\ lines across a large part of the remnant.

The use of IFS is since recently allowing to measure the emission line fluxes of an increasingly large number of star-forming galaxies both locally and at high redshift. The main goal of this study is to review the most widely used empirical oxygen calibrations, O3N2 and N2, by using new direct abundance measurements. We pay special attention to the expected uncertainty of these calibrations as a function of the index value or abundance derived and the presence of possible systematic offsets. This is possible thanks to the analysis of the most ambitious compilation of Te-based HII regions to date. This new dataset compiles the Te-based abundances of 603 HII regions extracted from the literature but also includes new measurements from the CALIFA survey. Besides providing new and improved empirical calibrations for the gas abundance, we also present here a comparison between our revisited calibrations with a total of 3423 additional CALIFA HII complexes with abundances derived using the ONS calibration by Pilyugin et al. (2010). The combined analysis of Te-based and ONS abundances allows us to derive their most accurate calibration to date for both the O3N2 and N2 single-ratio indicators, in terms of all statistical significance, quality and coverage of the space of parameters. In particular, we infer that these indicators show shallower abundance dependencies and statistically-significant offsets compared to those of Pettini and Pagel (2004), Nagao et al. (2006) and P\'erez-Montero and Contini (2009). The O3N2 and N2 indicators can be empirically applied to derive oxygen abundances calibrations from either direct abundance determinations with random errors of 0.18 and 0.16, respectively, or from indirect ones (but based on a large amount of data) reaching an average precision of 0.08 and 0.09 dex (random) and 0.02 and 0.08 dex (systematic; compared to the direct estimations),respectively.

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.

We use deep integral field spectroscopy data from the CALIFA survey to study the warm interstellar medium (WIM) of 32 nearby early-type galaxies (ETGs). We propose a tentative subdivision of our sample ETGs into two groups, according to their Ha equivalent width (EW) and Lyman continuum (LyC) photon escape fraction (PLF). Type i ETGs show nearly constant EWs and a PLF~0, suggesting that photoionization by post-AGB stars is the main driver of their faint extranuclear nebular emission. Type ii ETGs are characterized by very low, outwardly increasing EWs, and a PLF as large as ~0.9 in their centers. Such properties point to a low, and inwardly decreasing WIM density and/or volume filling factor. We argue that, because of extensive LyC photon leakage, emission-line luminosities and EWs are reduced in type ii ETG nuclei by at least one order of magnitude. Consequently, the line weakness of these ETGs is by itself no compelling evidence for their containing merely "weak" (sub-Eddington accreting) active galactic nuclei (AGN). In fact, LyC photon escape, which has heretofore not been considered, may constitute a key element in understanding why many ETGs with prominent signatures of AGN activity in radio continuum and/or X-ray wavelengths show only faint emission lines and weak signatures of AGN activity in their optical spectra. The LyC photon escape, in conjunction with dilution of nuclear EWs by line-of-sight integration through a triaxial stellar host, can systematically impede detection of AGN in gas-poor galaxy spheroids through optical emission-line spectroscopy. We further find that type i and ii ETGs differ little (~0.4 dex) in their mean BPT line ratios, which in both cases are characteristic of LINERs. This potentially hints at a degeneracy of the projected, luminosity-weighted BPT ratios for the specific 3D properties of the WIM in ETGs. (abridged)

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] & [$\alpha$/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.

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.

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.

We present the results on the study of the global and local M-Z relation based on the first data available from the CALIFA survey (150 galaxies). This survey provides integral field spectroscopy of the complete optical extent of each galaxy (up to 2-3 effective radii), with enough resolution to separate individual HII regions and/or aggregations. Nearly $\sim$3000 individual HII regions have been detected. The spectra cover the wavelength range between [OII]3727 and [SII]6731, with a sufficient signal-to-noise to derive the oxygen abundance and star-formation rate associated with each region. In addition, we have computed the integrated and spatially resolved stellar masses (and surface densities), based on SDSS photometric data. We explore the relations between the stellar mass, oxygen abundance and star-formation rate using this dataset. We derive a tight relation between the integrated stellar mass and the gas-phase abundance, with a dispersion smaller than the one already reported in the literature ($\sigma_{\Delta{\rm log(O/H)}}=$0.07 dex). Indeed, this dispersion is only slightly larger than the typical error derived for our oxygen abundances. However, we do not find any secondary relation with the star-formation rate, other than the one induced due to the primary relation of this quantity with the stellar mass. We confirm the result using the $\sim$3000 individual HII regions, for the corresponding local relations. Our results agree with the scenario in which gas recycling in galaxies, both locally and globally, is much faster than other typical timescales, like that of gas accretion by inflow and/or metal loss due to outflows. In essence, late-type/disk dominated galaxies seem to be in a quasi-steady situation, with a behavior similar to the one expected from an instantaneous recycling/closed-box model.

This work investigates the effect of the aperture size on derived galaxy properties for which we have spatially-resolved optical spectra. We focus on some indicators of star formation activity and dust attenuation for spiral galaxies that have been widely used in previous work on galaxy evolution. We have used 104 spiral galaxies from the CALIFA survey for which 2D spectroscopy with complete spatial coverage is available. From the 3D cubes we have derived growth curves of the most conspicuous Balmer emission lines (Halpha, Hbeta) for circular apertures of different radii centered at the galaxy's nucleus after removing the underlying stellar continuum. We find that the Halpha flux (f(Halpha)) growth curve follows a well defined sequence with aperture radius showing low dispersion around the median value. From this analysis, we derive aperture corrections for galaxies in different magnitude and redshift intervals. Once stellar absorption is properly accounted for, the f(Halpha)/f(Hbeta) ratio growth curve shows a smooth decline, pointing towards the absence of differential dust attenuation as a function of radius. Aperture corrections as a function of the radius are provided in the interval [0.3,2.5]R_50. Finally, the Halpha equivalent width (EW(Halpha)) growth curve increases with the size of the aperture and shows a very large dispersion for small apertures. This large dispersion prevents the use of reliable aperture corrections for this quantity. In addition, this result suggests that separating star-forming and quiescent galaxies based on observed EW(Halpha) through small apertures is likely to result in low EW(Halpha) star-forming galaxies begin classified as quiescent.

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.

We present the first public data release of the CALIFA survey. It consists of science-grade optical datacubes for the first 100 of eventually 600 nearby (0.005<z<0.03) galaxies, obtained with the integral-field spectrograph PMAS/PPak mounted on the 3.5m telescope at the Calar Alto observatory. The galaxies in DR1 already cover a wide range of properties in color-magnitude space, morphological type, stellar mass, and gas ionization conditions. This offers the potential to tackle a variety of open questions in galaxy evolution using spatially resolved spectroscopy. Two different spectral setups are available for each galaxy, (i) a low-resolution V500 setup covering the nominal wavelength range 3745-7500A with a spectral resolution of 6.0A (FWHM), and (ii) a medium-resolution V1200 setup covering the nominal wavelength range 3650-4840A with a spectral resolution of 2.3A (FWHM). We present the characteristics and data structure of the CALIFA datasets that should be taken into account for scientific exploitation of the data, in particular the effects of vignetting, bad pixels and spatially correlated noise. The data quality test for all 100 galaxies showed that we reach a median limiting continuum sensitivity of 1.0x10^-18erg/s/cm^2/A/arcsec^2 at 5635A and 2.2x10^-18erg/s/cm^2/A/arcsec^2 at 4500A for the V500 and V1200 setup respectively, which corresponds to limiting r and g band surface brightnesses of 23.6mag/arcsec^2 and 23.4mag/arcsec^2, or an unresolved emission-line flux detection limit of roughly 1x10^-17erg/s/cm^2/arcsec^2 and 0.6x10^-17erg/s/cm^2/arcsec^2, respectively. The median spatial resolution is 3.7", and the absolute spectrophotometric calibration is better than 15% (1sigma). We also describe the available interfaces and tools that allow easy access to this first public CALIFA data.

The growth of galaxies is one of the key problems in understanding the structure and evolution of the universe and its constituents. Galaxies can grow their stellar mass by accretion of halo or intergalactic gas clouds, or by merging with smaller or similar mass galaxies. The gas available translates into a rate of star formation, which controls the generation of metals in the universe. The spatially resolved history of their stellar mass assembly has not been obtained so far for any given galaxy beyond the Local Group. Here we demonstrate how massive galaxies grow their stellar mass inside-out. We report the results from the analysis of the first 105 galaxies of the largest to date three-dimensional spectroscopic survey of galaxies in the local universe (CALIFA). We apply the fossil record method of stellar population spectral synthesis to recover the spatially and time resolved star formation history of each galaxy. We show, for the first time, that the signal of downsizing is spatially preserved, with both inner and outer regions growing faster for more massive galaxies. Further, we show that the relative growth rate of the spheroidal component, nucleus and inner galaxy, that happened 5-7 Gyr ago, shows a maximum at a critical stellar mass ~10^10 Msun. We also find that galaxies less massive than ~10^10 Msun show a transition to outside-in growth, thus connecting with results from resolved studies of the growth of low mass galaxies.

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.

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.

We place our sample of 18 Virgo dwarf early-type galaxies (dEs) on the V-K - velocity dispersion, Faber-Jackson, and Fundamental Plane (FP) scaling relations for massive early-type galaxies (Es). We use a generalized velocity dispersion, which includes rotation, to be able to compare the location of both rotationally and pressure supported dEs with those of early and late-type galaxies. We find that dEs seem to bend the Faber-Jackson relation of Es to lower velocity dispersions, being the link between Es and dwarf spheroidal galaxies (dSphs). Regarding the FP relation, we find that dEs are significantly offset with respect to massive hot stellar systems, and re-casting the FP into the so-called kappa-space suggests that this offset is related to dEs having a total mass-to-light ratio higher than Es but still significantly lower than dSph galaxies. Given a stellar mass-to-light ratio based on the measured line indices of dEs, the FP offset allows us to infer that the dark matter fraction within the half light radii of dEs is on average >~ 42% (uncertainties of 17% in the K band and 20% in the V band), fully consistent with an independent estimate in an earlier paper in this series. We also find that dEs in the size-luminosity relation in the near-infrared, like in the optical, are offset from early-type galaxies, but seem to be consistent with late-type galaxies. We thus conclude that the scaling relations show that dEs are different from Es, and that they further strengthen our previous findings that dEs are closer to and likely formed from late-type galaxies.

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.

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.

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.

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{\beta} index and metallicity [Z/H]. The equivalent correlations for the Fe5015 index, {\alpha}-enhancement [{\alpha}/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.

We combine the estimated metallicities [Fe/H], abundances [\alpha/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 [\alpha/Fe] vs. [Fe/H], we find that there is a smooth variation with [\alpha/Fe], a proxy for age. At the same time, the orbital properties of the old stars with [\alpha/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 \alpha-young 'thin-disk' stars a significant contribution to the \alpha-old 'thick-disk' metal-rich stars. However, the remaining \alpha-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.

Context. We carry out a detailed 2D study of the ionised gas in the local universe galaxy UGC9837. In nearby galaxies, like the galaxy in question here, the spatial distribution of the physical properties can be studied in detail, providing benchmarks for galaxy formation models. Aims. Our aim is to derive detailed and spatially resolved physical properties of the ionised gas of UGC 9837. In addition, we derive an integrated spectrum of the galaxy and study how varying spatial coverage affects the derived integrated properties. We also study how the same properties would be seen if the galaxy was placed at a higher redshift and observed as part of one of the high-z surveys. Methods. UGC9837 was observed using the PMAS PPAK integral field unit. The spectra are reduced and calibrated and the stellar and ionised components separated. Using strong emission line ratios of the ionised gas, the source of ionisation, the dust extinction, the star formation rate, the electron density and the oxygen abundance derived from a total integrated spectrum, central integrated spectrum and individual fibre spectra are studied. Finally, the same properties are studied in a spectrum whose spatial resolution is degraded to simulate high-z observations. Results. The spatial distribution of the ionised gas properties is consistent with inside-out growing scenario of galaxies. We also find that lack of spatial coverage would bias the results derived from the integrated spectrum leading, e.g., to an under-estimation of ionisation and over-estimation of metallicity, if only the centre of the galaxy was covered by the spectrum. Our simulation of high-z observations shows that part of the spatial information, such as dust and SFR distribution would be lost while shallower gradients in metallicity and ionisation strength would be detected.

We present here the Calar Alto Legacy Integral Field Area (CALIFA) survey, which has been designed to provide a first step in this direction.We summarize the survey goals and design, including sample selection and observational strategy.We also showcase the data taken during the first observing runs (June/July 2010) and outline the reduction pipeline, quality control schemes and general characteristics of the reduced data. This survey is obtaining spatially resolved spectroscopic information of a diameter selected sample of $\sim600$ galaxies in the Local Universe (0.005< z <0.03). CALIFA has been designed to allow the building of two-dimensional maps of the following quantities: (a) stellar populations: ages and metallicities; (b) ionized gas: distribution, excitation mechanism and chemical abundances; and (c) kinematic properties: both from stellar and ionized gas components. CALIFA uses the PPAK Integral Field Unit (IFU), with a hexagonal field-of-view of $\sim1.3\sq\arcmin'$, with a 100% covering factor by adopting a three-pointing dithering scheme. The optical wavelength range is covered from 3700 to 7000 {\AA}, using two overlapping setups (V500 and V1200), with different resolutions: R\sim850 and R\sim1650, respectively. CALIFA is a legacy survey, intended for the community. The reduced data will be released, once the quality has been guaranteed. The analyzed data fulfill the expectations of the original observing proposal, on the basis of a set of quality checks and exploratory analysis. We conclude from this first look at the data that CALIFA will be an important resource for archaeological studies of galaxies in the Local Universe.

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.

We have developed a spherically symmetric dynamical model of a dwarf spheroidal galaxy using the Schwarzschild method. This type of modelling yields constraints both on the total mass distribution (e.g. enclosed mass and scale radius) as well as on the orbital structure of the system modelled (e.g. velocity anisotropy). Therefore not only can we derive the dark matter content of these systems, but also explore possible formation scenarios. Here we present preliminary results for the Sculptor dSph. We find that the mass of Sculptor within 1kpc is 8.5\times10^(7\pm0.05) M\odot, its anisotropy profile is tangentially biased and slightly more isotropic near the center. For an NFW profile, the preferred concentration (~15) is compatible with cosmological models. Very cuspy density profiles (steeper than NFW) are strongly disfavoured for Sculptor.

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...

Grigori Kuzmin was a very gifted dynamicist and one of the towering figures in the distinguished history of the Tartu Observatory. He obtained a number of important results in relative isolation which were later rediscovered in the West. This work laid the foundation for further advances in the theory of stellar systems in dynamical equilibrium, thereby substantially increasing our understanding of galaxy dynamics.

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.

[Abridged] We present ground-based MDM V-band and Spitzer/IRAC 3.6um-band photometric observations of the 72 representative galaxies of the SAURON Survey. In combination with the SAURON stellar velocity dispersion measured within an effective radius (se), this allows us to explore the location of our galaxies in the main scaling relations. We investigate the dependence of these relations on our recent kinematical classification of early-type galaxies (i.e. Slow/Fast Rotators) and the stellar populations. Slow Rotator and Fast Rotator E/S0 galaxies do not populate distinct locations in the scaling relations, although Slow Rotators display a smaller intrinsic scatter. Surprisingly, extremely young objects do not display the bluest (V-[3.6]) colours in our sample, as is usually the case in optical colours. This can be understood in the context of the large contribution of TP-AGB stars to the infrared, even for young populations, resulting in a very tight (V-[3.6]) - se relation that in turn allows us to define a strong correlation between metallicity and velocity dispersion. Many Sa galaxies appear to follow the Fundamental Plane defined by E/S0 galaxies. Galaxies that appear offset from the relations correspond mostly to objects with extremely young populations, with signs of on-going, extended star formation. We correct for this effect in the Fundamental Plane, by replacing luminosity with stellar mass using an estimate of the stellar mass-to-light ratio, so that all galaxies are part of a tight, single relation. The new estimated coefficients are consistent in both photometric bands and suggest that differences in stellar populations account for about half of the observed tilt with respect to the virial prediction. After these corrections, the Slow Rotator family shows almost no intrinsic scatter around the best-fit Fundamental Plane.

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.

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.

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)

[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.

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.

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).

(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\'c 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.

We investigate the accuracy of mass determinations M_BH of supermassive black holes in galaxies using dynamical models of the stellar kinematics. We compare 10 of our M_BH measurements, using integral-field OASIS kinematics, to published values. For a sample of 25 galaxies we confront our new M_BH derived using two modeling methods on the same OASIS data.

We combine SAURON integral field data of a representative sample of local early-type, red sequence galaxies with Spitzer/IRAC imaging in order to investigate the presence of trace star formation in these systems. With the Spitzer data, we identify galaxies hosting low-level star formation, as traced by PAH emission, with measured star formation rates that compare well to those estimated from other tracers. This star formation proceeds according to established scaling relations with molecular gas content, in surface density regimes characteristic of disk galaxies and circumnuclear starbursts. We find that star formation in early-type galaxies happens exclusively in fast-rotating systems and occurs in two distinct modes. In the first, star formation is a diffuse process, corresponding to widespread young stellar populations and high molecular gas content. The equal presence of co- and counter-rotating components in these systems strongly implies an external origin for the star-forming gas, and we argue that these star formation events may be the final stages of (mostly minor) mergers that build up the bulges of red sequence lenticulars. In the second mode of star formation, the process is concentrated into well-defined disk or ring morphologies, outside of which the host galaxies exhibit uniformly evolved stellar populations. This implies that these star formation events represent rejuvenations within previously quiescent stellar systems. Evidence for earlier star formation events similar to these in all fast rotating early-type galaxies suggests that this mode of star formation may be common to all such galaxies, with a duty cycle of roughly 1/10, and likely contributes to the embedded, co-rotating inner stellar disks ubiquitous in this population.

We present a preliminary analysis of the HI properties of early-type galaxies in the ATLAS3D sample. Using WSRT data for ~100 galaxies outside the Virgo cluster and data from the Alfalfa project for galaxies inside Virgo, we discuss the dependence of HI properties on environment. We detect HI in about half of the galaxies outside Virgo. For these systems, the HI morphology and kinematics change as a function of environment, going from regular, rotating systems around isolated galaxies to progressively more disturbed structures for galaxies with neighbours or in groups. In denser environment, inside Virgo, nearly none of the galaxies contains HI. We discuss future work in this field which will be enabled by next-generation, pre-SKA radio instruments. We present a simulated Apertif HI observation of an ATLAS3D early-type galaxy, showing how its appearance and detection level vary as a function of redshift.

We study the origin of the scaling relations of early-type galaxies (ETGs) by constructing detailed models of the stellar dynamics for the K-band selected, volume-limited ATLAS3D sample of 263 nearby ETGs, spanning a large range of masses and stellar velocity dispersions (60 < sigma < 350 km/s).

We investigate how well the intrinsic shape of early-type galaxies can be recovered when both photometric and two-dimensional stellar kinematic observations are available. We simulate these observations with galaxy models that are representative of observed oblate fast-rotator to triaxial slow-rotator early-type galaxies. By fitting realistic triaxial dynamical models to these simulated observations, we recover the intrinsic shape (and mass-to-light ratio), without making additional (ad-hoc) assumptions on the orientation. For (near) axisymmetric galaxies the dynamical modelling can strongly exclude triaxiality, but the regular kinematics do not further tighten the constraint on the intrinsic flattening significantly, so that the inclination is nearly unconstrained above the photometric lower limit even with two-dimensional stellar kinematics. Triaxial galaxies can have additional complexity in both the observed photometry and kinematics, such as twists and (central) kinematically decoupled components, which allows the intrinsic shape to be accurately recovered. For galaxies that are very round or show no significant rotation, recovery of the shape is degenerate, unless additional constraints such as from a thin disk are available.

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.

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.

[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.

We present the results of an investigation of the local escape velocity (Vesc) - line strength index relationship for 48 early type galaxies from the SAURON sample, the first such study based on a large sample of galaxies with both detailed integral field observations and extensive dynamical modelling. Values of Vesc are computed using Multi Gaussian Expansion (MGE) photometric fitting and axisymmetric, anisotropic Jeans' dynamical modelling simultaneously on HST and ground-based images. We determine line strengths and escape velocities at multiple radii within each galaxy, allowing an investigation of the correlation within individual galaxies as well as amongst galaxies. We find a tight correlation between Vesc and the line-strength indices. For Mgb we find that this correlation exists not only between different galaxies but also inside individual galaxies - it is both a local and global correlation. The relation within individual galaxies has the same slope and offset as the global relation to a good level of agreement, though there is significant intrinsic scatter in the local gradients. We transform our line strength index measurements to the single stellar population (SSP) equivalent age (t), metallicity ([Z/H]) and enhancement ([$\alpha$/Fe]) and carry out a principal component analysis of our SSP and Vesc data. We find that in this four-dimensional parameter space the galaxies in our sample are to a good approximation confined to a plane, given by $\log \mathrm({V}_{\mathrm{esc}}/500\mathrm{km/s}) = 0.85 \mathrm{[Z/H]} + 0.43 \log (\mathrm{t}/\mathrm{Gyrs})$ - 0.20. It is surprising that it is a combination of age and metallicity that is conserved; this may indicate a 'conspiracy' between age and metallicity or a weakness in the SSP models.

We use the integral-field spectrograph SAURON to measure the stellar line-of-sight velocity distribution and absorption line strengths out to four effective radii (Re) in the early-type galaxies NGC 3379 and NGC 821. With our newly developed observing technique we can now probe these faint regions in galaxies that were previously not accessible with traditional long-slit spectroscopy. We make optimal use of the large field-of-view and high throughput of the spectrograph: by adding the signal of all ~1400 lenslets into one spectrum, we obtain sufficient signal-to-noise in a few hours of observing time to reliably measure the absorption line kinematics and line strengths out to large radius. We find that the line strength gradients previously observed within 1 Re remain constant out to at least 4 Re, which puts constraints on the merger histories of these galaxies. The stellar halo populations are old and metal-poor. By constructing orbit-based Schwarzschild dynamical models we find that dark matter is necessary to explain the observed kinematics in NGC 3379 and NGC 821, with 30 - 50 per cent of the total matter being dark within 4 Re. The radial anisotropy in our best-fit halo models is less than in our models without halo, due to differences in orbital structure. The halo also has an effect on the Mgb - Vesc relation: its slope is steeper when a dark matter halo is added to the model.

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.

In order to attract and retain excellent researchers and diverse individuals in astrophysics, we recommend action be taken in several key areas impacting young scientists: (1) Maintain balance between large collaborations and individual projects through distribution of funding; encourage public releases of observational and simulation data for use by a broader community. (2) Improve the involvement of women, particularly at leading institutions. (3) Address the critical shortage of child care options and design reasonable profession-wide parental leave policies. (4) Streamline the job application and hiring process. We summarize our reasons for bringing these areas to the attention of the committee, and we suggest several practical steps that can be taken to address them.

Whereas considerable effort has been afforded in understanding the properties of galaxies, a full physical picture, connecting their baryonic and dark-matter content, super-massive black holes, and (metric) theories of gravity, is still ill-defined. Strong gravitational lensing furnishes a powerful method to probe gravity in the central regions of galaxies. It can (1) provide a unique detection-channel of dark-matter substructure beyond the local galaxy group, (2) constrain dark-matter physics, complementary to direct-detection experiments, as well as metric theories of gravity, (3) probe central super-massive black holes, and (4) provide crucial insight into galaxy formation processes from the dark matter point of view, independently of the nature and state of dark matter. To seriously address the above questions, a considerable increase in the number of strong gravitational-lens systems is required. In the timeframe 2010-2020, a staged approach with radio (e.g. EVLA, e-MERLIN, LOFAR, SKA phase-I) and optical (e.g. LSST and JDEM) instruments can provide 10^(2-4) new lenses, and up to 10^(4-6) new lens systems from SKA/LSST/JDEM all-sky surveys around ~2020. Follow-up imaging of (radio) lenses is necessary with moderate ground/space-based optical-IR telescopes and with 30-50m telescopes for spectroscopy (e.g. TMT, GMT, ELT). To answer these fundamental questions through strong gravitational lensing, a strong investment in large radio and optical-IR facilities is therefore critical in the coming decade. In particular, only large-scale radio lens surveys (e.g. with SKA) provide the large numbers of high-resolution and high-fidelity images of lenses needed for SMBH and flux-ratio anomaly studies.

The full spatio-chemo-dynamical structure of galaxies of all types and environments at low redshift provides a critical accompaniment to observations of galaxy formation at high redshift. The next decade brings the observational opportunity to strongly constrain nearby galaxies' histories of star formation and assembly, especially in the spheroids that comprise the large majority of the stellar mass in the Universe but have until now been difficult to study. In order to constrain the pathways to building up the spheroidal "red-sequence", various standard techniques in photometry and spectroscopy, particularly with resolved tracer populations like globular clusters and planetary nebulae, can be scaled up to comprehensive surveys as improved wide-field instrumentation is increasingly available. At the same time, progress in adaptive optics on giant telescopes could for the first time permit deep, resolved photometric and spectroscopic analysis of large samples of individual stars in these systems, thereby revolutionizing galaxy studies. Strong theoretical support is needed in order to understand the new observational constraints via detailed modeling and self-consistent simulations of star and galaxy formation throughout cosmic time.

We describe a new technique to measure stellar kinematics and line-strengths at large radii in nearby galaxies. Using the integral-field spectrograph SAURON as a 'photon-collector', we obtain spectra out to four effective radii (Re) in the early-type galaxy NGC 3379. By fitting orbit-based models to the extracted stellar velocity profile, we find that ~40% of the total mass within 5 Re is dark. The measured absorption line-strengths reveal a radial gradient with constant slope out to 4 Re.

[Abridged] We analysed two-dimensional maps of 48 early-type galaxies obtained with the SAURON and OASIS integral-field spectrographs using kinemetry, a generalisation of surface photometry to the higher order moments of the line-of-sight velocity distribution (LOSVD). In the SAURON sample, we find that 31% of early-type galaxies are single component systems. 91% of the multi-components systems have two kinematic subcomponents, the rest having three. In addition, 29% of galaxies have kinematically decoupled components, nuclear components with significant kinematic twists. We find that the velocity maps of fast rotators closely resemble those of inclined disks, except in the transition regions between kinematic subcomponents. In terms of E/S0 classification, this means that 74% of Es and 92% of S0s have components with disk-like kinematics. For the majority of fast rotators, the kinematic axial ratios are equal to or less than their photometric axial ratios, contrary to what is predicted with isotropic Jeans models viewed at different inclinations. The position angles of fast rotators are constant, while they vary abruptly in slow rotators. Velocity dispersion maps of face-on galaxies have shapes similar to the distribution of light. We constructed local (bin-by-bin) h3 - V/sigma and h4 - V/sigma diagrams from SAURON observations. We confirm the classical anti-correlation of h3 and V\sigma, but we also find that h3 is almost zero in some objects or even weakly correlated with V/sigma. The distribution of h4 for fast and slow rotators is mildly positive on average. The difference between slow and fast rotators is traceable throughout all moments of the LOSVD, with evidence for different intrinsic shapes and orbital contents and, hence, likely different evolutionary paths.

We construct axisymmetric and triaxial galaxy models with a phase-space distribution function that depends on linear combinations of the three exact integrals of motion for a separable potential. These Abel models, first introduced by Dejonghe & Laurent and subsequently extended by Mathieu & Dejonghe, are the axisymmetric and triaxial generalisations of the well-known spherical Osipkov-Merritt models. We show that the density and higher order velocity moments, as well as the line-of-sight velocity distribution (LOSVD) of these models can be calculated efficiently and that they capture much of the rich internal dynamics of early-type galaxies. We build a triaxial and oblate axisymmetric galaxy model with projected kinematics that mimic the two-dimensional kinematic observations that are obtained with integral-field spectrographs such as SAURON. We fit the simulated observations with axisymmetric and triaxial dynamical models constructed with our numerical implementation of Schwarzschild's orbit-superposition method. We find that Schwarzschild's method is able to recover the internal dynamics and three-integral distribution function of realistic models of early-type galaxies.

We present a flexible and efficient method to construct triaxial dynamical models of galaxies with a central black hole, using Schwarzschild's orbital superposition approach. Our method is general and can deal with realistic luminosity distributions, which project to surface brightness distributions that may show position angle twists and ellipticity variations. The models are fit to measurements of the full line-of-sight velocity distribution (wherever available). We verify that our method is able to reproduce theoretical predictions of a three-integral triaxial Abel model. In a companion paper (van de Ven, de Zeeuw & van den Bosch), we demonstrate that the method recovers the phase-space distribution function. We apply our method to two-dimensional observations of the E3 galaxy NGC 4365, obtained with the integral-field spectrograph SAURON, and study its internal structure, showing that the observed kinematically decoupled core is not physically distinct from the main body and the inner region is close to oblate axisymmetric.

We present HI observations of the elliptical galaxy NGC 2974, obtained with the Very Large Array. These observations reveal that the previously detected HI disc in this galaxy (Kim et al. 1988) is in fact a ring. By studying the harmonic expansion of the velocity field along the ring, we constrain the elongation of the halo and find that the underlying gravitational potential is consistent with an axisymmetric shape. We construct mass models of NGC 2974 by combining the HI rotation curve with the central kinematics of the ionised gas, obtained with the integral-field spectrograph SAURON. We introduce a new way of correcting the observed velocities of the ionised gas for asymmetric drift, and hereby disentangle the random motions of the gas caused by gravitational interaction from those caused by turbulence. To reproduce the observed flat rotation curve of the HI gas, we need to include a dark halo in our mass models. A pseudo-isothermal sphere provides the best model to fit our data, but we also tested an NFW halo and Modified Newtonian Dynamics (MOND), which fit the data marginally worse. The mass-to-light ratio M/L_I increases in NGC 2974 from 4.3 (M/L_I)sun at one effective radius to 8.5 (M/L_I)sun at 5 Re. This increase of M/L already suggests the presence of dark matter: we find that within 5 Re at least 55 per cent of the total mass is dark.

For a density that is not too sharply peaked towards the center, the local tidal field becomes compressive in all three directions. Available gas can then collapse and form a cluster of stars in the center, including or even being dominated by a central black hole. We show that for a wide range of (deprojected) Sersic profiles in a spherical potential, the tidal forces are compressive within a region which encloses most of the corresponding light of observed nuclear clusters in both late-type and early-type galaxies. In such models, tidal forces become disruptive nearly everywhere for relatively large Sersic indices n >= 3.5. We also show that the mass of a central massive object (CMO) required to remove all radial compressive tidal forces scales linearly with the mass of the host galaxy. If CMOs formed in (progenitor) galaxies with n ~ 1, we predict a mass fraction of ~ 0.1-0.5%, consistent with observations of nuclear clusters and super-massive black holes. While we find that tidal compression possibly drives the formation of CMOs in galaxies, beyond the central regions and on larger scales in clusters disruptive tidal forces might contribute to prevent gas from cooling.

We discuss some recent integral field spectroscopy using the SAURON instrument of a sample consisting of 24 early-type spirals, part of the SAURON Survey, and 18 late-type spirals. Using 2-dimensional maps of their stellar radial velocity, velocity dispersion, and absorption line strength, it is now much easier to understand the nature of nearby galactic bulges. We discuss a few highlights of this work, and point out some new ideas about the formation of galactic bulges.

Using the unique dataset obtained within the course of the SAURON project, a radically new view of the structure, dynamics and stellar populations of early-type galaxies has emerged. We show that galaxies come in two broad flavours (slow and fast rotators), depending on whether or not they exhibit clear large-scale rotation, as indicated via a robust measure of the specific angular momentum of baryons. This property is also linked with other physical characteristics of early-type galaxies, such as: the presence of dynamically decoupled cores, orbital structure and anisotropy, stellar populations and dark matter content. I here report on the observed link between this baryonic angular momentum and a mass sequence, and how this uniquely relates to the building of the red sequence via dissipative/dissipationless mergers and secular evolution.

Supermassive black holes are a key element in our understanding of how galaxies form. Most of the progress in this very active field of research is based on just ~30 determinations of black hole mass, accumulated over the past decade. We illustrate how integral-field spectroscopy, and in particular our OASIS modeling effort, can help improve the current situation.

Motivated by recent progress in the study of early-type galaxies owing to technological advances, the launch of new space telescopes and large ground-based surveys, we attempt a short review of our current understanding of the recent star-formation activity in such intriguing galactic systems.

In this proceeding we make use of the two-dimensional stellar kinematics of a representative sample of E and S0 galaxies obtained with the SAURON integral-field spectrograph to reveal that early-type galaxies appear in two broad flavours, depending on whether they exhibit clear large-scale rotation or not. We measure the level of rotation via a new parameter LambdaR and use it as a basis for a new kinematic classification that separates early-type galaxies into slow and fast rotators. With the aid of broad-band imaging we will reinforce this finding by comparing our kinematic results to the photometric properties of these two classes.

In this proceeding we look at the relationship between the photometric nuclear properties of early-type galaxies from Hubble Space Telescope imaging and their overall kinematics as observed with the SAURON integral-field spectrograph. We compare the inner slope of their photometric profiles and the Slow/Fast rotator classes, defined by the amplitude of a newly defined LambdaR parameter, to show that slow rotators tend to be more massive systems and display shallower inner profiles and fast rotators steper ones. It is important to remark, however, that there is not a one-to-one relationship between the two photometric and kinematic groups.

2D stellar kinematics of 48 representative E and S0 galaxies obtained with the SAURON IFS reveal that early-type galaxies appear in two broad flavours, depending on whether they exhibit clear large-scale rotation or not. We define a new parameter LambdaR which involves luminosity weighted averages over the full 2D kinematic field, as a proxy to quantify the observed projected stellar angular momentum per unit mass. We use it as a basis for a new kinematic classification: early-type galaxies are separated into slow and fast rotators (SRs, FRs), depending on whether they have LambdaR values within their effective radius Re below or above 0.1, respectively. SRs and FRs are shown to be physically distinct classes of galaxies, a result which cannot simply be the consequence of a biased viewing angle. FRs tend to be relatively low luminosity galaxies. SRs tend to be brighter and more massive galaxies, but are still spread over a wide range of absolute magnitude. 3 slow rotators of our sample, among the most massive ones, are consistent with zero rotation. Remarkably, all other SRs contain a large kpc-scale KDC. All FRs show well aligned photometric and kinemetric axes, and small velocity twists, in contrast with most SRs which exhibit significant misalignments and velocity twists. In a companion paper (Paper X), we also show that FRs and SRs are distinct classes in terms of their orbital distribution. We suggest that gas is a key ingredient in the formation and evolution of FRs, and that the slowest rotators are the extreme evolutionary end point reached deep in gravitational potential wells where dissipationless mergers had a major role in the evolution, and for which most of the baryonic angular momentum was expelled outwards. (abridged)

We present absorption line strength maps of a sample of 24 representative early-type spiral galaxies, mostly of type Sa, obtained as part of the SAURON survey of nearby galaxies using our custom-built integral-field spectrograph. Using high-quality spectra, spatially binned to a constant signal-to-noise, we measure several key age, metallicity and abundance ratio sensitive indices from the Lick/IDS system over a contiguous two-dimensional field including bulge and inner disc. We present maps of H beta, Fe 5015 and Mg b, for each galaxy The absorption line maps show that many galaxies contain some younger populations (<= 1 Gyr), distributed in small or large inner discs, or in circumnuclear star forming rings. In many cases these young stars are formed in circumnuclear mini-starbursts, which are dominating the light in the centres of some of the early-type spirals. These mini-starburst cause a considerable scatter in index-index diagrams such as Mg b- H beta and Mg b -Fe 5015, more than is measured for early-type galaxies. We find that the central regions of Sa galaxies display a wide range in ages, even within the galaxies. 50% of the sample show velocity dispersion drops in their centres. All of the galaxies of our sample lie on or below the Mg b- $\sigma$ relation for elliptical galaxies in the Coma cluster, and above the H beta absorption line - $\sigma$ relation for elliptical galaxies. If those relations are considered to be relations for the oldest local galaxies we see that our sample of spirals has a considerable scatter in age, with the largest scatter at the lowest $\sigma$. This is in disagreement with highly inclined samples, in which generally only old stellar populations are found in the central regions. All this can be understood if ... (see paper for rest of abstract)

We present absorption line-strength maps for a sample of 18 Sb-Sd galaxies observed using the integral-field spectrograph SAURON. The SAURON spectral range allows the measurement of the Lick/IDS indices Hbeta, Fe5015 and Mgb, which can be used to estimate the stellar population parameters. We present here the two-dimensional line-strength maps for each galaxy. From the maps, we learn that late-type spiral galaxies tend to have high Hbeta and low Fe5015 and Mgb values, and that the Hbeta index has often a positive gradient over the field, while the metal indices peak in the central region. We investigate the relations between the central line-strength indices and their correlations with morphological type and central velocity dispersion, and compare the observed behaviour with that for ellipticals, lenticulars and early-type spirals from the SAURON survey. We find that our galaxies lie below the Mg - sigma relation determined for elliptical galaxies and that the indices show a clear trend with morphological type. From the line-strength maps we calculate age, metallicity and abundance ratio maps and discuss the results from a one-SSP approach and from a two-SSP approach. Late-type galaxies are generally younger and more metal poor than ellipticals and have abundance ratios closer to solar values. We also explore a continuous star formation scenario, and try to recover the star formation history using the evolutionary models of Bruzual & Charlot (2003), assuming constant or exponentially declining star formation rate (SFR). We find a correlation between the e-folding time-scale tau of the starburst and the central velocity dispersion: more massive galaxies tend to have shorter tau, suggesting that the star formation happened long ago and has now basically ended, while for smaller objects with larger values of tau it is still active.

We analyse the orbital distribution of E and S0 galaxies using SAURON integral-field stellar kinematics. We construct the anisotropy (V/sigma,epsilon) diagram, for the 48 E/S0 galaxies from the SAURON survey. For a subsample of 24 galaxies consistent with axisymmetry, we use three-integral axisymmetric Schwarzschild dynamical models to recover the detailed orbital distribution and we find good agreement with the anisotropy derived from the (V/sigma,epsilon) diagram. In the companion Paper IX we show that the early-type galaxies can be subdivided into two classes of systems with or without a significant amount of specific stellar angular momentum. Here we show that the two classes have different distributions on the (V/sigma,epsilon) diagram. The slow rotators are more common among the most massive systems and are generally classified as E from photometry alone. Those in our sample tend to be fairly round (epsilon<0.3), but can have significant kinematical misalignments, indicating that as a class they are moderately triaxial, and span a range of anisotropies (delta<0.3). The fast rotators are generally fainter and are classified either E or S0. They can appear quite flattened (epsilon<0.7), do not show significant kinematical misalignments (unless barred or interacting), indicating they are nearly axisymmetric, and span an even larger range of anisotropies (delta<0.5). These results are confirmed when we extend our analysis to 18 additional E/S0 galaxies observed with SAURON. All these results support the idea that fast rotators are nearly oblate and contain disk-like components. The role of gas must have been important for their formation. The slow rotators are weakly triaxial. Current collisionless merger models seem unable to explain their detailed observed properties. (Abridged)

We discuss SAURON absorption line strength maps of a sample of 24 early-type spirals, mostly Sa. From the Lick indices H$\beta$, Mg b and Fe 5015 we derive SSP-ages and metallicities. By comparing the scaling relations of Mg b and H beta and central velocity dispersion with the same relation for the edge-on sample of Falcon-Barroso et al. (2002) we derive a picture in which the central regions of Sa galaxies contain at least 2 components: one (or more) thin, disc-like component, often containing recent star formation, and another, elliptical-like component, consisting of old stars and rotating more slowly, dominating the light above the plane. If one defines a bulge to be the component responsible for the light in excess of the outer exponential disc, then many Sa-bulges are dominated by a thin, disc-like component containing recent star formation.

By detecting ionised-gas emission in 75% of the cases, the SAURON integral-field spectroscopic survey has further demonstrated that early-type galaxies often display nebular emission. Furthermore, the SAURON data have shown that such emission comes with an intriguing variety of morphologies, kinematic behaviours and line ratios. Perhaps most puzzling was the finding that round and slowly rotating objects generally display uncorrelated stellar and gaseous angular momenta, consistent with an external origin for the gas, whereas flatter and fast rotating galaxies host preferentially co-rotating gas and stars, suggesting internal production of gas. Alternatively, a bias against the internal production of ionised gas and against the acquisition of retrograde material may be present in these two kinds of objects, respectively. In light of the different content of hot gas in these systems, with slowly rotating objects being the only systems capable of hosting massive X-ray halos, we suggest that a varying importance of evaporation of warm gas in the hot interstellar medium can contribute to explain the difference in the relative behaviour of gas and stars in these two kinds of objects. Namely, whereas in X-ray bright and slowly rotating galaxies stellar-loss material would quickly evaporate in the hot medium, in X-ray faint and fast rotating objects such material would be allowed to lose angular momentum and settle in a disk, which could also obstruct the subsequent acquisition of retrograde gas. Evidence for a connection between warm and hot gas phases, presumably driven by heat conduction, is presented for four slowly rotating galaxies with CHANDRA observations.

We give two examples of spiral galaxies that show non-circular gas motions in the inner kiloparsecs, from SAURON integral field spectroscopy. We use harmonic decomposition of the velocity field of the ionized gas to study the underlying mass distribution, employing linear theory. The higher order harmonic terms and the main kinematic features of the observed data are consistent with an analytically constructed simple bar model. We also present maps of a number of strong absorption lines in M 100, derive simple stellar populations and correlate them with features in the gas kinematics.

We derive the distribution of maser stars in the inner Milky Way (MW) based on an analysis of lV-diagrams (lVd) for two samples of maser stars: 771 OH/IR stars and 363 SiO-maser stars. They are all close to the plane of the MW and have long. from -45 to +45deg. The two lVds are compared and found to be very similar. They also compare well with the lVd of interstellar CO, but there are significant differences in detail between the stellar lVds and that of the ISM. Based on the qualitative discussion we divide the lVds into seven areas. In each area we compare the number of stars observed with those predicted by an assumed set of orbits in a galactic potential. This potential is axially symmetric but a weak rotating bar has been added. We conclude that the maser stars move on almost circular orbits outside of about 3.5 kpc, but that the orbits become more and more elongated when one goes deep inside our MW. We find a strong effect of the Corotation (CR) resonance (res) at 3.3 kpc, we see a small but noticeable effect of the Outer Lindblad res at 5 kpc and no effect of the Inner Lindblad res at r=0.8 kpc. We find a set of 6 groups of orbits that together predict counts in agreement with the counts of stars observed. We then calculate the trajectory of each orbit and so find the distribution of the maser stars in the plane of the MWG. This distribution has two new (but not unexpected) features. The first is a bar-like distribution within 2 kpc from the GC outlined. These orbits explain the high-vel stars near l=0deg in the forbidden and the permitted quadrants. The second feature are two "croissant"-like voids in the distribution close to the CR radius (3.3 kpc), which are the consequence of the presence of the CR res. We find excellent agreement with an earlier reconstruction by Sevenster (1999).

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.

We have developed an orbit-based method for constructing triaxial models of elliptical galaxies, which fit their observed surface brightness and kinematics. We have tested the method against analytical models with general distribution functions. Here, we present models that fit integral-field SAURON observations of NGC3379 and NGC4365.

We investigate with unprecedented accuracy the correlations between the dynamical mass-to-light ratio M/L and other global observables of E and S0 galaxies. We construct two-integral Jeans and three-integral Schwarzschild dynamical models for a sample of 25 E/S0 galaxies with SAURON integral-field stellar kinematics. We find a tight correlation of the form (M/L)=(3.80+/-0.14)*(sigma/200 km/s)^(0.84+/-0.07) between the dynamical M/L (in the I-band) and the luminosity-weighted second moment (sigma) of the line-of-sight velocity-distribution within Re. The observed rms scatter in M/L for our sample is 18%, while the inferred intrinsic scatter is ~13%. The (M/L)-sigma relation can be included in the remarkable series of tight correlations between sigma and other galaxy global observables. The comparison of the observed correlations with the predictions of the Fundamental Plane (FP), and with simple virial estimates, shows that the `tilt' of the FP of early-type galaxies, is due to a real M/L variation, while structural and orbital non-homology have a negligible effect. The virial mass is a reliable estimator of the mass in the central regions of galaxies. The best-fitting virial relation has the form (M/L)_vir=(5.0+/-0.1)*Re*sigma^2/(L*G). The comparison of the dynamical M/L with the (M/L)_pop inferred from the analysis of the stellar population, indicates a median dark matter fraction in early-type galaxies of ~30% of the total mass inside one Re. (Abridged)

We investigate the well-known correlations between the dynamical mass-to-light ratio M/L and other global observables of elliptical (E) and lenticular (S0) galaxies. We construct two-integral Jeans and three-integral Schwarzschild dynamical models for a sample of 25 E/S0 galaxies with SAURON integral-field stellar kinematics to about one effective (half-light) radius R_e. The comparison of the dynamical M/L with the (M/L)_pop inferred from the analysis of the stellar population, indicates that dark matter in early-type galaxies contributes ~30% of the total mass inside one R_e, in agreement with previous studies, with significant variations from galaxy to galaxy. Our results suggest a variation in M/L at constant (M/L)_pop, which seems to be linked to the galaxy dynamics. We speculate that fast rotating galaxies have lower dark matter fractions than the slow rotating and generally more massive ones.

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.

We use integral-field observation of the stellar kinematics obtained with SAURON in combination with Schwarzschild dynamical models to revisit our understanding of the classic (V/sigma)-epsilon anisotropy diagram of early-type galaxies.

We have used the ARGUS multi-object spectrometer at the CTIO 4m Blanco telescope to obtain 2756 radial velocity measurements for 1966 individual stars in the globular cluster omega Centauri brighter than blue photographic magnitude of about 16.5. Of these, 1589 stars are cluster members. A comparison with two independent radial velocity studies, carried out by Suntzeff & Kraft and by Mayor et al., demonstrates that the median error of our measurements is below 2 km/s for the stars brighter than B-magnitude 15, which constitute the bulk of the sample. The observed velocity dispersion decreases from about 15 km/s in the inner few arcmin to about 6 km/s at a radius of 25 arcmin. The cluster shows significant rotation, with a maximum amplitude of about 6 km/s in the radial zone between 6 and 10 arcmin. In a companion paper by van de Ven et al., we correct these radial velocities for the perspective rotation caused by the space motion of the cluster, and combine them with the internal proper motions of nearly 8000 cluster members measured by van Leeuwen et al., to construct a detailed dynamical model of omega Centauri and to measure its distance.

We determine the dynamical distance D, inclination i, mass-to-light ratio M/L and the intrinsic orbital structure of the globular cluster omega Cen, by fitting axisymmetric dynamical models to the ground-based proper motions of van Leeuwen et al. and line-of-sight velocities from four independent data-sets. We correct the observed velocities for perspective rotation caused by the space motion of the cluster, and show that the residual solid-body rotation component in the proper motions can be taken out without any modelling other than assuming axisymmetry. This also provides a tight constraint on D tan i. Application of our axisymmetric implementation of Schwarzschild's orbit superposition method to omega Cen reveals no dynamical evidence for a significant radial dependence of M/L. The best-fit dynamical model has a stellar V-band mass-to-light ratio M/L_V = 2.5 +/- 0.1 M_sun/L_sun and an inclination i = 50 +/- 4 degrees, which corresponds to an average intrinsic axial ratio of 0.78 +/- 0.03. The best-fit dynamical distance D = 4.8 +/- 0.3 kpc (distance modulus 13.75 +/- 0.13 mag) is significantly larger than obtained by means of simple spherical or constant-anisotropy axisymmetric dynamical models, and is consistent with the canonical value 5.0 +/- 0.2 kpc obtained by photometric methods. The total mass of the cluster is (2.5 +/- 0.3) x 10^6 M_sun. The best-fit model is close to isotropic inside a radius of about 10 arcmin and becomes increasingly tangentially anisotropic in the outer region, which displays significant mean rotation. This phase-space structure may well be caused by the effects of the tidal field of the Milky Way. The cluster contains a separate disk-like component in the radial range between 1 and 3 arcmin, contributing about 4% to the total mass.

We analyse SAURON kinematic maps of the inner kpc of the early-type (Sa) barred spiral galaxy NGC 5448. 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 displays 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, which is derived using linear theory. 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.

We have developed a method to construct realistic triaxial dynamical models for elliptical galaxies, allowing us to derive best-fitting parameters, such as the mass-to-light ratio and the black hole mass, and to study the orbital structure. We use triaxial theoretical Abel models to investigate the robustness of the method.

We analyse the Fundamental Plane (FP) of 26 strong gravitational lens galaxies with redshifts up to z~1, using tabulated data from Kochanek et al. (2000) and Rusin et al. (2003). The lens galaxies effectively form a mass-selected sample of early-type galaxies in environments of relatively low density. We analyse the FP and its evolution in the restframe Johnson B and Gunn r bands. Assuming that early-type galaxies are a homologous family, the FP then provides a direct measurement of the M/L ratio evolution. If we assume that the M/L ratios of field early-type galaxies evolve as power-laws, we find for the lens galaxies an evolution rate dlog(M/L)/dz of -0.62 +/-0.13 in restframe B, corresponding to a mean stellar formation redshift of <z*>=1.8 (-0.5/+1.4). This formation epoch does not differ significantly from that of cluster galaxies <z*cl>=2.0 (-0.2/+0.3), obtained by van Dokkum & Franx (2001) after correction for maximum progenitor bias. However, if we impose the constraint that lens and cluster galaxies that are of the same age have equal M/L ratios and we do not correct for progenitor bias, we find that the stellar populations of the lens galaxies are 10-15 % younger. We find that both the M/L ratios as well as the restframe colors of the lens galaxies show significant scatter. About half of the lens galaxies are consistent with an old cluster-like stellar populations, but the other galaxies are bluer and best fit by single burst models with younger stellar formation redshifts as low as z*~1. Moreover, the scatter in color is correlated with the scatter in M/L ratio. We interpret this as evidence of a significant age spread among the stellar populations of lens galaxies, whereas the ages of the stellar populations of the cluster galaxies are well approximated by a single formation epoch.

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 and six independent moments. By assuming that the potential is triaxial and of separable Staeckel form, the mixed moments vanish in confocal ellipsoidal coordinates. Consequently, 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. These equations were first derived by Lynden-Bell, over 40 years ago, but have resisted solution by standard methods. We present the general solution here. We consider the two-dimensional limiting cases first. We solve their Jeans equations by a new method which superposes singular solutions. The singular solutions, which are new, are standard Riemann-Green functions. The two-dimensional solutions are applied to non-axisymmetric discs, oblate and prolate spheroids, and also to the scale-free triaxial limit. We then extend the method of singular solutions to the triaxial case, and obtain a full solution. The general solution can be expressed in terms of complete (hyper)elliptic integrals which can be evaluated in a straightforward way, and provides the full set of second moments which can support a triaxial density distribution in a separable triaxial potential. (abridged)

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.

We have measured the two-dimensional kinematics and line-strength distributions of 72 representative nearby early-type galaxies, out to approximately one effective radius, with our panoramic integral-field spectrograph SAURON. The resulting maps reveal a rich variety in kinematical structures and linestrength distributions, indicating that early-type galaxies are more complex systems than often assumed. We are building detailed dynamical models for these galaxies, to derive their intrinsic shape and dynamical structure, and to determine the mass of the supermassive central black hole. Here we focus on two examples, the compact elliptical M32 and the E3 galaxy NGC4365. These objects represent two extreme cases: M32 has very regular kinematics which can be represented accurately by an axisymmetric model in which all stars rotate around the short axis, while NGC4365 is a triaxial galaxy with a prominent kinematically decoupled core, with an inner core that rotates about an axis that is nearly perpendicular to the rotation axis of the main body of the galaxy. Our dynamical models for these objects demonstrate that two-dimensional observations are essential for deriving the intrinsic orbital structure and dark matter content of galaxies.