In the last years, Astroinformatics has become a well defined paradigm for many fields of Astronomy. In this work we demonstrate the potential of a multidisciplinary approach to identify globular clusters (GCs) in the Fornax cluster of galaxies taking advantage of multi-band photometry produced by the VLT Survey Telescope using automatic self-adaptive methodologies. The data analyzed in this work consist of deep, multi-band, partially overlapping images centered on the core of the Fornax cluster. In this work we use a Neural-Gas model, a pure clustering machine learning methodology, to approach the GC detection, while a novel feature selection method ($\Phi$LAB) is exploited to perform the parameter space analysis and optimization. We demonstrate that the use of an Astroinformatics based methodology is able to provide GC samples that are comparable, in terms of purity and completeness with those obtained using single band HST data (Brescia et al. 2012) and two approaches based respectively on a morpho-photometric (Cantiello et al. 2018b) and a PCA analysis (D'Abrusco et al. 2015) using the same data discussed in this work.

The formation of Low mass X-ray binaries (LMXB) is favored within dense stellar systems such as Globular Clusters (GCs). The connection between LMXB and Globular Clusters has been extensively studied in the literature, but these studies have always been restricted to the innermost regions of galaxies. We present a study of LMXB in GCs within the central 1.5 deg^2 of the Fornax cluster with the aim of confirming the existence of a population of LMXB in intra-cluster GCs and understand if their properties are related to the host GCs, to the environment or/and to different formation channels.

We present results from an ongoing investigation using the Kilo Degree Survey (KiDS) on the VLT Survey Telescope (VST) to provide a census of ultra-compact massive galaxies (UCMGs), defined as galaxies with stellar masses $M_{\rm \star} > 8 \times 10^{10} \rm M_{\odot}$ and effective radii $R_{\rm e} < 1.5\,\rm kpc$. UCMGs, which are expected to have undergone very few merger events, provide a unique view on the accretion history of the most massive galaxies in the Universe. Over an effective sky area of nearly 330 square degrees, we select UCMG candidates from KiDS multi-colour images, which provide high quality structural parameters, photometric redshifts and stellar masses. Our sample of $\sim 1000$ photometrically selected UCMGs at $z < 0.5$ represents the largest sample of UCMG candidates assembled to date over the largest sky area. In this paper we present the first effort to obtain their redshifts using different facilities, starting with first results for 28 candidates with redshifts $z < 0.5$, obtained at NTT and TNG telescopes. We confirmed, as bona fide UCMGs, 19 out of the 28 candidates with new redshifts. A further 46 UCMG candidates are confirmed with literature spectroscopic redshifts (35 at $z < 0.5$), bringing the final cumulative sample of spectroscopically-confirmed lower-z UCMGs to 54 galaxies, which is the largest sample at redshifts below $0.5$. We use these spectroscopic redshifts to quantify systematic errors in our photometric selection, and use these to correct our UCMG number counts. We finally compare the results to independent datasets and simulations.

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.

Here we present positions and radial velocities for over 4000 globular clusters (GCs) in 27 nearby early-type galaxies from the SLUGGS survey. The SLUGGS survey is designed to be representative of elliptical and lenticular galaxies in the stellar mass range 10 $<$ log M$_{\ast}$/M$_{\odot}$ $<$ 11.7. The data have been obtained over many years, mostly using the very stable multi-object spectrograph DEIMOS on the Keck II 10m telescope. Radial velocities are measured using the calcium triplet lines with a velocity accuracy of $\pm$ 10-15 km/s. We use phase space diagrams (i.e. velocity--position diagrams) to identify contaminants such as foreground stars and background galaxies, and to show that the contribution of GCs from neighboring galaxies is generally insignificant. Likely ultra-compact dwarfs are tabulated separately. We find that the mean velocity of the GC system is close to that of the host galaxy systemic velocity, indicating that the GC system is in overall dynamical equilibrium within the galaxy potential. We also find that the GC system velocity dispersion scales with host galaxy stellar mass in a similar manner to the Faber-Jackson relation for the stellar velocity dispersion. Publication of these GC radial velocity catalogs should enable further studies in many areas, such as GC system substructure, kinematics, and host galaxy mass measurements.

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.

Integral field unit spectrographs allow the 2D exploration of the kinematics and stellar populations of galaxies, although they are generally restricted to small fields-of-view. Using the large field-of-view of the DEIMOS multislit spectrograph on Keck and our Stellar Kinematics using Multiple Slits (SKiMS) technique, we are able to extract sky-subtracted stellar light spectra to large galactocentric radii. Here we present a new DEIMOS mask design named SuperSKiMS that explores large spatial scales without sacrificing high spatial sampling. We simulate a set of observations with such a mask design on the nearby galaxy NGC 1023, showing that the kinematic and metallicity measurements can reach radii where the galaxy surface brightness is several orders of magnitude fainter than the sky. Such a technique is also able to reproduce the kinematic and metallicity 2D distributions obtained from literature integral field spectroscopy in the innermost galaxy regions. In particular, we use the simulated NGC 1023 kinematics to model its total mass distribution to large radii, obtaining comparable results with those from published integral field unit observation. Finally, from new spectra of NGC 1023 we obtain stellar 2D kinematics and metallicity distributions that show good agreement with integral field spectroscopy results in the overlapping regions. In particular, we do not find a significant offset between our SKiMS and the ATLAS$^\rm{3D}$ stellar velocity dispersion at the same spatial locations.

We study mass distributions within and beyond 5~effective radii ($R_{\rm e}$) in 23 early-type galaxies from the SLUGGS survey, using their globular cluster (GC) kinematic data. The data are obtained with Keck/DEIMOS spectrograph, and consist of line-of-sight velocities for ~$3500$ GCs, measured with a high precision of ~15 $\rm km\ s^{-1}$ per GC and extending out to $~13 R_{\rm e}$. We obtain the mass distribution in each galaxy using the tracer mass estimator of Watkins et al. and account for kinematic substructures, rotation of the GC systems and galaxy flattening in our mass estimates. The observed scatter between our mass estimates and results from the literature is less than 0.2 dex. The dark matter fraction within $5R_{\rm e}$ ($f_{\rm DM}$) increases from ~$0.6$ to ~$0.8$ for low- and high-mass galaxies, respectively, with some intermediate-mass galaxies ($M_*{\sim}10^{11}M_\odot$) having low $f_{\rm DM}\sim0.3$, which appears at odds with predictions from simple galaxy models. We show that these results are independent of the adopted orbital anisotropy, stellar mass-to-light ratio, and the assumed slope of the gravitational potential. However, the low $f_{\rm DM}$ in the ~$10^{11}M_\odot$ galaxies agrees with the cosmological simulations of Wu et al. where the pristine dark matter distribution has been modified by baryons during the galaxy assembly process. We find hints that these $M_*\sim10^{11}M_\odot$ galaxies with low $f_{\rm DM}$ have very diffuse dark matter haloes, implying that they assembled late. Beyond $5R_{\rm e}$, the $M/L$ gradients are steeper in the more massive galaxies and shallower in both low and intermediate mass 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.

Ultra diffuse galaxies (UDGs) have the sizes of giant galaxies but the luminosities of dwarfs. A key to understanding their origins comes from their total masses, but their low surface brightnesses ($\mu(V) \geq$ 25.0) generally prohibit dynamical studies. Here we report the first such measurements for a UDG (VCC~1287 in the Virgo cluster), based on its globular cluster system dynamics and size. From 7 GCs we measure a mean systemic velocity $v_{\rm sys}$ = 1071$^{+14}_{-15}$ km/s, thereby confirming a Virgo-cluster association. We measure a velocity dispersion of 33$^{+16}_{-10}$ km/s within 8.1 kpc, corresponding to an enclosed mass of $(4.5 \pm 2.8)\times10^{9}$ $M_{\odot}$ and a $g$-band mass-to-light ratio of $(M/L)_g = 106^{+126}_{-54}$. From the cumulative mass curve, along with the GC numbers, we estimate a virial mass of $\sim8\times10^{10}$ $M_{\odot}$, yielding a dark-to-stellar mass fraction of $\sim3000$. We show that this UDG is an outlier in $M_{\rm star} - M_{\rm halo}$ relations, suggesting extreme stochasticity in relatively massive star-forming halos in clusters. Finally, we discuss how counting GCs offers an efficient route to determining virial masses for UDGs.

[Abrigded] We have started a new deep, multi-imaging survey of the Fornax cluster, dubbed Fornax Deep Survey (FDS), at the VLT Survey Telescope. In this paper we present the deep photometry inside two square degrees around the bright galaxy NGC1399 in the core of the cluster. We found a very extended and diffuse envelope surrounding the luminous galaxy NGC1399: we map the surface brightness out to 33 arcmin (~ 192 kpc) from the galaxy center and down to about 31 mag/arcsec^2 in the g band. The deep photometry allows us to detect a faint stellar bridge in the intracluster region between NGC1399 and NGC1387. By analyzing the integrated colors of this feature, we argue that it could be due to the ongoing interaction between the two galaxies, where the outer envelope of NGC1387 on its east side is stripped away. By fitting the light profile, we found that it exists a physical break radius in the total light distribution at R=10 arcmin (~58 kpc) that sets the transition region between the bright central galaxy and the outer exponential stellar halo. We discuss the main implications of this work on the build-up of the stellar halo at the center of the Fornax cluster. By comparing with the numerical simulations of the stellar halo formation for the most massive BCGs, we find that the observed stellar halo mass fraction is consistent with a halo formed through the multiple accretion of progenitors with a stellar mass in the range 10^8 - 10^11 M_sun. This might suggest that the halo of NGC1399 has also gone through a major merging event. The absence of a significant number of luminous stellar streams and tidal tails out to 192 kpc suggests that the epoch of this strong interaction goes back to an early formation epoch. Therefore, differently from the Virgo cluster, the extended stellar halo around NGC1399 is characterised by a more diffuse and well-mixed component, including the ICL.

A strong correlation exists between the total mass of a globular cluster (GC) system and the virial halo mass of the host galaxy. However, the total halo mass in this correlation is a statistical measure conducted on spatial scales that are some ten times that of a typical GC system. Here we investigate the connection between GC systems and galaxy's dark matter on comparable spatial scales, using dynamical masses measured on a galaxy-by-galaxy basis. Our sample consists of 17 well-studied massive (stellar mass $\sim$10$^{11}$ M$_{\odot}$) early-type galaxies from the SLUGGS survey. We find the strongest correlation to be that of the blue (metal-poor) GC subpopulation and the dark matter content. This correlation implies that the dark matter mass of a galaxy can be estimated to within a factor of two from careful imaging of its GC system. The ratio of the GC system mass to that of the enclosed dark matter is nearly constant. We also find a strong correlation between the fraction of blue GCs and the fraction of enclosed dark matter, so that a typical galaxy with a blue GC fraction of 60 per cent has a dark matter fraction of 86 per cent over similar spatial scales. Both halo growth and removal (via tidal stripping) may play some role in shaping this trend. In the context of the two-phase model for galaxy formation, we find galaxies with the highest fractions of accreted stars to have higher dark matter fractions for a given fraction of blue GCs.

Early-type (E and S0) galaxies may have assembled via a variety of different evolutionary pathways. Here we investigate these pathways by comparing the stellar kinematic properties of 24 early-type galaxies from the SLUGGS survey with the hydrodynamical simulations of Naab et al. (2014). In particular, we use the kinematics of starlight up to 4 effective radii (R$_e$) as diagnostics of galaxy inner and outer regions, and assign each galaxy to one of six Naab et al. assembly classes. The majority of our galaxies (14/24) have kinematic characteristics that indicate an assembly history dominated by gradual gas dissipation and accretion of many gas-rich minor mergers. Three galaxies, all S0s, indicate that they have experienced gas-rich major mergers in their more recent past. One additional elliptical galaxy is tentatively associated with a gas-rich merger which results in a remnant galaxy with low angular momentum. Pathways dominated by gas-poor (major or minor) mergers dominate the mass growth of six galaxies. Most SLUGGS galaxies appear to have grown in mass (and size) via the accretion of stars and gas from minor mergers, with late major mergers playing a much smaller role. We find that the fraction of accreted stars correlates with the stellar mean age and metallicity gradient, but not with the slope of the total mass density profile. We briefly mention future observational and modelling approaches that will enhance our ability to accurately reconstruct the assembly histories of individual present day galaxies.

Due to longer dynamical timescales, the outskirts of early-type galaxies retain the footprint of their formation and assembly. Under the popular two-phase galaxy formation scenario, an initial in-situ phase of star formation is followed by minor merging and accretion of ex-situ stars leading to the expectation of observable transitions in the kinematics and stellar populations on large scales. However, observing the faint galactic outskirts is challenging, often leaving the transition unexplored. The large scale, spatially-resolved stellar kinematic data from the SAGES Legacy Unifying Galaxies and GlobularS (SLUGGS) survey are ideal for detecting kinematic transitions. We present kinematic maps out to 2.6 effective radii on average, kinemetry profiles, measurement of kinematic twists and misalignments, and the average outer intrinsic shape of 25 SLUGGS galaxies. We find good overall agreement in the kinematic maps and kinemetry radial profiles with literature. We are able to confirm significant radial modulations in rotational versus pressure support of galaxies with radius so that the central and outer rotational properties may be quite different. We also test the suggestion that galaxies may be more triaxial in their outskirts and find that while fast rotating galaxies were already shown to be axisymmetric in their inner regions, we are unable to rule out triaxiality in their outskirts. We compare our derived outer kinematic information to model predictions from a two-phase galaxy formation scenario. We find that the theoretical range of local outer angular momentum agrees well with our observations, but that radial modulations are much smaller than predicted.

Globular clusters (GCs) can be considered discrete, long-lived, dynamical tracers that retain crucial information about the assembly history of their parent galaxy. In this paper, we present a new catalogue of GC velocities and colours for the lenticular galaxy NGC 1023, we study their kinematics and spatial distribution, in comparison with the underlying stellar kinematics and surface brightness profile, and we test a new method for studying GC properties. Specifically, we decompose the galaxy light into its spheroid (assumed to represent the bulge + halo components) and disk components and use it to assign to each GC a probability of belonging to one of the two components. Then we model the galaxy kinematics, assuming a disk and spheroidal component, using planetary nebulae (PNe) and integrated stellar light. We use this kinematic model and the probability previously obtained from the photometry to recalculate for each GC its likelihood of being associated with the disk, the spheroid, or neither. We find that the reddest GCs are likely to be associated with the disk, as found for faint fuzzies in this same galaxy, suggesting that the disk of this S0 galaxy originated at z ~ 2. The majority of blue GCs are found likely to be associated with the spheroidal (hot) component. The method also allows us to identify objects that are unlikely to be in equilibrium with the system. In NGC1023 some of the rejected GCs form a substructure in phase space that is connected with NGC 1023 companion galaxy.

We present an integral field study of the internal structure, kinematics and stellar population of the almost edge-on, intermediate luminosity ($L_ {*}$) elliptical galaxy NGC 4697. We build extended 2-dimensional (2D) maps of the stellar kinematics and line-strengths of the galaxy up to $\sim 0.7 $ effective radii (R$_{eff}$) using a mosaic of 8 VIMOS (VIsible Multi-Objects Spectrograph on the VLT) integral-field unit pointings. We find clear evidence for a rotation-supported structure along the major axis from the 2D kinematical maps, confirming the previous classification of this system as a `fast-rotator'. We study the correlations between the third and fourth Gauss-Hermite moments of the line-of-sight velocity distribution (LOSVD) $h_3$ and $h_4$ with the rotation parameter ($V/\sigma$), and compare our findings to hydrodynamical simulations. We find remarkable similarities to predictions from gas-rich mergers. Based on photometry, we perform a bulge/disk decomposition and study the stellar population properties of the two components. The bulge and the disk show different stellar populations, with the stars in the bulge being older (age$_{\rm bulge}=13.5^{+1.4}_{-1.4}$ Gyr, age$_{\rm disk}=10.5^{+1.6}_{-2.0}$Gyr) and more metal-poor ($\mathrm{[M/H]_{bulge}} = -0.17^{+0.12}_{-0.1}$, $\mathrm{[M/H]_{disk}}=-0.03^{+0.02}_{-0.1}$). The evidence of a later-formed, more metal-rich disk embedded in an older, more metal-poor bulge, together with the LOSVD structure, supports a mass assembly scenario dominated by gas-rich minor mergers and possibly with a late gas-rich major merger that left a previously rapidly rotating system unchanged. The bulge and the disk do not show signs of different stellar Initial Mass Function slopes, and both match well with a Milky Way-like IMF.

The outer halo regions of early-type galaxies carry key information about their past accretion history. However, spectroscopically probing the stellar component at such galactocentric radii is still challenging. Using Keck/DEIMOS, we have been able to measure the metallicities of the stellar and globular cluster components in 12 early-type galaxies out to more than $10~\rm{R_{e}}$. We find similar metallicity gradients for the metal-poor and metal-rich globular cluster subpopulations, suggesting a common formation process for the two subpopulations. This is in conflict with most current theoretical predictions, where the metal-poor globular clusters are thought to be purely accreted and metal-rich globular clusters mostly formed in-situ. Moreover, we find that the globular cluster metallicity gradients show a trend with galaxy mass, being steeper in lower-mass galaxies than in higher-mass galaxies. This is similar to what we find for the outermost galaxy stars and suggests a more active accretion history, with a larger role played by major mergers, in the most massive galaxies. This conclusion is qualitatively consistent with expectations from two-phase galaxy assembly models. Finally, we find that the small difference in metallicity between galaxy stars and metal-rich globular clusters at $1~\rm{R_{e}}$ may correlate with galaxy mass. The origin of this difference is not currently clear.

We perform in-depth dynamical modelling of the luminous and dark matter (DM) content of the elliptical galaxy NGC 1407. Our strategy consists of solving the spherical Jeans equations for three independent dynamical tracers: stars, blue GCs and red GCs in a self-consistent manner. We adopt a maximum-likelihood Markov-Chain Monte Carlo fitting technique in the attempt to constrain the inner slope of the DM density profile (the cusp/core problem), and the stellar initial mass function (IMF) of the galaxy. We find the inner logarithmic slope of the DM density profiles to be $\gamma = 0.6\pm0.4$, which is consistent with either a DM cusp ($\gamma = 1$) or with a DM core $(\gamma = 0)$. Our findings are consistent with a Salpeter IMF, and marginally consistent with a Kroupa IMF. We infer tangential orbits for the blue GCs, and radial anisotropy for red GCs and stars. The modelling results are consistent with the virial mass--concentration relation predicted by $\Lambda$CDM simulations. The virial mass of NGC 1407 is $\log$ $M_{\rm vir} = 13.3 \pm 0.2 M_{\odot}$, whereas the stellar mass is $\log M_* = 11.8 \pm 0.1 M_{\odot}$. The overall uncertainties on the mass of NGC 1407 are only 5 per cent at the projected stellar effective radius. We attribute the disagreement between our results and previous X-ray results to the gas not being in hydrostatic equilibrium in the central regions of the galaxy. The halo of NGC 1407 is found be DM dominated, with a dynamical mass-to-light ratio of $M/L=260_{-100} ^{+174} M_{\odot}/L_{\odot, B}$. However, this value can be larger up to a factor of 3 depending on the assumed prior on the DM scale radius.

We study the total mass-density profile for a sample of 14 fast-rotator early-type galaxies (stellar masses $10.2<\log M_\ast/M_\odot<11.7$). We combine observations from the SLUGGS and Atlas3D surveys to map out the stellar kinematics in two-dimensions, out to a median radius for the sample of four half-light radii $R_e$ (or 10 kpc), and a maximum radius of 2.0-6.2 $R_e$ (or 4-21 kpc). We use axisymmetric dynamical models based on the Jeans equations, which allow for a spatially varying anisotropy, and employ quite general profiles for the dark halos, and in particular do not place any restriction on the profile slope. This is made possible by the availability of spatially extended two-dimensional kinematics. We find that our relatively simple models provide a remarkably good description of the observed kinematics. The resulting total density profiles are well described by a nearly-isothermal power law $\rho_{\rm tot}(r)\propto r^{-\gamma}$ from $R_e$/10 to at least 4$R_e$, the largest average deviation being 11%. The average logarithmic slope is $\langle\gamma\rangle=2.19\pm0.03$ with observed rms scatter of just $\sigma_\gamma=0.11$. This scatter out to large radii, where dark matter dominates, is as small as previously reported by lensing studies around $r\approx R_e/2$, where the stars dominate. Our bulge-halo conspiracy places much tighter constraints on galaxy formation models. It illustrates the power of two-dimensional stellar kinematics observations at large radii. It would now be important to test the generality of our results for different galaxy types and larger samples.

We present new wide-field photometry and spectroscopy of the globular clusters (GCs) around NGC 4649 (M60), the third brightest galaxy in the Virgo cluster. Imaging of NGC 4649 was assembled from a recently-obtained HST/ACS mosaic, and new Subaru/Suprime-Cam and archival CFHT/MegaCam data. About 1200 sources were followed up spectroscopically using combined observations from three multi-object spectrographs: Keck/DEIMOS, Gemini/GMOS and MMT/Hectospec. We confirm 431 unique GCs belonging to NGC 4649, a factor of 3.5 larger than previous datasets and with a factor of 3 improvement in velocity precision. We confirm significant GC colour bimodality and find that the red GCs are more centrally concentrated, while the blue GCs are more spatially extended. We infer negative GC colour gradients in the innermost 20 kpc and flat gradients out to large radii. Rotation is detected along the galaxy major axis for all tracers: blue GCs, red GCs, galaxy stars and planetary nebulae. We compare the observed properties of NGC 4649 with galaxy formation models. We find that formation via a major merger between two gas-poor galaxies, followed by satellite accretion, can consistently reproduce the observations of NGC 4649 at different radii. We find no strong evidence to support an interaction between NGC 4649 and the neighbouring spiral galaxy NGC 4647. We identify interesting GC kinematic features in our data, such as counter-rotating subgroups and bumpy kinematic profiles, which encode more clues about the formation history of NGC 4649.

As part of the SLUGGS survey, we stack 1137 Keck DEIMOS spectra of globular clusters from 10 galaxies to study their stellar populations in detail. The stacked spectra have median signal to noise ratios of $\sim 90$ \AA$^{-1}$. Besides the calcium triplet, we study weaker sodium, magnesium, titanium and iron lines as well as the H$\alpha$ and higher order Paschen hydrogen lines. In general, the stacked spectra are consistent with old ages and a Milky Way-like initial mass function. However, we see different metal line index strengths at fixed colour and magnitude, and differences in the calcium triplet--colour relation from galaxy to galaxy. We interpret this as strong evidence for variations in the globular cluster colour--metallicity relation between galaxies. Two possible explanations for the colour--metallicity relation variations are that the average ages of globular clusters vary from galaxy to galaxy or that the average abundances of light elements (i.e. He, C, N and O) differ between galaxies. Stacking spectra by magnitude, we see that the colours become redder and metal line indices stronger with brighter magnitudes. These trends are consistent with the previously reported `blue tilts' being mass--metallicity relations.

We introduce and provide the scientific motivation for a wide-field photometric and spectroscopic chemodynamical survey of nearby early-type galaxies (ETGs) and their globular cluster (GC) systems. The SLUGGS (SAGES Legacy Unifying Globulars and GalaxieS) survey is being carried out primarily with Subaru/Suprime-Cam and Keck/DEIMOS. The former provides deep gri imaging over a 900 arcmin^2 field-of-view to characterize GC and host galaxy colors and spatial distributions, and to identify spectroscopic targets. The NIR CaII triplet provides GC line-of-sight velocities and metallicities out to typically ~8 R_e, and to ~15 R_e in some cases. New techniques to extract integrated stellar kinematics and metallicities to large radii (~2-3 R_e) are used in concert with GC data to create two-dimensional velocity and metallicity maps for comparison with simulations of galaxy formation. The advantages of SLUGGS compared with other, complementary, 2D-chemodynamical surveys are its superior velocity resolution, radial extent, and multiple halo tracers. We describe the sample of 25 nearby ETGs, the selection criteria for galaxies and GCs, the observing strategies, the data reduction techniques, and modeling methods. The survey observations are nearly complete and more than 30 papers have so far been published using SLUGGS data. Here we summarize some initial results, including signatures of two-phase galaxy assembly, evidence for GC metallicity bimodality, and a novel framework for the formation of extended star clusters and ultracompact dwarfs. An integrated overview of current chemodynamical constraints on GC systems points to separate, in-situ formation modes at high redshifts for metal-poor and metal-rich GCs.

We present stellar kinematics of 22 nearby early-type galaxies (ETGs), based on two-dimensional (2D) absorption line stellar spectroscopy out to ~2-4 R_e (effective radii), as part of the ongoing SLUGGS Survey. The galaxies span a factor of 20 in intrinsic luminosity, as well as a full range of environment and ETG morphology. Our data consist of good velocity resolution (sigma_inst ~ 25 km/s) integrated stellar-light spectra extracted from the individual slitlets of custom made Keck/DEIMOS slitmasks. We extract stellar kinematics measurements (V, sigma, h_3, and h_4) for each galaxy. Combining with literature values from smaller radii, we present 2D spatially resolved maps of the large-scale kinematic structure in each galaxy. We find that the kinematic homogeneity found inside 1 R_e often breaks down at larger radii, where a variety of kinematic behaviors are observed. While central slow rotators remain slowly rotating in their halos, central fast rotators show more diversity, ranging from rapidly increasing to rapidly declining specific angular momentum profiles in the outer regions. There are indications that the outer trends depend on morphological type, raising questions about the proposed unification of the elliptical and lenticular (S0) galaxy families in the ATLAS^3D survey. Several galaxies in our sample show multiple lines of evidence for distinct disk components embedded in more slowly rotating spheroids, and we suggest a joint photometric-kinematic approach for robust bulge-disk decomposition. Our observational results appear generally consistent with a picture of two-phase (in-situ plus accretion) galaxy formation.

The previously clear division between small galaxies and massive star clusters is now occupied by objects called ultra compact dwarfs (UCDs) and compact ellipticals (cEs). Here we combine a sample of UCDs and cEs with velocity dispersions from the AIMSS project with literature data to explore their dynamical-to-stellar mass ratios. We confirm that the mass ratios of many UCDs in the stellar mass range 10$^6$ -- 10$^9$ M$_{\odot}$ are systematically higher than those for globular clusters which have mass ratios near unity. However, at the very highest masses in our sample, i.e. 10$^9$ -- 10$^{10}$ M$_{\odot}$, we find that cE galaxies also have mass ratios of close to unity, indicating their central regions are mostly composed of stars. Suggested explanations for the elevated mass ratios of UCDs have included a variable IMF, a central black hole, and the presence of dark matter. Here we present another possible explanation, i.e. tidal stripping. Under various assumptions, we find that the apparent variation in the mass ratio with stellar mass and stellar density can be qualitatively reproduced by published tidal stripping simulations of a dwarf elliptical galaxy. In the early stages of the stripping process the galaxy is unlikely to be in virial equilibrium. At late stages, the final remnant resembles the properties of $\sim$10$^7$ M$_{\odot}$ UCDs. Finally, we discuss the need for more detailed realistic modelling of tidal stripping over a wider range of parameter space, and observations to further test the stripping hypothesis.

We describe the structural and kinematic properties of the first compact stellar systems discovered by the AIMSS project. These spectroscopically confirmed objects have sizes ($\sim$6 $<$ R$_{\rm e}$ [pc] $<$ 500) and masses ($\sim$2$\times$10$^{6}$ $<$ M$_*$/M$_\odot$ $<$ 6$\times$10$^{9}$) spanning the range of massive globular clusters (GCs), ultra compact dwarfs (UCDs) and compact elliptical galaxies (cEs), completely filling the gap between star clusters and galaxies. Several objects are close analogues to the prototypical cE, M32. These objects, which are more massive than previously discovered UCDs of the same size, further call into question the existence of a tight mass-size trend for compact stellar systems, while simultaneously strengthening the case for a universal "zone of avoidance" for dynamically hot stellar systems in the mass-size plane. Overall, we argue that there are two classes of compact stellar systems: 1) massive star clusters and 2) a population closely related to galaxies. Our data provide indications for a further division of the galaxy-type UCD/cE population into two groups, one population that we associate with objects formed by the stripping of nucleated dwarf galaxies, and a second population that formed through the stripping of bulged galaxies or are lower-mass analogues of classical ellipticals. We find compact stellar systems around galaxies in low to high density environments, demonstrating that the physical processes responsible for forming them do not only operate in the densest clusters.

Faint fuzzies are a relatively new class of star cluster, first found in the nearby S0 galaxy NGC 1023 by Larsen & Brodie using WFPC2 images from the Hubble Space Telescope (HST). Here we investigate the star cluster system of NGC 1023 using an eight pointing mosaic of ACS images from HST. We identify blue and red normal globular clusters (two of which are particularly luminous and dense) and two ultra compact dwarf objects (with effective radius ~10 pc). With our more complete spatial coverage, we also find 81 red and 27 blue faint fuzzies (FFs). We confirm the association of the red FFs with the disk of NGC 1023, consistent with them being long-lived open clusters. Half of the blue FFs appear to be associated with the dwarf satellite galaxy NGC 1023A (which was largely absent from the original HST/WFPC2 coverage), while the remainder are spatially coincident with the densest HI gas that surrounds NGC 1023. The blue FFs have colours that are consistent with young (few 100 Myr old) star clusters that formed during the most recent interaction between NGC 1023 and its satellite NGC 1023A.

We study the mass and anisotropy distribution of the giant elliptical galaxy NGC 5846 using stars, as well as the red and blue globular cluster (GC) subpopulations. We break degeneracies in the dynamical models by taking advantage of the different phase space distributions of the two GC subpopulations to unambiguously constrain the mass of the galaxy and the anisotropy of the GC system. Red GCs show the same spatial distribution and behaviour as the starlight, whereas blue GCs have a shallower density profile, a larger velocity dispersion and a lower kurtosis, all of which suggest a different orbital distribution. We use a dispersion-kurtosis Jeans analysis and find that the solutions of separate analyses for the two GC subpopulations overlap in the halo parameter space. The solution converges on a massive dark matter halo, consistent with expectations from $\Lambda$CDM and WMAP7 cosmology in terms of virial mass ($\log M_{DM} \sim13.3 M_{sun}$) and concentration ($c_{vir}\sim8$). This is the first such analysis that solves the dynamics of the different GC subpopulations in a self-consistent manner. Our method improves the uncertainties on the halo parameter determination by a factor of two and opens new avenues for the use of elliptical galaxy dynamics as tests of predictions from cosmological simulations. The implied stellar mass-to-light ratio derived from the dynamical modelling is fully consistent with a Salpeter initial mass function (IMF) and rules out a bottom light IMF. The different GC subpopulations show markedly distinct orbital distributions at large radii, with red GCs having an anisotropy parameter $\beta\sim0.4$ outside $\sim3R_e$, and the blue GCs having $\beta\sim0.15$ at the same radii, while centrally ($\sim1R_e$) they are both isotropic. We discuss the implications of our findings within the two-phase formation scenario for early-type galaxies.

We report the discovery of a remarkable ultra-compact dwarf galaxy around the massive Virgo elliptical galaxy NGC 4649 (M60), which we term M60-UCD1. With a dynamical mass of 2.0 x 10^8 M_sun but a half-light radius of only ~ 24 pc, M60-UCD1 is more massive than any ultra-compact dwarfs of comparable size, and is arguably the densest galaxy known in the local universe. It has a two-component structure well-fit by a sum of Sersic functions, with an elliptical, compact (r_h=14 pc; n ~ 3.3) inner component and a round, exponential, extended (r_h=49 pc) outer component. Chandra data reveal a variable central X-ray source with L_X ~ 10^38 erg/s that could be an active galactic nucleus associated with a massive black hole or a low-mass X-ray binary. Analysis of optical spectroscopy shows the object to be old (~> 10 Gyr) and of solar metallicity, with elevated [Mg/Fe] and strongly enhanced [N/Fe] that indicates light element self-enrichment; such self-enrichment may be generically present in dense stellar systems. The velocity dispersion (~ 70 km/s) and resulting dynamical mass-to-light ratio (M/L_V=4.9 +/- 0.7) are consistent with---but slightly higher than---expectations for an old, metal-rich stellar population with a Kroupa initial mass function. The presence of a massive black hole or a mild increase in low-mass stars or stellar remnants is therefore also consistent with this M/L_V. The stellar density of the galaxy is so high that no dynamical signature of dark matter is expected. However, the properties of M60-UCD1 suggest an origin in the tidal stripping of a nucleated galaxy with M_B ~ -18 to -19.

It is not understood whether long-lived star clusters possess a continuous range of sizes and masses (and hence densities), or if rather, they should be considered as distinct types with different origins. Utilizing the Hubble Space Telescope (HST) to measure sizes, and long exposures on the Keck 10m telescope to obtain distances, we have discovered the first confirmed star clusters that lie within a previously claimed size-luminosity gap dubbed the `avoidance zone' by Hwang et al (2011). The existence of these star clusters extends the range of sizes, masses and densities for star clusters, and argues against current formation models that predict well-defined size-mass relationships (such as stripped nuclei, giant globular clusters or merged star clusters). The red colours of these gap objects suggests that they are not a new class of object but are related to Faint Fuzzies observed in nearby lenticular galaxies. We also report a number of low luminosity UCDs with sizes of up to 50 pc. Future, statistically complete, studies will be encouraged now that it is known that star clusters possess a continuous range of structural properties.

Understanding whether the bulge or the halo provides the primary link to the growth of supermassive black holes has strong implications for galaxy evolution and supermassive black hole formation itself. In this paper, we approach this issue by investigating extragalactic globular cluster (GC) systems, which can be used to probe the physics of both the bulge and the halo of the host galaxy. We study the relation between the supermassive black hole masses M_BH and the globular cluster system velocity dispersions sigma_GC using an updated and improved sample of 21 galaxies. We exploit the dichotomy of globular cluster system colours, to test if the blue and red globular clusters correlate differently with black hole mass. This may be expected if they trace the potentially different formation history of the halo and of the bulge of the host galaxy respectively. We find that M_BH correlates with the total GC system velocity dispersion, although not as strongly as claimed by recent work of Sadoun & Colin. We also examine the M_BH - sigma_GC relation for barred/bar-less and core/non-core galaxies, finding no significant difference, and for the first time we quantify the impact of radial gradients in the GC system velocity dispersion profile on the M_BH - sigma_GC relation. We additionally predict M_BH in 13 galaxies, including dwarf elliptical galaxies and the cD galaxy NGC 3311. We conclude that our current results cannot discriminate between the bulge/halo scenario. Although there is a hint that the red GC velocity dispersion might correlate better with M_BH than the blue GC velocity dispersion, the number statistics are still too low to be certain.

The origins of S0 galaxies remain obscure, with various mechanisms proposed for their formation, likely depending on environment. These mechanisms would imprint different signatures in the galaxies' stellar kinematics out to large radii, offering a method for distinguishing between them. We aim to study a sample of six S0 galaxies from a range of environments, and use planetary nebulae (PNe) as tracers of their stellar populations out to very large radii, to determine their kinematics in order to understand their origins. Using a special-purpose instrument, the Planetary Nebula Spectrograph, we observe and extract PNe catalogues for these six systems*. We show that the PNe have the same spatial distribution as the starlight, that the numbers of them are consistent with what would be expected in a comparable old stellar population in elliptical galaxies, and that their kinematics join smoothly onto those derived at smaller radii from conventional spectroscopy. The high-quality kinematic observations presented here form an excellent set for studying the detailed kinematics of S0 galaxies, in order to unravel their formation histories. We find that PNe are good tracers of stellar kinematics in these systems. We show that the recovered kinematics are largely dominated by rotational motion, although with significant random velocities in most cases.

We present a spectro-photometric survey of 2522 extragalactic globular clusters (GCs) around twelve early-type galaxies, nine of which have not been published previously. Combining space-based and multi-colour wide field ground-based imaging, with spectra from the Keck DEIMOS instrument, we obtain an average of 160 GC radial velocities per galaxy, with a high velocity precision of 15 km/s per GC. After studying the photometric properties of the GC systems, such as their spatial and colour distributions, we focus on the kinematics of metal-poor (blue) and metal-rich (red) GC subpopulations to an average distance of ~8 effective radii from the galaxy centre. Our results show that for some systems the bimodality in GC colour is also present in GC kinematics. The kinematics of the red GC subpopulations are strongly coupled with the host galaxy stellar kinematics. The blue GC subpopulations are more dominated by random motions, especially in the outer regions, and decoupled from the red GCs. Peculiar GC kinematic profiles are seen in some galaxies: the blue GCs in NGC 821 rotate along the galaxy minor axis, whereas the GC system of the lenticular galaxy NGC 7457 appears to be strongly rotation supported in the outer region. We supplement our galaxy sample with data from the literature and carry out a number of tests to study the kinematic differences between the two GC subpopulations. We confirm that the GC kinematics are coupled with the host galaxy properties and find that the velocity kurtosis and the slope of their velocity dispersion profiles is different between the two GC subpopulations in more massive galaxies.

There are only a few tracers available to probe the kinematics of individual early-type galaxies beyond one effective radius. Here we directly compare a sample of planetary nebulae (PNe), globular clusters (GCs) and galaxy starlight velocities out to ~4 effective radii, in the S0 galaxy NGC 2768. Using a bulge-to-disk decomposition of a K-band image we assign PNe and starlight to either the disk or the bulge. We show that the bulge PNe and bulge starlight follow the same radial density distribution as the red subpopulation of GCs, whereas the disk PNe and disk starlight are distinct components. We find good kinematic agreement between the three tracers to several effective radii (and with stellar data in the inner regions). Further support for the distinct nature of the two galaxy components come from our kinematic analysis. After separating the tracers into bulge and disk components we find the bulge to be a slowly rotating pressure-supported system, whereas the disk reveals a rapidly rising rotation curve with a declining velocity dispersion profile. The resulting V/sigma ratio for the disk resembles that of a spiral galaxy and hints at an origin for NGC 2768 as a transformed late-type galaxy. A two-component kinematic analysis for a sample of S0s will help to elucidate the nature of this class of galaxy.

Although the colour distribution of globular clusters in massive galaxies is well known to be bimodal, the spectroscopic metallicity distribution has been measured in only a few galaxies. After redefining the calcium triplet index-metallicity relation, we use our relation to derive the metallicity of 903 globular clusters in 11 early-type galaxies. This is the largest sample of spectroscopic globular cluster metallicities yet assembled. We compare these metallicities with those derived from Lick indices finding good agreement. In 6 of the 8 galaxies with sufficient numbers of high quality spectra we find bimodality in the spectroscopic metallicity distribution. Our results imply that most massive early-type galaxies have bimodal metallicity, as well as colour, distributions. This bimodality suggests that most massive galaxies early-type galaxies experienced two periods of star formation.

We present a comprehensive analysis of the spatial, kinematic, and chemical properties of stars and globular clusters (GCs) in the `ordinary' elliptical galaxy NGC 4494 using data from the Keck and Subaru telescopes. We derive galaxy surface brightness and colour profiles out to large galactocentric radii. We compare the latter to metallicities derived using the near-infrared Calcium Triplet. We obtain stellar kinematics out to ~3.5 effective radii. The latter appear flattened or elongated beyond ~1.8 effective radii in contrast to the relatively round photometric isophotes. In fact, NGC 4494 may be a flattened galaxy, possibly even an S0, seen at an inclination of ~45 degrees. We publish a catalogue of 431 GC candidates brighter than i=24 based on the photometry, of which 109 are confirmed spectroscopically and 54 have measured spectroscopic metallicities. We also report the discovery of 3 spectroscopically confirmed ultra-compact dwarfs around NGC 4494 with measured metallicities of -0.4<[Fe/H]<-0.3. Based on their properties, we conclude that they are simply bright GCs. The metal-poor globular clusters are found to be rotating with similar amplitude as the galaxy stars, while the metal-rich globular clusters show marginal rotation. We supplement our analysis with available literature data and results. Using model predictions of galaxy formation, and a suite of merger simulations, we find that many of the observational properties of NGC 4494 may be explained by formation in a relatively recent gas-rich major merger. Complete studies of individual galaxies incorporating a range of observational avenues and methods such as the one presented here will be an invaluable tool for constraining the fine details of galaxy formation models, especially at large galactocentric radii.