Source author record

Tom Broadhurst

Tom Broadhurst appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

80works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

80 published item(s)

preprint2022arXiv

JWST's PEARLS: A JWST/NIRCam view of ALMA sources

We report the results of James Webb Space Telescope/NIRCam observations of 19 (sub)millimeter (submm/mm) sources detected by the Atacama Large Millimeter Array (ALMA). The accurate ALMA positions allowed unambiguous identifications of their NIRCam counterparts. Taking gravitational lensing into account, these represent 16 distinct galaxies in three fields and constitute the largest sample of its kind to date. The counterparts' spectral energy distributions from rest-frame ultraviolet to near infrared provide photometric redshifts ($1<z<4.5$) and stellar masses ($M_*>10^{10.5}$ Msol), which are similar to sub-millimeter galaxy (SMG) hosts studied previously. However, our sample is fainter in submm/mm than the classic SMG samples are, and our sources exhibit a wider range of properties. They have dust-embedded star-formation rates as low as 10 Msol yr$^{-1}$, and the sources populate both the star-forming main sequence and the quiescent categories. The deep NIRCam data allow us to study the rest-frame near-IR morphologies. Excluding two multiply imaged systems and one quasar, the majority of the remaining sources are disk-like and show either little or no disturbance. This suggests that secular growth is a potential route for the assembly of high-mass disk galaxies. While a few hosts have large disks, the majority have small disks (median half-mass radius of 1.6 kpc). At this time, it is unclear whether this is due to the prevalence of small disks at these redshifts or some unknown selection effects of deep ALMA observations. A larger sample of ALMA sources with NIRCam observations will be able to address this question.

preprint2021arXiv

An analytic model for the sub-galactic matter power spectrum in fuzzy dark matter halos

Fuzzy dark matter (FDM), a scalar particle coupled to the gravitational field without self-interaction whose mass range is $m \sim 10^{-24} - 10^{-20}\ \rm{eV}$, is one of the promising alternative dark matter candidates to cold dark matter. The quantum interference pattern, which is a unique structure of FDM, can be seen in halos in cosmological FDM simulations. In this paper, we first provide an analytic model of the sub-galactic matter power spectrum originating from quantum clumps in FDM halos, in which the density distribution of the FDM is expressed by a superposition of quantum clumps whose size corresponds to the de Broglie wavelength of the FDM. These clumps are assumed to be distributed randomly such that the ensemble averaged density follows the halo profile such as the Navarro-Frenk-White profile. We then compare the convergence power spectrum projected along the line of sight around the Einstein radius, which is converted from the sub-galactic matter power spectrum, to that measured in the strong lens system SDSS J0252+0039. While we find that the current observation provides no useful constraint on the FDM mass, we show that future deep, high spatial resolution observations of strong lens systems can tightly constrain FDM with the mass around $10^{-22}\ \rm{eV}$.

preprint2020arXiv

A Soliton Solution for the Central Dark Masses in 47- Tuc Globular Cluster and Implications for the Axiverse

We offer a standing wave explanation for the rising proper motions of stars at the center of the globular cluster 47-Tucanae, amounting to $\simeq 0.44\%$ of the total mass. We show this can be explained as a solitonic core of dark matter composed of light bosons, $ m \geq 10^{-18} eV $, corresponding to $ \leq 0.27 pc$, as an alternative to a single black hole (BH) or a concentration of stellar BH remnants proposed recently. This is particularly important as having a concentrated stellar BH remnant with the above radii is very challenging without the heavy core since the three body encounters would prevent the BHs to be that concentrated. We propose this core develops from dark matter captured in the deep gravitational potential of this globular cluster as it orbits the dark halo of our galaxy. This boson may be evidence for a second light axion, additional to a lighter boson of $10^{-22} eV$, favored for the dominant dark matter implied by the large dark cores of dwarf spheroidal galaxies. The identification of two such light bosonic mass scales favors the generic string theory prediction of a wide, discrete mass spectrum of axionic scalar fields.

preprint2020arXiv

Dynamical Evidence of a Solitonic Core of $10^{9}M_\odot$ in the Milky Way

A wavelike solution for the non-relativistic universal dark matter (wave-DM) is rapidly gaining interest, following distinctive predictions of pioneering simulations of cosmic structure as an interference pattern of coherently oscillating bosons. A prominent solitonic standing wave is predicted at the center of every galaxy, representing the ground state, that has been identified with the wide, kpc scale dark cores of common dwarf-spheroidal galaxies, providing a boson mass of, $\simeq 10^{-22}$ eV. A denser soliton is predicted for Milky Way sized galaxies where momentum is higher, so the de Broglie scale of the soliton is smaller, $\simeq 100$ pc, of mass $\simeq 10^9 M_\odot$. Here we show the central motion of bulge stars in the Milky Way implies the presence of such a dark core, where the velocity dispersion rises inversely with radius to a maximum of $\simeq 130$ km/s, corresponding to an excess central mass of $\simeq 1.5\times 10^9 M_\odot$ within $\simeq 100$ pc, favouring a boson mass of $\simeq 10^{-22}$ eV. This quantitative agreement with such a unique and distinctive prediction is therefore strong evidence for a light bosonic solution to the long standing Dark Matter puzzle, such as the axions generic in String Theory.

preprint2020arXiv

Geometric Support for Dark Matter by an Unaligned Einstein Ring in Abell 3827

The non-detection of dark matter (DM) particles in increasingly stringent laboratory searches has encouraged alternative gravity theories where gravity is sourced only from visible matter. Here, we consider whether such theories can pass a two-dimensional test posed by gravitational lensing -- to reproduce a particularly detailed Einstein ring in the core of the galaxy cluster Abell 3827. We find that when we require the lensing mass distribution to strictly follow the shape (ellipticity and position angle) of the light distribution of cluster member galaxies, intracluster stars, and the X-ray emitting intracluster medium, we cannot reproduce the Einstein ring, despite allowing the mass-to-light ratios of these visible components to freely vary with radius to mimic alternative gravity theories. Alternatively, we show that the detailed features of the Einstein ring are accurately reproduced by allowing a smooth, freely oriented DM halo in the lens model, with relatively small contributions from the visible components at a level consistent with their observed brightnesses. This dominant DM component is constrained to have the same orientation as the light from the intracluster stars, indicating that the intracluster stars trace the gravitational potential of this component. The Einstein ring of Abell 3827 therefore presents a new challenge for alternative gravity theories: not only must such theories find agreement between the total lensing mass and visible mass, but they must also find agreement between the projected sky distribution of the lensing mass and that of the visible matter, a more stringent test than has hitherto been posed by lensing data.

preprint2020arXiv

Ghostly Galaxies as Solitons of Bose-Einstein Dark Matter

The large dark cores of common dwarf galaxies are unexplained by the standard heavy particle interpretation of dark matter. This puzzle is exacerbated by the discovery of a very large but barely visible, dark matter dominated galaxy Antlia II orbiting the Milky Way, uncovered by tracking star motions with the {\t Gaia} satellite. Although Antlia II has a low mass, its visible radius is more than double any known dwarf galaxy, with an unprecedentedly low density core. We show that Antlia II favors dark matter as a Bose-Einstein condensate, for which the ground state is a stable soliton with a core radius given by the de Broglie wavelength. The lower the galaxy mass, the larger the de Broglie wavelength, so the least massive galaxies should have the widest soliton cores of lowest density. An ultra-light boson of $m_ψ\sim 1.1 \times10^{-22}$ eV, accounts well for the large size and slowly moving stars within Antlia II, and agrees with boson mass estimates derived from the denser cores of more massive dwarf galaxies. For this very light boson, Antlia II is close to the lower limiting Jeans scale for galaxy formation permitted by the Uncertainty Principle, so other examples are expected but none significantly larger in size. This simple explanation for the puzzling dark cores of dwarf galaxies implies dark matter as an ultra-light boson, such as an axion generic in String Theory.

preprint2020arXiv

Interpreting LIGO/Virgo "Mass-Gap" events as lensed Neutron Star-Black Hole binaries

Gravitational lensing allows the detection of binary black holes (BBH) at cosmological distances with chirp masses that appear to be enhanced by $1+z$ in the range $1<z<4$, in good agreement with the reported BBH masses. We propose this effect also accounts for the puzzling "mass gap" events (MG) newly reported by LIGO/Virgo, as distant, lensed NSBH events with $1<z<4$. The fitted mass of the neutron star member becomes $(1+z)\times 1.4M_\odot$, and is therefore misclassified as a low mass black hole. In this way, we derive a redshift of $z\simeq 3.5$ and $z\simeq 1.0$ for two newly reported "mass asymmetric" events GW190412 \& GW190814, by interpreting them as lensed NSBH events, comprising a stellar mass black hole and neutron star. Over the past year an additional 31 BBH events and 5 MG events have been reported with high probability ($>95\%$), from which we infer a factor $\simeq 5$ higher intrinsic rate of NSBH events than BBH events, reflecting a higher proportion of neutron stars formed by early star formation. We predict a distinctive locus for lensed NSBH events in the observed binary mass plane, spanning $1<z<4$ with a narrow mass ratio, $q \simeq 0.2$, that can be readily tested when the waveform data are unlocked. All such events may show disrupted NS emission and are worthy of prompt follow-up as the high lensing magnification means EM detections are not prohibitive despite the high redshifts that we predict. Such lensed NSBH events provide an exciting prospect of directly charting the history of coalescing binaries via the cosmological redshift of their waveforms, determined relative to the characteristic mass of the neutron star member.

preprint2020arXiv

Multiple Images and Flux Ratio Anomaly of Fuzzy Gravitational Lenses

Extremely light bosonic wave dark matter ($ψ$DM) is an emerging dark matter candidate contesting the conventional cold dark matter paradigm and a model subject to intense scrutiny of late. This work for the first time reports testable salient features pertinent to gravitational lenses of $ψ$DM halos. $ψ$DM halos are distinctly filled with large-amplitude, small-scale density fluctuations with $δρ/ρ_{\rm halo}\sim 1$ in form of density granules. This halo yields ubiquitous flux ratio anomalies of a few tens of percent, as is typically found for lensed quasars, and may also produce rare hexad and octad images, for sources located in well-defined caustic zones. We have found new critical features appearing in the highly de-magnified lens center when the halo has sufficiently high surface density near a very compact massive core.

preprint2020arXiv

Multiple Ultralight Axionic Wave Dark Matter and Astronomical Structures

An ultralight scalar boson with mass $m_1 \simeq 10^{-22}$ eV is gaining credence as a Dark Matter (DM) candidate that explains the dark cores of dwarf galaxies as soliton waves. Such a boson is naturally interpreted as an axion generic in String Theory, with multiple light axions predicted in this context. We examine the possibility of soliton structures over a wide range of scales, accounting for galaxy core masses and the common presence of nuclear star clusters. We present a diagnostic soliton core mass-radius plot that provides a global view, indicating the existence of an additional axion with mass $m_2\simeq 10^{-20}$ eV, with the possibility of a third axion with mass $m_3 \gtrsim 0.5 \times 10^{-18}$ eV. We also argue that the relative mass densities measured for these axions are consistent with their cosmological production via the mis-alignment mechanism.

preprint2020arXiv

Soliton Random Walk and the Cluster-Stripping Problem in Ultralight Dark Matter

Simulations of ultralight, $\sim 10^{-22}\,\rm eV$, bosonic dark matter exhibit rich wave-like structure, including a soliton core within a surrounding halo that continuously self-interferes on the de Broglie scale. We show here that as an inherent consequence, the soliton undergoes a confined random walk at the base of the halo potential. This is significant for the fate of the ancient central star cluster in Eridanus II, as the agitated soliton gravitationally shakes the star cluster in and out of the soliton on a time scale of $\sim 100\,\rm Myr$, so complete tidal disruption of the star cluster can occur within $\sim 1\,\rm Gyr$. This destructive effect can be mitigated by tidal stripping of the halo of Eridanus II, thereby reducing the agitation, depending on its orbit around the Milky Way. Our simulations show the Milky Way tide affects the halo much more than the soliton, so the star cluster in Eridanus II can survive for over $5\,\rm Gyr$ within the soliton if it formed after significant halo stripping.

preprint2020arXiv

The BUFFALO HST Survey

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.

preprint2016arXiv

A Free-Form mass model of the Hubble Frontier Fields Cluster AS1063 (RXC J2248.7-4431) with over one hundred constraints

We derive a free-form mass distribution for the massive cluster AS1063 (z=0.348) using the completed optical imaging from the Hubble Frontier Fields programme. Based on a subset of 11 multiply lensed systems with spectroscopic redshift we produce a lens model that is accurate enough to unveil new multiply lensed systems, totalling over a 100 arclets, and to estimate their redshifts geometrically. Consistency is found between this precise model and that obtained using only the subset of lensed sources with spectroscopically measured redshifts. Although a relatively large elongation of the mass distribution is apparent relative to the X-ray map, no significant offset is found between the centroid of our mass distribution and that of the X-ray emission map, suggesting a relatively relaxed state for this cluster. For the well resolved lensed images we provide detailed model comparisons to illustrate the precision of our model and hence the reliability of our de-lensed sources. A clear linear structure is associated with one such source extending approx. 23 kpc in length, that could be an example of jet-induced star formation, at redshift z=3.1.

preprint2016arXiv

An Accurate Cluster Selection Function for the J-PAS Narrow-Band wide-field survey

The impending Javalambre Physics of the accelerating universe Astrophysical Survey (J-PAS) will be the first wide-field survey of $\gtrsim$ 8500 deg$^2$ to reach the `stage IV' category. Because of the redshift resolution afforded by 54 narrow-band filters, J-PAS is particularly suitable for cluster detection in the range z$<$1. The photometric redshift dispersion is estimated to be only $\sim 0.003$ with few outliers $\lesssim$ 4\% for galaxies brighter than $i\sim23$ AB, because of the sensitivity of narrow band imaging to absorption and emission lines. Here we evaluate the cluster selection function for J-PAS using N-body+semi-analytical realistic mock catalogues. We optimally detect clusters from this simulation with the Bayesian Cluster Finder, and we assess the completeness and purity of cluster detection against the mock data. The minimum halo mass threshold we find for detections of galaxy clusters and groups with both $>$80\% completeness and purity is $M_h \sim 5 \times 10^{13}M_{\odot}$ up to $z\sim 0.7$. We also model the optical observable, $M^*_{\rm CL}$-halo mass relation, finding a non-evolution with redshift and main scatter of $σ_{M^*_{\rm CL} | M_{\rm h}}\sim 0.14 \,dex$ down to a factor two lower in mass than other planned broad-band stage IV surveys, at least. For the $M_{\rm h} \sim 1 \times 10^{14}M_{\odot}$ Planck mass limit, J-PAS will arrive up to $z\sim 0.85$ with a $σ_{M^*_{\rm CL} | M_{\rm h}}\sim 0.12 \, dex$. Therefore J-PAS will provide the largest sample of clusters and groups up to $z\sim 0.8$ with a mass calibration accuracy comparable to X-ray data.

preprint2016arXiv

Contrasting Galaxy Formation from Quantum Wave Dark Matter, $ψ$DM, with $Λ$CDM, using Planck and Hubble Data

The newly established luminosity functions of high-z galaxies at $4 \lesssim z \lesssim 10$ can provide a stringent check on dark matter models that aim to explain the core properties of dwarf galaxies. The cores of dwarf spheroidal galaxies are understood to be too large to be accounted for by free streaming of warm dark matter without overly suppressing the formation of such galaxies. Here we demonstrate with cosmological simulations that wave dark matter, $ψ$DM, appropriate for light bosons such as axions, does not suffer this problem, given a boson mass of $m_ψ \ge 1.2 \times 10^{-22}{\,\rm eV}$ ($2σ$). In this case, the halo mass function is suppressed below $\sim 10^{10}{\,M_\odot}$ at a level that is consistent with the high-z luminosity functions, while simultaneously generating the kpc-scale cores in dwarf galaxies arising from the solitonic ground state in $ψ$DM. We demonstrate that the reionization history in this scenario is consistent with the Thomson optical depth recently reported by Planck, assuming a reasonable ionizing photon production rate. We predict that the luminosity function should turn over slowly around an intrinsic UV luminosity of $M_{\rm UV} \gtrsim -16$ at $z \gtrsim 4$. We also show that for galaxies magnified $\mathord{>}10\times$ in the Hubble Frontier Fields, $ψ$DM predicts an order of magnitude fewer detections than cold dark matter at $z \gtrsim 10$ down to $M_{\rm UV} \sim -15$, allowing us to distinguish between these very different interpretations for the observed coldness of dark matter.

preprint2016arXiv

Young Galaxy Candidates in the Hubble Frontier Fields - III. MACSJ0717.5+3745

In this paper we present the results of our search for and study of $z \gtrsim 6$ galaxy candidates behind the third Frontier Fields (FF) cluster, MACSJ0717.5+3745, and its parallel field, combining data from Hubble and Spitzer. We select 39 candidates using the Lyman Break technique, for which the clear non-detection in optical make the extreme mid-$z$ interlopers hypothesis unlikely. We also take benefit from $z \gtrsim 6$ samples selected using previous Frontier Fields datasets of Abell 2744 and MACS0416 to improve the constraints on the properties of very high-redshift objects. We compute the redshift and the physical properties, such emission lines properties, star formation rate, reddening, and stellar mass for all Frontier Fields objects from their spectral energy distribution using templates including nebular emission lines. We study the relationship between several physical properties and confirm the trend already observed in previous surveys for evolution of star formation rate with galaxy mass, and between the size and the UV luminosity of our candidates. The analysis of the evolution of the UV Luminosity Function with redshift seems more compatible with an evolution of density. Moreover, no robust $z\ge$8.5 object is selected behind the cluster field, and few $z$$\sim$9 candidates have been selected in the two previous datasets from this legacy survey, suggesting a strong evolution in the number density of galaxies between $z$$\sim$8 and 9. Thanks to the use of the lensing cluster, we study the evolution of the star formation rate density produced by galaxies with L$>$0.03L$^{\star}$, and confirm the strong decrease observed between $z$$\sim$8 and 9.

preprint2015arXiv

A Free-Form Prediction for the Reappearance of Supernova Refsdal in the Hubble Frontier Fields Cluster MACSJ1149.5+2223

The massive cluster MACSJ1149.5+2223(z=0.544) displays five very large lensed images of a well resolved spiral galaxy at $z_{\rm spect}=1.491$. It is within one of these images that the first example of a multiply-lensed supernova has been detected recently as part of the Grism Lens-Amplified Survey from Space. The depth of this data also reveals many HII regions within the lensed spiral galaxy which we identify between the five counter-images. Here we expand the capability of our free-form method to incorporate these HII regions locally, with other reliable lensed galaxies added for a global solution. This improved accuracy allows us to estimate when the Refsdal supernova will appear within the other lensed images of the spiral galaxy to an accuracy of $\sim$ 7\%. We predict this supernova will reappear in one of the counter-images (RA=11:49:36.025, DEC=+22:23:48.11, J2000) and on November 1$^{st}$ 2015 (with an estimated error of $\pm$ 25 days) it will be at the same phase as it was when it was originally discovered, offering a unique opportunity to study the early phases of this supernova and to examine the consistency of the mass model and the cosmological model that have an impact on the time delay prediction.

preprint2015arXiv

CLASH-X: A Comparison of Lensing and X-ray Techniques for Measuring the Mass Profiles of Galaxy Clusters

We present profiles of temperature (Tx), gas mass, and hydrostatic mass estimated from new and archival X-ray observations of CLASH clusters. We compare measurements derived from XMM and Chandra observations with one another and compare both to gravitational lensing mass profiles derived with CLASH HST and ground-based lensing data. Radial profiles of Chandra and XMM electron density and enclosed gas mass are nearly identical, indicating that differences in hydrostatic masses inferred from X-ray observations arise from differences in Tx measurements. Encouragingly, cluster Txs are consistent with one another at ~100-200 kpc radii but XMM Tx systematically decline relative to Chandra Tx at larger radii. The angular dependence of the discrepancy suggests additional investigation on systematics such as the XMM point spread function correction, vignetting and off-axis responses. We present the CLASH-X mass-profile comparisons in the form of cosmology-independent and redshift-independent circular-velocity profiles. Ratios of Chandra HSE mass profiles to CLASH lensing profiles show no obvious radial dependence in the 0.3-0.8 Mpc range. However, the mean mass biases inferred from the WL and SaWLens data are different. e.g., the weighted-mean value at 0.5 Mpc is <b> = 0.12 for the WL comparison and <b> = -0.11 for the SaWLens comparison. The ratios of XMM HSE mass profiles to CLASH lensing profiles show a pronounced radial dependence in the 0.3-1.0 Mpc range, with a weighted mean mass bias of value rising to <b>~0.3 at ~1 Mpc for the WL comparison and <b> of 0.25 for SaWLens comparison. The enclosed gas mass profiles from both Chandra and XMM rise to a value 1/8 times the total-mass profiles inferred from lensing at 0.5 Mpc and remain constant outside of that radius, suggesting that [8xMgas] profiles may be an excellent proxy for total-mass profiles at >0.5 Mpc in massive galaxy clusters.

preprint2015arXiv

Galaxy clusters and groups in the ALHAMBRA Survey

We present a catalogue of 348 galaxy clusters and groups with $0.2<z<1.2$ selected in the 2.78 $deg^2$ ALHAMBRA Survey. The high precision of our photometric redshifts, close to $1\%$, and the wide spread of the seven ALHAMBRA pointings ensure that this catalogue has better mass sensitivity and is less affected by cosmic variance than comparable samples. The detection has been carried out with the Bayesian Cluster Finder (BCF), whose performance has been checked in ALHAMBRA-like light-cone mock catalogues. Great care has been taken to ensure that the observable properties of the mocks photometry accurately correspond to those of real catalogues. From our simulations, we expect to detect galaxy clusters and groups with both $70\%$ completeness and purity down to dark matter halo masses of $M_h\sim3\times10^{13}\rm M_{\odot}$ for $z<0.85$. Cluster redshifts are expected to be recovered with $\sim0.6\%$ precision for $z<1$. We also expect to measure cluster masses with $σ_{M_h|M^*_{CL}}\sim0.25-0.35\, dex$ precision down to $\sim3\times10^{13}\rm M_{\odot}$, masses which are $50\%$ smaller than those reached by similar work. We have compared these detections with previous optical, spectroscopic and X-rays work, finding an excellent agreement with the rates reported from the simulations. We have also explored the overall properties of these detections such as the presence of a colour-magnitude relation, the evolution of the photometric blue fraction and the clustering of these sources in the different ALHAMBRA fields. Despite the small numbers, we observe tentative evidence that, for a fixed stellar mass, the environment is playing a crucial role at lower redshifts (z$<$0.5).

preprint2015arXiv

Hubble Space Telescope Combined Strong and Weak Lensing Analysis of the CLASH Sample: Mass and Magnification Models and Systematic Uncertainties

We present results from a comprehensive lensing analysis in HST data, of the complete CLASH cluster sample. We identify new multiple-images previously undiscovered allowing improved or first constraints on the cluster inner mass distributions and profiles. We combine these strong-lensing constraints with weak-lensing shape measurements within the HST FOV to jointly constrain the mass distributions. The analysis is performed in two different common parameterizations (one adopts light-traces-mass for both galaxies and dark matter while the other adopts an analytical, elliptical NFW form for the dark matter), to provide a better assessment of the underlying systematics - which is most important for deep, cluster-lensing surveys, especially when studying magnified high-redshift objects. We find that the typical (median), relative systematic differences throughout the central FOV are $\sim40\%$ in the (dimensionless) mass density, $κ$, and $\sim20\%$ in the magnification, $μ$. We show maps of these differences for each cluster, as well as the mass distributions, critical curves, and 2D integrated mass profiles. For the Einstein radii ($z_{s}=2$) we find that all typically agree within $10\%$ between the two models, and Einstein masses agree, typically, within $\sim15\%$. At larger radii, the total projected, 2D integrated mass profiles of the two models, within $r\sim2\arcmin$, differ by $\sim30\%$. Stacking the surface-density profiles of the sample from the two methods together, we obtain an average slope of $d\log (Σ)/d\log(r)\sim-0.64\pm0.1$, in the radial range [5,350] kpc. Lastly, we also characterize the behavior of the average magnification, surface density, and shear differences between the two models, as a function of both the radius from the center, and the best-fit values of these quantities.

preprint2015arXiv

Illuminating a Dark Lens : A Type Ia Supernova Magnified by the Frontier Fields Galaxy Cluster Abell 2744

SN HFF14Tom is a Type Ia Supernova (SN) discovered at z = 1.3457 +- 0.0001 behind the galaxy cluster Abell 2744 (z = 0.308). In a cosmology-independent analysis, we find that HFF14Tom is 0.77 +- 0.15 magnitudes brighter than unlensed Type Ia SNe at similar redshift, implying a lensing magnification of mu_obs = 2.03 +- 0.29. This observed magnification provides a rare opportunity for a direct empirical test of galaxy cluster lens models. Here we test 17 lens models, 13 of which were generated before the SN magnification was known, qualifying as pure "blind tests". The models are collectively fairly accurate: 8 of the models deliver median magnifications that are consistent with the measured mu to within 1-sigma. However, there is a subtle systematic bias: the significant disagreements all involve models overpredicting the magnification. We evaluate possible causes for this mild bias, and find no single physical or methodological explanation to account for it. We do find that model accuracy can be improved to some extent with stringent quality cuts on multiply-imaged systems, such as requiring that a large fraction have spectroscopic redshifts. In addition to testing model accuracies as we have done here, Type Ia SN magnifications could also be used as inputs for future lens models of Abell 2744 and other clusters, providing valuable constraints in regions where traditional strong- and weak-lensing information is unavailable.

preprint2015arXiv

Three-dimensional Multi-probe Analysis of the Galaxy Cluster A1689

We perform a 3D multi-probe analysis of the rich galaxy cluster A1689 by combining improved weak-lensing data from new BVRi'z' Subaru/Suprime-Cam observations with strong-lensing, X-ray, and Sunyaev-Zel'dovich effect (SZE) data sets. We reconstruct the projected matter distribution from a joint weak-lensing analysis of 2D shear and azimuthally integrated magnification constraints, the combination of which allows us to break the mass-sheet degeneracy. The resulting mass distribution reveals elongation with axis ratio ~0.7 in projection. When assuming a spherical halo, our full weak-lensing analysis yields a projected concentration of $c_{200c}^{2D}=8.9\pm 1.1$ ($c_{vir}^{2D}\sim 11$), consistent with and improved from earlier weak-lensing work. We find excellent consistency between weak and strong lensing in the region of overlap. In a parametric triaxial framework, we constrain the intrinsic structure and geometry of the matter and gas distributions, by combining weak/strong lensing and X-ray/SZE data with minimal geometric assumptions. We show that the data favor a triaxial geometry with minor-major axis ratio 0.39+/-0.15 and major axis closely aligned with the line of sight (22+/-10 deg). We obtain $M_{200c}=(1.2\pm 0.2)\times 10^{15} M_{\odot}/h$ and $c_{200c}=8.4\pm 1.3$, which overlaps with the $>1σ$ tail of the predicted distribution. The shape of the gas is rounder than the underlying matter but quite elongated with minor-major axis ratio 0.60+/-0.14. The gas mass fraction within 0.9Mpc is 10^{+3}_{-2}%. The thermal gas pressure contributes to ~60% of the equilibrium pressure, indicating a significant level of non-thermal pressure support. When compared to Planck's hydrostatic mass estimate, our lensing measurements yield a spherical mass ratio of $M_{Planck}/M_{GL}=0.70\pm 0.15$ and $0.58\pm 0.10$ with and without corrections for lensing projection effects, respectively.

preprint2015arXiv

Young Galaxy Candidates in the Hubble Frontier Fields. II. MACS,J0416-2403

We searched for z > 7 Lyman-break galaxies (LBGs) in the optical-to-mid-infrared Hubble Frontier Field and associated parallel field observations of the strong-lensing cluster MACS J0416-2403. We discovered 22 candidates, of which six lie at z > 9 and one lies at z > 10. Based on the Hubble and Spitzer photometry, all have secure photometric redshifts and a negligible probability of being at lower redshifts, according to their peak probability ratios, R. This substantial increase in the number of known high-redshift galaxies allows a solid determination of the luminosity function at z > 8. The number of high-z candidates in the parallel field is considerably higher than that in the Abell 2744 parallel field. Our candidates have median stellar masses of log(M_*) ~ 8.40^{+0.55}_{-0.31}~Msun, SFRs of ~ 1.6^{+0.5}_{-0.4} Msun yr^-1, and SFR-weighted ages of < 310^{+70}_{-140} Myr. Finally, we are able to put strong constraints on the z = 7,8,9 and 10 luminosity functions. One of the objects in the cluster field is a z ~ 10 candidate, with a magnification of mu ~ 20 +- 13. This object is likely the faintest z ~ 10 object known to date, allowing a first look into the extreme faint-end (L ~ 0.04L*) of the z ~ 10 luminosity function.

preprint2014arXiv

A Geometrically Supported $z\sim10$ Candidate Multiply-Imaged by the Hubble Frontier Fields Cluster Abell 2744

The deflection angles of lensed sources increase with their distance behind a given lens. We utilize this geometric effect to corroborate the $z_{phot}\simeq9.8$ photometric redshift estimate of a faint near-IR dropout, triply-imaged by the massive galaxy cluster Abell 2744 in deep Hubble Frontier Fields images. The multiple images of this source follow the same symmetry as other nearby sets of multiple images which bracket the critical curves and have well defined redshifts (up to $z_{spec}\simeq3.6$), but with larger deflection angles, indicating that this source must lie at a higher redshift. Similarly, our different parametric and non-parametric lens models all require this object be at $z\gtrsim4$, with at least 95\% confidence, thoroughly excluding the possibility of lower-redshift interlopers. To study the properties of this source we correct the two brighter images for their magnifications, leading to a SFR of $\sim0.3 M_{\odot}$/yr, a stellar mass of $\sim4\times10^{7} M_{\odot}$, and an age of $\lesssim220$ Myr (95\% confidence). The intrinsic apparent magnitude is 29.9 AB (F160W), and the rest-frame UV ($\sim1500 Å$) absolute magnitude is $M_{UV,AB}=-17.6$. This corresponds to $\sim0.1 L^{*}_{z=8}$ ($\sim0.2 L^{*}_{z=10}$, adopting $dM^{*}/dz\sim0.45$), making this candidate one of the least luminous galaxies discovered at $z\sim10$.

preprint2014arXiv

A Hydrodynamical Solution for the "Twin-Tailed" Colliding Galaxy Cluster "El Gordo"

The distinctive cometary X-ray morphology of the recently discovered massive galaxy cluster "El Gordo" (ACT-CT J0102-4915; z=0.87) indicates that an unusually high-speed collision is ongoing between two massive galaxy clusters. A bright X-ray "bullet" leads a "twin-tailed" wake, with the SZ centroid at the end of the Northern tail. We show how the physical properties of this system can be determined using our FLASH-based, N-body/hydrodynamic model, constrained by detailed X-ray, Sunyaev-Zel'dovich (SZ), and Hubble lensing and dynamical data. The X-ray morphology and the location of the two Dark Matter components and the SZ peak are accurately described by a simple binary collision viewed about 480 million years after the first core passage. We derive an impact parameter of ~300 kpc, and a relative initial infall velocity of ~2250 km/sec when separated by the sum of the two virial radii assuming an initial total mass of 2.15x10^(15) Msun and a mass ratio of 1.9. Our model demonstrates that tidally stretched gas accounts for the Northern X-ray tail along the collision axis between the mass peaks, and that the Southern tail lies off axis, comprising compressed and shock heated gas generated as the massive component plunges through the main cluster. The challenge for LCDM will be to find out if this physically extreme event can be plausibly accommodated when combined with the similarly massive, high infall velocity case of the "Bullet cluster" and other such cases being uncovered in the new SZ based surveys.

preprint2014arXiv

A Rigorous Free-form Lens Model of Abell 2744 to Meet the Hubble Frontier Fields Challenge

Hubble Frontier Fields (HFF) imaging of the most powerful lensing clusters provides access to the most magnified distant galaxies. The challenge is to construct lens models capable of describing these complex massive, merging clusters so that individual lensed systems can be reliably identified and their intrinsic properties accurately derived. We apply the free-form lensing method (WSLAP+) to A2744, providing a model independent map of the cluster mass, magnification, and geometric distance estimates to multiply-lensed sources. We solve simultaneously for a smooth cluster component on a pixel grid, together with local deflections by the cluster member galaxies. Combining model prediction with photometric redshift measurements, we correct and complete several systems recently claimed, and identify 4 new systems - totalling 65 images of 21 systems spanning a redshift range of 1.4<z<9.8. The reconstructed mass shows small enhancements in the directions where significant amounts of hot plasma can be seen in X-ray. We compare photometric redshifts with "geometric redshifts", finding a high level of self-consistency. We find excellent agreement between predicted and observed fluxes - with a best-fit slope of 0.999+-0.013 and an RMS of ~0.25 mag, demonstrating that our magnification correction of the lensed background galaxies is very reliable. Intriguingly, few multiply-lensed galaxies are detected beyond z~7.0, despite the high magnification and the limiting redshift of z~11.5 permitted by the HFF filters. With the additional HFF clusters we can better examine the plausibility of any pronounced high-z deficit, with potentially important implications for the reionization epoch and the nature of dark matter.

preprint2014arXiv

CLASH-VLT: Constraints on the Dark Matter Equation of State from Accurate Measurements of Galaxy Cluster Mass Profiles

A pressureless scenario for the Dark Matter (DM) fluid is a widely adopted hypothesis, despite the absence of a direct observational evidence. According to general relativity, the total mass-energy content of a system shapes the gravitational potential well, but different test particles perceive this potential in different ways depending on their properties. Cluster galaxy velocities, being $\ll$c, depend solely on the gravitational potential, whereas photon trajectories reflect the contributions from the gravitational potential plus a relativistic-pressure term that depends on the cluster mass. We exploit this phenomenon to constrain the Equation of State (EoS) parameter of the fluid, primarily DM, contained in galaxy clusters. We use the complementary information provided by the kinematic and lensing mass profiles of the galaxy cluster MACS 1206.2-0847 at $z=0.44$, as obtained in an extensive imaging and spectroscopic campaign within the CLASH survey. The unprecedented high quality of our data-set and the properties of this cluster are well suited to determine the EoS parameter of the cluster fluid. Since baryons contribute at most $15\%$ to the total mass in clusters and their pressure is negligible, the EoS parameter we derive describes the behavior of the DM fluid. We obtain the most stringent constraint on the DM EoS parameter to date, $w=(p_r+2\,p_t)/(3\,c^2ρ)=0.00\pm0.15\mathrm{(stat)}\pm0.08\mathrm{(syst)}$, averaged over the radial range $0.5\,\mathrm{Mpc}\leq$$r$$\leq$$r_{200}$, where $p_r$ and $p_t$ are the radial and tangential pressure, and $ρ$ is the density. We plan to further improve our constraint by applying the same procedure to all clusters from the ongoing CLASH-VLT program.

preprint2014arXiv

CLASH: Extreme Emission Line Galaxies and Their Implication on Selection of High-Redshift Galaxies

We utilize the CLASH (Cluster Lensing And Supernova survey with Hubble) observations of 25 clusters to search for extreme emission-line galaxies (EELGs). The selections are carried out in two central bands: F105W (Y105) and F125W (J125), as the flux of the central bands could be enhanced by the presence of [O III] 4959, 5007 at redshift of about 0.93-1.14 and 1.57-1.79, respectively. The multi-band observations help to constrain the equivalent widths of emission lines. Thanks to cluster lensing, we are able to identify 52 candidates down to an intrinsic limiting magnitude of 28.5 and to a rest-frame [O III] 4959,5007 equivalent width of about 3737 angstrom. Our samples include a number of EELGs at lower luminosities that are missed in other surveys, and the extremely high equivalent width can be only found in such faint galaxies. These EELGs can mimic the dropout feature similar to that of high redshift galaxies and contaminate the color-color selection of high redshift galaxies when the S/N ratio is limited or the band coverage is incomplete. We predict that the fraction of EELGs in the future high redshift galaxy selections cannot be neglected.

preprint2014arXiv

CLASH: Weak-Lensing Shear-and-Magnification Analysis of 20 Galaxy Clusters

We present a joint shear-and-magnification weak-lensing analysis of a sample of 16 X-ray-regular and 4 high-magnification galaxy clusters at 0.19<z<0.69 selected from the Cluster Lensing And Supernova survey with Hubble (CLASH). Our analysis uses wide-field multi-color imaging, taken primarily with Suprime-Cam on the Subaru Telescope. From a stacked shear-only analysis of the X-ray-selected subsample, we detect the ensemble-averaged lensing signal with a total signal-to-noise ratio of ~25 in the radial range of 200 to 3500kpc/h. The stacked tangential-shear signal is well described by a family of standard density profiles predicted for dark-matter-dominated halos in gravitational equilibrium, namely the Navarro-Frenk-White (NFW), truncated variants of NFW, and Einasto models. For the NFW model, we measure a mean concentration of $c_{200c}=4.01^{+0.35}_{-0.32}$ at $M_{200c}=1.34^{+0.10}_{-0.09} 10^{15}M_{\odot}$. We show this is in excellent agreement with Lambda cold-dark-matter (LCDM) predictions when the CLASH X-ray selection function and projection effects are taken into account. The best-fit Einasto shape parameter is $α_E=0.191^{+0.071}_{-0.068}$, which is consistent with the NFW-equivalent Einasto parameter of $\sim 0.18$. We reconstruct projected mass density profiles of all CLASH clusters from a joint likelihood analysis of shear-and-magnification data, and measure cluster masses at several characteristic radii. We also derive an ensemble-averaged total projected mass profile of the X-ray-selected subsample by stacking their individual mass profiles. The stacked total mass profile, constrained by the shear+magnification data, is shown to be consistent with our shear-based halo-model predictions including the effects of surrounding large-scale structure as a two-halo term, establishing further consistency in the context of the LCDM model.

preprint2014arXiv

Comparing gravitational redshifts of SDSS galaxy clusters with the magnification redshift enhancement of background BOSS galaxies

A clean measurement of the evolution of the galaxy cluster mass function can significantly improve our understanding of cosmology from the rapid growth of cluster masses below z < 0.5. Here we examine the consistency of cluster catalogues selected from the SDSS by applying two independent gravity-based methods using all available spectroscopic redshifts from the DR10 release. First, we detect a gravitational redshift related signal for 20,119 and 13,128 clusters with spectroscopic redshifts contained in the GMBCG and redMaPPer catalogues, respectively, at a level of $\sim - 10$ km s$^{-1}$. This we show is consistent with the magnitude expected using the richness-mass relations provided by the literature and after applying recently clarified relativistic and flux bias corrections. This signal is also consistent with the richest clusters in the larger catalogue of Wen et al. (2012), corresponding to $M_{200m} \gtrsim 2 \times 10^{14}\,\mathrm{M}_\odot\,h^{-1}$, however we find no significant detection of gravitational redshift signal for less riched clusters, which may be related to bulk motions from substructure and spurious cluster detections. Second, we find all three catalogues generate mass-dependent levels of lensing magnification bias, which enhances the mean redshift of flux-selected background galaxies from the BOSS survey. The magnitude of this lensing effect is generally consistent with the corresponding richness-mass relations advocated for the surveys. We conclude that all catalogues comprise a high proportion of reliable clusters, and that the GMBCG and redMaPPer cluster finder algorithms favor more relaxed clusters with a meaningful gravitational redshift signal, as anticipated by the red-sequence colour selection of the GMBCG and redMaPPer samples.

preprint2014arXiv

Consistent use of Type Ia supernovae highly magnified by galaxy clusters to constrain the cosmological parameters

We discuss how Type Ia supernovae (SNe) strongly magnified by foreground galaxy clusters should be self-consistently treated when used in samples fitted for the cosmological parameters. While the cluster lens magnification of a SN can be well constrained from sets of multiple images of various background galaxies with measured redshifts, its value is typically dependent on the fiducial set of cosmological parameters used to construct the mass model to begin with. In such cases, one should not naively demagnify the observed SN luminosity by the model magnification into the expected Hubble diagram, which would then create a bias, but take into account the cosmological parameters a-priori chosen to construct the mass model. We quantify the effect and find that a systematic error of typically a few percent, up to a few-dozen percent, per magnified SN, may be propagated onto a cosmological parameter fit, unless the cosmology assumed for the mass model is taken into account (the bias can be even larger if the SN is lying very near the critical curves). We also simulate how such a bias propagates onto the cosmological parameter fit using the Union2.1 sample, supplemented with strongly magnified SNe. The resulting bias on the deduced cosmological parameters is generally at the few percent level, if only few biased SNe are included, and increasing with the number of lensed SNe and their redshift. Samples containing magnified Type Ia SNe, e.g. from ongoing cluster surveys, should readily account for this possible bias.

preprint2014arXiv

Cosmic Structure as the Quantum Interference of a Coherent Dark Wave

The conventional cold, particle interpretation of dark matter (CDM) still lacks laboratory support and struggles with the basic properties of common dwarf galaxies, which have surprisingly uniform central masses and shallow density profiles. In contrast, galaxies predicted by CDM extend to much lower masses, with steeper, singular profiles. This tension motivates cold, wavelike dark matter ($ψ$DM) composed of a non-relativistic Bose-Einstein condensate, so the uncertainty principle counters gravity below a Jeans scale. Here we achieve the first cosmological simulations of this quantum state at unprecedentedly high resolution capable of resolving dwarf galaxies, with only one free parameter, $\bf{m_B}$, the boson mass. We demonstrate the large scale structure of this $ψ$DM simulation is indistinguishable from CDM, as desired, but differs radically inside galaxies. Connected filaments and collapsed haloes form a large interference network, with gravitationally self-bound solitonic cores inside every galaxy surrounded by extended haloes of fluctuating density granules. These results allow us to determine $\bf{m_B=(8.1^{+1.6}_{-1.7})\times 10^{-23}~eV}$ using stellar phase-space distributions in dwarf spheroidal galaxies. Denser, more massive solitons are predicted for Milky Way sized galaxies, providing a substantial seed to help explain early spheroid formation. Suppression of small structures means the onset of galaxy formation for $ψ$DM is substantially delayed relative to CDM, appearing at $\bf{z\lesssim 13}$ in our simulations.

preprint2014arXiv

Hubble Frontier Field Free-Form Mass Mapping of the Massive Multiple-Merging Cluster MACSJ0717.5+3745

We examine the latest data on the cluster MACSJ0717.5+3745 from the Hubble Frontier Fields campaign. The critically lensed area is the largest known of any lens and very irregular making it a challenge for parametric modelling. Using our Free-Form method we obtain an accurate solution, identify here many new sets of multiple images, doubling the number of constraints and improving the reconstruction of the dark matter distribution. Our reconstructed mass map shows several distinct central substructures with shallow density profiles, clarifying earlier work and defining well the relation between the dark matter distribution and the luminous and X-ray peaks within the critically lensed region. Using our free-form method, we are able to meaningfully subtract the mass contribution from cluster members to the deflection field to trace the smoothly distributed cluster dark matter distribution. We find 4 distinct concentrations, 3 of which are coincident with the luminous matter. The fourth peak has a significant offset from both the closest luminous and X-ray peaks. These findings, together with dynamical data from the motions of galaxies and gas will be important for uncovering the potentially important implications of this extremely massive and intriguing system.

preprint2014arXiv

The Orthogonally Aligned Dark Halo of an Edge-on Lensing Galaxy in the Hubble Frontier Fields: A Challenge for Modified Gravity

We examine a well resolved 8 arcsec lensed image that is symmetrically bent in the middle by an edge-on lenticular galaxy, in the Hubble Frontier Field (HFF) data of MACSJ0416.1-20403. The lengthy image is generated primarily by the large tangential shear of the cluster with a local secondary deflection by the member galaxy out to a limiting radius of about 18 kpc. The lensing lenticular galaxy is also well resolved and evidently lies nearly edge-on in projection. This fortuitous combination of a long arc intersecting an edge on galaxy provides us with an opportunity to place relatively strong constraints on the lensing effect of this galaxy. We can model the stellar lensing contribution using the observed pixels belonging to the galaxy, in 2D, and we add to this a standard parameterised dark halo component. Irrespective of the detailed choice of parameters we obtain a combined total mass of about 3E11 Msun. Depending on the dark halo parameters, the stellar contribution to this is limited to the range 5-15E10 Msun or 20-50 percent of the total mass, in good agreement with the independent stellar mass computed from the photometry of 5E10 Msun for a Chabrier IMF, or 8E10Msun for a Salpeter IMF. The major axis of the DM halo is constrained to be nearly orthogonal to the plane of the galaxy, within a range of about 15 degrees, and with an ellipticity e=015 corresponding to an axis ratio a/c=0.54. We show that these conclusions are very weakly dependent on the model of the cluster, or the additional influence of neighbouring galaxies or the properties of the lensed source. Alternative theories of gravity where the radial dependence is modified to avoid the need for DM are challenged by this finding since generically these must be tied to the baryonic component which here is a stellar disk oriented nearly orthogonally to the lensed image deflection.

preprint2014arXiv

Three Gravitationally Lensed Supernovae Behind CLASH Galaxy Clusters

We report observations of three gravitationally lensed supernovae (SNe) in the Cluster Lensing And Supernova survey with Hubble (CLASH) Multi-Cycle Treasury program. These objects, SN CLO12Car (z = 1.28), SN CLN12Did (z = 0.85), and SN CLA11Tib (z = 1.14), are located behind three different clusters, MACSJ1720.2+3536 (z = 0.391), RXJ1532.9+3021 (z = 0.345), and Abell 383 (z = 0.187), respectively. Each SN was detected in Hubble Space Telescope (HST) optical and infrared images. Based on photometric classification, we find that SNe CLO12Car and CLN12Did are likely to be Type Ia supernovae (SNe Ia), while the classification of SN CLA11Tib is inconclusive. Using multi-color light-curve fits to determine a standardized SN Ia luminosity distance, we infer that SN CLO12Car was approximately 1.0 +/- 0.2 mag brighter than field SNe Ia at a similar redshift and ascribe this to gravitational lens magnification. Similarly, SN CLN12Did is approximately 0.2 +/- 0.2 mag brighter than field SNe Ia. We derive independent estimates of the predicted magnification from CLASH strong+weak lensing maps of the clusters: 0.83 +/- 0.16 mag for SN CLO12Car, 0.28 +/- 0.08 mag for SN CLN12Did, and 0.43 +/- 0.11 mag for SN CLA11Tib. The two SNe Ia provide a new test of the cluster lens model predictions: we find that the magnifications based on the SN Ia brightness and those predicted by the lens maps are consistent. Our results herald the promise of future observations of samples of cluster-lensed SNe Ia (from the ground or space) to help illuminate the dark-matter distribution in clusters of galaxies, through the direct determination of absolute magnifications.

preprint2014arXiv

Understanding the Core-Halo Relation of Quantum Wave Dark Matter, $ψ$DM, from 3D Simulations

We examine the nonlinear structure of gravitationally collapsed objects that form in our simulations of wavelike cold dark matter ($ψ$DM), described by the Schrödinger-Poisson (SP) equation with a particle mass $\sim 10^{-22} {\rm eV}$. A distinct gravitationally self-bound solitonic core is found at the center of every halo, with a profile quite different from cores modeled in the warm or self-interacting dark matter scenarios. Furthermore, we show that each solitonic core is surrounded by an extended halo composed of large fluctuating dark matter granules which modulate the halo density on a scale comparable to the diameter of the solitonic core. The scaling symmetry of the SP equation and the uncertainty principle tightly relate the core mass to the halo specific energy, which, in the context of cosmological structure formation, leads to a simple scaling between core mass ($M_c$) and halo mass ($M_h$), $M_c \propto a^{-1/2} M_h^{1/3}$, where $a$ is the cosmic scale factor. We verify this scaling relation by (i) examining the internal structure of a statistical sample of virialized halos that form in our 3D cosmological simulations, and by (ii) merging multiple solitons to create individual virialized objects. Sufficient simulation resolution is achieved by adaptive mesh refinement and graphic processing units acceleration. From this scaling relation, present dwarf satellite galaxies are predicted to have kpc sized cores and a minimum mass of $\sim 10^8 {M_\odot}$, capable of solving the small-scale controversies in the cold dark matter model. Moreover, galaxies of $2\times10^{12} {M_\odot}$ at $z=8$ should have massive solitonic cores of $\sim 2\times10^9 {M_\odot}$ within $\sim 60 {\rm pc}$. Such cores can provide a favorable local environment for funneling the gas that leads to the prompt formation of early stellar spheroids and quasars.

preprint2013arXiv

CLASH: Complete Lensing Analysis of the Largest Cosmic Lens MACS J0717.5+3745 and Surrounding Structures

The galaxy cluster MACS J0717.5+3745 (z=0.55) is the largest known cosmic lens, with complex internal structures seen in deep X-ray, Sunyaev-Zel'dovich effect and dynamical observations. We perform a combined weak and strong lensing analysis with wide-field BVRi'z' Subaru/Suprime-Cam observations and 16-band Hubble Space Telescope observations taken as part of the Cluster Lensing And Supernova survey with Hubble (CLASH). We find consistent weak distortion and magnification measurements of background galaxies, and combine these signals to construct an optimally estimated radial mass profile of the cluster and its surrounding large-scale structure out to 5 Mpc/h. We find consistency between strong-lensing and weak-lensing in the region where these independent data overlap, <500 kpc/h. The two-dimensional weak-lensing map reveals a clear filamentary structure traced by distinct mass halos. We model the lensing shear field with 9 halos, including the main cluster, corresponding to mass peaks detected above 2.5σ_κ. The total mass of the cluster as determined by the different methods is M_{vir}=(2.8\pm0.4) \times 10^15 M_sun. Although this is the most massive cluster known at z>0.5, in terms of extreme value statistics we conclude that the mass of MACS J0717.5+3745 by itself is not in serious tension with LambdaCDM, representing only a ~2σ departure above the maximum simulated halo mass at this redshift.

preprint2013arXiv

Cluster Lensing Profiles Derived from a Redshift Enhancement of Magnified BOSS-Survey Galaxies

We report the first detection of a redshift-depth enhancement of background galaxies magnified by foreground clusters. Using 300,000 BOSS-Survey galaxies with accurate spectroscopic redshifts, we measure their mean redshift depth behind four large samples of optically selected clusters from the SDSS surveys, totalling 5,000-15,000 clusters. A clear trend of increasing mean redshift towards the cluster centers is found, averaged over each of the four cluster samples. In addition we find similar but noisier behaviour for an independent X-ray sample of 158 clusters lying in the foreground of the current BOSS sky area. By adopting the mass-richness relationships appropriate for each survey we compare our results with theoretical predictions for each of the four SDSS cluster catalogs. The radial form of this redshift enhancement is well fitted by a richness-to-mass weighted composite Navarro-Frenk-White profile with an effective mass ranging between M_200 ~ 1.4-1.8 10^14 M_sun for the optically detected cluster samples, and M_200 ~ 5.0 10^14 M_sun for the X-ray sample. This lensing detection helps to establish the credibility of these SDSS cluster surveys, and provides a normalization for their respective mass-richness relations. In the context of the upcoming bigBOSS, Subaru-PFS, and EUCLID-NISP spectroscopic surveys, this method represents an independent means of deriving the masses of cluster samples for examining the cosmological evolution, and provides a relatively clean consistency check of weak-lensing measurements, free from the systematic limitations of shear calibration.

preprint2013arXiv

Enabling Non-Parametric Strong Lensing Models to Derive Reliable Cluster Mass Distributions. WSLAP+

In the strong lensing regime non-parametric lens models struggle to achieve sufficient angular resolution for a meaningful derivation of the central cluster mass distribution. The problem lies mainly with cluster members which perturb lensed images and generate additional images, requiring high resolution modeling, even though we mainly wish to understand the relatively smooth cluster component. The required resolution is not achievable because the separation between lensed images is several times larger than the deflection angles by member galaxies, even for the deepest data. Here we bypass this limitation by incorporating a simple physical prior for member galaxies, using their observed positions and their luminosity scaled masses. This galaxy contribution is added to a relatively coarse Gaussian pixel grid for modeling the cluster mass distribution, extending our established WSLAP code (Diego et al. 2007). We test this new code with a simulation based on A1689, using the pixels belonging to multiply-lensed images and the observed member galaxies. Dealing with the cluster members this way leads to convergent solutions, without resorting to regularization, reproducing well the input cluster and substructures. We highlight the ability of this method to recover dark sub-components of the cluster, unrelated to member galaxies. Such anomalies can provide clues to the nature of invisible dark matter, but are hard to discover using parametrized models where substructures are defined by the visible data. With our increased resolution and stability we show, for the first time, that non-parametric models can be made sufficiently precise to locate multiply-lensed systems, thereby achieving fully self-consistent solutions without reliance on input systems from less objective means.

preprint2013arXiv

Galaxy halo truncation and Giant Arc Surface Brightness Reconstruction in the Cluster MACSJ1206.2-0847

In this work we analyze the mass distribution of MACSJ1206.2-0847, especially focusing on the halo properties of its cluster members. The cluster appears relaxed in its X-ray emission, but has significant amounts of intracluster light which is not centrally concentrated, suggesting that galaxy-scale interactions are still ongoing despite the overall relaxed state. The cluster lenses 12 background galaxies into multiple images and one galaxy at $z=1.033$ into a giant arc and its counterimage. The multiple image positions and the surface brightness distribution (SFB) of the arc which is bent around several cluster members are sensitive to the cluster galaxy halo properties. We model the cluster mass distribution with a NFW profile and the galaxy halos with two parameters for the mass normalization and extent of a reference halo assuming scalings with their observed NIR--light. We match the multiple image positions at an r.m.s. level of $0.85\arcsec$ and can reconstruct the SFB distribution of the arc in several filters to a remarkable accuracy based on this cluster model. The length scale where the enclosed galaxy halo mass is best constrained is about 5 effective radii -- a scale in between those accessible to dynamical and field strong lensing mass estimates on one hand and galaxy--galaxy weak lensing results on the other hand. The velocity dispersion and halo size of a galaxy with $m_{\rm 160W,AB}=19.2$ or $M_{\rm B,Vega}=-20.7$ are $σ=150 \rm kms^{-1}$ and $r\approx 26\pm 6 \rm kpc$, indicating that the halos of the cluster galaxies are tidally stripped. We also reconstruct the unlensed source (which is smaller by a factor of $\sim5.8$ in area), demonstrating the increase of morphological information due to lensing and conclude that this galaxy has likely star--forming spiral arms with a red (older) central component.

preprint2013arXiv

Tangential Velocity of the Dark Matter in the Bullet Cluster from Precise Lensed Image Redshifts

We show that the fast moving component of the "bullet cluster" (1E0657-56) can induce potentially resolvable redshift differences between multiply-lensed images of background galaxies. The moving cluster effect can be expressed as the scalar product of the lensing deflection angle with the tangential velocity of the mass components, and it is maximal for clusters colliding in the plane of the sky with velocities boosted by their mutual gravity. The bullet cluster is likely to be the best candidate for the first measurement of this effect due to the large collision velocity and because the lensing deflection and the cluster fields can be calculated in advance. We derive the deflection field using multiply-lensed background galaxies detected with the Hubble Space Telescope. The velocity field is modeled using self-consistent N-body/hydrodynamical simulations constrained by the observed X-ray and gravitational lensing features of this system. We predict that the triply-lensed images of systems "G" and "H" straddling the critical curve of the bullet component will show the largest frequency shifts up to ~0.5 km/sec. This is within the range of the Atacama Large Millimeter/sub-millimeter Array (ALMA) for molecular emission, and is near the resolution limit of the new generation high-throughput optical-IR spectrographs. A detection of this effect measures the tangential motion of the subclusters directly, thereby clarifying the tension with LCDM, which is inferred from gas motion less directly. This method may be extended to smaller redshift differences using the Ly-alpha forest towards QSOs lensed by more typical clusters of galaxies. More generally, the tangential component of the peculiar velocities of clusters derived by our method complements the radial component determined by the kinematic SZ effect, providing a full 3-dimensional description of velocities.

preprint2013arXiv

The Contribution of Halos with Different Mass Ratios to the Overall Growth of Cluster-Sized Halos

We provide a new observational test for a key prediction of the ΛCDM cosmological model: the contributions of mergers with different halo-to-main-cluster mass ratios to cluster-sized halo growth. We perform this test by dynamically analyzing seven galaxy clusters, spanning the redshift range $0.13 < z_c < 0.45$ and caustic mass range $0.4-1.5$ $10^{15} h_{0.73}^{-1}$ M$_{\odot}$, with an average of 293 spectroscopically-confirmed bound galaxies to each cluster. The large radial coverage (a few virial radii), which covers the whole infall region, with a high number of spectroscopically identified galaxies enables this new study. For each cluster, we identify bound galaxies. Out of these galaxies, we identify infalling and accreted halos and estimate their masses and their dynamical states. Using the estimated masses, we derive the contribution of different mass ratios to cluster-sized halo growth. For mass ratios between ~0.2 and ~0.7, we find a ~1 $σ$ agreement with ΛCDM expectations based on the Millennium simulations I and II. At low mass ratios, $\lesssim 0.2$, our derived contribution is underestimated since the detection efficiency decreases at low masses, $\sim 2 \times 10^{14}$ $h_{0.73}^{-1}$ M$_{\odot}$. At large mass ratios, $\gtrsim 0.7$, we do not detect halos probably because our sample, which was chosen to be quite X-ray relaxed, is biased against large mass ratios. Therefore, at large mass ratios, the derived contribution is also underestimated.

preprint2013arXiv

The Pre-merger Impact Velocity of the Binary Cluster A1750 from X-ray, Lensing and Hydrodynamical Simulations

Since the discovery of the "bullet cluster" several similar cases have been uncovered suggesting relative velocities well beyond the tail of high speed collisions predicted by the concordance LCDM model. However, quantifying such post-merger events with hydrodynamical models requires a wide coverage of possible initial conditions. Here we show that it is simpler to interpret pre-merger cases, such as A1750, where the gas between the colliding clusters is modestly affected, so that the initial conditions are clear. We analyze publicly available Chandra data confirming a significant increase in the projected X-ray temperature between the two cluster centers in A1750 consistent with our expectations for a merging cluster. We model this system with a self-consistent hydrodynamical simulation of dark matter and gas using the FLASH code. Our simulations reproduce well the X-ray data, and the measured redshift difference between the two clusters in the phase before the first core passage viewed at an intermediate projection angle. The deprojected initial relative velocity derived using our model is 1460 km/sec which is considerably higher than the predicted mean impact velocity for simulated massive haloes derived by recent LCDM cosmological simulations, but it is within the allowed range. Our simulations demonstrate that such systems can be identified using a multi-wavelength approach and numerical simulations, for which the statistical distribution of relative impact velocities may provide a definitive examination of a broad range of dark matter scenarios.

preprint2012arXiv

A Brightest Cluster Galaxy with an Extremely Large Flat Core

Hubble Space Telescope images of the galaxy cluster Abell 2261, obtained as part of the Cluster Lensing And Supernova survey with Hubble, show that the brightest galaxy in the cluster, A2261-BCG, has the largest core yet detected in any galaxy. The cusp radius of A2261-BCG is 3.2 kpc, twice as big as the next largest core known, and ~3x bigger than those typically seen in the most luminous BCGs. The morphology of the core in A2261-BCG is also unusual, having a flat or even slightly-depressed interior surface brightness profile, rather than the typical shallow cusp. This implies that the galaxy has a core with constant or even centrally decreasing stellar density. Interpretation of the core as an end product of the "scouring" action of a binary supermassive black hole implies a total black hole mass ~1E+10 M_sun from the extrapolation of most relationships between core structure and black hole mass. The core falls 1-sigma above the cusp-radius versus galaxy luminosity relation. Its large size in real terms, and the extremely large black hole mass required to generate it, raise the possibility that the core has been enlarged by additional processes, such as the ejection of the black holes that originally generated the core. The flat central stellar density profile is consistent with this hypothesis. The core is also displaced by 0.7 kpc from the center of the surrounding envelope, consistent with a local dynamical perturbation of the core.

preprint2012arXiv

A highly magnified candidate for a young galaxy seen when the Universe was 500 Myrs old

The early Universe at redshift z\sim6-11 marks the reionization of the intergalactic medium, following the formation of the first generation of stars. However, those young galaxies at a cosmic age of \lesssim 500 million years (Myr, at z \gtrsim 10) remain largely unexplored as they are at or beyond the sensitivity limits of current large telescopes. Gravitational lensing by galaxy clusters enables the detection of high-redshift galaxies that are fainter than what otherwise could be found in the deepest images of the sky. We report the discovery of an object found in the multi-band observations of the cluster MACS1149+22 that has a high probability of being a gravitationally magnified object from the early universe. The object is firmly detected (12 sigma) in the two reddest bands of HST/WFC3, and not detected below 1.2 μm, matching the characteristics of z\sim9 objects. We derive a robust photometric redshift of z = 9.6 \pm 0.2, corresponding to a cosmic age of 490 \pm 15Myr (i.e., 3.6% of the age of the Universe). The large number of bands used to derive the redshift estimate make it one of the most accurate estimates ever obtained for such a distant object. The significant magnification by cluster lensing (a factor of \sim15) allows us to analyze the object's ultra-violet and optical luminosity in its rest-frame, thus enabling us to constrain on its stellar mass, star-formation rate and age. If the galaxy is indeed at such a large redshift, then its age is less than 200 Myr (at the 95% confidence level), implying a formation redshift of zf \lesssim 14. The object is the first z>9 candidate that is bright enough for detailed spectroscopic studies with JWST, demonstrating the unique potential of galaxy cluster fields for finding highly magnified, intrinsically faint galaxies at the highest redshifts.

preprint2012arXiv

CLASH: Mass Distribution in and around MACS J1206.2-0847 from a Full Cluster Lensing Analysis

We derive an accurate mass distribution of the galaxy cluster MACS J1206.2-0847 (z=0.439) from a combined weak-lensing distortion, magnification, and strong-lensing analysis of wide-field Subaru BVRIz' imaging and our recent 16-band Hubble Space Telescope observations taken as part of the Cluster Lensing And Supernova survey with Hubble (CLASH) program. We find good agreement in the regions of overlap between several weak and strong lensing mass reconstructions using a wide variety of modeling methods, ensuring consistency. The Subaru data reveal the presence of a surrounding large scale structure with the major axis running approximately north-west south-east (NW-SE), aligned with the cluster and its brightest galaxy shapes, showing elongation with a \sim 2:1 axis ratio in the plane of the sky. Our full-lensing mass profile exhibits a shallow profile slope dlnΣ/dlnR\sim -1 at cluster outskirts (R>1Mpc/h), whereas the mass distribution excluding the NW-SE excess regions steepens further out, well described by the Navarro-Frenk-White form. Assuming a spherical halo, we obtain a virial mass M_{vir}=(1.1\pm 0.2\pm 0.1)\times 10^{15} M_{sun}/h and a halo concentration c_{vir} = 6.9\pm 1.0\pm 1.2 (\sim 5.7 when the central 50kpc/h is excluded), which falls in the range 4< <c> <7 of average c(M,z) predictions for relaxed clusters from recent Lambda cold dark matter simulations. Our full lensing results are found to be in agreement with X-ray mass measurements where the data overlap, and when combined with Chandra gas mass measurements, yield a cumulative gas mass fraction of 13.7^{+4.5}_{-3.0}% at 0.7Mpc/h (\approx 1.7r_{2500}), a typical value observed for high mass clusters.

preprint2012arXiv

CLASH: Precise New Constraints on the Mass Profile of Abell 2261

We precisely constrain the inner mass profile of Abell 2261 (z=0.225) for the first time and determine this cluster is not "over-concentrated" as found previously, implying a formation time in agreement with ΛCDM expectations. These results are based on strong lensing analyses of new 16-band HST imaging obtained as part of the Cluster Lensing and Supernova survey with Hubble (CLASH). Combining this with revised weak lensing analyses of Subaru wide field imaging with 5-band Subaru + KPNO photometry, we place tight new constraints on the halo virial mass M_vir = 2.2\pm0.2\times10^15 M\odot/h70 (within r \approx 3 Mpc/h70) and concentration c = 6.2 \pm 0.3 when assuming a spherical halo. This agrees broadly with average c(M,z) predictions from recent ΛCDM simulations which span 5 <~ <c> <~ 8. Our most significant systematic uncertainty is halo elongation along the line of sight. To estimate this, we also derive a mass profile based on archival Chandra X-ray observations and find it to be ~35% lower than our lensing-derived profile at r2500 ~ 600 kpc. Agreement can be achieved by a halo elongated with a ~2:1 axis ratio along our line of sight. For this elongated halo model, we find M_vir = 1.7\pm0.2\times10^15 M\odot/h70 and c_vir = 4.6\pm0.2, placing rough lower limits on these values. The need for halo elongation can be partially obviated by non-thermal pressure support and, perhaps entirely, by systematic errors in the X-ray mass measurements. We estimate the effect of background structures based on MMT/Hectospec spectroscopic redshifts and find these tend to lower Mvir further by ~7% and increase cvir by ~5%.

preprint2012arXiv

CLASH: Three Strongly Lensed Images of a Candidate z ~ 11 Galaxy

We present a candidate for the most distant galaxy known to date with a photometric redshift z = 10.7 +0.6 / -0.4 (95% confidence limits; with z < 9.5 galaxies of known types ruled out at 7.2-sigma). This J-dropout Lyman Break Galaxy, named MACS0647-JD, was discovered as part of the Cluster Lensing and Supernova survey with Hubble (CLASH). We observe three magnified images of this galaxy due to strong gravitational lensing by the galaxy cluster MACSJ0647.7+7015 at z = 0.591. The images are magnified by factors of ~8, 7, and 2, with the brighter two observed at ~26th magnitude AB (~0.15 uJy) in the WFC3/IR F160W filter (~1.4 - 1.7 um) where they are detected at >~ 12-sigma. All three images are also confidently detected at >~ 6-sigma in F140W (~1.2 - 1.6 um), dropping out of detection from 15 lower wavelength HST filters (~0.2 - 1.4 um), and lacking bright detections in Spitzer/IRAC 3.6um and 4.5um imaging (~3.2 - 5.0 um). We rule out a broad range of possible lower redshift interlopers, including some previously published as high redshift candidates. Our high redshift conclusion is more conservative than if we had neglected a Bayesian photometric redshift prior. Given CLASH observations of 17 high mass clusters to date, our discoveries of MACS0647-JD at z ~ 10.8 and MACS1149-JD1 at z ~ 9.6 are consistent with a lensed luminosity function extrapolated from lower redshifts. This would suggest that low luminosity galaxies could have reionized the universe. However given the significant uncertainties based on only two galaxies, we cannot yet rule out the sharp drop off in number counts at z >~ 10 suggested by field searches.

preprint2012arXiv

Cluster Lensing And Supernova survey with Hubble (CLASH): An Overview

The Cluster Lensing And Supernova survey with Hubble (CLASH) is a 524-orbit multi-cycle treasury program to use the gravitational lensing properties of 25 galaxy clusters to accurately constrain their mass distributions. The survey, described in detail in this paper, will definitively establish the degree of concentration of dark matter in the cluster cores, a key prediction of CDM. The CLASH cluster sample is larger and less biased than current samples of space-based imaging studies of clusters to similar depth, as we have minimized lensing-based selection that favors systems with overly dense cores. Specifically, twenty CLASH clusters are solely X-ray selected. The X-ray selected clusters are massive (kT > 5 keV; 5 - 30 x 10^14 M_solar) and, in most cases, dynamically relaxed. Five additional clusters are included for their lensing strength (Einstein radii > 35 arcsec at z_source = 2) to further quantify the lensing bias on concentration, to yield high resolution dark matter maps, and to optimize the likelihood of finding highly magnified high-redshift (z > 7) galaxies. The high magnification, in some cases, provides angular resolutions unobtainable with any current UVOIR facility and can yield z > 7 candidates bright enough for spectroscopic follow-up. A total of 16 broadband filters, spanning the near-UV to near-IR, are employed for each 20-orbit campaign on each cluster. These data are used to measure precise (sigma_phz < 0.02(1+z)) photometric redshifts for dozens of newly discovered multiply-lensed images per cluster. Observations of each cluster are spread over 8 epochs to enable a search, primarily in the parallel fields, for Type Ia supernovae at z > 1 to improve constraints on the time dependence of the dark energy equation of state and the evolution of such supernovae in an epoch when the universe is matter dominated.

preprint2012arXiv

Cluster Mass Profiles from a Bayesian Analysis of Weak Lensing Distortion and Magnification Measurements: Applications to Subaru Data

We directly construct model-independent mass profiles of galaxy clusters from combined weak-lensing distortion and magnification measurements within a Bayesian statistical framework,which allows for a full parameter-space extraction of the underlying signal. This method applies to the full range of radius outside the Einstein radius, and recovers the absolute mass normalization. We apply our method to deep Subaru imaging of five high-mass (>10^{15}M_{sun}) clusters, A1689, A1703, A370, Cl0024+17, and RXJ1347-11, to obtain accurate profiles to beyond the virial radius (r_{vir}). For each cluster the lens distortion and magnification data are shown to be consistent with each other, and the total signal-to-noise ratio of the combined measurements ranges from 13 to 24 per cluster. We form a model-independent mass profile from stacking the clusters, which is detected at 37σ out to R ~ 1.7r_{vir}. The projected logarithmic slope steepens from -1.01 \pm 0.09 at R ~ 0.1r_{vir} to -1.92 \pm 0.51 at R ~ 0.9r_{vir}. We also derive for each cluster inner strong-lensing based mass profiles from deep HST/ACS observations, which we show overlap well with the outer Subaru-based profiles and together are well described by a generalized form of the Navarro-Frenk-White profile, except for the ongoing merger RXJ1347-11, with modest variations in the central cusp slope (-dlnρ/dlnr < 0.9). The improvement here from adding the magnification measurements is significant, ~30% in terms of cluster mass profile measurements, compared with the lensing distortion signal.

preprint2012arXiv

Miscentring in Galaxy Clusters: Dark Matter to Brightest Cluster Galaxy Offsets in 10,000 SDSS Clusters

We characterise the typical offset between the Dark Matter (DM) projected centre and the Brightest Cluster Galaxy (BCG) in 10,000 SDSS clusters. To place constraints on the centre of DM, we use an automated strong-lensing analysis, mass-modelling technique which is based on the well-tested assumption that light traces mass. The cluster galaxies are modelled with a steep power-law, and the DM component is obtained by smoothing the galaxy distribution fitting a low-order 2D polynomial (via spline interpolation), while probing a whole range of polynomial degrees and galaxy power laws. We find that the offsets between the BCG and the peak of the smoothed light map representing the DM, Δ, are distributed equally around zero with no preferred direction, and are well described by a log-normal distribution with <log_{10}(Δ[h^{-1} Mpc])>=-1.895^{+0.003}_{-0.004}, and σ=0.501\pm0.004 (95% confidence levels), or <log_{10}(Δ[\arcsec])>=0.564\pm0.005, and σ=0.475\pm0.007. Some of the offsets originate in prior misidentifications of the BCG or other bright cluster members by the cluster finding algorithm, whose level we make an additional effort to assess, finding that ~10% of the clusters in the probed catalogue are likely to be misidentified, contributing to higher-end offsets in general agreement with previous studies. Our results constitute the first statistically-significant high-resolution distributions of DM-to-BCG offsets obtained in an observational analysis, and importantly show that there exists such a typical non-zero offset in the probed catalogue. The offsets show a weak positive correlation with redshift, so that higher separations are generally found for higher-z clusters in agreement with the hierarchical growth of structure, which in turn could help characterise the merger, relaxation and evolution history of clusters, in future studies. [ABRIDGED]

preprint2012arXiv

Probing ionizing radiation of L<~0.1L* star-forming galaxies at z>~3 with strong lensing

We show the effectiveness of strong lensing in the characterisation of Lyman continuum emission from faint L<~0.1L* star-forming galaxies at redshift >~ 3. Past observations of L>~L* galaxies at redshift >~3 have provided upper limits of the average escape fraction of ionising radiation of fesc~5%. Galaxies with relatively high fesc (>10%) seem to be particularly rare at these luminosities, there is therefore the need to explore fainter limits. Before the advent of giant ground based telescopes, one viable way to probe fesc down to 0.05-0.15L* is to exploit strong lensing magnification. This is investigated with Monte Carlo simulations that take into account the current observational capabilities. Adopting a lensing cross-section of 10 arcmin^2 within which the magnification is higher than 1 (achievable with about 4-5 galaxy clusters), with a U-band survey depth of 30(30.5) (AB, 1-sigma), it is possible to constrain fesc for z~3 star-forming galaxies down to 15(10)% at 3-sigma for L<0.15L* luminosities. This is particularly interesting if fesc increases at fainter luminosities, as predicted from various HI reionization scenarios and radiation transfer modelling. Ongoing observational programs on galaxy clusters are discussed and offer positive prospects for the future, even though from space the HST/WFC3 instrument represents the only option we have to investigate details of the spatial distribution of the Lyman continuum emission arising from z~2-4 galaxies.

preprint2012arXiv

Profiles of Dark Matter Velocity Anisotropy in Simulated Clusters

We report statistical results for dark matter (DM) velocity anisotropy, β, from a sample of some 6000 cluster-size halos (at redshift zero) identified in a ΛCDM hydrodynamical adaptive mesh refinement simulation performed with the Enzo code. These include profiles of β in clusters with different masses, relaxation states, and at several redshifts, modeled both as spherical and triaxial DM configurations. Specifically, although we find a large scatter in the DM velocity anisotropy profiles of different halos (across elliptical shells extending to at least ~$1.5 r_{vir}$), universal patterns are found when these are averaged over halo mass, redshift, and relaxation stage. These are characterized by a very small velocity anisotropy at the halo center, increasing outward to about 0.27 and leveling off at about $0.2 r_{vir}$. Indirect measurements of the DM velocity anisotropy fall on the upper end of the theoretically expected range. Though measured indirectly, the estimations are derived by using two different surrogate measurements - X-ray and galaxy dynamics. Current estimates of the DM velocity anisotropy are based on very small cluster sample. Increasing this sample will allow testing theoretical predictions, including the speculation that the decay of DM particles results in a large velocity boost. We also find, in accord with previous works, that halos are triaxial and likely to be more prolate when unrelaxed, whereas relaxed halos are more likely to be oblate. Our analysis does not indicate that there is significant correlation (found in some previous studies) between the radial density slope, γ, and β at large radii, $0.3 r_{vir} < r < r_{vir}$.

preprint2012arXiv

Spatially-resolved HST Grism Spectroscopy of a Lensed Emission Line Galaxy at z~1

We take advantage of gravitational lensing amplification by Abell 1689 (z=0.187) to undertake the first space-based census of emission line galaxies (ELGs) in the field of a massive lensing cluster. Forty-three ELGs are identified to a flux of i_775=27.3 via slitless grism spectroscopy. One ELG (at z=0.7895) is very bright owing to lensing magnification by a factor of ~4.5. Several Balmer emission lines detected from ground-based follow-up spectroscopy signal the onset of a major starburst for this low-mass galaxy (M_* = 2 x 10^9 solar masses) with a high specific star formation rate (~20 /Gyr). From the blue emission lines we measure a gas-phase oxygen abundance consistent with solar (12+log(O/H)=8.8 +/- 0.2). We break the continuous line-emitting region of this giant arc into seven ~1kpc bins (intrinsic size) and measure a variety of metallicity dependent line ratios. A weak trend of increasing metal fraction is seen toward the dynamical center of the galaxy. Interestingly, the metal line ratios in a region offset from the center by ~1kpc have a placement on the blue HII region excitation diagram with f([OIII])/f(Hbeta) and f([NeIII])/f(Hbeta) that can be fit by an AGN. This asymmetrical AGN-like behavior is interpreted as a product of shocks in the direction of the galaxy's extended tail, possibly instigated by a recent galaxy interaction.

preprint2012arXiv

The Universal Einstein Radius Distribution from 10,000 SDSS Clusters

We present results from strong-lens modelling of 10,000 SDSS clusters, to establish the universal distribution of Einstein radii. Detailed lensing analyses have shown that the inner mass distribution of clusters can be accurately modelled by assuming light traces mass, successfully uncovering large numbers of multiple-images. Approximate critical curves and the effective Einstein radius of each cluster can therefore be readily calculated, from the distribution of member galaxies and scaled by their luminosities. We use a subsample of 10 well-studied clusters covered by both SDSS and HST to calibrate and test this method, and show that an accurate determination of the Einstein radius and mass can be achieved by this approach "blindly", in an automated way, and without requiring multiple images as input. We present the results of the first 10,000 clusters analysed in the range $0.1<z<0.55$, and compare them to theoretical expectations. We find that for this all-sky representative sample the Einstein radius distribution is log-normal in shape, with $< Log(θ_{e}\arcsec)>=0.73^{+0.02}_{-0.03}$, $σ=0.316^{+0.004}_{-0.002}$, and with higher abundance of large $θ_{e}$ clusters than predicted by $Λ$CDM. We visually inspect each of the clusters with $θ_{e}>40 \arcsec$ ($z_{s}=2$) and find that $\sim20%$ are boosted by various projection effects detailed here, remaining with $\sim40$ real giant-lens candidates, with a maximum of $θ_{e}=69\pm12 \arcsec$ ($z_{s}=2$) for the most massive candidate, in agreement with semi-analytic calculations. The results of this work should be verified further when an extended calibration sample is available.

preprint2011arXiv

A Precise Cluster Mass Profile Averaged from the Highest-Quality Lensing Data

We outline our methods for obtaining high precision mass profiles, combining independent weak-lensing distortion, magnification, and strong-lensing measurements. For massive clusters the strong and weak lensing regimes contribute equal logarithmic coverage of the radial profile. The utility of high-quality data is limited by the cosmic noise from large scale structure along the line of sight. This noise is overcome when stacking clusters, as too are the effects of cluster asphericity and substructure, permitting a stringent test of theoretical models. We derive a mean radial mass profile of four similar mass clusters of high-quality HST and Subaru images, in the range R=40kpc/h to 2800kpc/h, where the inner radial boundary is sufficiently large to avoid smoothing from miscentering effects. The stacked mass profile is detected at 58-sigma significance over the entire radial range, with the contribution from the cosmic noise included. We show that the projected mass profile has a continuously steepening gradient out to beyond the virial radius, in remarkably good agreement with the standard Navarro-Frenk-White form predicted for the family of CDM-dominated halos in gravitational equilibrium. The central slope is constrained to lie in the range, -dlnρ/dln{r}=0.89^{+0.27}_{-0.39}. The mean concentration is c_{vir}=7.68^{+0.42}_{-0.40} (at a mean virial mass 1.54^{+0.11}_{-0.10}\times 10^{15} M_{sun}/h), which is high for relaxed, high-mass clusters, but consistent with LCDM when a sizable projection bias estimated from N-body simulations is considered. This possible tension will be more definitively explored with new cluster surveys, such as CLASH, LoCuSS, Subaru HSC, and XXM-XXL, to construct the c-M relation over a wider mass range.

preprint2011arXiv

A Weak Lensing Detection of the Cosmological Distance-Redshift Relation Behind Three Massive Clusters

The amplitude of weak lensing should increase with source distance, rising steeply behind a lens and saturating at high redshift, providing a model-independent means of measuring cosmic geometry. We measure the amplitude of weak lensing with redshift for three massive clusters, A370 (z=0.375), ZwCl0024+17 (z=0.395) and RXJ1347-11 (z=0.451), using deep, three-colour Subaru imaging. We define the depth of lensed populations with reference to the COSMOS and GOODS fields, providing a consistency check of photo-z estimates over a wide range of redshift and magnitude. The predicted distance-redshift relation is followed well for the deepest dataset, A370, for a wide range of cosmologies, and is consistent with less accurate data for the other two clusters. Scaling this result to a new survey of ~25 massive clusters should provide a useful cosmological constraint on w, complementing existing techniques, with distance measurements covering the untested redshift range, 1<z<5.

preprint2011arXiv

Triaxiality and non-thermal gas pressure in Abell 1689

Clusters of galaxies are uniquely important cosmological probes of the evolution of the large scale structure, whose diagnostic power depends quite significantly on the ability to reliably determine their masses. Clusters are typically modeled as spherical systems whose intracluster gas is in strict hydrostatic equilibrium (i.e., the equilibrium gas pressure is provided entirely by thermal pressure), with the gravitational field dominated by dark matter, assumptions that are only rough approximations. In fact, numerical simulations indicate that galaxy clusters are typically triaxial, rather than spherical, and that turbulent gas motions (induced during hierarchical merger events) provide an appreciable pressure component. Extending our previous work, we present results of a joint analysis of X-ray, weak and strong lensing measurements of Abell 1689. The quality of the data allows us to determine both the triaxial shape of the cluster and the level of non-thermal pressure that is required if the intracluster gas is in hydrostatic equilibrium. We find that the dark matter axis ratios are 1.24 +/- 0.13 and 2.02 +/- 0.01 on the plane of the sky and along the line of sight, respectively, and that about 20% of the pressure is non-thermal. Our treatment demonstrates that the dynamical properties of clusters can be determined in a (mostly) bias-free way, enhancing the use of clusters as more precise cosmological probes.

preprint2010arXiv

A Wide Area Survey for High-redshift Massive Galaxies. II. Near-infrared Spectroscopy of BzK-selected Massive Star-forming Galaxies

Results are presented from NIR spectroscopy of a sample of BzK-selected, massive star-forming galaxies (sBzKs) at 1.5<z<2.3 that were obtained with OHS/CISCO at Subaru and with SINFONI at VLT. Among the 28 sBzKs observed, Ha emission was detected in 14 objects, and for 11 of them the [NII]6583 was also measured. Multiwavelength photometry was also used to derive stellar masses and extinction parameters, whereas Ha and [NII] have allowed us to estimate SFR, metallicities, ionization mechanisms, and dynamical masses. In order to enforce agreement between SFRs from Ha with those derived from rest-frame UV and MIR, additional obscuration for the emission lines (that originate in HII regions) was required compared to the extinction derived from the slope of the UV continuum. We have also derived the stellar mass-metallicity relation, as well as the relation between stellar mass and specific SFR, and compared them to the results in other studies. At a given stellar mass, the sBzKs appear to have been already enriched to metallicities close to those of local star-forming galaxies of similar mass. The sBzKs presented here tend to have higher metallicities compared to those of UV-selected galaxies, indicating that NIR selected galaxies tend to be a chemically more evolved population. The sBzKs show specific SFRs that are systematically higher, by up to ~2 orders of magnitude, compared to those of local galaxies of the same mass. The empirical correlations between stellar mass and metallicity, and stellar mass and specific SFR are then compared with those of population synthesis models constructed either with the simple closed-box assumption, or within an infall scenario. Within the assumptions that are built-in such models, it appears that a short timescale for the star-formation (~100 Myr) and large initial gas mass appear to be required if one wants to reproduce both relations simultaneously.

preprint2010arXiv

Detailed Cluster Mass and Light profiles of A1703, A370 and RXJ1347-11 from Deep Subaru Imaging

Weak lensing work can be badly compromised by unlensed foreground and cluster members which dilute the true lensing signal. We show how the lensing amplitude in multi-colour space can be harnessed to securely separate cluster members from the foreground and background populations for three massive clusters, A1703 (z=0.258), A370 (z=0.375) and RXJ1347-11 (z=0.451) imaged with Subaru. The luminosity functions of these clusters when corrected for dilution, show similar faint-end slopes, α~= -1.0, with no marked faint-end upturn to our limit of M_R ~= -15.0, and only a mild radial gradient. In each case, the radial profile of the M/L ratio peaks at intermediate radius, ~=0.2r_{vir}, at a level of 300-500(M/L_R)_\odot, and then falls steadily towards ~100(M/L_R)_{\odot} at the virial radius, similar to the mean field level. This behaviour is likely due to the relative paucity of central late-type galaxies, whereas for the E/S0-sequence only a mild radial decline in M/L is found for each cluster. We discuss this behaviour in the context of detailed simulations where predictions for tidal stripping may now be tested accurately with observations.

preprint2010arXiv

Finding high-redshift dark stars with the James Webb Space Telescope

The first stars in the history of the Universe are likely to form in the dense central regions of 10^5-10^6 Msolar cold dark matter halos at z=10-50. The annihilation of dark matter particles in these environments may lead to the formation of so-called dark stars, which are predicted to be cooler, larger, more massive and potentially more long-lived than conventional population III stars. Here, we investigate the prospects of detecting high-redshift dark stars with the upcoming James Webb Space Telescope (JWST). We find that dark stars at z>6 are intrinsically too faint to be detected by JWST. However, by exploiting foreground galaxy clusters as gravitational telescopes, certain varieties of cool (Teff < 30000 K) dark stars should be within reach at redshifts up to z=10. If the lifetimes of dark stars are sufficiently long, many such objects may also congregate inside the first galaxies. We demonstrate that this could give rise to peculiar features in the integrated spectra of galaxies at high redshifts, provided that dark stars make up at least 1 percent of the total stellar mass in such objects.

preprint2010arXiv

Full Lensing Analysis of Abell 1703: Comparison of Independent Lens-Modelling Techniques

The inner mass-profile of the relaxed cluster Abell 1703 is analysed by two very different strong-lensing techniques applied to deep ACS and WFC3 imaging. Our parametric method has the accuracy required to reproduce the many sets of multiple images, based on the assumption that mass approximately traces light. We test this assumption with a fully non-parametric, adaptive grid method, with no knowledge of the galaxy distribution. Differences between the methods are seen on fine scales due to member galaxies which must be included in models designed to search for lensed images, but on the larger scale the general distribution of dark matter is in good agreement, with very similar radial mass profiles. We add undiluted weak-lensing measurements from deep multi-colour Subaru imaging to obtain a fully model-independent mass profile out to the virial radius and beyond. Consistency is found in the region of overlap between the weak and strong lensing, and the full mass profile is well-described by an NFW model of a concentration parameter, $c_{\rm vir}\simeq 7.15\pm0.5$ (and $M_{vir}\simeq 1.22\pm0.15 \times 10^{15}M_{\odot}/h$). Abell 1703 lies above the standard $c$--$M$ relation predicted for the standard $Λ$CDM model, similar to other massive relaxed clusters with accurately determined lensing-based profiles.

preprint2010arXiv

Quantifying the collisionless nature of dark matter and galaxies in A1689

We use extensive measurements of the cluster A1689 to assess the expected similarity in the dynamics of galaxies and dark matter (DM) in their motion as collisionless `particles' in the cluster gravitational potential. To do so we derive the radial profile of the specific kinetic energy of the cluster galaxies from the Jeans equation and observational data. Assuming that the specific kinetic energies of galaxies and DM are roughly equal, we obtain the mean value of the DM velocity anisotropy parameter, and the DM density profile. Since this deduced profile has a scale radius that is higher than inferred from lensing observations, we tested the validity of the assumption by repeating the analysis using results of simulations for the profile of the DM velocity anisotropy. Results of both analyses indicate a significant difference between the kinematics of galaxies and DM within $r \lesssim 0.3r_{\rm vir}$. This finding is reflected also in the shape of the galaxy number density profile, which flattens markedly with respect to the steadily rising DM profile at small radii. Thus, $r \sim 0.3r_{\rm vir}$ seems to be a transition region interior to which collisional effects significantly modify the dynamical properties of the galaxy population with respect to those of DM in A1689

preprint2010arXiv

Strong-Lensing Analysis of a Complete Sample of 12 MACS Clusters at z>0.5: Mass Models and Einstein Radii

We present the results of a strong-lensing analysis of a complete sample of 12 very luminous X-ray clusters at $z>0.5$ using HST/ACS images. Our modelling technique has uncovered some of the largest known critical curves outlined by many accurately-predicted sets of multiple images. The distribution of Einstein radii has a median value of $\simeq28\arcsec$ (for a source redshift of $z_{s}\sim2$), twice as large as other lower-$z$ samples, and extends to $55\arcsec$ for MACS J0717.5+3745, with an impressive enclosed Einstein mass of $7.4\times10^{14} M_{\odot}$. We find that 9 clusters cover a very large area ($>2.5 \sq \arcmin$) of high magnification ($μ> \times10$) for a source redshift of $z_{s}\sim8$, providing primary targets for accessing the first stars and galaxies. We compare our results with theoretical predictions of the standard $Λ$CDM model which we show systematically fall short of our measured Einstein radii by a factor of $\simeq1.4$, after accounting for the effect of lensing projection. Nevertheless, a revised analysis once arc redshifts become available, and similar analyses of larger samples, are needed in order to establish more precisely the level of discrepancy with $Λ$CDM predictions.

preprint2010arXiv

Strong-Lensing Analysis of MS 1358.4+6245: New Multiple Images and Implications for the Well-Resolved z=4.92 Galaxy

We present a strong-lensing analysis of the galaxy cluster MS 1358.4+6245 ($z=0.33$), in deep 6-band ACS/HST imaging. In addition to the well-studied system at $z=4.92$, our modelling method uncovers 19 new multiply-lensed images so that a total of 23 images and their redshifts are used to accurately constrain the inner mass distribution. We derive a relatively shallow inner mass profile, $d\log Σ/d\log r\simeq -0.33 \pm0.05$ ($r<200$ kpc), with a much higher magnification than estimated previously by models constrained only by the $z=4.92$ system. Using these many new images we can apply a non-parametric adaptive-grid method, which also yields a shallow mass profile without prior assumptions, strengthening our conclusions. The total magnification of the $z_s=4.92$ galaxy is high, about a $\sim100\times$ over its four images, so that the inferred source size, luminosity and star-formation rate are about $\sim5\times$ smaller than previous estimates, corresponding to a dwarf-sized galaxy of radius $\simeq1$ kpc. A detailed image of the interior morphology of the source is generated with a high effective resolution of only $\simeq$50 pc, thanks to the high magnification and to the declining angular diameter distance above $z\sim1.5$ for the standard cosmology, so that this image apparently represents the best resolved object known at high redshift.

preprint2010arXiv

The Highest Resolution Mass Map of Galaxy Cluster Substructure To Date Without Assuming Light Traces Mass: LensPerfect Analysis of Abell 1689

We present a strong lensing mass model of Abell 1689 which resolves substructures ~25 kpc across (including about ten individual galaxy subhalos) within the central ~400 kpc diameter. We achieve this resolution by perfectly reproducing the observed (strongly lensed) input positions of 168 multiple images of 55 knots residing within 135 images of 42 galaxies. Our model makes no assumptions about light tracing mass, yet we reproduce the brightest visible structures with some slight deviations. A1689 remains one of the strongest known lenses on the sky, with an Einstein radius of RE = 47.0" +/- 1.2" (143 +3/-4 kpc) for a lensed source at zs = 2. We find a single NFW or Sersic prole yields a good fit simultaneously (with only slight tension) to both our strong lensing (SL) mass model and published weak lensing (WL) measurements at larger radius (out to the virial radius). According to this NFW fit, A1689 has a mass of Mvir = 2.0 +0.5/-0.3 x 10^15 Msun / h70 (M200 = 1.8 +0.4/-0.3 x 10^15 Msun / h70) within the virial radius rvir = 3.0 +/- 0.2 Mpc / h70 (r200 = 2.4 +0.1/-0.2 Mpc / h70), and a central concentration cvir = 11.5 +1.5/-1.4 (c200 = 9.2 +/- 1.2). Our SL model prefers slightly higher concentrations than previous SL models, bringing our SL+WL constraints in line with other recent derivations. Our results support those of previous studies which find A1689 has either an anomalously large concentration or significant extra mass along the line of sight (perhaps in part due to triaxiality). If clusters are generally found to have higher concentrations than realized in simulations, this could indicate they formed earlier, perhaps as a result of early dark energy.

preprint2010arXiv

The Mass Structure of the Galaxy Cluster Cl0024+1654 from a Full Lensing Analysis of Joint Subaru and ACS/NIC3 Observations

We derive an accurate mass distribution of the rich galaxy cluster Cl0024+1654 (z=0.395) based on deep Subaru BR_{c}z' imaging and our recent comprehensive strong lensing analysis of HST/ACS/NIC3 observations. We obtain the weak lensing distortion and magnification of undilted samples of red and blue background galaxies by carefully combining all color and positional information. Unlike previous work, the weak and strong lensing are in excellent agreement where the data overlap. The joint mass profile continuously steepens out to the virial radius with only a minor contribution \sim 10% in the mass from known subcluster at a projected distance of \sim 700kpc/h. The projected mass distribution for the entire cluster is well fitted with a single Navarro-Frenk-White model with a virial mass, M_{vir} = (1.2 \pm 0.2) \times 10^{15} M_{sun}/h, and a concentration, c_{vir} = 9.2^{+1.4}_{-1.2}. This model fit is fully consistent with the depletion of the red background counts, providing independent confirmation. Careful examination and interpretation of X-ray and dynamical data strongly suggest that this cluster system is in a post collision state, which we show is consistent with our well-defined mass profile for a major merger occurring along the line of sight, viewed approximately 2-3Gyr after impact when the gravitational potential has had time to relax in the center, before the gas has recovered and before the outskirts are fully virialized. Finally, our full lensing analysis provides a model-independent constraint of M_{2D}(<r_{vir}) = (1.4 \pm 0.3) \times 10^{15} M_{sun}/h for the projected mass of the whole system, including any currently unbound material beyond the virial radius, which can constrain the sum of the two pre-merger cluster masses when designing simulations to explore this system.

preprint2009arXiv

Dynamical Study of A1689 from Wide-Field VLT/VIMOS Spectroscopy: Mass Profile, Concentration Parameter, and Velocity Anisotropy

We examine the dynamics structure of the rich cluster A1689, combining VLT/VIMOS spectroscopy with Subaru/Suprime-Cam imaging. The radial velocity distribution of $\sim 500$ cluster members is bounded by a pair of clearly defined velocity caustics, with a maximum amplitude of $\sim|4000|$ km/s at $\simeq$ 300 h$^{-1}$ kpc, beyond which the amplitude steadily declines, approaching zero velocity at a limiting radius of $\sim$ 2 h$^{-1}$ Mpc. We derive the 3D velocity anisotropy and galaxy number density profiles using a model-independent method to solve the Jeans equation, simultaneously incorporating the observed velocity dispersion profile, the galaxy counts from deep Subaru imaging, and our previously derived cluster mass profile from a joint lensing and X-ray analysis. The velocity anisotropy is found to be predominantly radial at large radius, becoming increasingly tangential towards the center, in accord with expectations. We also analyze the galaxy data independently of our previous analysis using two different methods: The first is based on a solution of the Jeans equation assuming an NFW form for the mass distribution, whereas in the second method the caustic amplitude is used to determine the escape velocity. The cluster virial mass derived by both of these dynamical methods is in good agreement with results from our earlier lensing and X-ray analysis. We also confirm the high NFW concentration parameter, with results from both methods combined to yield $c_{\rm vir}>13$ (1$σ$). The inferred virial radius is consistent with the limiting radius where the caustics approach zero velocity and where the counts of cluster members drop off, suggesting that infall onto A1689 is currently not significant.

preprint2009arXiv

Mass and Hot Baryons in Massive Galaxy Clusters from Subaru Weak Lensing and AMiBA SZE Observations

We present a multiwavelength analysis of a sample of four hot (T_X>8keV) X-ray galaxy clusters (A1689, A2261, A2142, and A2390) using joint AMiBA Sunyaev-Zel'dovich effect (SZE) and Subaru weak lensing observations, combined with published X-ray temperatures, to examine the distribution of mass and the intracluster medium (ICM) in massive cluster environments. Our observations show that A2261 is very similar to A1689 in terms of lensing properties. Many tangential arcs are visible around A2261, with an effective Einstein radius \sim 40 arcsec (at z \sim 1.5), which when combined with our weak lensing measurements implies a mass profile well fitted by an NFW model with a high concentration c_{vir} \sim 10, similar to A1689 and to other massive clusters. The cluster A2142 shows complex mass substructure, and displays a shallower profile (c_{vir} \sim 5), consistent with detailed X-ray observations which imply recent interaction. The AMiBA map of A2142 exhibits an SZE feature associated with mass substructure lying ahead of the sharp north-west edge of the X-ray core suggesting a pressure increase in the ICM. For A2390 we obtain highly elliptical mass and ICM distributions at all radii, consistent with other X-ray and strong lensing work. Our cluster gas fraction measurements, free from the hydrostatic equilibrium assumption, are overall in good agreement with published X-ray and SZE observations, with the sample-averaged gas fraction of <f_{gas}> = 0.133 \pm 0.027, for our sample <M_{vir}> = (1.2 \pm 0.1) \times 10^{15} M_{sun} h^{-1}. When compared to the cosmic baryon fraction f_b = Ω_b/Ω_m constrained by the WMAP 5-year data, this indicates <f_{gas}>/f_b = 0.78 \pm 0.16, i.e., (22 \pm 16)% of the baryons are missing from the hot phase of clusters.

preprint2009arXiv

New Multiply-Lensed Galaxies Identified in ACS/NIC3 Observations of Cl0024+1654 Using an Improved Mass Model

We present an improved strong-lensing analysis of Cl0024+1654 ($z$=0.39) using deep HST/ACS/NIC3 images, based on 33 multiply-lensed images of 11 background galaxies. These are found with a model that assumes mass approximately traces light, with a low order expansion to allow for flexibility on large scales. The model is constrained initially by the well known 5-image system ($z$=1.675) and refined as new multiply-lensed systems are identified using the model. Photometric redshifts of these new systems are then used to constrain better the mass profile by adopting the standard cosmological relation between redshift and lensing distance. Our model requires only 6 free parameters to describe well all positional and redshift data. The resulting inner mass profile has a slope of $d\log M/d\log r\simeq -0.55$, consistent with new weak-lensing measurements where the data overlap, at $r\simeq200$ kpc/$h_{70}$. The combined profile is well fitted by a high concentration NFW mass profile, $C_{\rm vir}\sim 8.6\pm1.6$, similar to other well studied clusters, but larger than predicted with standard $Λ$CDM. A well defined radial critical curve is generated by the model and is clearly observed at $r \simeq 12\arcsec$, outlined by elongated images pointing towards the centre of mass. The relative fluxes of the multiply-lensed images are found to agree well with the modelled magnifications, providing an independent consistency check.

preprint2009arXiv

The Largest Gravitational Lens: MACS J0717.5+3745 (z=0.546)

We identify 13 sets of multiply-lensed galaxies around MACS J0717.5+3745 ($z=0.546$), outlining a very large tangential critical curve of major axis $\sim2.8\arcmin$, filling the field of HST/ACS. The equivalent circular Einstein radius is $θ_{e}= 55 \pm 3\arcsec$ (at an estimated source redshift of $z_{s}\sim2.5$), corresponding to $r_e\simeq 350\pm 20 kpc$ at the cluster redshift, nearly three times greater than that of A1689 ($r_e\simeq 140 kpc$ for $z_{s}=2.5$). The mass enclosed by this critical curve is very large, $7.4\pm 0.5 \times 10^{14}M_{\odot}$ and only weakly model dependent, with a relatively shallow mass profile within $r<250 kpc$, reflecting the unrelaxed appearance of this cluster. This shallow profile generates a much higher level of magnification than the well known relaxed lensing clusters of higher concentration, so that the area of sky exceeding a magnification of $>10\times$, is $\simeq 3.5\sq\arcmin$ for sources with $z\simeq 8$, making MACS J0717.5+3745 a compelling target for accessing faint objects at high redshift. We calculate that only one such cluster, with $θ_{e}\ge 55\arcsec$, is predicted within $\sim 10^7$ Universes with $z\ge 0.55$, corresponding to a virial mass $\ge 3\times 10^{15} M_{\odot}$, for the standard $ΛCDM$ (WMAP5 parameters with $2σ$ uncertainties).

preprint2008arXiv

Comparison of Cluster Lensing Profiles with Lambda CDM Predictions

We derive lens distortion and magnification profiles of four well known clusters observed with Subaru. Each cluster is very well fitted by the general form predicted for Cold Dark Matter (CDM) dominated halos, with good consistency found between the independent distortion and magnification measurements. The inferred level of mass concentration is surprisingly high, 8<c_{vir}<15 (<c_{vir}> = 10.4 \pm 0.9), compared to the relatively shallow profiles predicted by the Lambda CDM model, c_{vir}=5.1 \pm 1.1 (for <M_{vir}> =1.25\times 10^{15}M_{\odot}/h). This represents a 4sigma discrepancy, and includes the relatively modest effects of projection bias and profile evolution derived from N-body simulations, which oppose each other with little residual effect. In the context of CDM based cosmologies, this discrepancy implies clusters collapse earlier (z\geq 1) than predicted (z<0.5), when the Universe was correspondingly denser.

preprint2007arXiv

Using Weak Lensing Dilution to Improve Measurements of the Luminous and Dark Matter in A1689

The E/SO sequence of a cluster defines a boundary redward of which a reliable weak lensing signal can be obtained from background galaxies, uncontaminated by cluster members. For bluer colors, both background and cluster members are present, reducing the distortion signal by the proportion of unlensed cluster members. In deep Subaru and HST/ACS images of A1689 the tangential distortion of galaxies with bluer colors falls rapidly toward the cluster center relative to the lensing signal of the red background. We use this dilution effect to derive the cluster light profile and luminosity function to large radius, with the advantage that no subtraction of far-field background counts is required. The light profile declines smoothly to the limit of the data, r<2Mpc/h, with a constant slope, dlog(L)/dlog(r)=-1.12+-0.06, unlike the lensing mass profile which steepens continuously with radius, so that M/L peaks at an intermediate radius, ~100kpc/h. A flatter behavior is found for the more physically meaningful ratio of dark-matter to stellar-matter, when accounting for the color-mass relation of cluster members. The cluster luminosity function has a flat slope, alpha=-1.05+-0.07, independent of radius and with no faint upturn to M_i'<-12. We establish that the very bluest objects are negligibly contaminated by the cluster V-i'<0.2, because their distortion profile rises towards the center following the red background, but offset higher by ~20%. This larger amplitude is consistent with the greater estimated depth of the faint blue galaxies, z~=2.0 compared to z~=0.85 for the red background, a purely geometric effect related to cosmological parameters. Finally, we improve upon our earlier mass profile by combining both the red and blue background populations, clearly excluding low concentration CDM profiles.

preprint2006arXiv

The effect of FIR emission from SDSS galaxies on the SFD Galactic extinction map

We compare the most successful and widely used map of Galactic dust extinction, provided by Schlegel, Finkbeiner & Davis (1998; hereafter SFD), to the galaxy number counts in the Sloan Digital Sky Survey (SDSS) photometric/spectroscopic DR4 sample. We divide the SDSS survey area into 69 disjoint subregions according to the dust extinction provided by SFD and compare the surface number density of galaxies in each subregion. As expected, the galaxy surface number density decreases with increasing extinction but only for SFD extinction values above about 0.1 to 0.2 magnitudes (depending on the band). At lower values of the SFD extinction, we find that the sky surface density of galaxies increases with increasing extinction, precisely the opposite of the effect expected from Galactic dust. We suggest that the far infrared (FIR) brightness of the sky in regions of true low dust extinction is significantly ``contaminated'' by the FIR emission from background galaxies. We show that such an explanation is both qualitatively and quantitatively consistent with the available data. Based on this interpretation we conclude that systematic errors in the SFD extinction map due to extragalactic FIR emission are quite small, of order hundredths of a magnitude, but nevertheless statistically detectable. (Abridged)

preprint2005arXiv

The ALHAMBRA Survey: For a systematic Study of Cosmic Evolution

ALHAMBRA is a project to gather data to sample a fraction of the Universe with enough precision to follow the evolution of its content and properties with z, a kind of Cosmic Tomography. It is defined as a 4 square degrees photometric survey with 20 contiguous, equal width, medium band filters covering from 3500 to 9700 A, plus the JHKs NIR bands. It was optimized to get (for a fixed observing time) the maximum number of objects with accurate classification and redshift and to be sensitive to relatively faint emission features. The observations will be carried out with the 3.5m telescope in Calar Alto (Spain). We intend to reach the limit AB=25 in all the optical filters from the bluest to 8300 A, and from 24.7 to 23.4 for the remainder. The expected limit in the NIR is fixed at Ks=20, H=21, J=22. The homogeneous and contiguous spectral coverage will result in several hundred thousand objects with accurate SED identification and z-values. The accuracy of the survey will allow to make detailed studies, including galactic, without the need for any further follow-up. Indeed, it will provide exciting targets for large telescopes, the GTC in particular.

preprint2004arXiv

Probing Halos of Galaxies at Very Large Radii using Background QSOs

(ABRIDGED) Gaseous halos of nine nearby galaxies (with redshifts cz < 6000 km/s) were probed at large galactocentric radii using background quasars observed with HST GHRS and STIS. The projected quasar-galaxy separations range from 55 to 387 h^{-1}_{75} kpc. Lyαabsorption lines were successfully detected in the spectra of five quasars, at impact parameters of up to 170 h^{-1}_{75} kpc from the center of the nearby galaxy, and in each case at wavelengths consistent with the galaxy's redshift. Our observations include the lowest redshift Lyαlines detected to date. HI velocity fields were obtained at the VLA for three of the galaxies in our sample (and in one case was available from the literature), to derive their rotation curves. When comparing the inner rotation curves of the galaxies with the velocity at large radius provided by the Lyαline it is apparent that it is very difficult to explain the observed Lyαvelocity as due to gas in an extended rotating disk. In most cases one would need to invoke large warps in the outer gas disks and also thick gas disks in order to reconcile the observed velocities with the predicted ones. In light of these results we conclude that in a typical galaxy there is no longer detectable atomic gas corotating in an extended disk at radii > 35 α^{-1}, where α^{-1} is the stellar disk exponential scale-length. The cosmic web is the most likely origin for the detected Lyαlines.

preprint2004arXiv

The Surprisingly Steep Mass Profile of Abell 1689, from a Lensing Analysis of Subaru Images

Subaru observations of A1689 (z=0.183) are used to derive an accurate, model-independent mass profile for the entire cluster, r<2 Mpc/h, by combining magnification bias and distortion measurements. The projected mass profile steepens quickly with increasing radius, falling away to zero at r~1.0 Mpc/h, well short of the anticipated virial radius. Our profile accurately matches onto the inner profile, r<200 kpc/h, derived from deep HST/ACS images. The combined ACS and Subaru information is well fitted by an NFW profile with virial mass, (1.93 \pm 0.20)10^15 M_sun, and surprisingly high concentration, c_vir=13.7^{+1.4}_{-1.1}, significantly larger than theoretically expected (c_vir~4), corresponding to a relatively steep overall profile. A slightly better fit is achieved with a steep power-law model that has its 2D logarithmic slope -3 and core radius theta_c~1.7' (r_c~210 kpc/h), whereas an isothermal profile is strongly rejected. These results are based on a reliable sample of background galaxies selected to be redder than the cluster E/S0 sequence. By including the faint blue galaxy population a much smaller distortion signal is found, demonstrating that blue cluster members significantly dilute the true signal for r~400 kpc/h. This contamination is likely to affect most weak lensing results to date.

preprint1998arXiv

Cloning Hubble Deep Fields: A Model-Independent Measurement of Galaxy Evolution

We present a model-independent method of quantifying galaxy evolution in high- resolution images, which we apply to the Hubble Deep Field (HDF). Our procedure is to k-correct the pixels belonging to the images of a complete set of bright galaxies and then to replicate each galaxy image to higher redshift by the product of its space density, 1/V_{max}, and the cosmological volume. The set of bright galaxies is itself selected from the HDF because presently the HDF provides the highest quality UV images of a redshift-complete sample of galaxies (31 galaxies with I<21.9, \bar{z}=0.5, and for which V/V_{max} is spread fairly). These galaxies are bright enough to permit accurate pixel-by-pixel k-corrections into the restframe UV (\sim 2000 A). We match the shot noise, spatial sampling and PSF smoothing of the HDF, resulting in entirely empirical and parameter free ``no-evolution'' deep fields of galaxies for direct comparison with the HDF. We obtain the following results. Faint HDF galaxies (I>24) are much smaller, more numerous, and less regular than our ``no-evolution'' extrapolation, for any relevant geometry. A higher proportion of HDF galaxies ``dropout'' in both U and B, indicating that some galaxies were brighter at higher redshifts than our ``cloned'' z\sim0.5 population. By simple image transformations we demonstrate that bolometric luminosity evolution generates galaxies which are too large and the contribution of any evolving dwarf population is uninterestingly small. A plausible fit is provided by `mass-conserving' density-evolution, consistent with hierarchical growth of small-scale structure. Finally, we show the potential for improvement using the Advanced Camera, with its superior UV and optical performance.

preprint1996arXiv

Autofib Redshift Survey: II -- The Evolution of the Galaxy Luminosity Function by Spectral Type

We determine the evolution of the galaxy luminosity function (LF) as a function of spectral type using the Autofib redshift survey, a compendium of over 1700 galaxy redshifts in various magnitude-limited samples spanning b_J=11.5-24.0. To carry out this analysis we have developed a cross-correlation technique which classifies faint galaxy spectra into one of six types based on local galaxy templates. Tests and simulations show that this technique yields classifications correct to within one type for more than 90% of the galaxies in our sample. We have also developed extensions of the step-wise maximum likelihood method and the STY parametric method for estimating LFs which are applicable to recovering an evolving LF from multiple samples. We find that: (i) The spectra and LF of E/S0 galaxies show no appreciable evolution out to at least z ~ 0.5. (ii) Early-type spirals show modest evolution, characterised by a gradual steepening of the faint end of their LF with redshift. (iii) Out to z ~ 0.5, the overall evolution of the galaxy population is dominated by changes seen in late-type spirals. The characteristic luminosity (L^*) of these galaxies appears to brighten with redshift and there are signs of strong density evolution (a rapid increase in ϕ^*). These effects appear to be luminosity dependent so that the LF steepens at higher redshift. These trends are accompanied by a steep increase in the median [OII] equivalent width, implying a rapid increase in the star-formation rate with redshift at fixed luminosity---a given star-formation rate is found at higher redshift in galaxies of higher luminosity. We find that these conclusions are robust with respect to spectral classification errors and the luminosity function estimator.

preprint1995arXiv

Autofib Redshift Survey - I. Evolution of the Galaxy Luminosity Function

We present a detailed determination of the restframe B-band galaxy luminosity function (LF) as a function of redshift and star formation activity from z=0 to z=0.75. The dataset used for this purpose is a combined sample of over 1700 redshifts spanning a wide range in apparent magnitude, 11.5<B<24.0, which we term the Autofib Redshift Survey. The sample includes various earlier magnitude-limited surveys constructed by our team as well as a new survey of 1026 redshifts measured for galaxies at intermediate magnitudes. The large range in apparent magnitude sampled allows us to investigate both the nature of the LF at low redshift (z<0.1) and possible evolution in its shape to z=0.75. We find that earlier bright surveys have underestimated the absolute normalisation of the LF. Because the shape of the local LF does not change with the survey apparent magnitude limit, it seems unlikely that the local deficiency arises from an underestimated population of low luminosity galaxies. Furthermore, surface brightness losses cannot be significant unless they conspire to retain the LF shape over a variety of detection thresholds. Our data directly demonstrates that the B-band LF evolves with redshift. This evolution is best represented as a steepening of the faint-end slope of the LF, from alpha=-1.1 at low redshift to alpha=-1.5 at z=0.5. Using [OII] emission as an indicator of star formation activity, we show that the LF of quiescent galaxies has remained largely unchanged since z=0.5, whereas the luminosity density of star-forming galaxies has declined by nearly a factor of 2. The steepening of the overall LF with lookback time is a direct consequence of the increasing space density of blue star-forming galaxies at moderate redshifts.

preprint1995arXiv

Mass Distributions of Clusters Using Gravitational Magnification

Lensing in the context of rich clusters is normally quantified from small image distortions, yielding a relative mass distribution in the limit of weak lensing. Here we show the magnification effect of lensing can also be mapped over a cluster, resulting in absolute mass determinations for the weak limit. Furthermore, given both magnification and distortion measurements, the mass distribution may be constrained in the strong regime. Methods for obtaining the magnification using spectroscopic and/or photometric information are discussed, for object detection within a fixed isophote or to a given flux limit. A map of the magnification around A1689 is constructed from the observed depletion of background red galaxy counts.