Accurate measurement of the cluster mass function is a crucial element in efforts to constrain structure formation models, the normalisation of the matter power spectrum and the cosmological matter density, and the nature and evolution of dark energy. Large weak lensing surveys of ~20,000 galaxy clusters and groups will be key tools in the observational pursuit of that goal. These weak lensing studies often proceed by stacking the lensing signals of many clusters and groups binned by mass-correlated observables such as richness and luminosity; typically such analyses ignore the triaxial structure of dark matter halos on the assumption that the averaging of many shear signals within each mass bin makes its effects (as large as factors of two in mass model parameter estimates from individual clusters) negligible. We test this assumption and find that triaxiality can bias 3D virial mass estimates compared to those for a spherical population by a few percent if suboptimal mass estimators are used. This bias affects not only direct lensing constraints on the mass function but can also affect the scatter and normalization of the mass-observable relations derived from lensing that are so crucial to constraining the cluster mass function with large samples. However, we demonstrate that a careful choice of mass estimator can remove the bias very effectively if the lensing signals from a sufficient number of triaxial halos are averaged together, and further quantify that sufficient number for adequate shape averaging. We thus show that by choosing observable bins to contain an adequate number of halos and by utilizing a carefully chosen 3D mass estimator stacked weak-lensing analyses can give unbiased constraints on the triaxial mass function.

Measuring the 3D distribution of mass on galaxy cluster scales is a crucial test of the LCDM model, providing constraints on the nature of dark matter. Recent work investigating mass distributions of individual galaxy clusters (e.g. Abell 1689) using weak and strong gravitational lensing has revealed potential inconsistencies between the predictions of structure formation models relating halo mass to concentration and those relationships as measured in massive clusters. However, such analyses employ simple spherical halo models while a growing body of work indicates that triaxial 3D halo structure is both common and important in parameter estimates. We recently introduced a Markov Chain Monte Carlo (MCMC) method to fit fully triaxial models to weak lensing data that gives parameter and error estimates that fully incorporate the true shape uncertainty present in nature. In this paper we apply that method to weak lensing data obtained with the ESO/MPG Wide-Field Imager for galaxy clusters A1689, A1835, and A2204, under a range of Bayesian priors derived from theory and from independent X-ray and strong lensing observations. For Abell 1689, using a simple strong lensing prior we find marginalized mean parameter values M_200 = (0.83 +- 0.16)x10^15 M_solar/h and C=12.2 +- 6.7, which are marginally consistent with the mass-concentration relation predicted in LCDM. The large error contours that accompany our triaxial parameter estimates more accurately represent the true extent of our limited knowledge of the structure of galaxy cluster lenses, and make clear the importance of combining many constraints from other theoretical, lensing (strong, flexion), or other observational (X-ray, SZ, dynamical) data to confidently measure cluster mass profiles. (Abridged)

Measuring the 3D distribution of mass on galaxy cluster scales is a crucial test of the LCDM model, providing constraints on the behaviour of dark matter. Recent work investigating mass distributions of individual galaxy clusters (e.g. Abell 1689) using weak and strong gravitational lensing has revealed potential inconsistencies between the predictions of structure formation models relating halo mass to concentration and those relationships as measured in massive clusters. However, such analyses employ simple spherical halo models while a growing body of work indicates that triaxial 3D halo structure is both common and important in parameter estimates. The very strong assumptions about the symmetry of the lensing halo implied with circular or perturbative elliptical NFW models are not physically motivated and lead to incorrect parameter estimates with significantly underestimated error bars. We here introduce a Markov Chain Monte Carlo (MCMC) method to fit fully triaxial models to weak lensing data that gives parameter and error estimates that fully incorporate the true uncertainty present in nature. Applying the MCMC triaxial fitting method to a population of NFW triaxial lenses drawn from the shape distribution of structure formation simulations, we find that including triaxiality cannot explain a population of massive, highly concentrated clusters within the framework of LCDM, but easily explains rare cases of apparently massive, highly concentrated, very efficient lensing clusters. Our MCMC triaxial NFW fitting method is easily expandable to include constraints from additional data types, and its application returns model parameters and errors that more accurately capture the true (and limited) extent of our knowledge of the structure of galaxy cluster lenses. (abridged)

The Hubble constant can be constrained using the time delays between multiple images of gravitationally lensed sources. In some notable cases, typical lensing analyses assuming isothermal galaxy density profiles produce low values for the Hubble constant, inconsistent with the result of the HST Key Project (72 +- 8 km/s/Mpc). Possible systematics in the values of the Hubble constant derived from galaxy lensing systems can result from a number of factors, e.g. neglect of environmental effects, assumption of isothermality, or contamination by line-of-sight structures. One additional potentially important factor is the triaxial structure of the lensing galaxy halo; most lens models account for halo shape simply by perturbing the projected spherical lensing potential, an approximation that is often necessary but that is inadequate at the levels of triaxiality predicted in the CDM paradigm. To quantify the potential error introduced by this assumption in estimates of the Hubble parameter, we strongly lens a distant galaxy through a sample of triaxial softened isothermal halos and use an MCMC method to constrain the lensing halo profile and the Hubble parameter from the resulting multiple image systems. We explore the major degeneracies between the Hubble parameter and several parameters of the lensing model, finding that without a way to accurately break these degeneracies accurate estimates of the Hubble parameter are not possible. Crucially, we find that triaxiality does not significantly bias estimates of the Hubble constant, and offer an analytic explanation for this behaviour in the case of isothermal profiles. Neglected triaxial halo shape cannot contribute to the low Hubble constant values derived in a number of galaxy lens systems.

Measuring the distribution of mass on galaxy cluster scales is a crucial test of the Lambda CDM model, providing constraints on the nature of dark matter. Recent work investigating mass distributions of individual galaxy clusters using gravitational lensing has illuminated potential inconsistencies between the predictions of structure formation models relating halo mass to concentration and those relationships as measured in massive clusters. However, such analyses typically employ only simple spherical halo models with canonical NFW slopes, while the halos formed in simulations show a range of more complex features. Here we investigate the impact of such expected deviations from the canonical NFW halo profile on mass and parameter estimation using weak gravitational lensing on massive cluster scales. The most important deviation is halo triaxiality because it is impossible even with fiducial weak lensing data to fully resolve the three-dimensional structure of the halo due to lensing's sensitivity only to projected mass. Significant elongation of the halo along the line of sight can cause the mass and concentration to be overestimated by as much as 50% and by a factor of 2, respectively, while foreshortening has the opposite effect. Additionally, triaxial halos in certain orientations are much better lenses than their spherical counterparts of the same mass, indicating that clusters chosen for study because of evident lensing are likely to be drawn from the high-triaxiality end of the halo shape distribution. While the effects of triaxiality alone may not be enough to fully explain the very high concentrations reported for some clusters, they go a long way in easing the tensions between observations and the predictions of the cold dark matter paradigm. (Abridged)

The distribution of mass on galaxy cluster scales is an important test of structure formation scenarios, providing constraints on the nature of dark matter itself. Several techniques have been used to probe the mass distributions of clusters, sometimes yielding results which are discrepant, or at odds with clusters formed in simulations - for example giving NFW concentration parameters much higher than expected in the standard CDM model. In addition, the velocity fields of some well studied galaxy clusters reveal the presence of several structures close to the line-of-sight, often not dynamically bound to the cluster itself. We investigate what impact such neighbouring but unbound massive structures would have on the determination of cluster profiles using weak gravitational lensing. Depending on its concentration and mass ratio to the primary halo, one secondary halo close to the line-of-sight can cause the estimated NFW concentration parameter to be significantly higher than that of the primary halo, and also cause the estimated mass to be biased high. Although it is difficult to envisage how this mechanism alone could yield concentrations as high as reported for some clusters, multiple haloes close to the line-of-sight, such as in the case of Abell 1689, can substantially increase the concentration parameter estimate. Together with the fact that clusters are triaxial, and that including baryonic physics also leads to an increase in the concentration of a dark matter halo, the tension between observations and the standard CDM model is eased. If the alignment with the secondary structure is imprecise, then the estimated concentration parameter can also be even lower than that of the primary halo, reinforcing the importance of identifying structures in cluster fields.