We present a detailed abundance analysis of the three brightest member stars at the top of the giant branch of the ultra-faint dwarf galaxy Grus~II. All stars exhibit a higher than expected $\mathrm{[Mg/Ca]}$ ratio compared to metal-poor stars in other ultra-faint dwarf galaxies and in the Milky Way halo. Nucleosynthesis in high mass ($\geqslant 20$M$_\odot$) core-collapse supernovae has been shown to create this signature. The abundances of this small sample (3) stars suggest the chemical enrichment of Grus~II could have occurred through substantial high-mass stellar evolution and is consistent with the framework of a top-heavy initial mass function. However, with only three stars it can not be ruled out that the abundance pattern is the result of a stochastic chemical enrichment at early times in the galaxy. The most metal-rich of the three stars also possesses a small enhancement in rapid neutron-capture ($r$-process) elements. The abundance pattern of the $r$-process elements in this star matches the scaled $r$-process pattern of the solar system and $r$-process enhanced stars in other dwarf galaxies and in the Milky Way halo, hinting at a common origin for these elements across a range of environments. All current proposed astrophysical sites of $r$-process element production are associated with high-mass stars, thus the possible top-heavy initial mass function of Grus~II would increase the likelihood of any of these events occurring. The time delay between the $\alpha$ and $r$-process element enrichment of the galaxy favors a neutron star merger as the origin of the $r$-process elements in Grus~II.

We perform a joint analysis of the counts and weak lensing signal of redMaPPer clusters selected from the Dark Energy Survey (DES) Year 1 dataset. Our analysis uses the same shear and source photometric redshifts estimates as were used in the DES combined probes analysis. Our analysis results in surprisingly low values for $S_8 =\sigma_8(\Omega_{\rm m}/0.3)^{0.5}= 0.65\pm 0.04$, driven by a low matter density parameter, $\Omega_{\rm m}=0.179^{+0.031}_{-0.038}$, with $\sigma_8-\Omega_{\rm m}$ posteriors in $2.4\sigma$ tension with the DES Y1 3x2pt results, and in $5.6\sigma$ with the Planck CMB analysis. These results include the impact of post-unblinding changes to the analysis, which did not improve the level of consistency with other data sets compared to the results obtained at the unblinding. The fact that multiple cosmological probes (supernovae, baryon acoustic oscillations, cosmic shear, galaxy clustering and CMB anisotropies), and other galaxy cluster analyses all favor significantly higher matter densities suggests the presence of systematic errors in the data or an incomplete modeling of the relevant physics. Cross checks with X-ray and microwave data, as well as independent constraints on the observable--mass relation from SZ selected clusters, suggest that the discrepancy resides in our modeling of the weak lensing signal rather than the cluster abundance. Repeating our analysis using a higher richness threshold ($\lambda \ge 30$) significantly reduces the tension with other probes, and points to one or more richness-dependent effects not captured by our model.

We present a simulated cosmology analysis using the second and third moments of the weak lensing mass (convergence) maps. The analysis is geared towards the third year (Y3) data from the Dark Energy Survey (DES), but the methodology can be applied to other weak lensing data sets. The second moment, or variances, of the convergence as a function of smoothing scale contains information similar to standard shear 2-point statistics. The third moment, or the skewness, contains additional non-Gaussian information. We present the formalism for obtaining the convergence maps from the measured shear and for obtaining the second and third moments of these maps given partial sky coverage. We estimate the covariance matrix from a large suite of numerical simulations. We test our pipeline through a simulated likelihood analyses varying 5 cosmological parameters and 10 nuisance parameters. Our forecast shows that the combination of second and third moments provides a 1.5 percent constraint on $S_8 \equiv \sigma_8 (\Omega_{\rm m}/0.3)^{0.5}$ for DES Y3 data. This is 20 percent better than an analysis using a simulated DES Y3 shear 2-point statistics, owing to the non-Gaussian information captured by the inclusion of higher-order statistics. The methodology developed here can be applied to current and future weak lensing datasets to make use of information from non-Gaussianity in the cosmic density field and to improve constraints on cosmological parameters.

The Tully-Fisher relation (TFR) links the stellar mass of a disk galaxy, $M_{\rm str}$, to its rotation speed: it is well approximated by a power law, shows little scatter, and evolves weakly with redshift. The relation has been interpreted as reflecting the mass-velocity scaling ($M\propto V^3$) of dark matter halos, but this interpretation has been called into question by abundance-matching (AM) models, which predict the galaxy-halo mass relation to be non-monotonic and rapidy evolving. We study the TFR of luminous spirals and its relation to AM using the EAGLE set of $\Lambda$CDM cosmological simulations. Matching both relations requires disk sizes to satisfy constraints given by the concentration of halos and their response to galaxy assembly. EAGLE galaxies approximately match these constraints and show a tight mass-velocity scaling that compares favourably with the observed TFR. The TFR is degenerate to changes in galaxy formation efficiency and the mass-size relation; simulations that fail to match the galaxy stellar mass function may fit the observed TFR if galaxies follow a different mass-size relation. The small scatter in the simulated TFR results because, at fixed halo mass, galaxy mass and rotation speed correlate strongly, scattering galaxies along the main relation. EAGLE galaxies evolve with lookback time following approximately the prescriptions of AM models and the observed mass-size relation of bright spirals, leading to a weak TFR evolution consistent with observation out to $z=1$. $\Lambda$CDM models that match both the abundance and size of galaxies as a function of stellar mass have no difficulty reproducing the observed TFR and its evolution.

The scaling of disk galaxy rotation velocity with baryonic mass (the "Baryonic Tully-Fisher" relation, BTF) has long confounded galaxy formation models. It is steeper than the M ~ V^3 scaling relating halo virial masses and circular velocities and its zero point implies that galaxies comprise a very small fraction of available baryons. Such low galaxy formation efficiencies may in principle be explained by winds driven by evolving stars, but the tightness of the BTF relation argues against the substantial scatter expected from such vigorous feedback mechanism. We use the APOSTLE/EAGLE simulations to show that the BTF relation is well reproduced in LCDM simulations that match the size and number of galaxies as a function of stellar mass. In such models, galaxy rotation velocities are proportional to halo virial velocity and the steep velocity-mass dependence results from the decline in galaxy formation efficiency with decreasing halo mass needed to reconcile the CDM halo mass function with the galaxy luminosity function. Despite the strong feedback, the scatter in the simulated BTF is smaller than observed, even when considering all simulated galaxies and not just rotationally-supported ones. The simulations predict that the BTF should become increasingly steep at the faint end, although the velocity scatter at fixed mass should remain small. Observed galaxies with rotation speeds below ~40 km/s seem to deviate from this prediction. We discuss observational biases and modeling uncertainties that may help to explain this disagreement in the context of LCDM models of dwarf galaxy formation.

The cold dark matter halo mass function is much steeper than the galaxy stellar mass function on galactic and subgalactic scales. This difference is usually reconciled by assuming that the galaxy formation efficiency drops sharply with decreasing halo mass, so that virtually no dwarf galaxies form in halos less massive than ~ 10^10 M_sun. In turn, this implies that, at any given radius, the dark mass enclosed by a galaxy must exceed a certain minimum. We use rotation curves of dwarf galaxies compiled from the literature to explore whether their enclosed mass is consistent with these constraints. We find that almost one half of the dwarfs in our sample with stellar mass between 10^6-10^7 M_sun are at odds with this restriction: either they live in halos with masses substantially below 10^10 M_sun or there is a mechanism capable of reducing the dark mass enclosed by some of the faintest dwarfs. Neither possibility is easily accommodated within the standard LCDM scenario. Extending galaxy formation to halos well below 10^10 M_sun would lead to large numbers of dwarf galaxies in excess of current estimates; at the same time, the extremely low stellar mass of the systems involved makes it unlikely that baryonic effects can reduce their dark matter content. Resolving this challenge seems to require new insights into dwarf galaxy formation, or perhaps a radical revision of the prevailing paradigm.