Search

We present an analysis of the spatial distribution of globular cluster (GC) systems of 118 nearby early-type galaxies in the Next Generation Virgo Cluster Survey (NGVS) and Mass Assembly of early-Type GaLAxies with their fine Structures (MATLAS) survey programs, which both used MegaCam on the Canada-France-Hawaii Telescope. We describe the procedure used to select GC candidates and fit the spatial distributions of GCs to a two-dimensional S\'ersic function, which provides effective radii (half number radii) and S\'ersic indices, and estimate background contamination by adding a constant term to the S'ersic function. In cases where a neighboring galaxy affects the estimation of the GC spatial distribution in the target galaxy, we fit two 2D S\'ersic functions, simultaneously. We also investigate the color distributions of GCs in our sample by using Gaussian Mixture Modeling. For GC systems with bimodal color distributions, we divide the GCs into blue and red subgroups and fit their respective spatial distributions with S\'ersic functions. Finally, we measure the total number of GCs based on our fitted S\'ersic function, and calculate the GC specific frequency., Comment: 31 pages, 10 Figures (115 additional figures are available in the ApJS online Journal), 2 Tables, accepted for publication in ApJS

The Virgo Cluster is the nearest substantial cluster of galaxies to the Milky Way and a cornerstone of the extragalactic distance scale. Here, we present JWST/NIRCam observations that simultaneously cover the cores and halos of ten galaxies in and around the Virgo Cluster and are designed to perform simultaneous measurements of the tip of the red giant branch (TRGB) and surface brightness fluctuations (SBF). Seven of the targets are within the Virgo Cluster and where we are able to resolve some of the cluster's substructure, while an additional three provide important constraints on Virgo infall. The seven galaxies within Virgo itself all have SBF measurements from the Advanced Camera for Surveys Virgo Cluster Survey (ACSVCS). After adjusting the ACSVCS measurements for the offset from our new JWST TRGB measurements, we determine a distance to the Virgo Cluster of d $=$ 16.17 $\pm$ 0.25 (stat) $\pm$ 0.47 (sys) Mpc. The work presented here is part of a larger program to develop a Population II distance scale through the TRGB and SBF that is completely independent of the prominent Cepheid + Type Ia supernova ladder. The TRGB distances to the galaxies presented here, when combined with future SBF measurements, will provide a crucial step forward for determining whether or not systematic errors can explain the well-known "Hubble tension'' or if there is significant evidence for cracks in the $\Lambda$CDM model., Comment: submitted to ApJ. comments welcome and appreciated. program name has changed from Paper I's "The Population II Extragalactic Distance Scale" to "The TRGB-SBF Project"

Differences between the local value of the Hubble constant measured via the distance ladder versus the value inferred from the cosmic microwave background with the assumption of the standard $\Lambda$CDM model have reached over 5$\sigma$ significance. To determine if this discrepancy is due to new physics or more mundane systematic errors, it is essential to remove as many sources of systematic uncertainty as possible by developing high-precision distance ladders that are independent of the traditional Cepheid and Type Ia supernovae route. Here we present JWST observations of three early-type Fornax Cluster galaxies, the first of fourteen observations from a Cycle 2 JWST program. Our modest integration times allow us to measure highly precise tip of the red giant branch (TRGB) distances, and will also be used to perform measurements of Surface Brightness Fluctuations (SBF). From these three galaxies, we determine an average TRGB distance modulus to the Fornax Cluster of $\mu$ = 31.424 $\pm$ 0.077 mag, or D = 19.3 $\pm$ 0.7 Mpc. With eleven more scheduled observations in nearby elliptical galaxies, our program will allow us set the zero point of the SBF scale to better than 2$\%$ for more distant measurements, charting a path towards a high-precision measurement of $H_{0}$ that is independent of the traditional Cepheid-SN Ia distance ladder., Comment: 14 pages, 9 figures, accepted to ApJ

The surface brightness fluctuation (SBF) method is a robust and efficient way of measuring distances to galaxies containing evolved stellar populations. Although many recent applications of the method have used space-based imaging, SBF remains a powerful technique for ground-based telescopes. Deep, wide-field imaging surveys with subarsecond seeing enable SBF measurements for numerous nearby galaxies. Using a preliminary calibration, Cantiello et al. (2018) presented SBF distances for 89 bright, mainly early-type galaxies observed in the Next Generation Virgo Cluster Survey (NGVS). Here, we present a refined calibration and SBF distances for 278 galaxies extending several magnitudes fainter than in previous work. The derived distances have uncertainties of 5-12\% depending on the properties of the individual galaxies, and our sample is more than three times larger than any previous SBF study of this region. Virgo has a famously complex structure with numerous subclusters, clouds and groups; we associate individual galaxies with the various substructures and map their three-dimensional spatial distribution. Curiously, subcluster A, centered around M87, appears to have two peaks in distance: the main peak at $\sim$16.5 Mpc and a smaller one at $\sim$19.4 Mpc. Subclusters B and C have distances of $\sim$15.8 Mpc. The W and W' groups form a filament-like structure, extending more than 15~Mpc behind the cluster with a commensurate velocity increase of $\sim$1000 \kms\ along its length. These measurements are a valuable resource for future studies of the relationship between galaxy properties and local environment within a dynamic and evolving region., Comment: 30 pages, 15 figures, Acccepted for publication on the ApJ

Nuclear star clusters (NSCs) are dense star clusters located at the centre of galaxies spanning a wide range of masses and morphologies. Analysing NSC occupation statistics in different environments provides an invaluable window into investigating early conditions of high-density star formation and mass assembly in clusters and group galaxies. We use HST/ACS deep imaging to obtain a catalogue of dwarf galaxies in two galaxy clusters in the Shapley Supercluster: the central cluster Abell 3558 and the northern Abell 1736a. The Shapley region is an ideal laboratory to study nucleation as it stands as the highest mass concentration in the nearby Universe. We investigate the NSC occurrence in quiescent dwarf galaxies as faint as $M_{I} = -10$ mag and compare it with all other environments where nucleation data is available. We use galaxy cluster/group halo mass as a proxy for the environment and employ a Bayesian logistic regression framework to model the nucleation fraction ($f_{n}$) as a function of galaxy luminosity and environment. We find a notably high $f_n$ in Abell 3558: at $M_{I} \approx -13.1$ mag, half the galaxies in the cluster host NSCs. This is higher than in the Virgo and Fornax clusters but comparable to the Coma Cluster. On the other hand, the $f_n$ in Abell 1736a is relatively lower, comparable to groups in the Local Volume. We find that the probability of nucleation varies with galaxy luminosity remarkably similarly in galaxy clusters. These results reinforce previous findings of the important role of the environment in NSC formation/growth., Comment: 20 pages, 10 figures, accepted for publication in the MNRAS

We study the size and structure of globular clusters (GC) systems of 118 early-type galaxies from the NGVS, MATLAS, and ACSVCS surveys. Fitting S\'ersic profiles, we investigate the relationship between effective radii of GC systems ($R_{e, \rm gc}$) and galaxy properties. GC systems are 2--4 times more extended than host galaxies across the entire stellar mass range of our sample ($10^{8.3} < M_* < 10^{11.6}~M_{\odot}$). The relationship between $R_{e, \rm gc}$ and galaxy stellar mass exhibits a characteristic "knee" at a stellar mass of $M_p \simeq 10^{10.8}$, similar to galaxy $R_e$--stellar mass relationship. We present a new characterization of the traditional blue and red GC color sub-populations, describing them with respect to host galaxy $(g'-i')$ color ($\Delta_{gi}$): GCs with similar colors to their hosts have a "red" $\Delta_{gi}$, and those significantly bluer GCs have a "blue" $\Delta_{gi}$. The GC populations with red $\Delta_{gi}$, even in dwarf galaxies, are twice as extended as the stars, suggesting that formation or survival mechanisms favor the outer regions. We find a tight correlation between $R_{e, \rm gc}$ and the total number of GCs, with intrinsic scatter $\lesssim 0.1$ dex spanning two and three orders of magnitude in size and number, respectively. This holds for both red and blue subpopulations, albeit with different slopes. Assuming that $N_{GC, Total}$ correlates with $M_{200}$, we find that the red GC systems have effective radii of roughly 1-5\% $R_{\rm 200}$, while the blue GC systems in massive galaxies can have sizes as large as $\sim$10\% $R_{\rm 200}$. Environmental dependence on $R_{e, \rm gc}$ is also found, with lower density environments exhibiting more extended GC systems at fixed mass., Comment: 28 pages, 18 Figures, 3 tables, accepted for publication in ApJ

Systematic studies have revealed hundreds of ultra-compact dwarf galaxies (UCDs) in the nearby Universe. With half-light radii $r_h$ of approximately 10-100 parsecs and stellar masses $M_*$ $\approx$ $10^6-10^8$ solar masses, UCDs are among the densest known stellar systems. Although similar in appearance to massive globular clusters, the detection of extended stellar envelopes, complex star formation histories, elevated mass-to-light ratio, and supermassive black holes suggest that some UCDs are remnant nuclear star clusters of tidally-stripped dwarf galaxies, or even ancient compact galaxies. However, only a few objects have been found in the transient stage of tidal stripping, and this assumed evolutionary path has never been fully traced by observations. Here we show that 106 galaxies in the Virgo cluster have morphologies that are intermediate between normal, nucleated dwarf galaxies and single-component UCDs, revealing a continuum that fully maps this morphological transition, and fills the `size gap' between star clusters and galaxies. Their spatial distribution and redder color are also consistent with stripped satellite galaxies on their first few pericentric passages around massive galaxies. The `ultra-diffuse' tidal features around several of these galaxies directly show how UCDs are forming through tidal stripping, and that this evolutionary path can include an early phase as a nucleated ultra-diffuse galaxy (UDG). These UCDs represent substantial visible fossil remnants of ancient dwarf galaxies in galaxy clusters, and more low-mass remnants probably remain to be found., Comment: Published in Nature. Accepted on September 15

Dwarf galaxies are important tracers of small-scale cosmological structure, yet much of our knowledge about these systems comes from the limited sample of dwarf galaxies within the Local Group. To make a comprehensive inventory of dwarf populations in the local Universe, we require effective methods for deriving distance estimates for large numbers of faint, low surface brightness objects. Here we test the surface brightness fluctuation (SBF) method, traditionally applied to brighter early-type galaxies, on a sample of 20 nearby dwarf galaxies detected in the COSMOS field. These objects are partially resolved in HST ACS images, and have confirmed redshift distances in the range 17-130 Mpc. We discuss the many model choices required in applying the SBF method, and explore how these affect the final distance estimates. Amongst other variations on the method, when applying the SBF method, we alter the standard equation to include a term accounting for the power spectrum of the background, greatly improving our results. For the most robust modelling choices, we find a roughly Gaussian SBF signal that correlates linearly with distance out to distances of 50-100 Mpc, but with only a fraction of the power expected. At larger distances, there is excess power relative to that predicted, probably from undetected point sources. Overall, obtaining accurate SBF distances to faint, irregular galaxies remains challenging, but may yet prove possible with the inclusion of more information about galaxy properties and point source populations, and the use of more advanced techniques., Comment: 18 pages, 26 figures, accepted by MNRAS

The Surface Brightness Fluctuation (SBF) method is a powerful tool for determining distances to early-type galaxies. The method measures the intrinsic variance in a galaxy's surface brightness distribution to determine its distance with an accuracy of about 5%. Here, we discuss the mathematical formalism behind the SBF technique, its calibration, and the practicalities of how measurements are performed. We review the various sources of uncertainties that affect the method and discuss how they can be minimized or controlled through careful observations and data analysis. The SBF technique has already been successfully applied to a large number of galaxies and used for deriving accurate constraints on the Hubble-Lema\^itre constant $H_0$. An approved JWST program will greatly reduce the systematic uncertainties by establishing a firm zero-point calibration using tip of the red giant branch (TRGB) distances. We summarize the existing results and discuss the excellent potential of the SBF method for improving the current constraints on $H_0$., Comment: 22 pages - 2nd version with minor revisions

We present Hubble Space Telescope photometry in optical (F475X) and near-infrared (F110W) bands of the globular cluster (GC) systems of the inner halos of a sample of 15 massive elliptical galaxies. The targets are selected from the volume-limited MASSIVE survey, and chosen to sample a range of environments from sparsely populated groups to BCGs in dense clusters. We also present a quantitative model of the relation between (F475X - F110W) colour and cluster metallicity [M/H], using simulated GCs. Because much of the GC population in such galaxies is built up through accretion, the metallicity distribution of the GC systems might be expected to vary with galaxy environment. The photometry is used to create a completeness-corrected metallicity distribution for each galaxy in the sample, and to fit a double Gaussian curve to each histogram in order to model the two standard red and blue subpopulations. Finally, the properties of the GC metallicity distribution are correlated against galaxy environment. We find that almost no GCS properties and host galaxy environmental properties are correlated, with the exception of a weak but consistent correlation between blue fraction and nth-nearest neighbour surface density. The results suggest that the systemic properties of the GCS, at least in the inner to mid-halo regions, are influenced more strongly by the local environment at early times, rather than by the environmental properties we see today., Comment: 21 pages, 18 figures. Accepted for publication in ApJ

The extragalactic distance scale is fundamental to our understanding of astrophysics and cosmology. In recent years, the surface brightness fluctuation (SBF) method, applied in the near-IR, has proven especially powerful for measuring galaxy distances, first with HST and now with a new JWST program to calibrate the method directly from the tip of the red giant branch (TRGB). So far, however, the distances from space have been gathered slowly, one or two at a time. With the Roman Space Telescope, we have the opportunity to measure uniformly high-quality SBF distances to thousands of galaxies out to hundreds of Mpc. The impact of these data on cosmology and galaxy studies depends on the specifics of the survey, including the filter selection, exposure depth, and (especially) the sky coverage. While the baseline HLWAS survey in four filters plus the grism would yield useful data, the impact would be limited by the relatively small area. A more optimal approach would concentrate on the most efficient passband (F146), adopt an exposure time sufficient to measure good quality distances well out into the Hubble flow, and then maximize the sky coverage within the total time constraints. Grism observations over the same area can provide the needed information on redshifts and spectral energy distributions for compact sources, while colors for larger objects can be obtained from lower resolution surveys. The proposed plan will enable accurate determination of the physical properties of thousands of nearby galaxies, an independent measure of the Hubble constant $H_0$ with negligible statistical error, and competitive constraints on $S_8{\,=\,}\sigma_8(\Omega_m/0.3)^{0.5}$. The resulting data set will be a phenomenal resource for a wide range of studies in astrophysics and cosmology., Comment: 5 pages; submitted to the call for Roman Core Community Survey white papers (revised version corrects a couple references)

The Nancy Grace Roman Telescope's High Latitude Wide Area Survey will have a number of synergies with the Vera Rubin Observatory's Legacy Survey of Space and Time (LSST), particularly for extragalactic star clusters. Understanding the nature of star clusters and star cluster systems are key topics in many areas of astronomy, chief among them stellar evolution, high energy astrophysics, galaxy assembly/dark matter, the extragalactic distance scale, and cosmology. One of the challenges will be disentangling the age/metallicity degeneracy because young ($\sim$Myr) metal-rich clusters have similar SEDs to old ($\sim$Gyr) metal-poor clusters. Rubin will provide homogeneous, $ugrizy$ photometric coverage, and measurements in the red Roman filters will help break the age-metallicity and age-extinction degeneracies, providing the first globular cluster samples that cover wide areas while essentially free of contamination from Milky Way stars. Roman's excellent spatial resolution will also allow measurements of cluster sizes. We advocate for observations of a large sample of galaxies with a range of properties and morphologies in the Rubin/LSST footprint matching the depth of the LSST Wide-Fast-Deep field $i$ band limit (26.3 mag), and recommend adding the F213 filter to the survey., Comment: white paper submitted for Roman CCS input

The Joint Observatories Kavli Science Forum in Chile was organised in a hybrid mode with the aim of encouraging collaborations, not only with the Chilean institutions but also between the different observing facilities based in Chile. The meeting featured scientific talks showing results obtained with the astronomical facilities based in Chile, but significant time was also dedicated to round-table discussions on Life Balance, Diversity-Equity-Inclusion, and the Road Ahead (i.e., the future of those Chile-based facilities)., Comment: 2 pages, 2 figures, Conference Report

We present wide-field, deep $K$-band photometry of 98 luminous early-type galaxies (ETGs) from the MASSIVE survey based on observations taken with the WIRCam instrument on the Canada-France-Hawaii Telescope. Using these images, we extract accurate total $K$-band luminosities ($L_K$) and half-light radii ($R_e$) for this sample of galaxies. We use these new values to explore the size-luminosity and Faber-Jackson relations for massive ETGs. Within this volume-limited sample, we find clear evidence for curvature in both relations, indicating that the most luminous galaxies tend to have larger sizes and smaller velocity dispersions than expected from a simple power-law fit to less luminous galaxies. Our measured relations are qualitatively consistent with the most massive elliptical galaxies forming largely through dissipationless mergers. When the sample is separated into fast and slow rotators, we find the slow rotators to exhibit similar changes in slope with increasing $L_K$, suggesting that low-mass and high-mass slow rotators have different formation histories. The curvatures in the $R_e-L_K$ and $\sigma-L_K$ relations cancel, leading to a relation between dynamical mass and luminosity that is well described by a single power-law: $R_e\sigma^2 \propto {L_K}^b$ with $b \approx 1.2$. This is consistent with the tilt of the fundamental plane observed in lower mass elliptical galaxies., Comment: 17 pages, 7 figures, accepted for publication in MNRAS

With Cosmicflows-4, distances are compiled for 55,877 galaxies gathered into 38,065 groups. Eight methodologies are employed, with the largest numbers coming from the correlations between the photometric and kinematic properties of spiral galaxies (TF) and elliptical galaxies (FP). Supernovae that arise from degenerate progenitors (type Ia Sne) are an important overlapping component. Smaller contributions come from distance estimates from the surface brightness fluctuations of elliptical galaxies and the luminosities and expansion rates of core collapse supernovae (SNII). Cepheid period-luminosity relation and tip of the red giant branch observations founded on local stellar parallax measurements along with the geometric maser distance to NGC 4258 provide the absolute scaling of distances. The assembly of galaxies into groups is an important feature of the study in facilitating overlaps between methodologies. Merging between multiple contributions within a methodology and between methodologies is carried out with Bayesian Markov chain Monte Carlo procedures. The final assembly of distances is compatible with a value of the Hubble constant of $H_0=74.6$ km s$^{-1}$ Mpc$^{-1}$ with the small statistical error of $\pm 0.8$ km s$^{-1}$ Mpc$^{-1}$ but a large potential systematic error of ~3 km s$^{-1}$ Mpc$^{-1}$. Peculiar velocities can be inferred from the measured distances. The interpretation of the field of peculiar velocities is complex because of large errors on individual components and invites analyses beyond the scope of this study., Comment: 38 pages, 24 figures. catalogs available at edd.ifa.hawaii.edu. Revised version, accepted to ApJ

The Gemini North Adaptive Optics (GNAO) facility is the upcoming AO facility for Gemini North providing a state-of-the-art AO system for surveys and time domain science in the era of JWST and Rubin operations. GNAO will be optimized to feed the Gemini infrared Multi Object Spectrograph (GIRMOS). While GIRMOS is the primary science driver for defining the capabilities of GNAO, any instrument operating with an f/32 beam can be deployed using GNAO. The GNAO project includes the development of a new laser guide star facility which will consist of four side-launched laser beams supporting the two primary AO modes of GNAO: a wide-field mode providing an improved image quality over natural seeing for a 2-arcminute circular field-of-view and a narrow-field mode providing near diffraction-limited performance over a 20x20 arcsecond square field-of-view. The GNAO wide field mode will enable GIRMOS's multi-IFU configuration in which the science beam to each individual IFU will be additionally corrected using multi-object AO within GIRMOS. The GNAO narrow field mode will feed the GIRMOS tiled IFU configuration in which all IFUs are combined into a "super"-IFU in the center of the field. GNAO also includes the development of a new Real Time Controller, a new GNAO Facility System Controller and finally the development of a new AO Bench. We present in this paper an overview of the GNAO facility and provide a status update of each product., Comment: SPIE conference 2022 Montreal

We take advantage of exquisitely deep optical imaging data from HST/ACS in the F475W ($g_{F475W}$) and F606W ($V_{F606W}$) bands, to study the properties of the globular cluster (GC) population in the intermediate mass lenticular galaxy PGC 087327, in the Hydra I galaxy cluster. We inspect the photometric (magnitudes and color) and morphometric (compactness, elongation, etc.) properties of sources lying in an area of $\sim19\times19$ kpc centered on PGC 087327, and compare them with four neighbouring fields over the same HST/ACS mosaic. This allowed us to identify a list of GC candidates and to inspect their properties using a background decontamination method. Relative to the four comparison fields, PGC 087327 shows a robust overdensity of GCs, $N_{GC}=82\pm9$. At the estimated magnitude of the galaxy, this number implies a specific frequency of $S_N=1.8\pm0.7$. In spite of the short wavelength interval available with the $g_{F475W}$ and $V_{F606W}$ passbands, the color distribution shows a clear bimodality with a blue peak at $\langle g_{F475W}{-}V_{F606W} \rangle =0.47\pm0.05$ mag and a red peak at $\langle g_{F475W}{-}V_{F606W}\rangle =0.62\pm0.03$ mag. We also observe the typical steeper slope of the radial distribution of red GCs relative to blue ones. Thanks to the unique depth of the available data, we characterize the GC luminosity function (GCLF) well beyond the expected GCLF turn-over. We find $g^{TOM}_{F475W} = 26.54\pm0.10$ mag and $V^{TOM}_{F606W} = 26.08 \pm 0.09$ mag, which after calibration yields a distance of $D_{GCLF} = 56.7 \pm 4.3(statistical) \pm 5.2(systematic)$ Mpc., Comment: Accepted for publication in A&A

We compare infrared surface brightness fluctuation (IR SBF) distances measured in galaxies that have hosted type Ia supernovae (SNIa) to distances estimated from SNIa light curve fits. We show that the properties of SNIa found in IR SBF hosts are very different from those exploding in Cepheid calibrators, therefore, this is a direct test of systematic uncertainties on estimation of the Hubble constant (Ho) using supernovae. The IR SBF results from Jensen et al. (2021; arXiv:2105.08299) provide a large and uniformly measured sample of IR SBF distances which we directly compare with distances to 25 SNIa host galaxies. We divide the Hubble flow SNIa into sub-samples that best match the divergent supernova properties seen in the IR SBF hosts and Cepheid hosts. We further divide the SNIa into a sample with light curve widths and host masses that are congruent to those found in the SBF-calibrated hosts. We refit the light curve stretch and color correlations with luminosity, and use these revised parameters to calibrate the Hubble flow supernovae with IR SBF calibrators. Relative to the Hubble flow, the average calibrator distance moduli vary by 0.03mag depending on the SNIa subsamples examined and this adds a 1.8% systematic uncertainty to our Hubble constant estimate. Based on the IR SBF calibrators, Ho=74.6$\pm$0.9(stat)$\pm$ 2.7(syst) km/s/Mpc, which is consistent with the Hubble constant derived from supernovae calibrated from Cepheid variables. We conclude that IR SBF provides reliable calibration of SNIa with a precision comparable to Cepheid calibrators, and with a significant saving in telescope time., Comment: 18 pages, 10 figures, 4 tables, accepted for publication in the Astrophysical Journal

We present a study of the stellar populations of globular clusters (GCs) in the Virgo Cluster core with a homogeneous spectroscopic catalog of 692 GCs within a major axis distance $R_{\rm maj} = $ 840 kpc from M87. We investigate radial and azimuthal variations in the mean age, total metallicity, [Fe/H], and $\alpha$-element abundance, of blue (metal-poor) and red (metal-rich) GCs using their co-added spectra. We find that the blue GCs have a steep radial gradient in [Z/H] within $R_{\rm maj} =$ 165 kpc, with roughly equal contributions from [Fe/H] and [$\alpha$/Fe], and flat gradients beyond. By contrast, the red GCs show a much shallower gradient in [Z/H], which is entirely driven by [Fe/H]. We use GC-tagged Illustris simulations to demonstrate an accretion scenario where more massive satellites (with more metal- and $\alpha$-rich GCs) sink further into the central galaxy than less massive ones, and where the gradient flattening occurs because of the low GC occupation fraction of low-mass dwarfs disrupted at larger distances. The dense environment around M87 may also cause the steep [$\alpha$/Fe] gradient of the blue GCs, mirroring what is seen in the dwarf galaxy population. The progenitors of red GCs have a narrower mass range than those of blue GCs, which makes their gradients shallower. We also explore spatial inhomogeneity in GC abundances, finding that the red GCs to the northwest of M87 are slightly more metal-rich. Future observations of GC stellar population gradients will be useful diagnostics of halo merger histories., Comment: 22 pages, 8 figures, 4 tables, accepted for publication in ApJ

In this paper we will list a few important goals that need to be addressed in the next decade, also taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant $H_0$, the $\sigma_8$--$S_8$ tension, and other less statistically significant anomalies. While these discordances can still be in part the result of systematic errors, their persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the necessity for new physics or generalisations beyond the standard model. In this paper, we focus on the $5.0\,\sigma$ tension between the {\it Planck} CMB estimate of the Hubble constant $H_0$ and the SH0ES collaboration measurements. After showing the $H_0$ evaluations made from different teams using different methods and geometric calibrations, we list a few interesting new physics models that could alleviate this tension and discuss how the next decade's experiments will be crucial. Moreover, we focus on the tension of the {\it Planck} CMB data with weak lensing measurements and redshift surveys, about the value of the matter energy density $\Omega_m$, and the amplitude or rate of the growth of structure ($\sigma_8,f\sigma_8$). We list a few interesting models proposed for alleviating this tension, and we discuss the importance of trying to fit a full array of data with a single model and not just one parameter at a time. Additionally, we present a wide range of other less discussed anomalies at a statistical significance level lower than the $H_0$--$S_8$ tensions which may also constitute hints towards new physics, and we discuss possible generic theoretical approaches that can collectively explain the non-standard nature of these signals.[Abridged], Comment: Contribution to Snowmass 2021. 224 pages, 27 figures. Accepted for publication in JHEAp

The detection of the accelerated expansion of the Universe has been one of the major breakthroughs in modern cosmology. Several cosmological probes (CMB, SNe Ia, BAO) have been studied in depth to better understand the nature of the mechanism driving this acceleration, and they are being currently pushed to their limits, obtaining remarkable constraints that allowed us to shape the standard cosmological model. In parallel to that, however, the percent precision achieved has recently revealed apparent tensions between measurements obtained from different methods. These are either indicating some unaccounted systematic effects, or are pointing toward new physics. Following the development of CMB, SNe, and BAO cosmology, it is critical to extend our selection of cosmological probes. Novel probes can be exploited to validate results, control or mitigate systematic effects, and, most importantly, to increase the accuracy and robustness of our results. This review is meant to provide a state-of-art benchmark of the latest advances in emerging beyond-standard cosmological probes. We present how several different methods can become a key resource for observational cosmology. In particular, we review cosmic chronometers, quasars, gamma-ray bursts, standard sirens, lensing time-delay with galaxies and clusters, cosmic voids, neutral hydrogen intensity mapping, surface brightness fluctuations, stellar ages of the oldest objects, secular redshift drift, and clustering of standard candles. The review describes the method, systematics, and results of each probe in a homogeneous way, giving the reader a clear picture of the available innovative methods that have been introduced in recent years and how to apply them. The review also discusses the potential synergies and complementarities between the various probes, exploring how they will contribute to the future of modern cosmology., Comment: Invited review article for Living Reviews in Relativity. Accepted version. 191 pages, 54 figures, 11 tables

We present a stellar dynamical mass measurement of a newly detected supermassive black hole (SMBH) at the center of the fast-rotating, massive elliptical galaxy NGC 2693 as part of the MASSIVE survey. We combine high signal-to-noise integral field spectroscopy (IFS) from the Gemini Multi-Object Spectrograph (GMOS) with wide-field data from the Mitchell Spectrograph at McDonald Observatory to extract and model stellar kinematics of NGC 2693 from the central $\sim 150$ pc out to $\sim2.5$ effective radii. Observations from Hubble Space Telescope (HST) WFC3 are used to determine the stellar light distribution. We perform fully triaxial Schwarzschild orbit modeling using the latest TriOS code and a Bayesian search in 6-D galaxy model parameter space to determine NGC 2693's SMBH mass ($M_\text{BH}$), stellar mass-to-light ratio, dark matter content, and intrinsic shape. We find $M_\text{BH} = \left(1.7\pm 0.4\right)\times 10^{9}\ M_\odot$ and a triaxial intrinsic shape with axis ratios $p=b/a=0.902 \pm 0.009$ and $q=c/a=0.721^{+0.011}_{-0.010}$, triaxiality parameter $T = 0.39 \pm 0.04$. In comparison, the best-fit orbit model in the axisymmetric limit and (cylindrical) Jeans anisotropic model of NGC 2693 prefer $M_\text{BH} = \left(2.4\pm 0.6\right)\times 10^{9}\ M_\odot$ and $M_\text{BH} = \left(2.9\pm 0.3\right)\times 10^{9}\ M_\odot$, respectively. Neither model can account for the non-axisymmetric stellar velocity features present in the IFS data., Comment: 18 pages, 9 figures, 2 tables, accepted to ApJ 2022/02/24

The stellar initial mass function (IMF) is a fundamental property in the measurement of stellar masses and galaxy star formation histories. In this work we focus on the most massive galaxies in the nearby universe $\log(M_{\star}/M_{\odot})>11.2$. We obtain high quality Magellan/LDSS-3 long slit spectroscopy with a wide wavelength coverage of $0.4\mu{\rm m}-1.01\mu{\rm m}$ for 41 early-type galaxies (ETGs) in the MASSIVE survey, and derive high S/N spectra within an aperture of $R_{\rm e}/8$. Using detailed stellar synthesis models, we constrain the elemental abundances and stellar IMF of each galaxy through full spectral modeling. All the ETGs in our sample have an IMF that is steeper than a Milky Way (Kroupa) IMF. The best-fit IMF mismatch parameter, $\alpha_{\rm IMF}=(M/L)/(M/L)_{\rm MW}$, ranges from 1.12 to 3.05, with an average of $\langle \alpha_{\rm IMF} \rangle=1.84$, suggesting that on average, the IMF is more bottom-heavy than Salpeter. Comparing the estimated stellar mass with the dynamical mass, we find that most galaxies have stellar masses smaller than their dynamical masses within the $1\sigma$ uncertainty. We complement our sample with lower-mass galaxies from the literature, and confirm that $\log(\alpha_{\rm IMF})$ is positively correlated with $\log(\sigma)$, $\log(M_{\star})$, and $\log(M_{\rm dyn})$. The IMF in the centers of more massive ETGs is more bottom-heavy. In addition, we find that $\log(\alpha_{\rm IMF})$ is positively correlated with both [Mg/Fe] and the estimated total metallicity [Z/H]. We find suggestive evidence that the effective stellar surface density $\Sigma_{\rm Kroupa}$ might be responsible for the variation of $\alpha_{\rm IMF}$. We conclude that $\sigma$, [Mg/Fe] and [Z/H] are the primary drivers of the global stellar IMF variation., Comment: 23 pages, 11 figures, submitted to ApJ

We measured high-quality surface brightness fluctuation (SBF) distances for a sample of 63 massive early-type galaxies using the WFC3/IR camera on the Hubble Space Telescope. The median uncertainty on the SBF distance measurements is 0.085 mag, or 3.9% in distance. Achieving this precision at distances of 50 to 100 Mpc required significant improvements to the SBF calibration and data analysis procedures for WFC3/IR data. Forty-two of the galaxies are from the MASSIVE Galaxy Survey, a complete sample of massive galaxies within ~100 Mpc; the SBF distances for these will be used to improve the estimates of the stellar and central supermassive black hole masses in these galaxies. Twenty-four of the galaxies are Type Ia supernova hosts, useful for calibrating SN Ia distances for early-type galaxies and exploring possible systematic trends in the peak luminosities. Our results demonstrate that the SBF method is a powerful and versatile technique for measuring distances to galaxies with evolved stellar populations out to 100 Mpc and constraining the local value of the Hubble constant., Comment: Accepted for publication in Astrophysical Journal Supplement Series; 22 pages, 7 figures, with 61 additional figures to be published as an online figure set

We use deep high resolution \textit{HST/ACS} imaging of two fields in the core of the Coma cluster to investigate the occurrence of nuclear star clusters (NSCs) in quiescent dwarf galaxies as faint as $M_{I} = -10$ mag. We employ a hierarchical Bayesian logistic regression framework to model the faint end of the nucleation fraction ($f_{n}$) as a function of both galaxy luminosity and environment. We find that $f_n$ is remarkably high in Coma: at $M_{I} \approx -13$ mag half of the cluster dwarfs still host prominent NSCs. Comparison with dwarf systems in nearby clusters and groups shows that, within the uncertainties, the rate at which the probability of nucleation varies with galaxy luminosity is nearly universal. On the other hand, the fraction of nucleated galaxies at fixed luminosity does exhibit an environmental dependence. More massive environments feature higher nucleation fractions and fainter values of the half-nucleation luminosity, which roughly scales with host halo virial mass as $L_{I,f_{n50}} \propto \mathcal{M}_{200}^{-0.2}$. Our results reinforce the role of galaxy luminosity/mass as a major driver of the efficiency of NSC formation and also indicate a clear secondary dependence on the environment, hence paving the way to more refined theoretical models., Comment: 15 pages, 11 figures, accepted to MNRAS

We present a measurement of the Hubble constant $H_0$ from surface brightness fluctuation (SBF) distances for 63 bright, mainly early-type galaxies out to 100 Mpc observed with the Wide Field Camera 3 Infrared Channel (WFC3/IR) on the Hubble Space Telescope (HST). The sample is drawn from several independent HST imaging programs using the F110W bandpass of WFC3/IR. The majority of galaxies are in the 50 to 80 Mpc range and come from the MASSIVE galaxy survey. The median statistical uncertainty on individual distance measurements is 4%. We construct the Hubble diagram with these IR SBF distances and constrain $H_0$ using {four} different treatments of the galaxy velocities. For the SBF zero point calibration, we use both the existing tie to Cepheid variables, updated for consistency with the latest determination of the distance to the Large Magellanic Cloud from detached eclipsing binaries, and a new tie to the tip of the red giant branch (TRGB) calibrated from the maser distance to NGC4258. These two SBF calibrations are consistent with each other and with theoretical predictions from stellar population models. From a weighted average of the Cepheid and TRGB calibrations, we derive $H_0=73.3{\,\pm\,}0.7{\,\pm\,}2.4$ km/s/Mpc, where the error bars reflect the statistical and systematic uncertainties. This result accords well with recent measurements of $H_0$ from Type~Ia supernovae, time delays in multiply lensed quasars, and water masers. The systematic uncertainty could be reduced to below 2% by calibrating the SBF method with precision TRGB distances for a statistical sample of massive early-type galaxies out to the Virgo cluster measured with the James Webb Space Telescope., Comment: Replaced with version accepted to ApJ. 17 pages, including detailed Appendix on the Tip of the Red Giant Branch calibration; have added an additional peculiar velocity model to the analysis

We present a study of ultra compact dwarf (UCD) galaxies in the Virgo cluster based mainly on imaging from the Next Generation Virgo Cluster Survey (NGVS). Using $\sim$100 deg$^{2}$ of $u^*giz$ imaging, we have identified more than 600 candidate UCDs, from the core of Virgo out to its virial radius. Candidates have been selected through a combination of magnitudes, ellipticities, colors, surface brightnesses, half-light radii and, when available, radial velocities. Candidates were also visually validated from deep NGVS images. Subsamples of varying completeness and purity have been defined to explore the properties of UCDs and compare to those of globular clusters and the nuclei of dwarf galaxies with the aim of delineating the nature and origins of UCDs. From a surface density map, we find the UCDs to be mostly concentrated within Virgo's main subclusters, around its brightest galaxies. We identify several subsamples of UCDs -- i.e., the brightest, largest, and those with the most pronounced and/or asymmetric envelopes -- that could hold clues to the origin of UCDs and possible evolutionary links with dwarf nuclei. We find some evidence for such a connection from the existence of diffuse envelopes around some UCDs, and comparisons of radial distributions of UCDs and nucleated galaxies within the cluster., Comment: 70 pages, 24 figures, 6 tables, accepted for publication in ApJS

We use high spatial resolution stellar velocity maps from the Gemini GMOS integral-field spectrograph (IFS) and wide-field velocity maps from the McDonald Mitchell IFS to study the stellar velocity profiles and kinematic misalignments from $\sim 200$ pc to $\sim 20$ kpc in 20 early-type galaxies with stellar mass $M_* > 10^{11.7} M_\odot$ in the MASSIVE survey. While 80% of the galaxies have low spins ($\lambda < 0.1$) and low rotational velocities ($< 50$ km/s) in both the central region and the main body, we find a diverse range of velocity features and misalignment angles. For the 18 galaxies with measurable central kinematic axes, 10 have well aligned kinematic axis and photometric major axis, and the other 8 galaxies have misalignment angles that are distributed quite evenly from $15^\circ$ to the maximal value of $90^\circ$. There is a strong correlation between central kinematic misalignment and galaxy spin, where all 4 galaxies with significant spins have well aligned kinematic and photometric axes, but only 43% of the low-spin galaxies are well aligned. The central and main-body kinematic axes within a galaxy are not always aligned. When the two kinematic axes are aligned ($\sim 60$% of the cases), they are either also aligned with the photometric major axis or orthogonal to it. We find 13 galaxies to also exhibit noticeable local kinematic twists, and one galaxy to have a counter-rotating core. A diverse assembly history consisting of multiple gas-poor mergers of a variety of progenitor mass ratios and orbits is likely to be needed to account for the predominance of low spins and the wide range of central and main-body velocity features reported here for local massive ETGs., Comment: Accepted for publication in The Astrophysical Journal. 15 pages, 7 figures

We present a new stellar dynamical mass measurement ($M_\mathrm{BH}$) of the supermassive black hole (SMBH) in NGC 1453, a fast-rotating massive elliptical galaxy in the MASSIVE survey. We measure stellar kinematics in 135 spatial bins in the central 1.5 kpc by 2 kpc region of the galaxy using high signal-to-noise ($S/N \sim 130$) spectra from the Gemini-North GMOS integral-field spectrograph (IFS). Combining with wide-field IFS kinematics out to $\sim 3$ effective radii and stellar light distributions from Hubble Space Telescope (HST) WFC3 images, we perform Schwarzschild orbit-based mass modeling in the axisymmetric limit to constrain the mass components in NGC 1453. The best-fit black hole mass is $M_\mathrm{BH} =(2.9 \pm 0.4) 10^9 M_\odot$; the mass models without a central black hole are excluded at the $8.7\sigma$ level. The NGC 1453 black hole lies within the intrinsic scatter of the SMBH and galaxy scaling relations, unlike three other galaxies hosting $\gtrsim 10^{10} M_\odot$ SMBHs in the MASSIVE sample. The high-$S/N$ GMOS spectra enable us to determine 8 moments of the Gauss-Hermite expansion of the line-of-sight velocity distributions (LOSVDs), which are used as constraints in the orbit modeling. The stellar orbits in the mass models are further constrained to produce negligible $h_9$ through $h_{12}$ to minimize spurious behavior in the LOSVDs. We show that truncating the series at $h_4$, as was often done in prior work, leads to a much weaker constraint on the inferred $M_\mathrm{BH}$ for NGC 1453. Furthermore, we discuss precautions and modifications that are needed to achieve axisymmetry in triaxial orbit codes that use the Schwarzschild method to sample the start space of stellar orbits in triaxial gravitational potentials., Comment: Accepted for publication in The Astrophysical Journal. 22 pages, 12 figures; References updated

We describe the goals and first results of a Program for Imaging of the PERseus cluster of galaxies (PIPER). The first phase of the program builds on imaging of fields obtained with the Hubble Space Telescope (HST) ACS/WFC and WFC3/UVIS cameras. Our PIPER target fields with HST include major early-type galaxies including the active central giant NGC 1275; known Ultra-Diffuse Galaxies; and the Intracluster Medium. The resulting two-color photometry in F475W and F814W reaches deep enough to resolve and measure the globular cluster (GC) populations in the Perseus member galaxies. Here we present initial results for eight pairs of outer fields that confirm the presence of Intergalactic GCs (IGCs) in fields as distant as 740 kpc from the Perseus center (40\% of the virial radius of the cluster). Roughly 90% of these IGCs are identifiably blue (metal-poor) but there is a clear trace of a red (metal-rich) component as well, even at these very remote distances., Comment: In press for The Astrophysical Journal

As part of the National Science Foundation funded "Gemini in the Era of MultiMessenger Astronomy" (GEMMA) program, Gemini Observatory is developing GNAO, a widefield adaptive optics (AO) facility for Gemini-North on Maunakea, the only 8m-class open-access telescope available to the US astronomers in the northern hemisphere. GNAO will provide the user community with a queue-operated Multi-Conjugate AO (MCAO) system, enabling a wide range of innovative solar system, Galactic, and extragalactic science with a particular focus on synergies with JWST in the area of time-domain astronomy. The GNAO effort builds on institutional investment and experience with the more limited block-scheduled Gemini Multi-Conjugate System (GeMS), commissioned at Gemini South in 2013. The project involves close partnerships with the community through the recently established Gemini AO Working Group and the GNAO Science Team, as well as external instrument teams. The modular design of GNAO will enable a planned upgrade to a Ground Layer AO (GLAO) mode when combined with an Adaptive Secondary Mirror (ASM). By enhancing the natural seeing by an expected factor of two, GLAO will vastly improve Gemini North's observing efficiency for seeing-limited instruments and strengthen its survey capabilities for multi-messenger astronomy., Comment: 12 pages, 6 figures, ASTRO 2020 Decadal Survey

We use observations from the GEMINI-N/GMOS integral-field spectrograph (IFS) to obtain spatially resolved stellar kinematics of the central $\sim 1$ kpc of 20 early-type galaxies (ETGs) with stellar masses greater than $10^{11.7} M_\odot$ in the MASSIVE survey. Together with observations from the wide-field Mitchell IFS at McDonald Observatory in our earlier work, we obtain unprecedentedly detailed kinematic maps of local massive ETGs, covering a scale of $\sim 0.1-30$ kpc. The high ($\sim 120$) signal-to-noise of the GMOS spectra enable us to obtain two-dimensional maps of the line-of-sight velocity, velocity dispersion $\sigma$, as well as the skewness $h_3$ and kurtosis $h_4$ of the stellar velocity distributions. All but one galaxy in the sample have $\sigma(R)$ profiles that increase towards the center, whereas the slope of $\sigma(R)$ at one effective radius ($R_e$) can be of either sign. The $h_4$ is generally positive, with 14 of the 20 galaxies having positive $h_4$ within the GMOS aperture and 18 having positive $h_4$ within $1 R_e$. The positive $h_4$ and rising $\sigma(R)$ towards small radii are indicative of a central black hole and velocity anisotropy. We demonstrate the constraining power of the data on the mass distributions in ETGs by applying Jeans anisotropic modeling (JAM) to NGC~1453, the most regular fast rotator in the sample. Despite the limitations of JAM, we obtain a clear $\chi^2$ minimum in black hole mass, stellar mass-to-light ratio, velocity anisotropy parameters, and the circular velocity of the dark matter halo., Comment: Accepted to ApJ

The globular cluster (GC) specific frequency ($S_N$), defined as the number of GCs per unit galactic luminosity, represents the efficiency of GC formation (and survival) compared to field stars. Despite the naive expectation that star cluster formation should scale directly with star formation, this efficiency varies widely across galaxies. To explore this variation we measure the z-band GC specific frequency ($S_{N,z}$) for 43 early-type galaxies (ETGs) from the Hubble Space Telescope (HST)/Advanced Camera for Surveys (ACS) Fornax Cluster Survey. Combined with the homogenous measurements of $S_{N,z}$ in 100 ETGs from the HST/ACS Virgo Cluster Survey from Peng et al. (2008), we investigate the dependence of $S_{N,z}$ on mass and environment over a range of galaxy properties. We find that $S_{N,z}$ behaves similarly in the two galaxy clusters, despite the clusters' order-of-magnitude difference in mass density. The $S_{N,z}$ is low in intermediate-mass ETGs ($-20

In this paper we study the molecular gas content of a representative sample of 67 of the most massive early-type galaxies in the local universe, drawn uniformly from the MASSIVE survey. We present new IRAM-30m telescope observations of 30 of these galaxies, allowing us to probe the molecular gas content of the entire sample to a fixed molecular-to-stellar mass fraction of 0.1%. The total detection rate in this representative sample is 25$^{+5.9}_{-4.4}$%, and by combining the MASSIVE and ATLAS$^{\rm 3D}$ molecular gas surveys we find a joint detection rate of 22.4$^{+2.4}_{-2.1}$%. This detection rate seems to be independent of galaxy mass, size, position on the fundamental plane, and local environment. We show here for the first time that true slow rotators can host molecular gas reservoirs, but the rate at which they do so is significantly lower than for fast-rotators. Objects with a higher velocity dispersion at fixed mass (a higher kinematic bulge fraction) are less likely to have detectable molecular gas, and where gas does exist, have lower molecular gas fractions. In addition, satellite galaxies in dense environments have $\approx$0.6 dex lower molecular gas-to-stellar mass ratios than isolated objects. In order to interpret these results we created a toy model, which we use to constrain the origin of the gas in these systems. We are able to derive an independent estimate of the gas-rich merger rate in the low-redshift universe. These gas rich mergers appear to dominate the supply of gas to ETGs, but stellar mass loss, hot halo cooling and transformation of spiral galaxies also play a secondary role., Comment: 19 pages, 12 figures, accepted by MNRAS

We measure the stellar populations as a function of radius for 90 early-type galaxies (ETGs) in the MASSIVE survey, a volume-limited integral-field spectroscopic (IFS) galaxy survey targeting all northern-sky ETGs with absolute K-band magnitude M_K < -25.3 mag, or stellar mass M* 4x10^11 M_sun, within 108 Mpc. We are able to measure reliable stellar population parameters for individual galaxies out to 10-20 kpc (1-3 R_e) depending on the galaxy. Focusing on ~R_e (~10 kpc), we find significant correlations between the abundance ratios, sigma, and M* at large radius, but we also find that the abundance ratios saturate in the highest-mass bin. We see a strong correlation between the kurtosis of the line of sight velocity distribution (h4) and the stellar population parameters beyond R_e. Galaxies with higher radial anisotropy appear to be older, with metal-poorer stars and enhanced [alpha/Fe]. We suggest that the higher radial anisotropy may derive from more accretion of small satellites. Finally, we see some evidence for correlations between environmental metrics (measured locally and on >5 Mpc scales) and the stellar populations, as expected if satellites are quenched earlier in denser environments., Comment: 17 pages, 7 figures, Appendix not included here due to size constraints. Posted after responding to referee's comments

The conclusion that the diffuse galaxy NGC1052-DF2 is deficient in dark matter depends critically on the distance, as noted in the original publication and subsequent works. The results of our preliminary assessment of the SBF distance have already been referenced on the web. Here we present the details of a more careful analysis, supply some context, and comment on the remaining uncertainty., Comment: 2 pages, accepted by Research Notes of the American Astronomical Society; text is a bit laconic because of the RNAAS word limit

Intra-cluster (IC) populations are expected to be a natural result of the hierarchical assembly of clusters, yet their low space densities make them difficult to detect and study. We present the first definitive kinematic detection of an IC population of globular clusters (GCs) in the Virgo cluster, around the central galaxy, M87. This study focuses on the Virgo core for which the combination of NGVS photometry and follow-up spectroscopy allows us to reject foreground star contamination and explore GC kinematics over the full Virgo dynamical range. The GC kinematics changes gradually with galactocentric distance, decreasing in mean velocity and increasing in velocity dispersion, eventually becoming indistinguishable from the kinematics of Virgo dwarf galaxies at $\mathrm{R>320\, kpc}$. By kinematically tagging M87 halo and intra-cluster GCs we find that 1) the M87 halo has a smaller fraction ($52\pm3\%$) of blue clusters with respect to the IC counterpart ($77\pm10\%$), 2) the $(g'-r')_{0}$ vs $(i'-z')_{0}$ color-color diagrams reveal a galaxy population that is redder than the IC population that may be due to a different composition in chemical abundance and progenitor mass, and 3) the ICGC distribution is shallower and more extended than the M87 GCs, yet still centrally concentrated. The ICGC specific frequency, $S_{N,\mathrm{ICL}}=10.2\pm4.8$, is consistent with what is observed for the population of quenched, low-mass galaxies within 1~Mpc from the cluster's center. The IC population at Virgo's center is thus consistent with being an accreted component from low-mass galaxies tidally stripped or disrupted through interactions, with a total mass of $\mathrm{M_{ICL,tot}=10.8\pm0.1\times10^{11}M_{\odot}}$., Comment: 31 pages, 13 figures, 1 table, accepted for publication in ApJ

Ultra-diffuse galaxies (UDGs) are unusual galaxies with low luminosities, similar to classical dwarf galaxies, but sizes up to $\sim\!5$ larger than expected for their mass. Some UDGs have large populations of globular clusters (GCs), something unexpected in galaxies with such low stellar density and mass. We have carried out a comprehensive study of GCs in both UDGs and classical dwarf galaxies at comparable stellar masses using HST observations of the Coma cluster. We present new imaging for 33 Dragonfly UDGs with the largest effective radii ($>2$ kpc), and additionally include 15 UDGs and 54 classical dwarf galaxies from the HST/ACS Coma Treasury Survey and the literature. Out of a total of 48 UDGs, 27 have statistically significant GC systems, and 11 have candidate nuclear star clusters. The GC specific frequency ($S_N$) varies dramatically, with the mean $S_N$ being higher for UDGs than for classical dwarfs. At constant stellar mass, galaxies with larger sizes (or lower surface brightnesses) have higher $S_N$, with the trend being stronger at higher stellar mass. At lower stellar masses, UDGs tend to have higher $S_N$ when closer to the center of the cluster, i.e., in denser environments. The fraction of UDGs with a nuclear star cluster also depends on environment, varying from $\sim\!40$\% in the cluster core, where it is slightly lower than the nucleation fraction of classical dwarfs, to $\lesssim20\%$ in the outskirts. Collectively, we observe an unmistakable diversity in the abundance of GCs, and this may point to multiple formation routes., Comment: 15 pages, 10 figures, 2 tables, accepted for publication in ApJ

Substructure in globular cluster (GC) populations around large galaxies is expected in galaxy formation scenarios that involve accretion or merger events, and it has been searched for using direct associations between GCs and structure in the diffuse galaxy light, or with GC kinematics. Here, we present a search for candidate substructures in the GC population around the Virgo cD galaxy M87 through the analysis of the spatial distribution of the GC colors.~The study is based on a sample of $\sim\!1800$ bright GCs with high-quality $u,g,r,i,z,K_s$ photometry, selected to ensure a low contamination by foreground stars or background galaxies.~The spectral energy distributions of the GCs are associated with formal estimates of age and metallicity, which are representative of its position in a 4-D color-space relative to standard single stellar population models.~Dividing the sample into broad bins based on the relative formal ages, we observe inhomogeneities which reveal signatures of GC substructures.~The most significant of these is a spatial overdensity of GCs with relatively young age labels, of diameter $\sim\!0.1$\,deg ($\sim\!30\,$kpc), located to the south of M87.~The significance of this detection is larger than about 5$\sigma$ after accounting for estimates of random and systematic errors.~Surprisingly, no large Virgo galaxy is present in this area, that could potentially host these GCs.~But candidate substructures in the M87 halo with equally elusive hosts have been described based on kinematic studies in the past.~The number of GC spectra available around M87 is currently insufficient to clarify the nature of the new candidate substructure., Comment: 11 pages, 6 figures, Accepted for publication in ApJ

We use spatially resolved two-dimensional stellar velocity maps over a $107"\times 107"$ field of view to investigate the kinematic features of 90 early-type galaxies above stellar mass $10^{11.5}M_\odot$ in the MASSIVE survey. We measure the misalignment angle $\Psi$ between the kinematic and photometric axes and identify local features such as velocity twists and kinematically distinct components. We find 46% of the sample to be well aligned ($\Psi < 15^{\circ}$), 33% misaligned, and 21% without detectable rotation (non-rotators). Only 24% of the sample are fast rotators, the majority of which (91%) are aligned, whereas 57% of the slow rotators are misaligned with a nearly flat distribution of $\Psi$ from $15^{\circ}$ to $90^{\circ}$. 11 galaxies have $\Psi \gtrsim 60^{\circ}$ and thus exhibit minor-axis ("prolate") rotation in which the rotation is preferentially around the photometric major axis. Kinematic misalignments occur more frequently for lower galaxy spin or denser galaxy environments. Using the observed misalignment and ellipticity distributions, we infer the intrinsic shape distribution of our sample and find that MASSIVE slow rotators are consistent with being mildly triaxial, with mean axis ratios of $b/a=0.88$ and $c/a=0.65$. In terms of local kinematic features, 51% of the sample exhibit kinematic twists of larger than $20^{\circ}$, and 2 galaxies have kinematically distinct components. The frequency of misalignment and the broad distribution of $\Psi$ reported here suggest that the most massive early-type galaxies are mildly triaxial, and that formation processes resulting in kinematically misaligned slow rotators such as gas-poor mergers occur frequently in this mass range., Comment: Accepted to MNRAS

We present near-infrared observations of 35 of the most massive early-type galaxies in the local universe. The observations were made using the infrared channel of the Hubble Space Telescope Wide Field Camera 3 in the F110W (1.1 $\mu$m) filter. We measured surface brightness profiles and elliptical isophotal fit parameters from the nuclear regions out to a radius of ~10 kpc in most cases. We find that 37% (13) of the galaxies in our sample have isophotal position angle rotations greater than 20 degrees over the radial range imaged by WFC3/IR, which is often due to the presence of neighbors or multiple nuclei. Most galaxies in our sample are significantly rounder near the center than in the outer regions. This sample contains six fast rotators and 28 slow rotators. We find that all fast rotators are either disky or show no measurable deviation from purely elliptical isophotes. Among slow rotators, significantly disky and boxy galaxies occur with nearly equal frequency. The galaxies in our sample often exhibit changing isophotal shapes, sometimes showing both significantly disky and boxy isophotes at different radii. The fact that parameters vary widely between galaxies and within individual galaxies is evidence that these massive galaxies have complicated formation histories, and some of them have experienced recent mergers and have not fully relaxed. These data demonstrate the value of high spatial resolution IR imaging of galaxies and provide measurements necessary for determining stellar masses, dynamics, and black hole masses in high mass galaxies., Comment: 31 pages, 39 figures, accepted for publication in the Astrophysical Journal

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(

We measure the radial profiles of the stellar velocity dispersions, $\sigma(R)$, for 90 early-type galaxies (ETGs) in the MASSIVE survey, a volume-limited integral-field spectroscopic (IFS) galaxy survey targeting all northern-sky ETGs with absolute $K$-band magnitude $M_K < -25.3$ mag, or stellar mass $M_* > 4 \times 10^{11} M_\odot$, within 108 Mpc. Our wide-field 107" $\times$ 107" IFS data cover radii as large as 40 kpc, for which we quantify separately the inner (2 kpc) and outer (20 kpc) logarithmic slopes $\gamma_{\rm inner}$ and $\gamma_{\rm outer}$ of $\sigma(R)$. While $\gamma_{\rm inner}$ is mostly negative, of the 56 galaxies with sufficient radial coverage to determine $\gamma_{\rm outer}$ we find 36% to have rising outer dispersion profiles, 30% to be flat within the uncertainties, and 34% to be falling. The fraction of galaxies with rising outer profiles increases with $M_*$ and in denser galaxy environment, with 10 of the 11 most massive galaxies in our sample having flat or rising dispersion profiles. The strongest environmental correlations are with local density and halo mass, but a weaker correlation with large-scale density also exists. The average $\gamma_{\rm outer}$ is similar for brightest group galaxies, satellites, and isolated galaxies in our sample. We find a clear positive correlation between the gradients of the outer dispersion profile and the gradients of the velocity kurtosis $h_4$. Altogether, our kinematic results suggest that the increasing fraction of rising dispersion profiles in the most massive ETGs are caused (at least in part) by variations in the total mass profiles rather than in the velocity anisotropy alone., Comment: Accepted/in press, MNRAS

Large samples of globular clusters (GC) with precise multi-wavelength photometry are becoming increasingly available and can be used to constrain the formation history of galaxies. We present the results of an analysis of Milky Way (MW) and Virgo core GCs based on five optical-near-infrared colors and ten synthetic stellar population models. For the MW GCs, the models tend to agree on photometric ages and metallicities, with values similar to those obtained with previous studies. When used with Virgo core GCs, for which photometry is provided by the Next Generation Virgo cluster Survey (NGVS), the same models generically return younger ages. This is a consequence of the systematic differences observed between the locus occupied by Virgo core GCs and models in panchromatic color space. Only extreme fine-tuning of the adjustable parameters available to us can make the majority of the best-fit ages old. Although we cannot exclude that the formation history of the Virgo core may lead to more conspicuous populations of relatively young GCs than in other environments, we emphasize that the intrinsic properties of the Virgo GCs are likely to differ systematically from those assumed in the models. Thus, the large wavelength coverage and photometric quality of modern GC samples, such as used here, is not by itself sufficient to better constrain the GC formation histories. Models matching the environment-dependent characteristics of GCs in multi-dimensional color space are needed to improve the situation., Comment: 18 pages, 12 figures, Accepted for publication in ApJ

We analyse the environmental properties of 370 local early-type galaxies (ETGs) in the MASSIVE and ATLAS3D surveys, two complementary volume-limited integral-field spectroscopic (IFS) galaxy surveys spanning absolute $K$-band magnitude $-21.5 > M_K > -26.6$, or stellar mass $8 \times 10^{9} < M_* < 2 \times 10^{12} M_\odot$. We find these galaxies to reside in a diverse range of environments measured by four methods: group membership (whether a galaxy is a brightest group/cluster galaxy, satellite, or isolated), halo mass, large-scale mass density (measured over a few Mpc), and local mass density (measured within the $N$th neighbour). The spatially resolved IFS stellar kinematics provide robust measurements of the spin parameter $\lambda_e$ and enable us to examine the relationship among $\lambda_e$, $M_*$, and galaxy environment. We find a strong correlation between $\lambda_e$ and $M_*$, where the average $\lambda_e$ decreases from $\sim 0.4$ to below 0.1 with increasing mass, and the fraction of slow rotators $f_{\rm slow}$ increases from $\sim 10$% to 90%. We show for the first time that at fixed $M_*$, there are almost no trends between galaxy spin and environment; the apparent kinematic morphology-density relation for ETGs is therefore primarily driven by $M_*$ and is accounted for by the joint correlations between $M_*$ and spin, and between $M_*$ and environment. A possible exception is that the increased $f_{\rm slow}$ at high local density is slightly more than expected based only on these joint correlations. Our results suggest that the physical processes responsible for building up the present-day stellar masses of massive galaxies are also very efficient at reducing their spin, in any environment., Comment: Accepted to MNRAS

We present the first systematic investigation of the existence, spatial distribution, and kinematics of warm ionized gas as traced by the [O II] 3727AA emission line in 74 of the most massive galaxies in the local Universe. All of our galaxies have deep integral field spectroscopy from the volume- and magnitude-limited MASSIVE survey of early-type galaxies with stellar mass log(M_*/M_sun) > 11.5 (M_K < -25.3 mag) and distance D < 108 Mpc. Of the 74 galaxies in our sample, we detect warm ionized gas in 28, which yields a global detection fraction of 38\pm6% down to a typical [O II] equivalent width limit of 2AA. MASSIVE fast rotators are more likely to have gas than MASSIVE slow rotators with detection fractions of 80\pm10% and 28\pm6%, respectively. The spatial extents span a wide range of radii (0.6 - 18.2 kpc; 0.1 - 4R_e), and the gas morphologies are diverse, with 17/28 = 61\pm9% being centrally concentrated, 8/28 = 29\pm9% exhibiting clear rotation out to several kpc, and 3/28 = 11\pm6% being extended but patchy. Three out of four fast rotators show kinematic alignment between the stars and gas, whereas the two slow rotators with robust kinematic measurements available exhibit kinematic misalignment. Our inferred warm ionized gas masses are roughly ~10^5M_sun. The emission line ratios and radial equivalent width profiles are generally consistent with excitation of the gas by the old underlying stellar population. We explore different gas origin scenarios for MASSIVE galaxies and find that a variety of physical processes are likely at play, including internal gas recycling, cooling out of the hot gaseous halo, and gas acquired via mergers., Comment: Re-submitted to ApJ after addressing referee report; main body is 14 pages with 6 figures; with appendices, total of 24 pages and 14 figures; comments welcome

The total mass M_GCS in the globular cluster (GC) system of a galaxy is empirically a near-constant fraction of the total mass M_h = M_bary + M_dark of the galaxy, across a range of 10^5 in galaxy mass. This trend is radically unlike the strongly nonlinear behavior of total stellar mass M_star versus M_h. We discuss extensions of this trend to two more extreme situations: (a) entire clusters of galaxies, and (b) the Ultra-Diffuse Galaxies (UDGs) recently discovered in Coma and elsewhere. Our calibration of the ratio \eta_M = M_GCS / M_h from normal galaxies, accounting for new revisions in the adopted mass-to-light ratio for GCs, now gives \eta_M = 2.9 \times 10^{-5} as the mean absolute mass fraction. We find that the same ratio appears valid for galaxy clusters and UDGs. Estimates of \eta_M in the four clusters we examine tend to be slightly higher than for individual galaxies, butmore data and better constraints on the mean GC mass in such systems are needed to determine if this difference is significant. We use the constancy of \eta_M to estimate total masses for several individual cases; for example, the total mass of the Milky Way is calculated to be M_h = 1.1 \times 10^{12} M_sun. Physical explanations for the uniformity of \eta_M are still descriptive, but point to a picture in which massive, dense star clusters in their formation stages were relatively immune to the feedback that more strongly influenced lower-density regions where most stars form., Comment: In press for Astrophysical Journal