This page contains a collection of results used to validate the behavior of Galacticus in calculations relevant to strong lensing measurements.

PonosV is a dark matter only N-body simulation run by Fiacconi et al. (2016), simulating a halo with a virial mass of \(1.2\times10^{13}\mathrm{M}_\odot\) at \(z=0\) using WMAP cosmological parameters. Fiacconi et al. (2016) measure the projected mass function of subhalos at \(z=0.7\) at a radius of \(0.02 R_\mathrm{vir}\), and fit to a functional form: $$ \frac{\mathrm{d}^2n}{\mathrm{d}\log M_\mathrm{sub} \mathrm{d}A} = \Sigma_\mathrm{sub} \left( \frac{M_\mathrm{sub}}{m_0} \right)^{-\alpha}, $$ where \(M_\mathrm{sub}\) is the bound mass of the subhalo, and \(m_0 = 10^8\mathrm{M}_\odot\) is a pivot mass. To compare to this, several realizations of a mass-matched halo are run using Galacticus. From these, the normalization, \(\Sigma_\mathrm{sub}\), and slope, \(\alpha\), are measured at \(z=0.7\) around the radius \(0.02 R_\mathrm{vir}\), and in the range \(10^8 < M_\mathrm{sub}/\mathrm{M}_\odot < 10^9\). Below, we report the percentage of realizations which result in \(\alpha\) and \(\Sigma_\mathrm{sub}\) above/below the value reported by Fiacconi et al. (2016) for PonosV. There is significant scatter in both \(\alpha\) and \(\Sigma_\mathrm{sub}\) - a model is considered successful if no less than 5% of realizations lie above/below the value measured in PonosV.