The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

We present a new galaxy cluster lens modeling approach, hybrid-Lenstool, that is implemented in the publicly available modeling software Lenstool. hybrid-Lenstool combines a parametric approach to model the core of the cluster, and a non-parametric (free-form) approach to model the outskirts. hybrid-Lenstool optimizes both strong- and weak-lensing constraints simultaneously (Joint-Fit), providing a self-consistent reconstruction of the cluster mass distribution on all scales. In order to demonstrate the capabilities of the new algorithm, we tested it on a simulated cluster. hybrid-Lenstool yields more accurate reconstructed mass distributions than the former Sequential-Fit approach where the parametric and the non-parametric models are optimized successively. Indeed, we show with the simulated cluster that the mass density profile reconstructed with a Sequential-Fit deviates form the input by $2-3\sigma$ at all scales while the Joint-Fit gives a profile that is within $1-1.5\sigma$ of the true value. This gain in accuracy is consequential for recovering mass distributions exploiting cluster lensing and therefore for all applications of clusters as cosmological probes. Finally we found that the Joint-Fit approach yields shallower slope of the inner density profile than the Sequential-Fit approach, thus revealing possible biases in previous lensing studies.

In the context of the BUFFALO (Beyond Ultra-deep Frontier Fields And Legacy Observations) survey, we present a new analysis of the merging galaxy cluster MACS\,J0416.1-2403 ($z = 0.397$) and its parallel field using the data collected by the Hubble Frontier Fields (HFF) campaign. In this work, we measure the surface mass density from a weak-lensing analysis, and characterise the overall matter distribution in both the cluster and parallel fields. The surface mass distribution derived for the parallel field shows clumpy overdensities connected by filament-like structures elongated in the direction of the cluster core. We also characterise the X-ray emission of the cluster, and compare it with the lensing mass distribution. We identify five substructures at the $>5\sigma$ level over the two fields, four of them being in the cluster one. Furthermore, three of them are located close to the edges of the field of view, and border issues can significantly hamper the determination of their physical parameters. Finally, we compare our results with the predicted subhalo distribution of one of the Hydrangea/C-EAGLE simulated cluster. Significant differences are obtained suggesting the simulated cluster is at a more advanced evolutionary state than MACS\,J0416.1-2403. Our results anticipate the upcoming BUFFALO observations that will link the two HFF fields, extending further the \emph{HST} coverage, and thus allowing a better characterisation of the reported substructures.

Satellite galaxies in galaxy clusters represent a significant fraction of the global galaxy population. Because of the unusual dense environment of clusters, their evolution is driven by different mechanisms than the ones affecting field or central galaxies. Understanding the different interactions they are subject to, and how they are influenced by them, is therefore an important step towards explaining the global picture of galaxy evolution. In this paper, we use the publicly-available high resolution hydrodynamical simulation Illustris-1 to study satellite galaxies in the three most massive host haloes (with masses $M_{200} > 10^{14}\,h^{-1}\rm{M}_{\odot}$) at $z=0$. We measure the Stellar-to-Halo Mass Relation (hereafter SHMR) of the galaxies, and find that for satellites it is shifted towards lower halo masses compared to the SHMR of central galaxies. We provide simple fitting functions for both the central and satellite SHMR. To explain the shift between the two, we follow the satellite galaxies since their time of accretion into the clusters, and quantify the impact of dark matter stripping and star formation. We find that subhaloes start losing their dark matter as soon as they get closer than $\sim 1.5\times R_{\rm{vir}}$ to the centre of their host, and that up to 80\% of their dark matter content gets stripped during infall. On the other hand, star formation quenching appears to be delayed, and galaxies continue to form stars for a few Gyr after accretion. The combination of these two effects impacts the ratio of stellar to dark matter mass which varies drastically during infall, from 0.03 to 0.3.

In Montero-Dorta et al. 2017, we show that luminous red galaxies (LRGs) from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) at $z\sim0.55$ can be divided into two groups based on their star formation histories. So-called fast-growing LRGs assemble $80\%$ of their stellar mass at $z\sim5$, whereas slow-growing LRGs reach the same evolutionary state at $z\sim1.5$. We further demonstrate that these two subpopulations present significantly different clustering properties on scales of $\sim1 - 30 \mathrm{Mpc}$. Here, we measure the mean halo mass of each subsample using the galaxy-galaxy lensing technique, in the $\sim190\deg^2$ overlap of the LRG catalogue and the CS82 and CFHTLenS shear catalogues. We show that fast- and slow-growing LRGs have similar lensing profiles, which implies that they live in haloes of similar mass: $\log\left(M_{\rm halo}^{\rm fast}/h^{-1}\mathrm{M}_{\odot}\right) = 12.85^{+0.16}_{-0.26}$ and $\log\left(M_{\rm halo}^{\rm slow}/h^{-1}\mathrm{M}_{\odot}\right) =12.92^{+0.16}_{-0.22}$. This result, combined with the clustering difference, suggests the existence of galaxy assembly bias, although the effect is too subtle to be definitively proven given the errors on our current weak-lensing measurement. We show that this can soon be achieved with upcoming surveys like DES.

In the hierarchical formation model, galaxy clusters grow by accretion of smaller groups or isolated galaxies. During the infall into the centre of a cluster, the properties of accreted galaxies change. In particular, both observations and numerical simulations suggest that its dark matter halo is stripped by the tidal forces of the host. We use galaxy-galaxy weak lensing to measure the average mass of dark matter haloes of satellite galaxies as a function of projected distance to the centre of the host, for different stellar mass bins. Assuming that the stellar component of the galaxy is less disrupted by tidal stripping, stellar mass can be used as a proxy of the infall mass. We study the stellar to halo mass relation of satellites as a function of the cluster-centric distance to measure tidal stripping. We use the shear catalogues of the DES science verification archive, the CFHTLenS and the CFHT Stripe 82 (CS82) surveys, and we select satellites from the redMaPPer catalogue of clusters. For galaxies located in the outskirts of clusters, we find a stellar to halo mass relation in good agreement with the theoretical expectations from \citet{moster2013} for central galaxies. In the centre of the cluster, we find that this relation is shifted to smaller halo mass for a given stellar mass. We interpret this finding as further evidence for tidal stripping of dark matter haloes in high density environments.

We analyze the spectra of 300,000 luminous red galaxies (LRGs) with stellar masses $M_* \gtrsim 10^{11} M_{\odot}$ from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). By studying their star-formation histories, we find two main evolutionary paths converging into the same quiescent galaxy population at $z\sim0.55$. Fast-growing LRGs assemble $80\%$ of their stellar mass very early on ($z\sim5$), whereas slow-growing LRGs reach the same evolutionary state at $z\sim1.5$. Further investigation reveals that their clustering properties on scales of $\sim$1-30 Mpc are, at a high level of significance, also different. Fast-growing LRGs are found to be more strongly clustered and reside in overall denser large-scale structure environments than slow-growing systems, for a given stellar-mass threshold. Our results imply a dependence of clustering on stellar-mass assembly history (naturally connected to the mass-formation history of the corresponding halos) for a homogeneous population of similar mass and color, which constitutes a strong observational evidence of galaxy assembly bias.

Current and future large redshift surveys, as the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey (SDSS-IV/eBOSS) or the Dark Energy Spectroscopic Instrument (DESI), will use emission-line galaxies (ELG) to probe cosmological models by mapping the large-scale structure of the Universe in the redshift range $0.6 < z < 1.7$. With current data, we explore the halo-galaxy connection by measuring three clustering properties of $g$-selected ELGs as matter tracers in the redshift range $0.6 < z < 1$: (i) the redshift-space two-point correlation function using spectroscopic redshifts from the BOSS ELG sample and VIPERS; (ii) the angular two-point correlation function on the footprint of the CFHT-LS; (iii) the galaxy-galaxy lensing signal around the ELGs using the CFHTLenS. We interpret these observations by mapping them onto the latest high-resolution MultiDark Planck N-body simulation, using a novel (Sub)Halo-Abundance Matching technique that accounts for the ELG incompleteness. ELGs at $z\sim0.8$ live in halos of $(1\pm 0.5)\times10^{12}\,h^{-1}$M$_{\odot}$ and 22.5$\pm2.5$% of them are satellites belonging to a larger halo. The halo occupation distribution of ELGs indicates that we are sampling the galaxies in which stars form in the most efficient way, according to their stellar-to-halo mass ratio.

A measurement by microlensing of the planetary mass function of planets with masses ranging from 5M_E to 10M_J and orbital radii from 0.5 to 10 AU was reported recently. A strategy for extending the mass range down to (1-3)M_E is proposed here. This entails monitoring the peaks of a few tens of microlensing events with moderately high magnifications with 1-2m class telescopes. Planets of a few Earth masses are found to produce deviations of ~ 5% to the peaks of microlensing light curves with durations ~ (0.7-3)hr in events with magnification ~ 100 if the projected separation of the planet lies in the annular region (0.85-1.2)r_E. Similar deviations are produced by Earth mass planets in the annular region (0.95-1.05)r_E. It is possible that sub-Earths could be detected very close to the Einstein ring if they are sufficiently abundant, and also planetary systems with more than one low mass planet.