We present a novel fitting formula for the halo concentration enhancement in chameleon $f(R)$ gravity relative to General Relativity (GR). The formula is derived by employing a large set of $N$-body simulations of the Hu-Sawicki $f(R)$ model which cover a wide range of model and cosmological parameters, resolutions and simulation box sizes. The complicated dependence of the concentration on halo mass $M$, redshift $z$, and the $f(R)$ and cosmological parameters can be combined into a simpler form that depends only on a rescaled mass $M/10^{p_2}$, with $p_2\equiv1.5\log_{10}\left[|{\bar{f}_R(z)}|/(1+z)\right]+21.64$ and $\bar{f}_R(z)$ the background scalar field at $z$, irrespective of the $f(R)$ model parameter. Our fitting formula can describe the concentration enhancement well for redshifts $z\leq3$, nearly 7 orders of magnitude in $M/10^{p_2}$ and five decades in halo mass. This is part of a series of works which aims to provide a general framework for self-consistent and unbiased tests of gravity using present and upcoming galaxy cluster surveys. The fitting formula, which is the first quantitative model for the concentration enhancement due to chameleon type modified gravity, is an important part in this framework and will allow continuous exploration of the parameter space. It can also be used to model other statistics such as the matter power spectrum.

We propose a new framework for testing gravity using cluster observations, which aims to provide an unbiased constraint on modified gravity models from Sunyaev Zel'dovich (SZ) and X-ray cluster counts and the cluster gas fraction, among other possible observables. Focusing on a popular $ f(R) $ model of gravity, we propose a novel procedure to recalibrate mass scaling relations from $ \Lambda $CDM to $ f(R) $ gravity for SZ and X-ray cluster observables. We find that the complicated modified gravity effects can be simply modelled as a dependence on a combination of the background scalar field and redshift, $ f_R(z)/(1+z) $, regardless of the $ f(R) $ model parameter. By employing a large suite of N-body simulations, we demonstrate that a theoretically derived tanh fitting formula is in excellent agreement with the dynamical mass enhancement of dark matter haloes for a large range of background field parameters and redshifts. Our framework is sufficiently flexible to allow for tests of other models and inclusion of further observables. The one-parameter description of the dynamical mass enhancement can have important implications on the theoretical modelling of observables and on practical tests of gravity.

We have traced the spatial distributions of intermediate-age and old stars in nine dwarf galaxies in the distant parts of the Local Group, using multi-epoch 3.6 and 4.5 micron data from the DUST in Nearby Galaxies with Spitzer (DUSTiNGS) survey. Using complementary optical imaging from the Hubble Space Telescope, we identify the tip of the red giant branch (TRGB) in the 3.6 micron photometry, separating thermally-pulsating asymptotic giant branch (TP-AGB) stars from the larger red giant branch (RGB) populations. Unlike the constant TRGB in the I-band, at 3.6 micron the TRGB magnitude varies by ~0.7 mag, making it unreliable as a distance indicator. The intermediate-age and old stars are well mixed in two-thirds of the sample with no evidence of a gradient in the ratio of the intermediate-age to old stellar populations outside the central ~1-2'. Variable AGB stars are detected in the outer extremities of the galaxies, indicating that chemical enrichment from these dust-producing stars may occur in the outer regions of galaxies with some frequency. Theories of structure formation in dwarf galaxies must account for the lack of radial gradients in intermediate-age populations and the presence of these stars in the outer extremities of dwarfs. Finally, we identify unique features in individual galaxies, such as extended tidal features in Sex A and Sag DIG and a central concentration of AGB stars in the inner regions of NGC 185 and NGC 147.