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.

The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) will soon start to scan thousands of square degrees of the northern extragalactic sky with a unique set of $56$ optical filters from a dedicated $2.55$m telescope, JST, at the Javalambre Astrophysical Observatory. Before the arrival of the final instrument (a 1.2 Gpixels, 4.2deg$^2$ field-of-view camera), the JST was equipped with an interim camera (JPAS-Pathfinder), composed of one CCD with a 0.3deg$^2$ field-of-view and resolution of 0.23 arcsec pixel$^{-1}$. To demonstrate the scientific potential of J-PAS, with the JPAS-Pathfinder camera we carried out a survey on the AEGIS field (along the Extended Groth Strip), dubbed miniJPAS. We observed a total of $\sim 1$ deg$^2$, with the $56$ J-PAS filters, which include $54$ narrow band (NB, $\rm{FWHM} \sim 145$Angstrom) and two broader filters extending to the UV and the near-infrared, complemented by the $u,g,r,i$ SDSS broad band (BB) filters. In this paper we present the miniJPAS data set, the details of the catalogues and data access, and illustrate the scientific potential of our multi-band data. The data surpass the target depths originally planned for J-PAS, reaching $\rm{mag}_{\rm {AB}}$ between $\sim 22$ and $23.5$ for the NB filters and up to $24$ for the BB filters ($5\sigma$ in a $3$~arcsec aperture). The miniJPAS primary catalogue contains more than $64,000$ sources extracted in the $r$ detection band with forced photometry in all other bands. We estimate the catalogue to be complete up to $r=23.6$ for point-like sources and up to $r=22.7$ for extended sources. Photometric redshifts reach subpercent precision for all sources up to $r=22.5$, and a precision of $\sim 0.3$% for about half of the sample. (Abridged)

We examine the massive colliding cluster El Gordo, one of the most massive clusters at high redshift. We use a free-form lensing reconstruction method that avoids making assumptions about the mass distribution. We use data from the RELICS program and identify new multiply lensed system candidates. The new set of constraints and free-form method provides a new independent mass estimate of this intriguing colliding cluster. Our results are found to be consistent with earlier parametric models, indirectly confirming the assumptions made in earlier work. By fitting a double gNFW profile to the lens model, and extrapolating to the virial radius, we infer a total mass for the cluster of $M_{200c}=(1.08^{+0.65}_{-0.12})\times10^{15}$M$_{\odot}$. We estimate the uncertainty in the mass due to errors in the photometric redshifts, and discuss the uncertainty in the inferred virial mass due to the extrapolation from the lens model. We also find in our lens map a mass overdensity corresponding to the large cometary tail of hot gas, reinforcing its interpretation as a large tidal feature predicted by hydrodynamical simulations that mimic El Gordo. Finally, we discuss the observed relation between the plasma and the mass map, finding that the peak in the projected mass map may be associated with a large concentration of colder gas, exhibiting possible star formation. El Gordo is one of the first clusters that will be observed with JWST, which is expected to unveil new high redshift lensed galaxies around this interesting cluster, and provide a more accurate estimation of its mass.

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.

Strong lensing of gravitational waves is more likely for distant sources but predicted event rates are highly uncertain with many astrophysical origins proposed. Here we open a new avenue to estimate the event rate of strongly lensed systems by exploring the amplitude of the stochastic gravitational wave background (SGWB). This method can provide a direct upper bound on the high redshift binary coalescing rates, which can be translated into an upper bound on the expected rate of strongly lensed systems. We show that from the ongoing analysis of the Laser Interferometer Gravitational-wave Observatory (LIGO)-Virgo and in the future from the LIGO-Virgo design sensitivity stringent bounds on the lensing event rate can be imposed using the SGWB signal. Combining measurements of loud gravitational wave events with an unresolved stochastic background detection will improve estimates of the numbers of lensed events at high redshift. The proposed method is going to play a crucial in understanding the population of lensed and unlensed systems from gravitational wave observations.

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.

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.

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_\psi \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.

Dark matter as a Bose-Einstein condensate, such as the axionic scalar field particles of String Theory, can explain the coldness of dark matter on large scales. Pioneering simulations in this context predict a rich wave-like structure, with a ground state soliton core in every galaxy surrounded by a halo of excited states that interfere on the de Broglie scale. This de Broglie scale is largest for low mass galaxies as momentum is lower, providing a simple explanation for the wide cores of dwarf spheroidal galaxies. Here we extend these "wave dark matter" ($\psi$DM) predictions to the newly discovered class of "Ultra Diffuse Galaxies" (UDG) that resemble dwarf spheroidal galaxies but with more extended stellar profiles. Currently the best studied example, DF44, has a uniform velocity dispersion of $\simeq 33$km/s, extending to a least 3 kpc, that we show is reproduced by our $\psi$DM simulations with a soliton radius of $\simeq 0.5$ kpc. In the $\psi$DM context, we show the relatively flat dispersion profile of DF44 lies between massive galaxies with compact dense solitons, as may be present in the Milky Way on a scale of 100pc and lower mass galaxies where the velocity dispersion declines centrally within a wide, low density soliton, like Antlia II, of radius 3Kpc. In contrast, standard CDM requires excessive tangential stellar motions predominate to counter the inherent central cusp which would otherwise enhance the central velocity dispersion.

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.

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.

The primary and secondary masses of the binary black holes (BBH) reported by LIGO/Virgo are correlated with a narrow dispersion that appears to increase in proportion to mass. The mean binary mass ratio $1.45\pm0.07$ we show is consistent with pairs drawn randomly from the mass distribution of black holes in our Galaxy. However, BBH masses are concentrated around $\simeq 30M_\odot$, whereas black holes in our Galaxy peak at $\simeq 10M_\odot$. This mass difference can be reconciled by gravitational lensing magnification which allows distant events to be detected with typically $z\simeq 2$, so the waveform is reduced in frequency by $1+z$, and hence the measured chirp masses appear 3 times larger than their intrinsic values. This redshift enhancement also accounts for the dispersion of primary and secondary masses, both of which should increase as $1+z$, thereby appearing to scale with mass, in agreement with the data. Thus the BBH component masses provide independent support for lensing, implying most high chirp mass events have intrinsic masses like the stellar mass black holes in our Galaxy, coalescing at $z>1$, with only two low mass BBH detections, of $\simeq 10M_\odot$ as expected for unlensed events in the local Universe, $z\simeq 0.1$. This lensing solution requires a rapidly declining BBH event rate below $z<1$, which together with the observed absence of BBH spin suggests most events originate within young globular clusters at $z>1$, via efficient binary capture of stellar mass black holes with randomly oriented spins.

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.

We carry out a detailed strong lensing analysis of a sub-sample of eight galaxy clusters of the CLASH survey, in the redshift range of $ z_{\rm cluster} = [0.23-0.59]$, using extensive spectroscopic information, primarily MUSE archival data complemented with CLASH-VLT redshift measurements. Different models are tested in each cluster depending on the complexity of its mass distribution and on the number of detected multiple images. Four clusters show more than five spectroscopically confirmed multiple image families. In this sample, we do not make use of families that are only photometrically identified, in order to reduce model degeneracies and systematics due to the potential misidentifications of some multiple images. We present spectroscopic confirmation of 27 multiply lensed sources, with no previous spectroscopic measurements, spanning over the redshift range of $z_{\rm src}=[0.7-6.1]$. Moreover, we confirm an average of $48$ galaxy members in the core of each cluster, thanks to the high efficiency and large field of view of MUSE. Despite having different properties (i.e., number of mass components, total mass, redshift, etc), the projected total mass and mass density profiles of all clusters have very similar shapes, when rescaled by independent measurements of $M_{200c}$ and $R_{200c}$. Specifically, we measure the mean value of the projected total mass of our cluster sample within 10 (20)% of $R_{200c}$ to be 0.13 (0.32) of $M_{200c}$, with a remarkably small scatter of 5 (6)%. Furthermore, the large number of high-z sources and the precise magnification maps derived in this work for four clusters add up to the sample of high-quality gravitational telescopes to be used to study the faint and distant Universe. The strong lensing models and the full redshift catalogues from MUSE are publicly available.

Our aim is to determine the distribution of stellar population parameters (extinction, age, metallicity, and star formation rate) of quiescent galaxies within the rest-frame stellar mass$-$colour and $UVJ$ colour$-$colour diagrams corrected for extinction up to $z\sim1$. These novel diagrams reduce the contamination in samples of quiescent galaxies owing to dust-reddened galaxies, and they provide useful constraints on stellar population parameters. We set constraints on the stellar population parameters of quiescent galaxies combining the ALHAMBRA multi-filter photo-spectra with our SED-fitting code MUFFIT, making use of composite stellar population models. The extinction obtained by MUFFIT allowed us to remove dusty star-forming (DSF) galaxies from the sample of red $UVJ$ galaxies. The distributions of stellar population parameters across these rest-frame diagrams are revealed after the dust correction and are fitted by the LOESS method to reduce uncertainty effects. Quiescent galaxy samples defined via classical $UVJ$ diagrams are typically contaminated by a $\sim20$% fraction of DSF galaxies. A significant part of the galaxies in the green valley are actually obscured star-forming galaxies ($\sim30-65$%). Consequently, the transition of galaxies from the blue cloud to the red sequence, and hence the related mechanisms for quenching, seems to be much more efficient and faster than previously reported. The rest-frame stellar mass$-$colour and $UVJ$ colour$-$colour diagrams are useful for constraining the age, metallicity, extinction, and star formation rate of quiescent galaxies by only their redshift, rest-frame colours, and/or stellar mass. Dust correction plays an important role in understanding how quiescent galaxies are distributed in these diagrams and is key to performing a pure selection of quiescent galaxies via intrinsic colours.

The clustering amplitude of 7143 clusters from the Sloan Digital Sky Survey (SDSS) is found to increase with cluster mass, closely agreeing with the Gaussian random field hypothesis for structure formation. The amplitude of the observed cluster correlation exceeds the predictions from pure cold dark matter (CDM) simulation by $\simeq 6\%$ for the standard Planck-based values of the cosmological parameters. We show that this excess can be naturally accounted for by free streaming of light neutrinos, which opposes gravitational growth, so clusters formed at fixed mass are fewer and hence more biased than for a pure CDM density field. An enhancement of the cluster bias by 7\% matches the observations, corresponding to a total neutrino mass, $m_{\nu} = 0.119 \pm 0.034$ eV at 67\% confidence level, for the standard relic neutrino density. If ongoing laboratory experiments favor a normal neutrino mass hierarchy then we may infer a somewhat larger total mass than the minimum oscillation based value, $\sum m_{\nu} \simeq 0.056eV$, with 90\% confidence. Much higher precision can be achieved by applying our method to a larger sample of more distant clusters with weak lensing derived masses.

We aim at constraining the stellar population properties of quiescent galaxies. These properties reveal how these galaxies evolved and assembled since $z\sim1$ up to the present time. Combining the ALHAMBRA multi-filter photo-spectra with the SED-fitting code MUFFIT, we build a complete catalogue of quiescent galaxies via the dust-corrected stellar mass vs colour diagram. This catalogue includes stellar population properties, such as age, metallicity, extinction, stellar mass and photometric redshift, retrieved from the analysis of composited populations based on two independent sets of SSP models. We develop and apply a novel methodology to provide, for the first time, the analytic probability distribution functions (PDFs) of mass-weighted age, metallicity, and extinction of quiescent galaxies as a function of redshift and stellar mass. We adopt different star formation histories to discard potential systematics in the analysis. The number density of quiescent galaxies is found to increase since $z\sim1$, with a more substantial variation at lower mass. Quiescent galaxies feature extinction $A_V<0.6$, with median values in the range $A_V = 0.15\mathrm{-}0.3$. At increasing stellar mass, quiescent galaxies are older and more metal rich since $z\sim1$. A detailed analysis of the PDFs reveals that the evolution of quiescent galaxies is not compatible with passive evolution and a slight decrease is hinted at median metallicity $0.1\mathrm{-}0.2$~dex. The intrinsic dispersion of the age and metallicity PDFs show a dependence with stellar mass and/or redshift. These results are consistent with both sets of SSP models and the alternative SFH assumptions explored. Consequently, the quiescent population must undergo an evolutive pathway including mergers and/or remnants of star formation to reconcile the observed trends, where the `progenitor' bias should also be taken into account.

We take advantage of new 100-hour Spitzer/IRAC observations available for MACS0647-JD, a strongly lensed $z\approx11$ galaxy candidate, to provide improved constraints on its physical properties. Probing the physical properties of galaxies at $z>8$ is challenging due to the inherent faintness of such sources and the limited wavelength coverage available. Thanks to the high $\approx$2-6$\times$ lensing magnification of the multiple images of MACS0647-JD, we can use the sensitive Spitzer/IRAC data to probe the rest-frame optical fluxes of MACS0647-JD and investigate its physical properties including the age and the stellar mass. In deriving Spitzer/IRAC fluxes for MACS0647-JD, great care is taken in coping with the impact of three bright ($\approx$8-16 mag) stars in our field to ensure robust results. Assuming a constant star formation rate, the age, stellar mass, and rest-frame UV slope we estimate for MACS0647-JD based on a stack of the photometry are log$_{10}$(age/yr) = 8.6$^{+0.1}_{-2.1}$, log$_{10}$(M$_{*}$/M$_{\odot}$) = 9.1$^{+0.2}_{-1.4}$, and $\beta = -$1.3$\pm$0.6, respectively. We compare our results with expectations from the EAGLE simulation and find that MACS0647-JD has properties consistent with corresponding to the most massive and rapidly star-forming galaxies in the simulation. We also find that its radius, 105$\pm$28 pc, is a factor of $\approx$2 smaller than the mean size in a separate simulation project DRAGONS. Interestingly enough, the observed size is similar to the small sizes seen in very low-luminosity $z\approx$6-10 galaxies behind lensing clusters.

Recent near-IR power-spectra and panchromatic Extragalactic Background Light measurements provide upper limits on the near-IR surface brightness (SB>31 mag/arcsec^2) that may come from Pop III stars and accretion disks around resulting stellar-mass black holes (BHs) in the epoch of First Light (z=7-17). Physical parameters for zero metallicity Pop III stars at z>7 can be estimated from MESA stellar evolution models through helium-depletion, and for BH accretion disks from quasar microlensing results and multicolor accretion models. Second-generation stars can form at higher multiplicity, so that BH accretion disks may be fed by Roche-lobe overflow from lower-mass companions in their AGB stage. The near-IR SB constraints can be used to calculate the number of caustic transits behind lensing clusters that JWST and the 25~39 m ground-based telescopes may detect for both Pop III stars and stellar mass BH accretion disks. Because Pop III stars and stellar mass BH accretion disks have sizes of a few x 10^-11 arcsec at z>7, typical caustic magnifications can be mu=10^4~10^5, with rise times of hours and decline times of < 1 year for cluster transverse velocities of v<1000 km/s. Microlensing by intracluster medium objects can modify transit magnifications, and lengthen visibility times. Depending on BH masses, accretion-disk radii and feeding efficiencies, stellar-mass BH accretion-disk caustic transits could outnumber those from Pop III stars. To observe Pop III caustic transits directly may require monitoring 3~30 lensing clusters to AB< 29 mag over a decade or more. Such a program must be started with JWST in Cycle 1, and -- depending on the role of microlensing in the Intra Cluster Light -- should be continued for decades with the GMT and TMT, where JWST and the ground-based telescopes each will play a unique and strongly complementary role.

Microlenses with typical stellar masses (a few ${\rm M}_{\odot}$) have traditionally been disregarded as potential sources of gravitational lensing effects at LIGO/Virgo frequencies, since the time delays are often much smaller than the inverse of the frequencies probed by LIGO/Virgo, resulting in negligible interference effects at LIGO/Virgo frequencies. While this is true for isolated microlenses in this mass regime, we show how, under certain circumstances and for realistic scenarios, a population of microlenses (for instance stars and remnants from a galaxy halo or from the intracluster medium) embedded in a macromodel potential (galaxy or cluster) can conspire together to produce time delays of order one millisecond which would produce significant interference distortions in the observed strains. At sufficiently large magnification factors (of several hundred), microlensing effects should be common in gravitationally lensed gravitational waves. We explore the regime where the predicted signal falls in the frequency range probed by LIGO/Virgo. We find that stellar mass microlenses, permeating the lens plane, and near critical curves, can introduce interference distortions in strongly lensed gravitational waves. For those lensed events with negative parity, (or saddle points, never studied before in the context of gravitational waves), and that take place near caustics of macromodels, they are more likely to produce measurable interference effects at LIGO/Virgo frequencies. This is the first study that explores the effect of a realistic population of microlenses, plus a macromodel, on strongly lensed gravitational waves.

Individual highly magnified stars have been recently discovered at lookback times of more than half the age of the Universe, in lensed galaxies that straddle the critical curves of massive galaxy clusters. Having confirmed their detectability, it is now important to carry out systematic searches for them in order to establish their frequency, and in turn learn about the statistical properties of high-redshift stars and of the granularity of matter in the foreground deflector. Here we report the discovery of a highly magnified star at redshift $z = 0.94$ in a strongly lensed arc behind a Hubble Frontier Field galaxy cluster, MACS J0416.1-2403, discovered as part of a systematic archival search. The bright transient (dubbed "Warhol") was discovered in Hubble Space Telescope data taken on 2014 September 15 and 16. This single image faded over a period of two weeks, and observations taken on 2014 September 1 show that the duration of the microlensing event was at most four weeks in total. The light curve may also exhibit slow changes over a period of years consistent with the level of microlensing expected from stars responsible for the intracluster light (ICL) of the cluster. Optical and infrared observations taken near peak brightness can be fit by a stellar spectrum with moderate host-galaxy extinction. A blue supergiant matches the measured spectral energy distribution near peak, implying a temporary magnification of at least several thousand. While the spectrum of an O-type star would also fit the transient's spectral energy distribution, extremely luminous O-type stars are much less common than blue supergiants. The short timescale of the event and the estimated effective temperature indicate that the lensed source is an extremely magnified star.

We identify a binary black hole (BBH) merger that appears to be multiply lensed by an intervening galaxy. The LIGO/Virgo events GW170809 and GW170814 have indistinguishable waveforms separated by 5 days, and overlap on the sky within the 90\% credible region. Their strain amplitudes are also similar, implying a modest relative magnification ratio, as expected for a pair of lensed gravitational waves. The phase of the two events is also consistent with being the same, adding more evidence in support of both events originating from the same BBH merger. The difference in the published inferred distances of each event can then be interpreted as following from their different magnifications. The observed chirp masses of both events are also similar, as expected for a pair of lensed events, with a common detected value of $29.1^{+1.3}_{-1.0}M_{\odot}$, lying at the peak of the observed distribution of chirp masses. We infer this case is a prototypical example of a lensed event that supports our lensing prediction \cite{Broadhurst2018} according to which, cosmologically distant, magnified BBH comprise most of the LIGO/Virgo events with chirp masses enhanced above $\simeq 15M_{\odot}$ by the cosmological expansion. From our predictions we estimate an intrinsic, unlensed, chirp mass of $\simeq 10-12 M_\odot$, with a source redshift in the range $0.9<z<2.5$. We also outline a joint analysis over all baseline permutations that can stringently test our lensing interpretation of these two events. More generally, lensed events effectively multiply the number of baseline permutations and motivates the use of more interferometers for round the clock coverage of all repeat events of a given source, in order to maximise the orbital details and sky localization of lensed BBH sources.

Intermediate-mass black holes are the missing link that connects stellar-mass to supermassive black holes and are key to understanding galaxy evolution. Gravitational waves, like photons, can be lensed, leading to discernable effects such as diffraction or repeated signals. We investigate the detectability of intermediate-mass black hole deflectors in the LIGO-Virgo detector network. In particular, we simulate gravitational waves with variable source distributions lensed by an astrophysical population of intermediate-mass black holes, and use standard LIGO tools to infer the properties of these lenses. We find detections of intermediate-mass black holes at $98\%$ confidence level over a wide range of binary and lens parameters. Therefore, we conclude that intermediate-mass black holes could be detected through lensing of gravitational waves in the LIGO-Virgo detector network.

At cosmological distances, gravitational lensing can provide a direct measurement of supermassive black hole (SMBH) masses irrespective of their luminosities. Here, we directly estimate the mass of a SMBH in the brightest cluster galaxy (BCG) of MACS J1149+2223.5 at $z=0.54$ through one of the multiply-lensed images of a background spiral galaxy at $z=1.49$ projected close to the BCG. In this particular image, an intrinsically compact region in one of the spiral arms is lensed into an arc that curves towards the BCG center. This arc has a radius of curvature of only $\sim$0."6, betraying the presence of a local compact deflector. Its curvature is most simply reproduced by a point-like object with a mass of $8.4^{+4.3}_{-1.8}\times10^{9}M_\odot$, similar to SMBH masses in local elliptical galaxies having comparable luminosities. The SMBH is noticeably offset by $4.4\pm0.3$ kpc from the BCG light centre, plausibly the result of a kick imparted $\sim2.0\times10^7$ years ago during the merger of two SMBHs, placing it just beyond the stellar core. A similar curvature can be produced by replacing the offset SMBH with a compact galaxy having a mass of $\sim2\times 10^{10}M_\odot$ within a cutoff radius of $<4$ kpc, and an unusually large $M/L>50(M/L)_\odot$ to make it undetectable in the deep Hubble Frontiers Fields image, at or close to the cluster redshift; such a lensing galaxy, however, perturbs the adjacent lensed images in an undesirable way.

Acting as powerful gravitational lenses, the strong lensing galaxy clusters of the deep Hubble Frontier Fields (HFF) program permit access to lower-luminosity galaxies lying at higher redshifts than hitherto possible. We analyzed the HFF to measure the volume density of Lyman-break galaxies at $z > 4.75$ by identifying a complete and reliable sample up to $z \simeq 10$. A marked deficit of such galaxies was uncovered in the highly magnified regions of the clusters relative to their outskirts, implying that the magnification of the sky area dominates over additional faint galaxies magnified above the flux limit. This negative magnification bias is consistent with a slow rollover at the faint end of the UV luminosity function, and indicates a preference for Bose-Einstein condensate dark matter with a light boson mass of $m_\mathrm{B} \simeq 10^{-22} \, \mathrm{eV}$ over standard cold dark matter. We emphasize that measuring the magnification bias requires no correction for multiply lensed images (with typically three or more images per source), whereas directly reconstructing the luminosity function will lead to an overestimate unless such images can be exhaustively matched up, especially at the faint end that is accessible only in the strongly lensed regions. In addition, we detected a distinctive downward transition in galaxy number density at $z \gtrsim 8$, which may be linked to the relatively late reionization reported by Planck. Our results suggests that JWST will likely peer into an "abyss" with essentially no galaxies detected in deep NIR imaging at $z > 10$.

We reconstruct the two-dimensional (2D) matter distributions in 20 high-mass galaxy clusters selected from the CLASH survey by using the new approach of performing a joint weak lensing analysis of 2D shear and azimuthally averaged magnification measurements. This combination allows for a complete analysis of the field, effectively breaking the mass-sheet degeneracy. In a Bayesian framework, we simultaneously constrain the mass profile and morphology of each individual cluster assuming an elliptical Navarro-Frenk-White halo characterized by the mass, concentration, projected axis ratio, and position angle of the projected major axis.. We find that spherical mass estimates of the clusters from azimuthally averaged weak-lensing measurements in previous work are in excellent agreement with our results from a full 2D analysis. Combining all 20 clusters in our sample, we detect the elliptical shape of weak-lensing halos at the $5\sigma$ significance level within a scale of 2Mpc$/h$. The median projected axis ratio is $0.67\pm 0.07$ at a virial mass of $M_\mathrm{vir}=(15.2\pm 2.8) \times 10^{14} M_\odot$, which is in agreement with theoretical predictions of the standard collisionless cold dark matter model. We also study misalignment statistics of the brightest cluster galaxy, X-ray, thermal Sunyaev-Zel'dovich effect, and strong-lensing morphologies with respect to the weak-lensing signal. Among the three baryonic tracers studied here, we find that the X-ray morphology is best aligned with the weak-lensing mass distribution, with a median misalignment angle of $21\pm 7$ degrees. We also conduct a stacked quadrupole shear analysis assuming that the X-ray major axis is aligned with that of the projected mass distribution. This yields a consistent axis ratio of $0.67\pm 0.10$, suggesting again a tight alignment between the intracluster gas and dark matter.

Gravitational waves can be focussed by the gravity of an intervening galaxy, just like light, thereby magnifying binary merging events in the far Universe. High magnification by galaxies is found to be responsible for the brightest sources detected in sky surveys, but the low angular resolution of LIGO/Virgo is insufficient to check this lensing possibility directly. Here we find that the first six binary black hole (BBH) merging events reported by LIGO/Virgo show clear evidence for lensing in the plane of observed mass and source distance. The four lowest frequency events follow an apparent locus in this plane, which we can reproduce by galaxy lensing, where the higher the magnification, the generally more distant the source so the wave train is stretched more by the Universal expansion, by factors of 2-4. This revises the reported BBH distances upwards by an order of magnitude, equal to the square root of the magnification. Furthermore, the reported black hole masses must be decreased by 2-4 to counter the larger stretch factor, since the orbital frequency is used to derive the black hole masses. This lowers the masses to 5-15 solar masses, well below the puzzlingly high values of 20-35 solar masses otherwise estimated, with the attraction of finding agreement in mass with black holes orbiting stars in our own Galaxy, thereby implying a stellar origin for the low frequency events in the far Universe. We also show that the other two BBH events of higher frequency detected by LIGO/VIRGO, lie well below the lensing locus, consistent with being nearby and unlensed. If this apparent division between local and distant lensed events is reinforced by new detections then the spins and masses of stellar black holes can be compared over a timespan of 10 billion years by LIGO/Virgo.

Galaxy-cluster gravitational lenses can magnify background galaxies by a total factor of up to ~50. Here we report an image of an individual star at redshift z=1.49 (dubbed "MACS J1149 Lensed Star 1 (LS1)") magnified by >2000. A separate image, detected briefly 0.26 arcseconds from LS1, is likely a counterimage of the first star demagnified for multiple years by a >~3 solar-mass object in the cluster. For reasonable assumptions about the lensing system, microlensing fluctuations in the stars' light curves can yield evidence about the mass function of intracluster stars and compact objects, including binary fractions and specific stellar evolution and supernova models. Dark-matter subhalos or massive compact objects may help to account for the two images' long-term brightness ratio.

A galaxy cluster acts as a cosmic telescope over background galaxies but also as a cosmic microscope of the lens imperfections. The diverging magnification of lensing caustics enhances the microlensing effect of substructure present within the lensing mass. Fine-scale structure can be accessed as a moving background source brightens and disappears when crossing these caustics. The recent recognition of a distant lensed star near the Einstein radius of the galaxy cluster MACSJ1149.5+2223 (Kelly et al. 2017) allows the rare opportunity to reach subsolar mass microlensing through a super-critical column of cluster matter. Here we compare these observations with high-resolution ray-tracing simulations that include stellar microlensing set by the observed intracluster starlight and also primordial black holes that may be responsible for the recently observed LIGO events. We explore different scenarios with microlenses from the intracluster medium and black holes, including primordial ones, and examine strategies to exploit these unique alignments. We find that the best constraints on the fraction of compact dark matter in the small-mass regime can be obtained in regions of the cluster where the intracluster medium plays a negligible role. This new lensing phenomenon should be widespread and can be detected within modest-redshift lensed galaxies so that the luminosity distance is not prohibitive for detecting individual magnified stars. Continuous {\it Hubble Space Telescope} monitoring of several such optimal arcs will be rewarded by an unprecedented mass spectrum of compact objects that can contribute to uncovering the nature of dark matter.

We have selected a sample of eleven massive clusters of galaxies observed by the Hubble Space Telescope in order to study the impact of the dynamical state on the IntraCluster Light (ICL) fraction, the ratio of total integrated ICL to the total galaxy member light. With the exception of the Bullet cluster, the sample is drawn from the Cluster Lensing and Supernova Survey and the Frontier Fields program, containing five relaxed and six merging clusters. The ICL fraction is calculated in three optical filters using the CHEFs IntraCluster Light Estimator, a robust and accurate algorithm free of a priori assumptions. We find that the ICL fraction in the three bands is, on average, higher for the merging clusters, ranging between $\sim7-23\%$, compared with the $\sim 2-11\%$ found for the relaxed systems. We observe a nearly constant value (within the error bars) in the ICL fraction of the regular clusters at the three wavelengths considered, which would indicate that the colors of the ICL and the cluster galaxies are, on average, coincident and, thus, their stellar populations. However, we find a higher ICL fraction in the F606W filter for the merging clusters, consistent with an excess of lower-metallicity/younger stars in the ICL, which could have migrated violently from the outskirts of the infalling galaxies during the merger event.

Context. Knowing the exact shape of the UV luminosity function of high-redshift galaxies is important in order to understand the star formation history of the early universe. However, the uncertainties, especially at the faint and bright ends of the LFs, are still significant. Aims. In this paper, we study the UV luminosity function of redshift z = 2.5 - 4.5 galaxies in 2.38 deg^2 of ALHAMBRA data with I <= 24. Thanks to the large area covered by ALHAMBRA, we particularly constrain the bright end of the luminosity function. We also calculate the cosmic variance and the corresponding bias values for our sample and derive their host dark matter halo masses. Methods. We use a novel methodology based on redshift and magnitude probability distribution functions (PDFs). This methodology robustly takes into account the uncertainties due to redshift and magnitude errors, shot noise and cosmic variance, and models the luminosity function in two dimensions (z; M_UV ). Results. We find an excess of bright ~ M*_UV galaxies as compared to the studies based on broad-band photometric data. However, our results agree well with the luminosity function of the magnitude-selected spectroscopic VVDS data. We measure high bias values, b ~ 8 - 10, that are compatible with the previous measurements considering the redshifts and magnitudes of our galaxies and further reinforce the real high-redshift nature of our bright galaxies. Conclusions. We call into question the shape of the luminosity function at its bright end; is it a double power-law as suggested by the recent broad-band photometric studies or rather a brighter Schechter function, as suggested by our multi-filter analysis and the spectroscopic VVDS data.

The recent discovery of fast transient events near critical curves of massive galaxy clusters, which are interpreted as highly magnified individual stars in giant arcs due to caustic crossing, opens up the possibility of using such microlensing events to constrain a range of dark matter models such as primordial black holes and scalar field dark matter. Based on a simple analytic model, we study lensing properties of a point mass lens embedded in a high magnification region, and derive the dependence of the peak brightness, microlensing time scales, and event rates on the mass of the point mass lens as well as the radius of a source star that is magnified. We find that the lens mass and source radius of the first event MACS J1149 Lensed Star 1 (LS1) are constrained, with the lens mass range of $0.1~M_\odot \lesssim M \lesssim 4\times 10^3M_\odot$ and the source radius range of $40~R_\odot \lesssim R \lesssim 260~R_\odot$. In the most plausible case with $M\approx 0.3~M_\odot$ and $R\approx 180~R_\odot$, the source star should have been magnified by a factor of $\approx 4300$ at the peak. The derived lens properties are fully consistent with the interpretation that MACS J1149 LS1 is a microlensing event produced by a star that contributes to the intra-cluster light. We argue that compact dark matter models with high fractional mass densities for the mass range $10^{-5}M_\odot \lesssim M\lesssim 10^2M_\odot$ are inconsistent with the observation of MACS J1149 LS1 because such models predict too low magnifications. Our work demonstrates a potential use of caustic crossing events in giant arcs to constrain compact dark matter.

We use hydrodynamical/N-body simulations to interpret the newly discovered Bullet-cluster-like merging cluster, ZwCl 0008.8+5215 (ZwCl 0008 hereafter), where a dramatic collision is apparent from multi-wavelength observations. We have been able to find a self-consistent solution for the radio, X-ray, and lensing phenomena by projecting an off-axis, binary cluster encounter viewed just after first core passage. A pair radio relics traces well the leading and trailing shock fronts that our simulation predict, providing constraints on the collision parameters. We can also account for the observed distinctive comet-like X-ray morphology and the positions of the X-ray peaks relative to the two lensing mass centroids and the two shock front locations. Relative to the Bullet cluster, the total mass is about 70% lower, ($1.2\pm0.1) \times 10^{15}$ Msun, with a correspondingly lower infall velocity, $1800\pm300$ km/s, and an impact parameter of $400\pm100$ kpc. As a result, the gas component of the infalling cluster is not trailing significantly behind the associated dark matter as in the case of the Bullet cluster. The degree of agreement we find between all the observables provides strong evidence that dark matter is effectively collisionless on large scales calling into question other claims and theories that advocate modified gravity.

Light Axionic Dark Matter, motivated by string theory, is increasingly favored for the "no-WIMP era". Galaxy formation is suppressed below a Jeans scale, of $\simeq 10^8 M_\odot$ by setting the axion mass to, $m_B \sim 10^{-22}$eV, and the large dark cores of dwarf galaxies are explained as solitons on the de-Broglie scale. This is persuasive, but detection of the inherent scalar field oscillation at the Compton frequency, $\omega_B= (2.5{\rm \, months})^{-1}(m_B/10^{-22}eV)$, would be definitive. By evolving the coupled Schr\"odinger-Poisson equation for a Bose-Einstein condensate, we predict the dark matter is fully modulated by de-Broglie interference, with a dense soliton core of size $\simeq 150pc$, at the Galactic center. The oscillating field pressure induces General Relativistic time dilation in proportion to the local dark matter density and pulsars within this dense core have detectably large timing residuals, of $\simeq 400nsec/(m_B/10^{-22}eV)$. This is encouraging as many new pulsars should be discovered near the Galactic center with planned radio surveys. More generally, over the whole Galaxy, differences in dark matter density between pairs of pulsars imprints a pairwise Galactocentric signature that can be distinguished from an isotropic gravitational wave background.

We present a technique that permits the analysis of stellar population gradients in a relatively low cost way compared to IFU surveys analyzing a vastly larger samples as well as out to larger radii. We developed a technique to analyze unresolved stellar populations of spatially resolved galaxies based on photometric multi-filter surveys. We derived spatially resolved stellar population properties and radial gradients by applying a Centroidal Voronoi Tesselation and performing a multi-color photometry SED fitting. This technique has been applied to a sample of 29 massive (M$_{\star}$ > 10$^{10.5}$ M$_{\odot}$), early-type galaxies at $z$ < 0.3 from the ALHAMBRA survey. We produced detailed 2D maps of stellar population properties (age, metallicity and extinction). Radial structures have been studied and luminosity-weighted and mass-weighted gradients have been derived out to 2 - 3.5 R$_\mathrm{eff}$. We find the gradients of early-type galaxies to be on average flat in age ($\nabla$log Age$_\mathrm{L}$ = 0.02 $\pm$ 0.06 dex/R$_\mathrm{eff})$ and negative in metallicity ($\nabla$[Fe/H]$_\mathrm{L}$ = - 0.09 $\pm$ 0.06 dex/R$_\mathrm{eff}$). Overall, the extinction gradients are flat ($\nabla$A$_\mathrm{v}$ = - 0.03 $\pm$ 0.09 mag/R$_\mathrm{eff}$ ) with a wide spread. These results are in agreement with previous studies that used standard long-slit spectroscopy as well as with the most recent integral field unit (IFU) studies. According to recent simulations, these results are consistent with a scenario where early-type galaxies were formed through major mergers and where their final gradients are driven by the older ages and higher metallicity of the accreted systems. We demonstrate the scientific potential of multi-filter photometry to explore the spatially resolved stellar populations of local galaxies and confirm previous spectroscopic trends from a complementary technique.

A massive galaxy cluster can serve as a magnifying glass for distant stellar populations, with strong gravitational lensing exposing details in the lensed background galaxies that would otherwise be undetectable. The MACS J0416.1-2403 cluster (hereafter MACS0416) is one of the most efficient lenses in the sky, and in 2014 it was observed with high-cadence imaging from the Hubble Space Telescope (HST). Here we describe two unusual transient events that appeared behind MACS0416 in a strongly lensed galaxy at redshift $z = 1.0054 \pm 0.0002$. These transients---designated HFF14Spo-NW and HFF14Spo-SE and collectively nicknamed "Spock"---were faster and fainter than any supernova (SN), but significantly more luminous than a classical nova. They reached peak luminosities of $\sim10^{41}$ erg s$^{-1}$ ($M_{\rm AB} < -14$ mag) in 5 rest-frame days, then faded below detectability in roughly the same time span. Models of the cluster lens suggest that these events may be spatially coincident at the source plane, but are most likely not temporally coincident. We find that HFF14Spo can be explained as a luminous blue variable (LBV), a recurrent nova (RN), or a pair of stellar microlensing events. To distinguish between these hypotheses will require a clarification of the positions of nearby critical curves, along with high-cadence monitoring of the field that could detect new transient episodes in the host galaxy.

We present accurate photometric redshifts for galaxies observed by the Cluster Lensing and Supernova survey with Hubble (CLASH). CLASH observed 25 massive galaxy cluster cores with the Hubble Space Telescope in 16 filters spanning 0.2 - 1.7 $\mu$m. Photometry in such crowded fields is challenging. Compared to our previously released catalogs, we make several improvements to the photometry, including smaller apertures, ICL subtraction, PSF matching, and empirically measured uncertainties. We further improve the Bayesian Photometric Redshift (BPZ) estimates by adding a redder elliptical template and by inflating the photometric uncertainties of the brightest galaxies. The resulting photometric redshift accuracies are dz/(1+z) $\sim$ 0.8\%, 1.0\%, and 2.0\% for galaxies with I-band F814W AB magnitudes $<$ 18, 20, and 23, respectively. These results are consistent with our expectations. They improve on our previously reported accuracies by a factor of 4 at the bright end and a factor of 2 at the faint end. Our new catalog includes 1257 spectroscopic redshifts, including 382 confirmed cluster members. We also provide stellar mass estimates. Finally, we include lensing magnification estimates of background galaxies based on our public lens models. Our new catalog of all 25 CLASH clusters is available via MAST. The analysis techniques developed here will be useful in other surveys of crowded fields, including the Frontier Fields and surveys carried out with J-PAS and JWST.

We construct a large, redshift complete sample of distant galaxy clusters by correlating Sloan Digital Sky Survey (SDSS) Data Release 12 (DR12) redshifts with clusters identified with the red-sequence Matched-filter Probabilistic Percolation (redMaPPer) algorithm. Our spectroscopic completeness is 97% for ~ 7,000 clusters within the redMaPPer selection limit, $z \leqslant$ 0.325, so that our cluster correlation functions are much more precise than earlier work and not suppressed by photometric redshifts. We derive an accurate power-law mass-richness relation from the observed abundance with respect to the mass function from Millennium XXL (MXXL) simulation, adjusted to the Planck weighted cosmology. The number density of clusters is found to decline by 20% over the range 0.1 $< z <$ 0.3, in good agreement with the evolution predicted by MXXL. Our projected three-dimensional correlation function scales with richness, $\lambda$, rising from $r_0=$ 14 h$^{-1}$ Mpc at $\lambda\simeq$ 25, to $r_0=$ 22 h$^{-1}$ Mpc at $\lambda\simeq$ 60, with a gradient that matches MXXL when applying our mass-richness relation, whereas the observed amplitude of the correlation function at $\left<z\right>=$ 0.24 exceeds the MXXL prediction by 20% at the $\simeq$ 2.5$\sigma$ level. This tension cannot be blamed on spurious, randomly located clusters as this would reduce the correlation amplitude. Full consistency between the correlation function and the abundances is achievable for the pre-Planck values of $\sigma_8=$ 0.9, $\Omega_m=$ 0.25, and h = 0.73, matching the improved distance ladder estimate of the Hubble constant.

Using Planck satellite data, we construct SZ gas pressure profiles for a large, volume-complete sample of optically selected clusters. We have defined a sample of over 8,000 redMaPPer clusters from the Sloan Digital Sky Survey (SDSS), within the volume-complete redshift region 0.100 < z < 0.325, for which we construct Sunyaev-Zel'dovich (SZ) effect maps by stacking Planck data over the full range of richness. Dividing the sample into richness bins we simultaneously solve for the mean cluster mass in each bin together with the corresponding radial pressure profile parameters, employing an MCMC analysis. These profiles are well detected over a much wider range of cluster mass and radius than previous work, showing a clear trend towards larger break radius with increasing cluster mass. Our SZ-based masses fall ~24% below the mass-richness relations from weak lensing, in a similar fashion as the "hydrostatic bias" related with X-ray derived masses. We correct for this bias to derive an optimal mass-richness relation finding a slope 1.22 +/- 0.04 and a pivot mass log(M_500/M_0)= 14.432 +/- 0.041, evaluated at a richness lambda=60. Finally, we derive a tight Y_500-M_500 relation over a wide range of cluster mass, with a power law slope equal to 1.72 +/- 0.07, that agrees well with the independent slope obtained by the Planck team with an SZ-selected cluster sample, but extends to lower masses with higher precision.

We present a light-traces-mass (LTM) strong-lensing model of the massive lensing cluster MACS J2135.2-0102 ($z$=0.33; hereafter MACS2135), known in part for hosting the Cosmic Eye galaxy lens. MACS2135 is also known to multiply-lens a $z=$2.3 sub-mm galaxy near the Brightest Cluster Galaxy (BCG), as well as a prominent, triply-imaged system at a large radius of $\sim$37" south of the BCG. We use the latest available Hubble imaging to construct an accurate lensing model for this cluster, identifying six new multiply-imaged systems with the guidance of our LTM method, so that we have roughly quadrupled the number of lensing constraints. We determine that MACS2135 is amongst the top lensing clusters known, comparable in size to the Hubble Frontier Fields. For a source at $z_{s}=2.32$, we find an effective Einstein radius of $\theta_{e}=27\pm3$", enclosing $1.12 \pm0.16 \times10^{14}$ $M_{\odot}$. We make our lens model, including mass and magnification maps, publicly available, in anticipation of searches for high-$z$ galaxies with the James Webb Space Telescope for which this cluster is a compelling target.

The colliding cluster, CIZA J2242.8+5301, displays a spectacular, almost 2 Mpc long shock front with a radio based Mach number M ~ 5, that is puzzlingly large compared with the X-ray estimate of M ~ 2.5. The extent to which the X-ray temperature jump is diluted by cooler unshocked gas projected through the cluster currently lacks quantification. Thus, here we apply our self-consistent N-body/hydro-dynamical code (based on FLASH) to model this binary cluster encounter. We can account for the location of the shock front and also the elongated X-ray emission by tidal stretching of the gas and dark matter between the two cluster centers. The required total mass is $8.9 \times 10^{14}$ Msun with a 1.3:1 mass ratio favoring the southern cluster component. The relative velocity we derive is $\simeq 2500$ km/s initially between the two main cluster components, with an impact parameter of 120 kpc. This solution implies that the shock temperature jump derived from the low angular resolution X-ray satellite SUZAKU is underestimated by a factor of two, due to cool gas in projection, bringing the observed X-ray and radio estimates into agreement. We propose that the complex southern relics in CIZA J2242.8+5301, have been broken up as the southerly moving "back" shocked gas impacts the gas still falling in along the collision axis. Finally, we use our model to generate Compton-y maps to estimate the reduction in radio flux caused by the thermal Sunyaev-Zel'dovich (SZ) effect. At 30 GHz, this amounts to $\Delta S_n = -0.072$ mJy/arcmin$^2$ and $\Delta S_s = -0.075$ mJy/arcmin$^2$ at the locations of the northern and southern shock fronts respectively. Our model estimate agrees with previous empirical estimates that have inferred the measured radio spectra can be significantly affected by the SZ effect, with implications for charged particle acceleration models of the radio relics.

We show how LIGO is expected to detect coalescing binary black holes at $z>1$, that are lensed by the intervening galaxy population. Gravitational magnification, $\mu$, strengthens gravitational wave signals by $\sqrt{\mu}$, without altering their frequencies, which if unrecognised leads to an underestimate of the event redshift and hence an overestimate of the binary mass. High magnifications can be reached for coalescing binaries because the region of intense gravitational wave emission during coalescence is so small ($\sim$100km), permitting very close projections between lensing caustics and gravitational-wave events. Our simulations incorporate accurate waveforms convolved with the LIGO power spectral density. Importantly, we include the detection dependence on sky position and orbital orientation, which for the LIGO configuration translates into a wide spread in observed redshifts and chirp masses. Currently we estimate a detectable rate of lensed events \rateEarly{}, that rises to \rateDesign{}, at LIGO's design sensitivity limit, depending on the high redshift rate of black hole coalescence.

We search for high-redshift dropout galaxies behind the Hubble Frontier Fields (HFF) galaxy cluster MACS J1149.5+2223, a powerful cosmic lens that has revealed a number of unique objects in its field. Using the deep images from the Hubble and Spitzer space telescopes, we find 11 galaxies at z>7 in the MACS J1149.5+2223 cluster field, and 11 in its parallel field. The high-redshift nature of the bright z~9.6 galaxy MACS1149-JD, previously reported by Zheng et al., is further supported by non-detection in the extremely deep optical images from the HFF campaign. With the new photometry, the best photometric redshift solution for MACS1149-JD reduces slightly to z=9.44 +/- 0.12. The young galaxy has an estimated stellar mass of (7 +/- 2)X10E8 Msun, and was formed at z=13.2 +1.9-1.6 when the universe was ~300 Myr old. Data available for the first four HFF clusters have already enabled us to find faint galaxies to an intrinsic magnitude of M(UV) ~ -15.5, approximately a factor of ten deeper than the parallel fields.

Our goal is to study the evolution of the $B-$band luminosity function (LF) since $z=1$ using ALHAMBRA data. We used the photometric redshift and the $I-$band selection magnitude probability distribution functions (PDFs) of those ALHAMBRA galaxies with $I\leq24$ mag to compute the posterior LF. We statistically studied quiescent and star-forming galaxies using the template information encoded in the PDFs. The LF covariance matrix in redshift-magnitude-galaxy type space was computed, including the cosmic variance. That was estimated from the intrinsic dispersion of the LF measurements in the 48 ALHAMBRA sub-fields. The uncertainty due to the photometric redshift prior is also included in our analysis. We modelled the LF with a redshift-dependent Schechter function affected by the same selection effects than the data. The measured ALHAMBRA LF at $0.2\leq z<1$ and the evolving Schechter parameters both for quiescent and star-forming galaxies agree with previous results in the literature. The estimated redshift evolution of $M_B^* \propto Qz$ is $Q_{\rm SF}=-1.03\pm0.08$ and $Q_{\rm Q}=-0.80\pm0.08$, and of $\log \phi^* \propto Pz$ is $P_{\rm SF}=-0.01\pm0.03$ and $P_{\rm Q}=-0.41\pm0.05$. The measured faint-end slopes are $\alpha_{\rm SF}=-1.29\pm0.02$ and $\alpha_{\rm Q}=-0.53\pm0.04$. We find a significant population of faint quiescent galaxies, modelled by a second Schechter function with slope $\beta=-1.31\pm0.11$. We find a factor $2.55\pm0.14$ decrease in the luminosity density $j_B$ of star-forming galaxies, and a factor $1.25\pm0.16$ increase in the $j_B$ of quiescent ones since $z=1$, confirming the continuous build-up of the quiescent population with cosmic time. The contribution of the faint quiescent population to $j_B$ increases from 3% at $z=1$ to 6% at $z=0$. The developed methodology will be applied to future multi-filter surveys such as J-PAS.

A multiply-lensed galaxy, MACS0647-JD, with a probable photometric redshift of $z\simeq 10.7^{+0.6}_{-0.4}$ is claimed to constitute one of the very earliest known galaxies, formed well before reionization was completed. However, spectral evidence that MACS0647-JD lies at high redshift has proven infeasible and so here we seek an independent lensing based "geometric redshift" derived from the angles between the three lensed images of MACS0647-JD, using our free-form mass model (WSLAP+) for the lensing cluster MACSJ0647.7+7015 (at $z=0.591$). Our lens model uses the 9 sets of multiple images, including those of MACS0647-JD, identified by the CLASH survey towards this cluster. We convincingly exclude the low redshift regime of $z<3$, for which convoluted critical curves are generated by our method, as the solution bends to accommodate the wide angles of MACS0647-JD for this low redshift. Instead, a best fit to all sets of lensed galaxy positions and redshifts provides a geometric redshift of $z\simeq 10.8^{+0.3}_{-0.4}$ for MACS0647-JD, strongly supporting the higher photometric redshift solution. Importantly, we find a tight linear relation between the relative brightnesses of all 9 sets of multiply lensed images and their relative magnifications as predicted by our model. This agreement provides a benchmark for the quality of the lens model, and establishes the robustness of our free-form lensing method for measuring model-independent geometric source distances and for deriving objective central cluster mass distributions. After correcting for its magnification the luminosity of MACS0647-JD remains relatively high at $M_{UV}=-19.4$, which is within a factor of a few in flux of some surprisingly luminous $z\simeq 10$--$11$ candidates discovered recently in Hubble black field surveys.

The original ALHAMBRA catalogue contained over 400,000 galaxies selected using a synthetic F814W image, to the magnitude limit AB(F814W)$\approx$24.5. Given the photometric redshift depth of the ALHAMBRA multiband data (<z>=0.86) and the approximately $I$-band selection, there is a noticeable bias against red objects at moderate redshift. We avoid this bias by creating a new catalogue selected in the $K_s$ band. This newly obtained catalogue is certainly shallower in terms of apparent magnitude, but deeper in terms of redshift, with a significant population of red objects at $z>1$. We select objects using the $K_s$ band images, which reach an approximate AB magnitude limit $K_s \approx 22$. We generate masks and derive completeness functions to characterize the sample. We have tested the quality of the photometry and photometric redshifts using both internal and external checks. Our final catalogue includes $\approx 95,000$ sources down to $K_s \approx 22$, with a significant tail towards high redshift. We have checked that there is a large sample of objects with spectral energy distributions that correspond to that of massive, passively evolving galaxies at $z > 1$, reaching as far as $z \approx 2.5$. We have tested the possibility of combining our data with deep infrared observations at longer wavelengths, particularly Spitzer IRAC data.

We derive a free-form mass distribution for the unrelaxed cluster A370 (z=0.375), using the latest Hubble Frontier Fields images and GLASS spectroscopy. Starting from a reliable set of 10 multiply lensed systems we produce a free-form lens model that identifies ~ 80 multiple-images. Good consistency is found between models using independent subsamples of these lensed systems, with detailed agreement for the well resolved arcs. The mass distribution has two very similar concentrations centred on the two prominent Brightest Cluster Galaxies (or BCGs), with mass profiles that are accurately constrained by a uniquely useful system of long radially lensed images centred on both BCGs. We show that the lensing mass profiles of these BCGs are mainly accounted for by their stellar mass profiles, with a modest contribution from dark matter within r<100 kpc of each BCG. This conclusion may favour a cooled cluster gas origin for BCGs, rather than via mergers of normal galaxies for which dark matter should dominate over stars. Growth via merging between BCGs is, however, consistent with this finding, so that stars still dominate over dark matter .

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.

We present VIMOS-VLT spectroscopy of the Frontier Fields cluster MACS~J0416.1-2403 (z=0.397). Taken as part of the CLASH-VLT survey, the large spectroscopic campaign provided more than 4000 reliable redshifts, including ~800 cluster member galaxies. The unprecedented sample of cluster members at this redshift allows us to perform a highly detailed dynamical and structural analysis of the cluster out to ~2.2$r_{200}$ (~4Mpc). Our analysis of substructures reveals a complex system composed of a main massive cluster ($M_{200}$~0.9$\times 10^{15} M_{\odot}$) presenting two major features: i) a bimodal velocity distribution, showing two central peaks separated by $\Delta V_{rf}$~1100 km s$^{-1}$ with comparable galaxy content and velocity dispersion, ii) a projected elongation of the main substructures along the NE-SW direction, with a prominent subclump ~600 kpc SW of the center and an isolated BCG approximately halfway between the center and the SW clump. We also detect a low mass structure at z~0.390, ~10' S of the cluster center, projected at ~3Mpc, with a relative line-of-sight velocity of $\Delta V_{rf}$~-1700 km s$^{-1}$. The cluster mass profile that we obtain through our dynamical analysis deviates significantly from the "universal" NFW, being best fit by a Softened Isothermal Sphere model instead. The mass profile measured from the galaxy dynamics is found to be in relatively good agreement with those obtained from strong and weak lensing, as well as with that from the X-rays, despite the clearly unrelaxed nature of the cluster. Our results reveal overall a complex dynamical state of this massive cluster and support the hypothesis that the two main subclusters are being observed in a pre-collisional phase, in line with recent findings from radio and deep X-ray data. With this article we also release the entire redshift catalog of 4386 sources in the field of this cluster.

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.

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, $\psi$DM, appropriate for light bosons such as axions, does not suffer this problem, given a boson mass of $m_{\psi} \ge 1.2 \times 10^{-22}{\,\rm eV}$ ($2\sigma$). 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 $\psi$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, $\psi$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.

We study the clustering of galaxies as a function of spectral type and redshift in the range $0.35 < z < 1.1$ using data from the Advanced Large Homogeneous Area Medium Band Redshift Astronomical (ALHAMBRA) survey. The data cover 2.381 deg$^2$ in 7 fields, after applying a detailed angular selection mask, with accurate photometric redshifts [$\sigma_z < 0.014(1+z)$] down to $I_{AB} < 24$. From this catalog we draw five fixed number density, redshift-limited bins. We estimate the clustering evolution for two different spectral populations selected using the ALHAMBRA-based photometric templates: quiescent and star-forming galaxies. For each sample, we measure the real-space clustering using the projected correlation function. Our calculations are performed over the range $[0.03,10.0] h^{-1}$ Mpc, allowing us to find a steeper trend for $r_p \lesssim 0.2 h^{-1}$ Mpc, which is especially clear for star-forming galaxies. Our analysis also shows a clear early differentiation in the clustering properties of both populations: star-forming galaxies show weaker clustering with evolution in the correlation length over the analysed redshift range, while quiescent galaxies show stronger clustering already at high redshifts, and no appreciable evolution. We also perform the bias calculation where similar segregation is found, but now it is among the quiescent galaxies where a growing evolution with redshift is clearer. These findings clearly corroborate the well known colour-density relation, confirming that quiescent galaxies are mainly located in dark matter halos that are more massive than those typically populated by star-forming galaxies.

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 $\sigma_{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 $\sigma_{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.

We perform a comprehensive study of the total mass distribution of the galaxy cluster RXCJ2248 ($z=0.348$) with a set of high-precision strong lensing models, which take advantage of extensive spectroscopic information on many multiply lensed systems. In the effort to understand and quantify inherent systematics in parametric strong lensing modelling, we explore a collection of 22 models where we use different samples of multiple image families, parametrizations of the mass distribution and cosmological parameters. As input information for the strong lensing models, we use the CLASH HST imaging data and spectroscopic follow-up observations, carried out with the VIMOS and MUSE spectrographs, to identify bona-fide multiple images. A total of 16 background sources, over the redshift range $1.0-6.1$, are multiply lensed into 47 images, 24 of which are spectroscopically confirmed and belong to 10 individual sources. The cluster total mass distribution and underlying cosmology in the models are optimized by matching the observed positions of the multiple images on the lens plane. We show that with a careful selection of a sample of spectroscopically confirmed multiple images, the best-fit model reproduces their observed positions with a rms of $0.3$ in a fixed flat $\Lambda$CDM cosmology, whereas the lack of spectroscopic information lead to biases in the values of the model parameters. Allowing cosmological parameters to vary together with the cluster parameters, we find (at $68\%$ confidence level) $\Omega_m=0.25^{+0.13}_{-0.16}$ and $w=-1.07^{+0.16}_{-0.42}$ for a flat $\Lambda$CDM model, and $\Omega_m=0.31^{+0.12}_{-0.13}$ and $\Omega_\Lambda=0.38^{+0.38}_{-0.27}$ for a universe with $w=-1$ and free curvature. Using toy models mimicking the overall configuration of RXCJ2248, we estimate the impact of the line of sight mass structure on the positional rms to be $0.3\pm 0.1$.(ABRIDGED)

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.

Supernova "Refsdal," multiply imaged by cluster MACSJ1149.5+2223, represents a rare opportunity to make a true blind test of model predictions in extragalactic astronomy, on a time scale that is short compared to a human lifetime. In order to take advantage of this event, we produced seven gravitational lens models with five independent methods, based on Hubble Space Telescope (HST) Hubble Frontier Field images, along with extensive spectroscopic follow-up observations by HST, the Very Large and the Keck Telescopes. We compare the model predictions and show that they agree reasonably well with the measured time delays and magnification ratios between the known images, even though these quantities were not used as input. This agreement is encouraging, considering that the models only provide statistical uncertainties, and do not include additional sources of uncertainties such as structure along the line of sight, cosmology, and the mass sheet degeneracy. We then present the model predictions for the other appearances of SN "Refsdal." A future image will reach its peak in the first half of 2016, while another image appeared between 1994 and 2004. The past image would have been too faint to be detected in existing archival images. The future image should be approximately one third as bright as the brightest known image (i.e., H_AB~25.7 mag at peak and H_AB~26.7 mag six months before peak), and thus detectable in single-orbit HST images. We will find out soon whether our predictions are correct.

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.

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.

The relative cosmic variance ($\sigma_v$) is a fundamental source of uncertainty in pencil-beam surveys and, as a particular case of count-in-cell statistics, can be used to estimate the bias between galaxies and their underlying dark-matter distribution. Our goal is to test the significance of the clustering information encoded in the $\sigma_v$ measured in the ALHAMBRA survey. We measure the cosmic variance of several galaxy populations selected with $B-$band luminosity at $0.35 \leq z < 1.05$ as the intrinsic dispersion in the number density distribution derived from the 48 ALHAMBRA subfields. We compare the observational $\sigma_v$ with the cosmic variance of the dark matter expected from the theory, $\sigma_{v,{\rm dm}}$. This provides an estimation of the galaxy bias $b$. The galaxy bias from the cosmic variance is in excellent agreement with the bias estimated by two-point correlation function analysis in ALHAMBRA. This holds for different redshift bins, for red and blue subsamples, and for several $B-$band luminosity selections. We find that $b$ increases with the $B-$band luminosity and the redshift, as expected from previous work. Moreover, red galaxies have a larger bias than blue galaxies, with a relative bias of $b_{\rm rel} = 1.4 \pm 0.2$. Our results demonstrate that the cosmic variance measured in ALHAMBRA is due to the clustering of galaxies and can be used to characterise the $\sigma_v$ affecting pencil-beam surveys. In addition, it can also be used to estimate the galaxy bias $b$ from a method independent of correlation functions.

We consistently analyse for the first time the impact of survey depth and spatial resolution on the most used morphological parameters for classifying galaxies through non-parametric methods: Abraham and Conselice-Bershady concentration indices, Gini, M20 moment of light, asymmetry, and smoothness. Three different non-local datasets are used, ALHAMBRA and SXDS (examples of deep ground-based surveys), and COSMOS (deep space-based survey). We used a sample of 3000 local, visually classified galaxies, measuring their morphological parameters at their real redshifts (z ~ 0). Then we simulated them to match the redshift and magnitude distributions of galaxies in the non-local surveys. The comparisons of the two sets allow to put constraints on the use of each parameter for morphological classification and evaluate the effectiveness of the commonly used morphological diagnostic diagrams. All analysed parameters suffer from biases related to spatial resolution and depth, the impact of the former being much stronger. When including asymmetry and smoothness in classification diagrams, the noise effects must be taken into account carefully, especially for ground-based surveys. M20 is significantly affected, changing both the shape and range of its distribution at all brightness levels.We suggest that diagnostic diagrams based on 2 - 3 parameters should be avoided when classifying galaxies in ground-based surveys, independently of their brightness; for COSMOS they should be avoided for galaxies fainter than F814 = 23.0. These results can be applied directly to surveys similar to ALHAMBRA, SXDS and COSMOS, and also can serve as an upper/lower limit for shallower/deeper ones.

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 $\sigma_{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).

We present MUFFIT, a new generic code optimized to retrieve the main stellar population parameters of galaxies in photometric multi-filter surveys, and we check its reliability and feasibility with real galaxy data from the ALHAMBRA survey. Making use of an error-weighted $\chi^2$-test, we compare the multi-filter fluxes of galaxies with the synthetic photometry of mixtures of two single stellar populations at different redshifts and extinctions, to provide through a Monte Carlo method the most likely range of stellar population parameters (mainly ages and metallicities), extinctions, redshifts, and stellar masses. To improve the diagnostic reliability, MUFFIT identifies and removes from the analysis those bands that are significantly affected by emission lines. We highlight that the retrieved age-metallicity locus for a sample of $z \le 0.22$ early-type galaxies in ALHAMBRA at different stellar mass bins are in very good agreement with the ones from SDSS spectroscopic diagnostics. Moreover, a one-to-one comparison between the redshifts, ages, metallicities, and stellar masses derived spectroscopically for SDSS and by MUFFIT for ALHAMBRA reveals good qualitative agreements in all the parameters. In addition, and using as input the results from photometric-redshift codes, MUFFIT improves the photometric-redshift accuracy by $\sim 10$-$20\%$, and it also detects nebular emissions in galaxies, providing physical information about their strengths. Our results show the potential of multi-filter galaxy data to conduct reliable stellar population studies with the appropiate analysis techniques, as MUFFIT.

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\sigma$ 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.

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.

We report on our first set of spectroscopic Hubble Space Telescope observations of the z~11 candidate galaxy strongly lensed by the MACSJ0647.7+7015 galaxy cluster. The three lensed images are faint and we show that these early slitless grism observations are of sufficient depth to investigate whether this high-redshift candidate, identified by its strong photometric break at ~1.5 micron, could possibly be an emission line galaxy at a much lower redshift. While such an interloper would imply the existence of a rather peculiar object, we show here that such strong emission lines would clearly have been detected. Comparing realistic, two-dimensional simulations to these new observations we would expect the necessary emission lines to be detected at >5 sigma while we see no evidence for such lines in the dispersed data of any of the three lensed images. We therefore exclude that this object could be a low redshift emission line interloper, which significantly increases the likelihood of this candidate being a bona fide z~11 galaxy.

Context. Most observational results on the high redshift restframe UV-bright galaxies are based on samples pinpointed using the so called dropout technique or Ly-alpha selection. However, the availability of multifilter data allows now replacing the dropout selections by direct methods based on photometric redshifts. In this paper we present the methodology to select and study the population of high redshift galaxies in the ALHAMBRA survey data. Aims. Our aim is to develop a less biased methodology than the traditional dropout technique to study the high redshift galaxies in ALHAMBRA and other multifilter data. Thanks to the wide area ALHAMBRA covers, we especially aim at contributing in the study of the brightest, less frequent, high redshift galaxies. Methods. The methodology is based on redshift probability distribution functions (zPDFs). It is shown how a clean galaxy sample can be obtained by selecting the galaxies with high integrated probability of being within a given redshift interval. However, reaching both a complete and clean sample with this method is challenging. Hence, a method to derive statistical properties by summing the zPDFs of all the galaxies in the redshift bin of interest is introduced. Results. Using this methodology we derive the galaxy rest frame UV number counts in five redshift bins centred at z=2.5, 3.0, 3.5, 4.0, and 4.5, being complete up to the limiting magnitude at m_UV(AB)=24. With the wide field ALHAMBRA data we especially contribute in the study of the brightest ends of these counts, sampling well the surface densities down to m_UV(AB)=21-22. Conclusions. We show that using the zPDFs it is easy to select a clean sample of high redshift galaxies. We also show that statistical analysis of the properties of galaxies is better done using a probabilistic approach, which takes into account both the incompleteness and contamination in a natural way.

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, $\kappa$, and $\sim20\%$ in the magnification, $\mu$. 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 (\Sigma)/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.

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.

We present a detailed mass reconstruction and a novel study on the substructure properties in the core of the CLASH and Frontier Fields galaxy cluster MACS J0416.1-2403. We show and employ our extensive spectroscopic data set taken with the VIMOS instrument as part of our CLASH-VLT program, to confirm spectroscopically 10 strong lensing systems and to select a sample of 175 plausible cluster members to a limiting stellar mass of log(M_*/M_Sun) ~ 8.6. We reproduce the measured positions of 30 multiple images with a remarkable median offset of only 0.3" by means of a comprehensive strong lensing model comprised of 2 cluster dark-matter halos, represented by cored elliptical pseudo-isothermal mass distributions, and the cluster member components. The latter have total mass-to-light ratios increasing with the galaxy HST/WFC3 near-IR (F160W) luminosities. The measurement of the total enclosed mass within the Einstein radius is accurate to ~5%, including systematic uncertainties. We emphasize that the use of multiple-image systems with spectroscopic redshifts and knowledge of cluster membership based on extensive spectroscopic information is key to constructing robust high-resolution mass maps. We also produce magnification maps over the central area that is covered with HST observations. We investigate the galaxy contribution, both in terms of total and stellar mass, to the total mass budget of the cluster. When compared with the outcomes of cosmological $N$-body simulations, our results point to a lack of massive subhalos in the inner regions of simulated clusters with total masses similar to that of MACS J0416.1-2403. Our findings of the location and shape of the cluster dark-matter halo density profiles and on the cluster substructures provide intriguing tests of the assumed collisionless, cold nature of dark matter and of the role played by baryons in the process of structure formation.

We examine the nonlinear structure of gravitationally collapsed objects that form in our simulations of wavelike cold dark matter ($\psi$DM), described by the Schr\"{o}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.

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.

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.

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.

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.

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 $\alpha_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.

Our goal is to develop and test a novel methodology to compute accurate close pair fractions with photometric redshifts. We improve the current methodologies to estimate the merger fraction f_m from photometric redshifts by (i) using the full probability distribution functions (PDFs) of the sources in redshift space, (ii) including the variation in the luminosity of the sources with z in both the selection of the samples and in the luminosity ratio constrain, and (iii) splitting individual PDFs into red and blue spectral templates to deal robustly with colour selections. We test the performance of our new methodology with the PDFs provided by the ALHAMBRA photometric survey. The merger fractions and rates from the ALHAMBRA survey are in excellent agreement with those from spectroscopic work, both for the general population and for red and blue galaxies. With the merger rate of bright (M_B <= -20 - 1.1z) galaxies evolving as (1+z)^n, the power-law index n is larger for blue galaxies (n = 2.7 +- 0.5) than for red galaxies (n = 1.3 +- 0.4), confirming previous results. Integrating the merger rate over cosmic time, we find that the average number of mergers per galaxy since z = 1 is N_m = 0.57 +- 0.05 for red galaxies and N_m = 0.26 +- 0.02 for blue galaxies. Our new methodology exploits statistically all the available information provided by photometric redshift codes and provides accurate measurements of the merger fraction by close pairs only using photometric redshifts. Current and future photometric surveys will benefit of this new methodology.

We utilise a two-color Lyman-Break selection criterion to search for z~9-10 galaxies over the first 19 clusters in the CLASH program. A systematic search yields three z~9-10 candidates. While we have already reported the most robust of these candidates, MACS1149-JD, two additional z~9 candidates are also found and have H_{160}-band magnitudes of ~26.2-26.9. A careful assessment of various sources of contamination suggests <~1 contaminants for our z~9-10 selection. To determine the implications of these search results for the LF and SFR density at z~9, we introduce a new differential approach to deriving these quantities in lensing fields. Our procedure is to derive the evolution by comparing the number of z~9-10 galaxy candidates found in CLASH with the number of galaxies in a slightly lower redshift sample (after correcting for the differences in selection volumes), here taken to be z~8. This procedure takes advantage of the fact that the relative volumes available for the z~8 and z~9-10 selections behind lensing clusters are not greatly dependent on the details of the lensing models. We find that the normalization of the UV LF at z~9 is just 0.28_{-0.20}^{+0.39}\times that at z~8, ~1.4_{-0.8}^{+3.0}x lower than extrapolating z~4-8 LF results. While consistent with the evolution in the UV LF seen at z~4-8, these results marginally favor a more rapid evolution at z>8. Compared to similar evolutionary findings from the HUDF, our result is less insensitive to large-scale structure uncertainties, given our many independent sightlines on the high-redshift universe.

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 \AA$) 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$.

We utilize 16 band Hubble Space Telescope (HST) observations of 18 lensing clusters obtained as part of the Cluster Lensing And Supernova survey with Hubble (CLASH) Multi-Cycle Treasury program to search for $z\sim6-8$ galaxies. We report the discovery of 204, 45, and 13 Lyman-break galaxy candidates at $z\sim6$, $z\sim7$, and $z\sim8$, respectively, identified from purely photometric redshift selections. This large sample, representing nearly an order of magnitude increase in the number of magnified star-forming galaxies at $z\sim 6-8$ presented to date, is unique in that we have observations in four WFC3/UVIS UV, seven ACS/WFC optical, and all five WFC3/IR broadband filters, which enable very accurate photometric redshift selections. We construct detailed lensing models for 17 of the 18 clusters to estimate object magnifications and to identify two new multiply lensed $z \gtrsim 6$ candidates. The median magnifications over the 17 clusters are 4, 4, and 5 for the $z\sim6$, $z\sim7$, and $z\sim8$ samples, respectively, over an average area of 4.5 arcmin$^2$ per cluster. We compare our observed number counts with expectations based on convolving "blank" field UV luminosity functions through our cluster lens models and find rough agreement down to $\sim27$ mag, where we begin to suffer significant incompleteness. In all three redshift bins, we find a higher number density at brighter observed magnitudes than the field predictions, empirically demonstrating for the first time the enhanced efficiency of lensing clusters over field surveys. Our number counts also are in general agreement with the lensed expectations from the cluster models, especially at $z\sim6$, where we have the best statistics.

We report the discovery of 24 Lyman-break candidates at 7<z<10.5, in the Hubble Frontier Fields (HFF) imaging data of Abell 2744 (z=0.308), plus Spiter/IRAC data and archival ACS data. The sample includes a triple image system with a photometric redshift of z~7.4. This high redshift is geometrically confirmed by our lens model corresponding to deflection angles that are 12% larger than the lower-redshift systems used to calibrate the lens model at z=2.019. The majority of our high-redshift candidates are not expected to be multiply lensed given their locations in the image plane and the brightness of foreground galaxies, but are magnified by factors of ~1.3-15, so that we are seeing further down the luminosity function than comparable deep field imaging. It is apparent that the redshift distribution of these sources does not smoothly extend over the full redshift range accessible at z<12, but appears to break above z=9. Nine candidates are clustered within a small region of 20" across, representing a potentially unprecedented concentration. Given the poor statistics, however, we must await similar constraints from the additional HFF clusters to properly examine this trend. The physical properties of our candidates are examined using the range of lens models developed for the HFF program by various groups including our own, for a better estimate of underlying systematics. Our spectral-energy-distribution fits for the brightest objects suggest stellar masses of ~1E9 solar mass, star-formation rates of ~4 solar mass per yr, and a typical formation redshift of z<19.

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.

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 ($\psi$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 $\psi$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 $\psi$DM is substantially delayed relative to CDM, appearing at $\bf{z\lesssim 13}$ in our simulations.

We present a free form mass reconstruction of the massive lensing cluster MACSJ0416.1-2403 using the latest Hubble Frontier Fields data. Our model independent method finds that the extended lensing pattern is generated by two elongated, closely projected clusters of similar mass. Our lens model identifies new lensed images with which we improve the accuracy of the dark matter distribution. We find that the bimodal mass distribution is nearly coincident with the bimodal X-ray emission, but with the two dark matter peaks lying closer together than the centroids of the X-ray emisison. We show this can be achieved if the collision has occurred close to the plane and such that the cores are deflected around each other. The projected mass profiles of both clusters are well constrained because of the many interior lensed images, leading to surprisingly flat mass profiles of both components in the region 15-100 kpc. We discuss the extent to which this may be generated by tidal forces in our dynamical model which are large during an encounter of this type as the cores "graze" each other. The relative velocity between the two cores is estimated to be about 1200 km/s and mostly along the line of sight so that our model is consistent with the relative redshift difference between the two cD galaxies (dz = 0.04).

The ALHAMBRA (Advance Large Homogeneous Area Medium Band Redshift Astronomical) survey has observed 8 different regions of the sky, including sections of the COSMOS, DEEP2, ELAIS, GOODS-N, SDSS and Groth fields using a new photometric system with 20 contiguous ~ $300\AA$ filters covering the optical range, combining them with deep $JHKs$ imaging. The observations, carried out with the Calar Alto 3.5m telescope using the wide field (0.25 sq. deg FOV) optical camera LAICA and the NIR instrument Omega-2000, correspond to ~700hrs on-target science images. The photometric system was designed to maximize the effective depth of the survey in terms of accurate spectral-type and photo-zs estimation along with the capability of identification of relatively faint emission lines. Here we present multicolor photometry and photo-zs for ~438k galaxies, detected in synthetic F814W images, complete down to I~24.5 AB, taking into account realistic noise estimates, and correcting by PSF and aperture effects with the ColorPro software. The photometric ZP have been calibrated using stellar transformation equations and refined internally, using a new technique based on the highly robust photometric redshifts measured for emission line galaxies. We calculate photometric redshifts with the BPZ2 code, which includes new empirically calibrated templates and priors. Our photo-zs have a precision of $dz/(1+z_s)=1%$ for I<22.5 and 1.4% for 22.5<I<24.5. Precisions of less than 0.5% are reached for the brighter spectroscopic sample, showing the potential of medium-band photometric surveys. The global $P(z)$ shows a mean redshift <z>=0.56 for I<22.5 AB and <z>=0.86 for I<24.5 AB. The data presented here covers an effective area of 2.79 sq. deg, split into 14 strips of 58.5'x15.5' and represents ~32 hrs of on-target.

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.

We study the clustering of galaxies as function of luminosity and redshift in the range $0.35 < z < 1.25$ using data from the Advanced Large Homogeneous Area Medium Band Redshift Astronomical (ALHAMBRA) survey. The ALHAMBRA data used in this work cover $2.38 \mathrm{deg}^2$ in 7 independent fields, after applying a detailed angular selection mask, with accurate photometric redshifts, $\sigma_z \lesssim 0.014 (1+z)$, down to $I_{\rm AB} < 24$. Given the depth of the survey, we select samples in $B$-band luminosity down to $L^{\rm th} \simeq 0.16 L^{*}$ at $z = 0.9$. We measure the real-space clustering using the projected correlation function, accounting for photometric redshifts uncertainties. We infer the galaxy bias, and study its evolution with luminosity. We study the effect of sample variance, and confirm earlier results that the COSMOS and ELAIS-N1 fields are dominated by the presence of large structures. For the intermediate and bright samples, $L^{\rm med} \gtrsim 0.6L^{*}$, we obtain a strong dependence of bias on luminosity, in agreement with previous results at similar redshift. We are able to extend this study to fainter luminosities, where we obtain an almost flat relation, similar to that observed at low redshift. Regarding the evolution of bias with redshift, our results suggest that the different galaxy populations studied reside in haloes covering a range in mass between $\log_{10}[M_{\rm h}/(h^{-1}\mathrm{M}_{\odot})] \gtrsim 11.5$ for samples with $L^{\rm med} \simeq 0.3 L^{*}$ and $\log_{10}[M_{\rm h}/(h^{-1}\mathrm{M}_{\odot})] \gtrsim 13.0$ for samples with $L^{\rm med} \simeq 2 L^{*}$, with typical occupation numbers in the range of $\sim 1 - 3$ galaxies per halo.

We present a new determination of the concentration-mass relation for galaxy clusters based on our comprehensive lensing analysis of 19 X-ray selected galaxy clusters from the Cluster Lensing and Supernova Survey with Hubble (CLASH). Our sample spans a redshift range between 0.19 and 0.89. We combine weak lensing constraints from the Hubble Space Telescope (HST) and from ground-based wide field data with strong lensing constraints from HST. The result are reconstructions of the surface-mass density for all CLASH clusters on multi-scale grids. Our derivation of NFW parameters yields virial masses between 0.53 x 10^15 and 1.76 x 10^15 M_sol/h and the halo concentrations are distributed around c_200c ~ 3.7 with a 1-sigma significant negative trend with cluster mass. We find an excellent 4% agreement between our measured concentrations and the expectation from numerical simulations after accounting for the CLASH selection function based on X-ray morphology. The simulations are analyzed in 2D to account for possible biases in the lensing reconstructions due to projection effects. The theoretical concentration-mass (c-M) relation from our X-ray selected set of simulated clusters and the c-M relation derived directly from the CLASH data agree at the 90% confidence level.

We present the results of a numerical study based on the analysis of the MUSIC-2 simulations, aimed at estimating the expected concentration-mass relation for the CLASH cluster sample. We study nearly 1400 halos simulated at high spatial and mass resolution, which were projected along many lines-of-sight each. We study the shape of both their density and surface-density profiles and fit them with a variety of radial functions, including the Navarro-Frenk-White, the generalised Navarro-Frenk-White, and the Einasto density profiles. We derive concentrations and masses from these fits and investigate their distributions as a function of redshift and halo relaxation. We use the X-ray image simulator X-MAS to produce simulated Chandra observations of the halos and we use them to identify objects resembling the X-ray morphologies and masses of the clusters in the CLASH X-ray selected sample. We also derive a concentration-mass relation for strong-lensing clusters. We find that the sample of simulated halos which resemble the X-ray morphology of the CLASH clusters is composed mainly by relaxed halos, but it also contains a significant fraction of un-relaxed systems. For such a sample we measure an average 2D concentration which is ~11% higher than found for the full sample of simulated halos. After accounting for projection and selection effects, the average NFW concentrations of CLASH clusters are expected to be intermediate between those predicted in 3D for relaxed and super-relaxed halos. Matching the simulations to the individual CLASH clusters on the basis of the X-ray morphology, we expect that the NFW concentrations recovered from the lensing analysis of the CLASH clusters are in the range [3-6], with an average value of 3.87 and a standard deviation of 0.61. Simulated halos with X-ray morphologies similar to those of the CLASH clusters are affected by a modest orientation bias.

The Javalambre-Physics of the Accelerated Universe Astrophysical Survey (J-PAS) is a narrow band, very wide field Cosmological Survey to be carried out from the Javalambre Observatory in Spain with a purpose-built, dedicated 2.5m telescope and a 4.7 sq.deg. camera with 1.2Gpix. Starting in late 2015, J-PAS will observe 8500sq.deg. of Northern Sky and measure $0.003(1+z)$ photo-z for $9\times10^7$ LRG and ELG galaxies plus several million QSOs, sampling an effective volume of $\sim 14$ Gpc$^3$ up to $z=1.3$ and becoming the first radial BAO experiment to reach Stage IV. J-PAS will detect $7\times 10^5$ galaxy clusters and groups, setting constrains on Dark Energy which rival those obtained from its BAO measurements. Thanks to the superb characteristics of the site (seeing ~0.7 arcsec), J-PAS is expected to obtain a deep, sub-arcsec image of the Northern sky, which combined with its unique photo-z precision will produce one of the most powerful cosmological lensing surveys before the arrival of Euclid. J-PAS unprecedented spectral time domain information will enable a self-contained SN survey that, without the need for external spectroscopic follow-up, will detect, classify and measure $\sigma_z\sim 0.5\%$ redshifts for $\sim 4000$ SNeIa and $\sim 900$ core-collapse SNe. The key to the J-PAS potential is its innovative approach: a contiguous system of 54 filters with $145\AA$ width, placed $100\AA$ apart over a multi-degree FoV is a powerful "redshift machine", with the survey speed of a 4000 multiplexing low resolution spectrograph, but many times cheaper and much faster to build. The J-PAS camera is equivalent to a 4.7 sq.deg. "IFU" and it will produce a time-resolved, 3D image of the Northern Sky with a very wide range of Astrophysical applications in Galaxy Evolution, the nearby Universe and the study of resolved stellar populations.

We aim at constraining the assembly history of clusters by studying the intra cluster light (ICL) properties, estimating its contribution to the fraction of baryons in stars, f*, and understanding possible systematics/bias using different ICL detection techniques. We developed an automated method, GALtoICL, based on the software GALAPAGOS to obtain a refined version of typical BCG+ICL maps. We applied this method to our test case MACS J1206.2-0847, a massive cluster located at z=0.44, that is part of the CLASH sample. Using deep multi-band SUBARU images, we extracted the surface brightness (SB) profile of the BCG+ICL and we studied the ICL morphology, color, and contribution to f* out to R500. We repeated the same analysis using a different definition of the ICL, SBlimit method, i.e. a SB cut-off level, to compare the results. The most peculiar feature of the ICL in MACS1206 is its asymmetric radial distribution, with an excess in the SE direction and extending towards the 2nd brightest cluster galaxy which is a Post Starburst galaxy. This suggests an interaction between the BCG and this galaxy that dates back to t <= 1.5 Gyr. The BCG+ICL stellar content is 8% of M_(*,500) and the (de-) projected baryon fraction in stars is f*=0.0177 (0.0116), in excellent agreement with recent results. The SBlimit method provides systematically higher ICL fractions and this effect is larger at lower SB limits. This is due to the light from the outer envelopes of member galaxies that contaminate the ICL. Though more time consuming, the GALtoICL method provides safer ICL detections that are almost free of this contamination. This is one of the few ICL study at redshift z > 0.3. At completion, the CLASH/VLT program will allow us to extend this analysis to a statistically significant cluster sample spanning a wide redshift range: 0.2<z<0.6.

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.

Hubble Space Telescope imaging of the galaxy cluster Abell 1689 has revealed an exceptional number of strongly lensed multiply-imaged galaxies, including high-redshift candidates. Previous studies have used this data to obtain the most detailed dark matter reconstructions of any galaxy cluster to date, resolving substructures ~25 kpc across. We examine Abell 1689 (hereafter, A1689) non-parametrically, combining strongly lensed images and weak distortions from wider field Subaru imaging, and we incorporate member galaxies to improve the lens solution. Strongly lensed galaxies are often locally affected by member galaxies, however, these perturbations cannot be recovered in grid based reconstructions because the lensing information is too sparse to resolve member galaxies. By adding luminosity-scaled member galaxy deflections to our smooth grid we can derive meaningful solutions with sufficient accuracy to permit the identification of our own strongly lensed images, so our model becomes self consistent. We identify 11 new multiply lensed system candidates and clarify previously ambiguous cases, in the deepest optical and NIR data to date from Hubble and Subaru. Our improved spatial resolution brings up new features not seen when the weak and strong lensing effects are used separately, including clumps and filamentary dark matter around the main halo. Our treatment means we can obtain an objective mass ratio between the cluster and galaxy components, for examining the extent of tidal stripping of the luminous member galaxies. We find a typical mass-to-light ratios of M/L_B = 21 inside the r<1 arcminute region that drops to M/L_B = 17 inside the r<40 arcsecond region. Our model independence means we can objectively evaluate the competitiveness of stacking cluster lenses for defining the geometric lensing-distance-redshift relation in a model independent way.

We present a strong lensing system in which a double source is imaged 5 times by 2 early-type galaxies. We take advantage in this target of the multi-band photometry obtained as part of the CLASH program, complemented by the spectroscopic data of the VLT/VIMOS and FORS2 follow-up campaign. We use a photometric redshift of 3.7 for the source and confirm spectroscopically the membership of the 2 lenses to the galaxy cluster MACS J1206.2-0847 at redshift 0.44. We exploit the excellent angular resolution of the HST/ACS images to model the 2 lenses in terms of singular isothermal sphere profiles and derive robust effective velocity dispersions of (97 +/- 3) and (240 +/- 6) km/s. The total mass distribution of the cluster is also well characterized by using only the local information contained in this lensing system, that is located at a projected distance of more than 300 kpc from the cluster luminosity center. According to our best-fitting lensing and composite stellar population models, the source is magnified by a total factor of 50 and has a luminous mass of about (1.0 +/- 0.5) x 10^{9} M_{Sun}. By combining the total and luminous mass estimates of the 2 lenses, we measure luminous over total mass fractions projected within the effective radii of 0.51 +/- 0.21 and 0.80 +/- 0.32. With these lenses we can extend the analysis of the mass properties of lens early-type galaxies by factors that are about 2 and 3 times smaller than previously done with regard to, respectively, velocity dispersion and luminous mass. The comparison of the total and luminous quantities of our lenses with those of astrophysical objects with different physical scales reveals the potential of studies of this kind for investigating the internal structure of galaxies. These studies, made possible thanks to the CLASH survey, will allow us to go beyond the current limits posed by the available lens samples in the field.

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\rho)=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 $\rho$ is the density. We plan to further improve our constraint by applying the same procedure to all clusters from the ongoing CLASH-VLT program.

Our goal is to estimate empirically, for the first time, the cosmic variance that affects merger fraction studies based on close pairs. We compute the merger fraction from photometric redshift close pairs with 10h^-1 kpc <= rp <= 50h^-1 kpc and Dv <= 500 km/s, and measure it in the 48 sub-fields of the ALHAMBRA survey. We study the distribution of the measured merger fractions, that follow a log-normal function, and estimate the cosmic variance sigma_v as the intrinsic dispersion of the observed distribution. We develop a maximum likelihood estimator to measure a reliable sigma_v and avoid the dispersion due to the observational errors (including the Poisson shot noise term). The cosmic variance of the merger fraction depends mainly on (i) the number density of the populations under study, both for the principal (n_1) and the companion (n_2) galaxy in the close pair, and (ii) the probed cosmic volume V_c. We find a significant dependence on neither the search radius used to define close companions, the redshift, nor the physical selection (luminosity or stellar mass) of the samples. We provide a parametrisation of the cosmic variance with n_1, n_2, and V_c, sigma_v = 0.48 n_1^{-0.54} V_c^{-0.48} (n_2/n_1)^{-0.37}. Thanks to this prescription, future merger fraction studies based on close pairs could account properly for the cosmic variance on their results.

The Cluster Lensing And Supernovae survey with Hubble (CLASH) is an Hubble Space Telescope (HST) Multi-Cycle Treasury program observing 25 massive galaxy clusters. CLASH observations are carried out in 16 bands from UV to NIR to derive accurate and reliable estimates of photometric redshifts. We present the CLASH photometric redshifts and study the photometric redshift accuracy of the arcs in more detail for the case of MACS1206.2-0847. We use the publicly available Le Phare and BPZ photometric redshift codes on 17 CLASH galaxy clusters. Using Le Phare code for objects with StoN>=10, we reach a precision of 3%(1+z) for the strong lensing arcs, which is reduced to 2.4%(1+z) after removing outliers. For galaxies in the cluster field the corresponding values are 4%(1+z) and 3%(1+z). Using mock galaxy catalogues, we show that 3%(1+z) precision is what one would expect from the CLASH photometry when taking into account extinction from dust, emission lines and the finite range of SEDs included in the photo-z template library. We study photo-z results for different aperture photometry and find that the SExtractor isophotal photometry works best. Le Phare and BPZ give similar photo-z results for the strong lensing arcs as well as galaxies of the cluster field. Results are improved when optimizing the photometric aperture shape showing an optimal aperture size around 1" radius giving results which are equivalent to isophotal photometry. Tailored photometry of the arcs improve the photo-z results.

We present a quintuply lensed z ~ 6 candidate discovered in the field of the galaxy cluster RXC J2248.7-4431 (z ~ 0.348) targeted within the Cluster Lensing and Supernova survey with Hubble (CLASH) and selected in the deep HST Frontier Fields survey. Thanks to the CLASH 16-band HST imaging, we identify the quintuply lensed z ~ 6 candidate as an optical dropout in the inner region of the cluster, the brightest image having magAB=24.81+-0.02 in the f105w filter. We perform a detailed photometric analysis to verify its high-z and lensed nature. We get as photometric redshift z_phot ~ 5.9, and given the extended nature and NIR colours of the lensed images, we rule out low-z early type and galactic star contaminants. We perform a strong lensing analysis of the cluster, using 13 families of multiple lensed images identified in the HST images. Our final best model predicts the high-z quintuply lensed system with a position accuracy of 0.8''. The magnifications of the five images are between 2.2 and 8.3, which leads to a delensed UV luminosity of L_1600 ~ 0.5L*_1600 at z=6. We also estimate the UV slope from the observed NIR colours, finding a steep beta=-2.89+-0.38. We use singular and composite stellar population SEDs to fit the photometry of the hiz candidate, and we conclude that it is a young (age <300 Myr) galaxy with mass of M ~ 10^8Msol, subsolar metallicity (Z<0.2Zsol) and low dust content (AV ~ 0.2-0.4).

We present VIsible Multi-Object Spectrograph (VIMOS) observations of a z 6 galaxy quintuply imaged by the Frontier Fields galaxy cluster RXC J2248.7-4431 (z=0.348). This sub-L^*, high-z galaxy has been recently discovered by Monna et al. (2013) using dropout techniques with the 16-band HST photometry acquired as part of the Cluster Lensing And Supernova survey with Hubble (CLASH). Obtained as part of the CLASH-VLT survey, the VIMOS medium-resolution spectra of this source show a very faint continuum between ~8700A and ~9300A and a prominent emission line at 8643A, which can be readily identified with Lyman-alpha at z=6.110. The emission line exhibits an asymmetric profile, with a more pronounced red wing. The rest-frame equivalent width of the line is EW=79+-10A. After correcting for magnification, the star-formation rate (SFR) estimated from the Lya line is SFR(Lya)=11 M_{sol}/yr and that estimated from the UV data is SFR(UV)=3 M_{sol}/yr. We estimate that the effective radius of the source is R_e<~0.4 kpc, which implies a star formation surface mass density Sigma_{SFR}>6 M_{sol}/yr/kpc^2 and, using the Kennicutt-Schmidt relation, a gas surface mass density Sigma_{gas}>10^3 M_{sol}/pc^2. Our results support the idea that this magnified, distant galaxy is a young and compact object with 0.4 L^* at z=6, with comparable amount of mass in gas and stars. Future follow-up observations with ALMA will provide valuable insight into the SFR and molecular gas content of this source. In the spirit of the Frontier Fields initiative, we also publish the redshifts of several multiply imaged sources and other background objects which will help improving the strong lensing model of this galaxy cluster.

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\sigma_\kappa. 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{\sigma} departure above the maximum simulated halo mass at this redshift.

We provide a new observational test for a key prediction of the \Lambda 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 $\sigma$ agreement with \Lambda 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.

ALHAMBRA is a photometric survey designed to trace the cosmic evolution and cosmic variance. It covers a large area of ~ 4 sq. deg in 8 fields, where 7 fields overlap with other surveys, allowing to have complementary data in other wavelengths. All observations were carried out in 20 continuous, medium band (30 nm width) optical and 3 near-infrared (JHK) bands, providing the precise measurements of photometric redshifts. In addition, morphological classification of galaxies is crucial for any kind of galaxy formation and cosmic evolution studies, providing the information about star formation histories, their environment and interactions, internal perturbations, etc. We present a morphological classification of > 40,000 galaxies in the ALHAMBRA survey. We associate to every galaxy a probability to be early-type using the automated Bayesian code galSVM. Despite of the spatial resolution of the ALHAMBRA images (~ 1 arcsec), for 22,051 galaxies we obtained the contamination by other type of less than 10%. Of those, 1,640 and 10,322 galaxies are classified as early- (down to redshifts ~ 0.5) and late-type (down to redshifts ~ 1.0), respectively, with magnitudes F613W < 22.0. In addition, for magnitude range 22.0 < F613W < 23.0 we classified other 10,089 late-type galaxies with redshifts < 1.3. We show that the classified objects populate the expected regions in the colour-mass and colour-magnitude planes. The presented dataset is especially attractive given the homogeneous multi-wavelength coverage available in the ALHAMBRA fields, and is intended to be used in a variety of scientific applications. The low-contamination catalogue (< 10%) is made publicly available with the present paper.

We use an unprecedented data-set of about 600 redshifts for cluster members, obtained as part of a VLT/VIMOS large programme, to constrain the mass profile of the z=0.44 cluster MACS J1206.2-0847 over the radial range 0-5 Mpc (0-2.5 virial radii) using the MAMPOSSt and Caustic methods. We then add external constraints from our previous gravitational lensing analysis. We invert the Jeans equation to obtain the velocity-anisotropy profiles of cluster members. With the mass-density and velocity-anisotropy profiles we then obtain the first determination of a cluster pseudo-phase-space density profile. The kinematics and lensing determinations of the cluster mass profile are in excellent agreement. This is very well fitted by a NFW model with mass M200=(1.4 +- 0.2) 10^15 Msun and concentration c200=6 +- 1, only slightly higher than theoretical expectations. Other mass profile models also provide acceptable fits to our data, of (slightly) lower (Burkert, Hernquist, and Softened Isothermal Sphere) or comparable (Einasto) quality than NFW. The velocity anisotropy profiles of the passive and star-forming cluster members are similar, close to isotropic near the center and increasingly radial outside. Passive cluster members follow extremely well the theoretical expectations for the pseudo-phase-space density profile and the relation between the slope of the mass-density profile and the velocity anisotropy. Star-forming cluster members show marginal deviations from theoretical expectations. This is the most accurate determination of a cluster mass profile out to a radius of 5 Mpc, and the only determination of the velocity-anisotropy and pseudo-phase-space density profiles of both passive and star-forming galaxies for an individual cluster [abridged]

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.

Growing observational evidence now indicates that nebular line emission has a significant impact on the rest-frame optical fluxes of z~5-7 galaxies observed with Spitzer. This line emission makes z~5-7 galaxies appear more massive, with lower specific star formation rates. However, corrections for this line emission have been very difficult to perform reliably due to huge uncertainties on the overall strength of such emission at z>~5.5. Here, we present the most direct observational evidence yet for ubiquitous high-EW [OIII]+Hbeta line emission in Lyman-break galaxies at z~7, while also presenting a strategy for an improved measurement of the sSFR at z~7. We accomplish this through the selection of bright galaxies in the narrow redshift window z~6.6-7.0 where the IRAC 4.5 micron flux provides a clean measurement of the stellar continuum light. Observed 4.5 micron fluxes in this window contrast with the 3.6 micron fluxes which are contaminated by the prominent [OIII]+Hbeta lines. To ensure a high S/N for our IRAC flux measurements, we consider only the brightest (H_{160}<26 mag) magnified galaxies we have identified in CLASH and other programs targeting galaxy clusters. Remarkably, the mean rest-frame optical color for our bright seven-source sample is very blue, [3.6]-[4.5]=-0.9+/-0.3. Such blue colors cannot be explained by the stellar continuum light and require that the rest-frame EW of [OIII]+Hbeta be greater than 637 Angstroms for the average source. The bluest four sources from our seven-source sample require an even more extreme EW of 1582 Angstroms. Our derived lower limit for the mean [OIII]+Hbeta EW could underestimate the true EW by ~2x based on a simple modeling of the redshift distribution of our sources. We can also set a robust lower limit of >~4 Gyr^-1 on the specific star formation rates based on the mean SED for our seven-source sample. (abridged)

We aim to illustrate the potentiality of the Advanced Large, Homogeneous Area, Medium-Band Redshift Astronomical (ALHAMBRA) survey to investigate the high redshift universe through the detection of quasi stellar objects (QSOs) at redshifts larger than 5. The search for z>5 QSOs candidates was done by fitting an extensive library of spectral energy distributions --including active and non-active galaxy templates as well as stars-- to the photometric database of the ALHAMBRA survey (composed of 20 optical medium-band plus the 3 broad-band JHKs filters). Our selection over ~1 square degree of ALHAMBRA data (~1/4 of the total area covered by the survey), combined with GTC/OSIRIS spectroscopy, has yielded the identification of an optically faint QSO at very high redshift (z = 5.41). The QSO has an absolute magnitude of ~-24 at the 1450{\AA} continuum, a bolometric luminosity of ~2x10^46 erg/s and an estimated black hole mass of ~10^8 Msolar. This QSO adds itself to a reduced number of known UV faint sources at these redshifts. The preliminary derived space density is compatible with the most recent determinations of the high-z QSO luminosity functions (QLF). This new detection shows how ALHAMBRA, as well as forthcoming well designed photometric surveys, can provide a wealth of information on the origin and early evolution of this kind of objects.

We report the PACS-100um/160um detections of a sample of 42 GALEX-selected and FIR-detected Lyman break galaxies (LBGs) at z ~ 1 located in the COSMOS field and analyze their ultra-violet (UV) to far-infrared (FIR) properties. The detection of these LBGs in the FIR indicates that they have a dust content high enough so that its emission can be directly detected. According to a spectral energy distribution (SED) fitting with stellar population templates to their UV-to-near-IR observed photometry, PACS-detected LBGs tend to be bigger, more massive, dustier, redder in the UV continuum, and UV-brighter than PACS-undetected LBGs. PACS-detected LBGs at z ~ 1 are mostly disk-like galaxies and are located over the green-valley and red sequence of the color-magnitude diagram of galaxies at their redshift. By using their UV and IR emission, we find that PACS-detected LBGs tend to be less dusty and have slightly higher total star-formation rates (SFRs) than other PACS-detected UV-selected galaxies within their same redshift range. As a consequence of the selection effect due to the depth of the FIR observations employed, all our PACS-detected LBGs are LIRGs. However, none of them are in the ULIRG regime, where the FIR observations are complete. The finding of ULIRGs-LBGs at higher redshifts suggests an evolution of the FIR emission of LBGs with cosmic time. In an IRX-$\beta$ diagram, PACS-detected LBGs at z ~ 1 tend to be located around the relation for local starburst similarly to other UV-selected PACS-detected galaxies at their same redshift. Consequently, the dust-correction factors obtained with their UV continuum slope allow to determine their total SFR, unlike at higher redshifts. However, the dust attenuation derived from UV to NIR SED fitting overestimates the total SFR for most of our PACS-detected LBGs in age-dependent way: the overestimation factor is higher in younger galaxies.

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.

Aims. We investigate the degree of improvement in dark energy constraints that can be achieved by extending Type Ia Supernova (SN Ia) samples to redshifts z > 1.5 with the Hubble Space Telescope (HST), particularly in the ongoing CANDELS and CLASH multi-cycle treasury programs. Methods. Using the popular CPL parametrization of the dark energy, w = w0 +wa(1-a), we generate mock SN Ia samples that can be projected out to higher redshifts. The synthetic datasets thus generated are fitted to the CPL model, and we evaluate the improvement that a high-z sample can add in terms of ameliorating the statistical and systematic uncertainties on cosmological parameters. Results. In an optimistic but still very achievable scenario, we find that extending the HST sample beyond CANDELS+CLASH to reach a total of 28 SN Ia at z > 1.0 could improve the uncertainty in the wa parameter by up to 21%. The corresponding improvement in the figure of merit (FoM) would be as high as 28%. Finally, we consider the use of high-redshift SN Ia samples to detect non-cosmological evolution in SN Ia luminosities with redshift, finding that such tests could be undertaken by future spacebased infrared surveys using the James Webb Space Telescope (JWST).

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 $\sigma=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.

We take advantage of the exceptional photometric coverage provided by the combination of GALEX data in the UV and the ALHAMBRA survey in the optical and near-IR to analyze the physical properties of a sample of 1225 GALEX-selected Lyman break galaxies (LBGs) at $0.8 \lesssim z \lesssim 1.2$ located in the COSMOS field. This is the largest sample of LBGs studied at that redshift range so far. According to a spectral energy distribution (SED) fitting with synthetic stellar population templates, we find that LBGs at $z \sim 1$ are mostly young galaxies with a median age of 341 Myr and have intermediate dust attenuation, $\ < E_s (B-V) \ > \sim 0.20$. Due to their selection criterion, LBGs at $z \sim 1$ are UV-bright galaxies and have high dust-corrected total SFR, with a median value of 16.9 $M_\odot {\rm yr}^{-1}$. Their median stellar mass is $\log{\left(M_*/M_\odot \right)} = 9.74$. We obtain that the dust-corrected total SFR of LBGs increases with stellar mass and the specific SFR is lower for more massive galaxies. Only 2% of the galaxies selected through the Lyman break criterion have an AGN nature. LBGs at $z \sim 1$ are mostly located over the blue cloud of the color-magnitude diagram of galaxies at their redshift, with only the oldest and/or the dustiest deviating towards the green valley and red sequence. Morphologically, 69% of LBGs are disk-like galaxies, with the fraction of interacting, compact, or irregular systems being much lower, below 12%. LBGs have a median effective radius of 2.5 kpc and bigger galaxies have higher total SFR and stellar mass. Comparing to their high-redshift analogues, we find evidence that LBGs at lower redshifts are bigger, redder in the UV continuum, and have a major presence of older stellar populations in their SEDs. However, we do not find significant difference in the distributions of stellar mass or dust attenuation.

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.

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.

We perform a strong-lensing analysis of the merging galaxy cluster MACS J0416.1-2403 (M0416; z=0.42) in recent CLASH/HST observations. We identify 70 new multiple images and candidates of 23 background sources in the range 0.7<z_{phot}<6.14 including two probable high-redshift dropouts, revealing a highly elongated lens with axis ratio ~5:1, and a major axis of ~100\arcsec (z_{s}~2). Compared to other well-studied clusters, M0416 shows an enhanced lensing efficiency. Although the critical area is not particularly large (~0.6 \square\arcmin; z_{s}~2), the number of multiple images, per critical area, is anomalously high. We calculate that the observed elongation boosts the number of multiple images, \emph{per critical area}, by a factor of ~2.5\times, due to the increased ratio of the caustic area relative to the critical area. Additionally, we find that the observed separation between the two main mass components enlarges the critical area by a factor of ~2. These geometrical effects can account for the high number (density) of multiple images observed. We find in numerical simulations, that only ~4% of the clusters (with M_{vir}>6 x 10^{14} h^{-1}M_{\odot}) exhibit as elongated critical curves as M0416.

Small galaxies consisting entirely of population III (pop III) stars may form at high redshifts, and could constitute one of the best probes of such stars. Here, we explore the prospects of detecting gravitationally lensed pop III galaxies behind the galaxy cluster J0717.5+3745 (J0717) with both the Hubble Space Telescope (HST) and the upcoming James Webb Space Telescope (JWST). By projecting simulated catalogs of pop III galaxies at z~7-15 through the J0717 magnification maps, we estimate the lensed number counts as a function of flux detection threshold. We find that the ongoing HST survey CLASH, targeting a total of 25 galaxy clusters including J0717, potentially could detect a small number of pop III galaxies if ~1% of the baryons in these systems have been converted into pop III stars. Using JWST exposures of J0717, this limit can be pushed to ~0.1% of the baryons. Ultra-deep JWST observations of unlensed fields are predicted to do somewhat worse, but will be able to probe pop III galaxies with luminosities intermediate between those detectable in HST/CLASH and in JWST observations of J0717. We also explain how current measurements of the galaxy luminosity function at z=7-10 can be used to constrain pop III galaxy models with very high star formation efficiencies (~10% of the baryons converted into pop III stars).

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.

We report statistical results for dark matter (DM) velocity anisotropy, \beta, from a sample of some 6000 cluster-size halos (at redshift zero) identified in a \Lambda CDM hydrodynamical adaptive mesh refinement simulation performed with the Enzo code. These include profiles of \beta\ 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, \gamma, and \beta\ at large radii, $0.3 r_{vir} < r < r_{vir}$.

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, \Delta, are distributed equally around zero with no preferred direction, and are well described by a log-normal distribution with <log_{10}(\Delta [h^{-1} Mpc])>=-1.895^{+0.003}_{-0.004}, and \sigma=0.501\pm0.004 (95% confidence levels), or <log_{10}(\Delta [\arcsec])>=0.564\pm0.005, and \sigma=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]

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.

We examine the inner mass distribution of the relaxed galaxy cluster Abell 383 in deep 16-band HST/ACS+WFC3 imaging taken as part of the CLASH multi-cycle treasury program. Our program is designed to study the dark matter distribution in 25 massive clusters, and balances depth with a wide wavelength coverage to better identify lensed systems and generate precise photometric redshifts. This information together with the predictive strength of our strong-lensing analysis method identifies 13 new multiply-lensed images and candidates, so that a total of 27 multiple-images of 9 systems are used to tightly constrain the inner mass profile, $d\log \Sigma/d\log r\simeq -0.6\pm 0.1$ (r<160 kpc). We find consistency with the standard distance-redshift relation for the full range spanned by the lensed images, 1.01<z<6.03, with the higher redshift sources deflected through larger angles as expected. The inner mass profile derived here is consistent with the results of our independent weak-lensing analysis of wide-field Subaru images, with good agreement in the region of overlap. The overall mass profile is well fitted by an NFW profile with M_{vir}=(5.37^{+0.70}_{-0.63}\pm 0.26) x 10^{14}M_{\odot}/h and a relatively high concentration, c_{vir}=8.77^{+0.44}_{-0.42}\pm 0.23, which lies above the standard c-M relation similar to other well-studied clusters. The critical radius of Abell 383 is modest by the standards of other lensing clusters, r_{E}\simeq16\pm2\arcsec (for z_s=2.55), so the relatively large number of lensed images uncovered here with precise photometric redshifts validates our imaging strategy for the CLASH survey. In total we aim to provide similarly high-quality lensing data for 25 clusters, 20 of which are X-ray selected relaxed clusters, enabling a precise determination of the representative mass profile free from lensing bias. (ABRIDGED)

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{\sigma} 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{\rho}/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.

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\Sigma/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.

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.

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.

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 {\mu}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.

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(\theta_{e}\arcsec)>=0.73^{+0.02}_{-0.03}$, $\sigma=0.316^{+0.004}_{-0.002}$, and with higher abundance of large $\theta_{e}$ clusters than predicted by $\Lambda$CDM. We visually inspect each of the clusters with $\theta_{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 $\theta_{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.

We characterize the ability of the ALHAMBRA survey to assign accurate photo-z's to BLAGN and QSOs based on their ALHAMBRA very-low-resolution optical-NIR spectroscopy. A sample of 170 spectroscopically identified BLAGN and QSOs have been used together with a library of templates (including SEDs from AGN, normal, starburst galaxies and stars) in order to fit the 23 photometric data points provided by ALHAMBRA in the optical and NIR (20 medium-band optical filters plus the standard JHKs). We find that the ALHAMBRA photometry is able to provide an accurate photo-z and spectral classification for ~88% of the spectroscopic sources over 2.5 deg^2 in different areas of the survey, all of them brighter than m678=23.5 (equivalent to r(SLOAN)~24.0). The derived photo-z accuracy is better than 1% and comparable to the most recent results in other cosmological fields. The fraction of outliers (~12%) is mainly caused by the larger photometric errors for the faintest sources and the intrinsic variability of the BLAGN/QSO population. A small fraction of outliers may have an incorrectly assigned spectroscopic redshift. The definition of the ALHAMBRA survey in terms of the number of filters, filter properties, area coverage and depth is able to provide photometric redshifts for BLAGN/QSOs with a precision similar to any previous survey that makes use of medium-band optical photometry. In agreement with previous literature results, our analysis also reveals that, in the 0<z<4 redshift interval, very accurate photo-z can be obtained without the use of near-IR broadband photometry at the expense of a slight increase of outliers. The NIR importance is expected to increase at higher redshifts (z>4). These results are relevant for the design of future optical follow-ups of surveys with a large fraction of BLAGN, as it is the case for X-rays or radio surveys.

We report the discovery of a z_{phot}=6.18^{+0.05}_{-0.07} (95% confidence level) dwarf galaxy, lensed into four images by the galaxy cluster MACS J0329.6-0211 (z_{l}=0.45). The galaxy is observed as a high-redshift dropout in HST/ACS/WFC3 CLASH and Spitzer/IRAC imaging. Its redshift is securely determined due to a clear detection of the Lyman-break in the 18-band photometry, making this galaxy one of the highest-redshift multiply-lensed objects known to date with an observed magnitude of F125W=24.00\pm0.04 AB mag for its highest-magnified image. We also present the first strong-lensing analysis of this cluster uncovering 15 additional multiply-imaged candidates of five lower-redshift sources spanning the range z_{s}~2-4. The mass model independently supports the high photometric redshift and reveals magnifications of 11.6^{+8.9}_{-4.1}, 17.6^{+6.2}_{-3.9}, 3.9^{+3.0}_{-1.7}, and 3.7^{+1.3}_{-0.2}, respectively, for the four images of the high-redshift galaxy. With this we construct a source image with a physical resolution of ~200 pc when the universe was ~0.9 Gyr old, where the z~6.2 galaxy occupies a source-plane area of approximately 2.2 kpc^{2}. Modeling the observed spectral energy distribution using population synthesis models, we find a demagnified stellar mass of ~10^{9} {M}_{sun}, subsolar metallicity (Z/Z_{sun}~0.5), low dust content (A_{V}~0.1 mag), a demagnified star formation rate (SFR) of ~3.2 {M}_{sun} yr^{-1}, and a specific SFR of ~3.4 Gyr^{-1}, all consistent with the properties of local dwarf galaxies.

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 {\Lambda}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 {\Lambda}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%.

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.

We report the discovery of seven strongly lensed Lyman break galaxy (LBG) candidates at z~7 detected in Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) imaging of Abell 1703. The brightest candidate, called A1703-zD1, has an observed (lensed) magnitude of 24.0 AB (26 sigma) in the WFC3/IR F160W band, making it 0.2 magnitudes brighter than the z_850-dropout candidate recently reported behind the Bullet Cluster and 0.7 magnitudes brighter than the previously brightest known z~7.6 galaxy, A1689-zD1. With a cluster magnification of ~9, this source has an intrinsic magnitude of H_160 = 26.4 AB, a strong z_850 - J_125 break of 1.7 magnitudes, and a photometric redshift of z~6.7. Additionally, we find six other bright LBG candidates with H_160 band magnitudes of 24.9-26.4, photometric redshifts z~6.4 - 8.8, and magnifications mu~3-40. Stellar population fits to the ACS, WFC3/IR, and \Spitzer/IRAC data for A1703-zD1 and A1703-zD4 yield stellar masses (0.7 - 3.0) x 10^{9} M_sun, stellar ages 5-180 Myr, and star-formation rates ~7.8 M_sun/yr, and low reddening with A_V <= 0.7. The source-plane reconstruction of the exceptionally bright candidate A1703-zD1 exhibits an extended structure, spanning ~4 kpc in the z~6.7 source plane, and shows three resolved star-forming knots of radius r~0.4 kpc.

We present a strong-lensing analysis of the galaxy cluster MACS J1206.2-0847 ($z$=0.44) using UV, Optical, and IR, HST/ACS/WFC3 data taken as part of the CLASH multi-cycle treasury program, with VLT/VIMOS spectroscopy for some of the multiply-lensed arcs. The CLASH observations, combined with our mass-model, allow us to identify 47 new multiply-lensed images of 12 distant sources. These images, along with the previously known arc, span the redshift range $1\la z\la5.5$, and thus enable us to derive a detailed mass distribution and to accurately constrain, for the first time, the inner mass-profile of this cluster. We find an inner profile slope of $d\log \Sigma/d\log \theta\simeq -0.55\pm 0.1$ (in the range [1\arcsec, 53\arcsec], or $5\la r \la300$ kpc), as commonly found for relaxed and well-concentrated clusters. Using the many systems uncovered here we derive credible critical curves and Einstein radii for different source redshifts. For a source at $z_{s}\simeq2.5$, the critical curve encloses a large area with an effective Einstein radius of $\theta_{E}=28\pm3\arcsec$, and a projected mass of $1.34\pm0.15\times10^{14} M_{\odot}$. From the current understanding of structure formation in concordance cosmology, these values are relatively high for clusters at $z\sim0.5$, so that detailed studies of the inner mass distribution of clusters such as MACS J1206.2-0847 can provide stringent tests of the $\Lambda$CDM paradigm.

We present a detailed strong lensing, weak lensing and X-ray analysis of Abell 2744 (z = 0.308), one of the most actively merging galaxy clusters known. It appears to have unleashed `dark', `ghost', `bullet' and `stripped' substructures, each ~10^14 solar masses. The phenomenology is complex and will present a challenge for numerical simulations to reproduce. With new, multiband HST imaging, we identify 34 strongly-lensed images of 11 galaxies around the massive Southern `core'. Combining this with weak lensing data from HST, VLT and Subaru, we produce the most detailed mass map of this cluster to date. We also perform an independent analysis of archival Chandra X-ray imaging. Our analyses support a recent claim that the Southern core and Northwestern substructure are post-merger and exhibit morphology similar to the Bullet Cluster viewed from an angle. From the separation between X-ray emitting gas and lensing mass in the Southern core, we derive a new and independent constraint on the self-interaction cross section of dark matter particles sigma/m <~ 3 \pm 1 cm^2 g^-1. In the Northwestern substructure, the gas, dark matter, and galaxy components have become separated by much larger distances. Most curiously, the `ghost' clump (primarily gas) leads the `dark' clump (primarily dark matter) by more than 150 kpc. We propose an enhanced `ram-pressure slingshot' scenario which may have yielded this reversal of components with such a large separation, but needs further confirmation by follow-up observations and numerical simulations. A secondary merger involves a second `bullet' clump in the North and an extremely `stripped' clump to the West. The latter appears to exhibit the largest separation between dark matter and X-ray emitting baryons detected to date in our sky.

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{\rho}/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.

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.

A long-standing problem of strong lensing by galaxy clusters regards the observed high rate of giant gravitational arcs as compared to the predictions in the framework of the "standard" cosmological model. Recently, few other inconsistencies between theoretical expectations and observations have been claimed which regard the large size of the Einstein rings and the high concentrations of few clusters with strong lensing features. All of these problems consistently indicate that observed galaxy clusters may be gravitational lenses stronger than expected. We use clusters extracted from the MareNostrum Universe to build up mock catalogs of galaxy clusters selected through their X-ray flux. We use these objects to estimate the probability distributions of lensing cross sections, Einstein rings, and concentrations for the sample of 12 MACS clusters at $z>0.5$ presented in Ebeling et al. (2007) and discussed in Zitrin et al. (2010). We find that simulated clusters produce $\sim 50%$ less arcs than observed clusters do. The medians of the distributions of the Einstein ring sizes differ by $\sim 25%$ between simulations and observations. We estimate that, due to cluster triaxiality and orientation biases affecting the lenses with the largest cross sections, the concentrations of the individual MACS clusters inferred from the lensing analysis should be up to a factor of $\sim 2$ larger than expected from the $\Lambda$CDM model. The arc statistics, the Einstein ring, and the concentration problems in strong lensing clusters are mitigated but not solved on the basis of our analysis. Nevertheless, due to the lack of redshifts for most of the multiple image systems used for modeling the MACS clusters, the results of this work will need to be verified with additional data. The upcoming CLASH program will provide an ideal sample for extending our comparison (abridged).

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.

Accurate mass determination of clusters of galaxies is crucial if they are to be used as cosmological probes. However, there are some discrepancies between cluster masses determined based on gravitational lensing, and X-ray observations assuming strict hydrostatic equilibirium (i.e., the equilibrium gas pressure is provided entirely by thermal pressure). Cosmological simulations suggest that turbulent gas motions remaining from hierarchical structure formation may provide a significant contribution to the equilibrium pressure in clusters. We analyze a sample of massive clusters of galaxies drawn from high resolution cosmological simulations, and find a significant contribution (20%-45%) from non-thermal pressure near the center of relaxed clusters, and, in accord with previous studies, a minimum contribution at about 0.1 Rvir, growing to about 30%-45% at the virial radius, Rvir. Our results strongly suggest that relaxed clusters should have significant non-thermal support in their core region. As an example, we test the validity of strict hydrostatic equilibirium in the well-studied massive galaxy cluster Abell 1689 using the latest high resolution gravitational lensing and X-ray observations. We find a contribution of about 40% from non-thermal pressure within the core region of A1689, suggesting an alternate explanation for the mass discrepancy: the strict hydrostatic equilibirium is not valid in this region.

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 \Sigma/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.

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 ($\mu > \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 $\Lambda$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 $\Lambda$CDM predictions.

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 $\Lambda$CDM model, similar to other massive relaxed clusters with accurately determined lensing-based profiles.

Ly-alpha emitters (LAEs) are seen everywhere in the redshift domain from local to z~7. Far-infrared (FIR) counterparts of LAEs at different epochs could provide direct clues on dust content, extinction, and spectral energy distribution (SED) for these galaxies. We search for FIR counterparts of LAEs that are optically detected in the GOODS-North field at redshift z~2.2 using data from the Herschel Space Telescope with the Photodetector Array Camera and Spectrometer (PACS). The LAE candidates were isolated via color-magnitude diagram using the medium-band photometry from the ALHAMBRA Survey, ancillary data on GOODS-North, and stellar population models. According to the fitting of these spectral synthesis models and FIR/optical diagnostics, most of them seem to be obscured galaxies whose spectra are AGN-dominated. From the analysis of the optical data, we have observed a fraction of AGN or composite over source total number of ~0.75 in the LAE population at z~2.2, which is marginally consistent with the fraction previously observed at z=2.25 and even at low redshift (0.2<z<0.45), but significantly different from the one observed at redshift ~3, which could be compatible either with a scenario of rapid change in the AGN fraction between the epochs involved or with a non detection of obscured AGN in other z=2-3 LAE samples due to lack of deep FIR observations. We found three robust FIR (PACS) counterparts at z~2.2 in GOODS-North. This demonstrates the possibility of finding dust emission in LAEs even at higher redshifts.

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

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.

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.

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.

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.

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, \alpha ~= -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.

This paper presents the characterization of the optical range of the ALHAMBRA photometric system, a 20 contiguous, equal-width, medium-band CCD system with wavelength coverage from 3500A to 9700A. The photometric description of the system is done by presenting the full response curve as a product of the filters, CCD and atmospheric transmission curves, and using some first and second order moments of this response function. We also introduce the set of standard stars that defines the system, formed by 31 classic spectrophotometric standard stars which have been used in the calibration of other known photometric systems, and 288 stars, flux calibrated homogeneously, from the Next Generation Spectral Library (NGSL). Based on the NGSL, we determine the transformation equations between Sloan Digital Sky Survey (SDSS) ugriz photometry and the ALHAMBRA photometric system, in order to establish some relations between both systems. Finally we develop and discuss a strategy to calculate the photometric zero points of the different pointings in the ALHAMBRA project.

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 $\theta_{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 $\theta_{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 $\Lambda CDM$ (WMAP5 parameters with $2\sigma$ uncertainties).

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 $\Lambda$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.

We identify the largest known lensed images of a single spiral galaxy, lying close to the centre of the distant cluster MACS J1149.5+2223 ($z=0.544$). These images cover a total area of $\simeq 150 \Box\arcsec$ and are magnified $\simeq 200$ times. Unusually, there is very little image distortion implying the central mass distribution is almost uniform over a wide area ($r\simeq200 kpc$) with a surface density equal to the critical density for lensing, corresponding to maximal lens magnification. Many fainter multiply-lensed galaxies are also uncovered by our model, outlining a very large tangential critical curve, of radius $r\simeq 170 kpc$, posing a potential challenge for the standard LCDM-Cosmology. Because of the uniform central mass distribution a particularly clean measurement of the mass of the brightest cluster galaxy is possible here, for which we infer stars contribute most of the mass within a limiting radius of $\simeq 30 kpc$, with a mass-to-light ratio of $M/L_{B}\simeq 4.5(M/L)_{\odot}$. This cluster with its uniform and central mass distribution acts analogously to a regular magnifying glass, converging light without distorting the images, resulting in the most powerful lens yet discovered for accessing the faint high-$z$ Universe.

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$\sigma$). 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.

We report on the discovery of three bright, strongly-lensed objects behind Abell 1703 and CL0024+16 from a dropout search over 25 square arcminutes of deep NICMOS data, with deep ACS optical coverage. They are undetected in the deep ACS images below 8500 A and have clear detections in the J and H bands. Fits to the ACS, NICMOS and IRAC data yield robust photometric redshifts in the range z~6-7 and largely rule out the possibility that they are low-redshift interlopers. All three objects are extended, and resolved into a pair of bright knots. The bright i-band dropout in Abell 1703 has an H-band AB magnitude of 23.9, which makes it one of the brightest known galaxy candidates at z>5.5. Our model fits suggest a young, massive galaxy only ~ 60 million years old with a mass of ~ 1E10 solar mass. The dropout galaxy candidates behind CL0024+16 are separated by 2.5" (~ 2 kpc in the source plane), and have H-band AB magnitudes of 25.0 and 25.6. Lensing models of CL0024+16 suggest that the objects have comparable intrinsic magnitudes of AB ~ 27.3, approximately one magnitude fainter than L* at z~6.5. Their similar redshifts, spectral energy distribution, and luminosities, coupled with their very close proximity on the sky, suggest that they are spatially associated, and plausibly are physically bound. Combining this sample with two previously reported, similarly magnified galaxy candidates at z~6-8, we find that complex systems with dual nuclei may be a common feature of high-redshift galaxies.

The ALHAMBRA survey aims to cover 4 square degrees using a system of 20 contiguous, equal width, medium-band filters spanning the range 3500 A to 9700 A plus the standard JHKs filters. Here we analyze deep near-IR number counts of one of our fields (ALH08) for which we have a relatively large area (0.5 square degrees) and faint photometry (J=22.4, H=21.3 and K=20.0 at the 50% of recovery efficiency for point-like sources). We find that the logarithmic gradient of the galaxy counts undergoes a distinct change to a flatter slope in each band: from 0.44 at [17.0, 18.5] to 0.34 at [19.5, 22.0] for the J band; for the H band 0.46 at [15.5, 18.0] to 0.36 at [19.0, 21.0], and in Ks the change is from 0.53 in the range [15.0, 17.0] to 0.33 in the interval [18.0, 20.0]. These observations together with faint optical counts are used to constrain models that include density and luminosity evolution of the local type-dependent luminosity functions. Our models imply a decline in the space density of evolved early-type galaxies with increasing redshift, such that only 30% - 50% of the bulk of the present day red-ellipticals was already in place at z~1.

We conduct a search for z>~7 dropout galaxies behind 11 massive lensing clusters using 21 arcmin**2 of deep HST NICMOS, ACS, and WFPC2 image data. In total, over this entire area, we find only one robust z~7 z-dropout candidate (previously reported around Abell 1689). Four less robust z-dropout and J-dropout candidates are also found. The nature of the four weaker candidates could not be precisely determined due to the limited depth of the available optical data, but detailed simulations suggest that all four could be low-redshift interlopers. We compare these numbers with what we might expect using the z~7 UV luminosity function (LF) determined from field searches. We predict 2.7 z~7 z-dropouts and 0.3 z~9 J-dropouts over our cluster search area, in reasonable agreement with our observational results, given the small numbers. The number of z>~7 candidates we find in the present search are much lower than has been reported in several previous studies of the prevalence of z>~7 galaxies behind lensing clusters. To understand these differences, we examined z>~7 candidates in other studies and conclude that only a small fraction are likely to be z>~7 galaxies. Our findings support models that show that gravitational lensing from clusters is of the most value for detecting galaxies at magnitudes brighter than L* (H<~ 27) where the LF is expected to be very steep. Use of these clusters to constrain the faint-end slope or determine the full LF is likely of less value due to the shallower effective slope measured for the LF at fainter magnitudes, as well as significant uncertainties introduced from modelling both the gravitational lensing and incompleteness.

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 = \Omega_b/\Omega_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.

In the next years, several cosmological surveys will rely on imaging data to estimate the redshift of galaxies, using traditional filter systems with 4-5 optical broad bands; narrower filters improve the spectral resolution, but strongly reduce the total system throughput. We explore how photometric redshift performance depends on the number of filters n_f, characterizing the survey depth through the fraction of galaxies with unambiguous redshift estimates. For a combination of total exposure time and telescope imaging area of 270 hrs m^2, 4-5 filter systems perform significantly worse, both in completeness depth and precision, than systems with n_f >= 8 filters. Our results suggest that for low n_f, the color-redshift degeneracies overwhelm the improvements in photometric depth, and that even at higher n_f, the effective photometric redshift depth decreases much more slowly with filter width than naively expected from the reduction in S/N. Adding near-IR observations improves the performance of low n_f systems, but still the system which maximizes the photometric redshift completeness is formed by 9 filters with logarithmically increasing bandwidth (constant resolution) and half-band overlap, reaching ~0.7 mag deeper, with 10% better redshift precision, than 4-5 filter systems. A system with 20 constant-width, non-overlapping filters reaches only ~0.1 mag shallower than 4-5 filter systems, but has a precision almost 3 times better, dz = 0.014(1+z) vs. dz = 0.042(1+z). We briefly discuss a practical implementation of such a photometric system: the ALHAMBRA survey.

The Einstein radius of a cluster provides a relatively model-independent measure of the mass density of a cluster within a projected radius of ~ 150 kpc, large enough to be relatively unaffected by gas physics. We show that the observed Einstein radii of four well-studied massive clusters, for which reliable virial masses are measured, lie well beyond the predicted distribution of Einstein radii in the standard LambdaCDM model. Based on large samples of numerically simulated cluster-sized objects with virial masses ~ 1e15 solar, the predicted Einstein radii are only 15-25'', a factor of two below the observed Einstein radii of these four clusters. This is because the predicted mass profile is too shallow to exceed the critical surface density for lensing at a sizable projected radius. After carefully accounting for measurement errors as well as the biases inherent in the selection of clusters and the projection of mass measured by lensing, we find that the theoretical predictions are excluded at a 4-sigma significance. Since most of the free parameters of the LambdaCDM model now rest on firm empirical ground, this discrepancy may point to an additional mechanism that promotes the collapse of clusters at an earlier time thereby enhancing their central mass density.

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.

We present the results of deep spectroscopy for the central region of the cluster lens SDSS J1004+4112 with the Subaru telescope. A secure detection of an emission line of the faint blue stellar object (component E) near the center of the brightest cluster galaxy (G1) confirms that it is the central fifth image of the lensed quasar system. In addition, we measure the stellar velocity dispersion of G1 to be sigma_* = 352+-13 km/s. We combine these results to obtain constraints on the mass M_BH of the putative black hole (BH) at the center of the inactive galaxy G1, and hence on the M_BH-sigma_* relation at the lens redshift z_l=0.68. From detailed mass modeling, we place an upper limit on the black hole mass, M_BH < 2.1x10^{10}M_sun at 1-sigma level (<3.1x10^{10}M_sun at 3-sigma), which is consistent with black hole masses expected from the local and redshift-evolved M_BH-sigma_* relations, M_BH~10^{9}-10^{10}M_sun.

We present a high signal-to-noise spectrum of a bright galaxy at z = 4.9 in 14 h of integration on VLT FORS2. This galaxy is extremely bright, i_850 = 23.10 +/- 0.01, and is strongly-lensed by the foreground massive galaxy cluster Abell 1689 (z=0.18). Stellar continuum is seen longward of the Ly-alpha emission line at ~7100 \AA, while intergalactic H I produces strong absorption shortward of Ly-alpha. Two transmission spikes at ~6800 Angstroms (A) and ~7040 A are also visible, along with other structures at shorter wavelengths. Although fainter than a QSO, the absence of a strong central ultraviolet flux source in this star forming galaxy enables a measurement of the H I flux transmission in the intergalactic medium (IGM) in the vicinity of a high redshift object. We find that the effective H I optical depth of the IGM is remarkably high within a large 14 Mpc (physical) region surrounding the galaxy compared to that seen towards QSOs at similar redshifts. Evidently, this high-redshift galaxy is located in a region of space where the amount of H I is much larger than that seen at similar epochs in the diffuse IGM. We argue that observations of high-redshift galaxies like this one provide unique insights on the nascent stages of baryonic large-scale structures that evolve into the filamentary cosmic web of galaxies and clusters of galaxies observed in the present universe.

(ABRIDGED) We describe the first results of the ALHAMBRA survey which provides cosmic tomography of the evolution of the contents of the Universe over most of Cosmic history. Our approach employs 20 contiguous, equal-width, medium-band filters covering from 3500 to 9700 A, plus the JHKs bands, to observe an area of 4 sqdeg on the sky. The optical photometric system has been designed to maximize the number of objects with accurate classification by SED and redshift, and to be sensitive to relatively faint emission lines. The observations are being carried out with the Calar Alto 3.5m telescope using the cameras LAICA and O-2000. The first data confirm that we are reaching the expected magnitude limits of AB<~25 mag in the optical filters from the blue to 8300 A, and from AB=24.7 to 23.4 for the redder ones. The limit in the NIR is (Vega) K_s~20, H~21, J~22. We expect to obtain accurate redshift values, Delta z/(1+z) <~ 0.03 for about 5x10^5 galaxies with I<~25 (60% complete), and z_med=0.74. This accuracy, together with the homogeneity of the selection function, will allow for the study of the redshift evolution of the large scale structure, the galaxy population and its evolution with redshift, the identification of clusters of galaxies, and many other studies, without the need for any further follow-up. It will also provide targets for detailed studies with 10m-class telescopes. Given its area, spectral coverage and its depth, apart from those main goals, the ALHAMBRA-Survey will also produce valuable data for galactic studies.

We present the first results from the ALHAMBRA survey. ALHAMBRA will cover a relatively wide area (4 square degrees) using a purposely-designed set of 20 medium-band filters, down to an homogeneous magnitude limit AB~25 in most of them, adding also deep near-infrared imaging in JHK. To this aim we are using the Calar Alto 3.5m telescope. A small area of the ALHAMBRA survey has already been observed through our complete filter set, and this allows for the first time to check all the steps of the survey, including the pipelines that have been designed for the project, the fulfilment of the data quality expectations, the calibration procedures, and the photometric redshift machinery for which ALHAMBRA has been optimised. We present here the basic results regarding the properties of the galaxy sample selected in a 15x15 square arcmin area of the ALHAMBRA-8 field, which includes approximately 10000 galaxies with precise photometric redshift measurements. In a first estimate, approximately 500 of them must be galaxies with z>2.

Using HST and Spitzer IRAC imaging, we report the discovery of a very bright strongly lensed Lyman break galaxy (LBG) candidate at z~7.6 in the field of the massive galaxy cluster Abell 1689. The galaxy candidate, which we refer to as A1689-zD1, shows a strong z-J break of at least 2.2 mag and is completely undetected (<1 sigma) in HST/ACS g, r, i, and z-band data. These properties, combined with the very blue J-H and H-[4.5] colors, are exactly the properties of an z~7.6 LBG and can only be reasonably fit by a star-forming galaxy at z=7.6 +/- 0.4. Attempts to reproduce these properties with a model galaxy at z<4 yield particularly poor fits. A1689-zD1 has an observed (lensed) magnitude of 24.7 AB (8 sigma) in the NICMOS H band and is ~1.3 mag brighter than the brightest-known z-dropout galaxy. When corrected for the cluster magnification of 9.3 at z~7.6, the candidate has an intrinsic magnitude of H=27.1 AB, or about an L* galaxy at z~7.6. The source-plane deprojection shows that the star formation is occurring in compact knots of size ~<300 pc. The best-fit stellar population synthesis models yield a median redshift of 7.6, stellar masses (1.6-3.9) x 10^9 M_sun, stellar ages 45-320 Myr, star-formation rates ~<7.6 M_sun/yr, and low reddening with A_V <= 0.3. These properties are generally similar to those of LBGs found at z~5-6. The inferred stellar ages suggest a formation redshift of z~8-10 (t~<0.63 Gyr). A1689-zD1 is the brightest observed, highly reliable z>7.0 galaxy candidate found to date.

[Abridged] We derive a projected 2D mass map of the well studied galaxy cluster A1689 based on an entropy-regularized maximum-likelihood combination of the lens magnification and distortion of red background galaxies registered in deep Subaru images. The known strong lensing information is also readily incorporated in this approach, represented as a central pixel with a mean surface density close to the critical value. We also utilize the distortion measurements to locally downweight the intrinsic clustering noise, which otherwise perturbs the depletion signal. The projected mass profile continuously steepens with radius and is well fitted by the NFW model, but with a surprising large concentration c_vir=13.4^{+5.3}_{-3.3}, lying far from the predicted value of c_vir ~ 5, corresponding to the measured virial mass, M_vir=(2.1\pm 0.2)\times 10^15M_{\odot}, posing a challenge to the standard assumptions defining the LCDM model. We examine the consistency of our results with estimates derived with the standard weak lensing estimators and by comparison with the inner mass profile obtained from strong lensing. All the reconstructions tested here imply a virial mass in the range, M_vir=(1.5-2.1)\times 10^15M_{\odot}, and the combined ACS and Subaru-2D mass reconstruction yields a tight constraint on the concentration parameter, c_vir=12.7 \pm 1(stat.) \pm 2.8(systematic) (c_200 ~ 10), improving upon the statistical accuracy of our earlier 1D analysis. Importantly, our best fitting profile properly reproduces the observed Einstein radius of 45 arcsec (z_s=1), in contrast to other weak lensing work, reporting lower concentration profiles, which underestimate the observed Einstein radius.

We present the first result of a survey for strong galaxy-galaxy lenses in Sloan Digital Sky Survey (SDSS) images. SDSS J082728.70+223256.4 was selected as a lensing candidate using selection criteria based on the color and positions of objects in the SDSS photometric catalog. Follow-up imaging and spectroscopy showed this object to be a lensing system. The lensing galaxy is an elliptical at z = 0.349 in a galaxy cluster. The lensed galaxy has the spectrum of a post-starburst galaxy at z = 0.766. The lensing galaxy has an estimated mass of $\sim 1.2 \times 10^{12} M_{\odot}$ and the corresponding mass to light ratio in the B-band is $\sim 26 M_{\odot}/L_{\odot}$ inside 1.1 effective radii of the lensing galaxy. Our study shows how catalogs drawn from multi-band surveys can be used to find strong galaxy-galaxy lenses having multiple lens images. Our strong lensing candidate selection based on photometry-only catalogs will be useful in future multi-band imaging surveys such as SNAP and LSST.

We present a new approach to gravitational lens massmap reconstruction. Our massmap solutions perfectly reproduce the positions, fluxes, and shears of all multiple images. And each massmap accurately recovers the underlying mass distribution to a resolution limited by the number of multiple images detected. We demonstrate our technique given a mock galaxy cluster similar to Abell 1689 which gravitationally lenses 19 mock background galaxies to produce 93 multiple images. We also explore cases in which far fewer multiple images are observed, such as four multiple images of a single galaxy. Massmap solutions are never unique, and our method makes it possible to explore an extremely flexible range of physical (and unphysical) solutions, all of which perfectly reproduce the data given. Each reconfiguration of the source galaxies produces a new massmap solution. An optimization routine is provided to find those source positions (and redshifts, within uncertainties) which produce the "most physical" massmap solution, according to a new figure of merit developed here. Our method imposes no assumptions about the slope of the radial profile nor mass following light. But unlike "non-parametric" grid-based methods, the number of free parameters we solve for is only as many as the number of observable constraints (or slightly greater if fluxes are constrained). For each set of source positions and redshifts, massmap solutions are obtained "instantly" via direct matrix inversion by smoothly interpolating the deflection field using a recently developed mathematical technique. Our LensPerfect software is straightforward and easy to use and is made publicly available via our website.

We identify a third image in the unique quasar lens SDSS J1029+2623, the second known quasar lens produced by a massive cluster of galaxies. The spectrum of the third image shows similar emission and absorption features, but has a redder continuum than the other two images which can be explained by differential extinction or microlensing. We also identify several lensed arcs. Our observations suggest a complicated structure of the lens cluster at z~0.6. We argue that the three lensed images are produced by a naked cusp on the basis of successful mass models, the distribution of cluster member galaxies, and the shapes and locations of the lensed arcs. Lensing by a naked cusp is quite rare among galaxy-scale lenses but is predicted to be common among large-separation lensed quasars. Thus the discovery can be viewed as support for an important theoretical prediction of the standard cold dark matter model.

We carry out a comprehensive joint analysis of high quality HST/ACS and Chandra measurements of A1689, from which we derive mass, temperature, X-ray emission and abundance profiles. The X-ray emission is smooth and symmetric, and the lensing mass is centrally concentrated indicating a relaxed cluster. Assuming hydrostatic equilibrium we deduce a 3D mass profile that agrees simultaneously with both the lensing and X-ray measurements. However, the projected temperature profile predicted with this 3D mass profile exceeds the observed temperature by ~30% at all radii, a level of discrepancy comparable to the level found for other relaxed clusters. This result may support recent suggestions from hydrodynamical simulations that denser, more X-ray luminous small-scale structure can bias observed temperature measurements downward at about the same (~30%) level. We determine the gas entropy at 0.1r_{vir} (where r_{vir} is the virial radius) to be ~800 keV cm^2, as expected for a high temperature cluster, but its profile at >0.1r_{vir} has a power-law form with index ~0.8, considerably shallower than the ~1.1 index advocated by theoretical studies and simulations. Moreover, if a constant entropy ''floor'' exists at all, then it is within a small region in the inner core, r<0.02r_{vir}, in accord with previous theoretical studies of massive clusters.

We present results from a weak lensing analysis of the galaxy cluster A1689 (z=0.183) based on deep wide-field imaging data taken with Suprime-Cam on Subaru telescope. A maximum entropy method has been used to reconstruct directly the projected mass distribution of A1689 from combined lensing distortion and magnification measurements of red background galaxies.The resulting mass distribution is clearly concentrated around the cD galaxy, and mass and light in the cluster are similarly distributed in terms of shape and orientation. The azimuthally-averaged mass profile from the two-dimensional reconstruction is in good agreement with the earlier results from the Subaru one-dimensional analysis of the weak lensing data, supporting the assumption of quasi-circular symmetry in the projected mass distribution of the cluster.

We present a comprehensive mass reconstruction of the rich galaxy cluster Cl 0024+17 at z~0.4 from ACS data, unifying both strong- and weak-lensing constraints. The weak-lensing signal from a dense distribution of background galaxies (~120 per square arcmin) across the cluster enables the derivation of a high-resolution parameter-free mass map. The strongly-lensed objects tightly constrain the mass structure of the cluster inner region on an absolute scale, breaking the mass-sheet degeneracy. The mass reconstruction of Cl 0024+17 obtained in such a way is remarkable. It reveals a ringlike dark matter substructure at r~75" surrounding a soft, dense core at r~50". We interpret this peculiar sub-structure as the result of a high-speed line-of-sight collision of two massive clusters 1-2 Gyr ago. Such an event is also indicated by the cluster velocity distribution. Our numerical simulation with purely collisionless particles demonstrates that such density ripples can arise by radially expanding, decelerating particles that originally comprised the pre-collision cores. Cl 0024+17 can be likened to the bullet cluster 1E0657-56, but viewed $along$ the collision axis at a much later epoch. In addition, we show that the long-standing mass discrepancy for Cl 0024+17 between X-ray and lensing can be resolved by treating the cluster X-ray emission as coming from a superposition of two X-ray systems. The cluster's unusual X-ray surface brightness profile that requires a two isothermal sphere description supports this hypothesis.

We present results of the HST Advanced Camera for Surveys spectroscopic ground-based redshift survey in the field of A1689. We measure 98 redshifts, increasing the number of spectroscopically confirmed objects by sixfold. We present two spectra from this catalog of the Sextet Arcs, images which arise from a strongly-lensed Lyman Break Galaxy (LBG) at a redshift of z=3.038. Gravitational lensing by the cluster magnifies its flux by a factor of ~16 and produces six separate images with a total r-band magnitude of r_625=21.7. The two spectra, each of which represents emission from different regions of the LBG, show H I and interstellar metal absorption lines at the systemic redshift. Significant variations are seen in Ly-alpha profile across a single galaxy, ranging from strong absorption to a combination of emission plus absorption. A spectrum of a third image close to the brightest arc shows Ly-alpha emission at the same redshift as the LBG, arising from either another spatially distinct region of the galaxy, or from a companion galaxy close to the LBG. Taken as a group, the Ly-alpha equivalent width in these three spectra decreases with increasing equivalent width of the strongest interstellar absorption lines. We discuss how these variations can be used to understand the physical conditions in the LBG. Intrinsically, this LBG is faint, ~0.1L*, and forming stars at a modest rate, ~4 solar masses per year. We also detect absorption line systems toward the Sextet Arcs at z=2.873 and z=2.534. The latter system is seen across two of our spectra.

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.

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)

(Abridged.) We present F435W (B), F606W (Broad V), and F814W (Broad I) coronagraphic images of the debris disk around Beta Pictoris obtained with HST's Advanced Camera for Surveys. We confirm that the previously reported warp in the inner disk is a distinct secondary disk inclined by ~5 deg from the main disk. The main disk's northeast extension is linear from 80 to 250 AU, but the southwest extension is distinctly bowed with an amplitude of ~1 AU over the same region. Both extensions of the secondary disk appear linear, but not collinear, from 80 to 150 AU. Within ~120 AU of the star, the main disk is ~50% thinner than previously reported. The surface-brightness profiles along the spine of the main disk are fitted with four distinct radial power laws between 40 and 250 AU, while those of the secondary disk between 80 and 150 AU are fitted with single power laws. These discrepancies suggest that the two disks have different grain compositions or size distributions. The F606W/F435W and F814W/F435W flux ratios of the composite disk are nonuniform and asymmetric about both projected axes of the disk. Within ~120 AU, the m_F435W-m_F606W and m_F435W-m_F814W colors along the spine of the main disk are ~10% and ~20% redder, respectively, than those of Beta Pic. These colors increasingly redden beyond ~120 AU, becoming 25% and 40% redder, respectively, than the star at 250 AU. We compare the observed red colors within ~120 AU with the simulated colors of non-icy grains having a radial number density ~r^-3 and different compositions, porosities, and minimum grain sizes. The observed colors are consistent with those of compact or moderately porous grains of astronomical silicate and/or graphite with sizes >0.15-0.20 um, but the colors are inconsistent with the blue colors expected from grains with porosities >90%.

We have combined deep BRIz' imaging over 2x940 arcmin^2 fields obtained with the Suprime-Cam on the Subaru telescope with JKs imaging with the SOFI camera at the New Technology Telescope to search for high-redshift massive galaxies. K-band selected galaxies have been identified over an area of ~920 arcmin^2 to K_Vega=19.2, of which 320 arcmin^2 are complete to K_Vega=20. The BzK selection technique was used to obtain complete samples of ~500 candidate massive star-forming galaxies (sBzKs) and ~160 candidate massive, passively-evolving galaxies (pBzKs), both at 1.4<z<2.5. With the (R-K)_Vega > 5 criterion we also identified ~850 extremely red objects (EROs). The surface density of sBzKs and pBzKs is found to 1.20+/-0.05 arcmin^{-2} and 0.38+/-0.03 arcmin^{-2}, respectively. Both sBzKs and pBzKs are strongly clustered, at a level at least comparable to that of EROs, with pBzKs appearing more clustered than sBzKs. We estimate the reddening, star formation rates (SFRs) and stellar masses (M_*) of the sBzKs, confirming that to K_Vega~20 median values are M_*~10^{11}M_sun, SFR 190M_sun yr^{-1}, and E(B-V)~0.44. The most massive sBzKs are also the most actively star-forming, an effect which can be seen as a manifestation of downsizing at early epochs. The space density of massive pBzKs at z~1.4-2 is 20%+/-7% that of similarly massive early-type galaxies at z~0, and similar to that of sBzKs of the same mass. We argue that star formation quenching in these sBzKs will result in nearly doubling the space density of massive early-type galaxies, thus matching their local density.

A five square arcminute region around the luminous radio-loud quasar SDSS J0836+0054 (z=5.8) hosts a wealth of associated galaxies, characterized by very red (1.3 < i_775 - z_{850} < 2.0) color. The surface density of these z~5.8 candidates is approximately six times higher than the number expected from deep ACS fields. This is one of the highest galaxy overdensities at high redshifts, which may develop into a group or cluster. We also find evidence for a substructure associated with one of the candidates. It has two very faint companion objects within two arcseconds, which are likely to merge. The finding supports the results of a recent simulation that luminous quasars at high redshifts lie on the most prominent dark-matter filaments and are surrounded by many fainter galaxies. The quasar activity from these regions may signal the buildup of a massive system.

We have identified three multiply imaged galaxies in Hubble Space Telescope images of the redshift z=0.68 cluster responsible for the large-separation quadruply lensed quasar, SDSS J1004+4112. Spectroscopic redshifts have been secured for two of these systems using the Keck I 10m telescope. The most distant lensed galaxy, at z=3.332, forms at least four images, and an Einstein ring encompassing 3.1 times more area than the Einstein ring of the lensed QSO images at z=1.74, due to the greater source distance. For a second multiply imaged galaxy, we identify Ly_alpha emission at a redshift of z=2.74. The cluster mass profile can be constrained from near the center of the brightest cluster galaxy, where we observe both a radial arc and the fifth image of the lensed quasar, to the Einstein radius of the highest redshift galaxy, ~110 kpc. Our preliminary modeling indicates that the mass approximates an elliptical body, with an average projected logarithmic gradient of ~-0.5. The system is potentially useful for a direct measurement of world models in a previously untested redshift range.

The steep mass profile of A1689 derived from recent detailed lensing observations is not readily reconciled with the low concentration halos predicted by the standard CDM model. However, halo triaxiality may act to bias the profile constraints derived assuming a spherically symmetric mass distribution, since lensing relates only to the projected mass distribution. A degree of halo triaxiality is inherent to the CDM structure formation, arising from the collision-less nature of the dark matter. Here we compare the CDM-based model predictions of triaxial halo with the precise lensing measurements of A1689 based on the ACS/HST and Subaru data, over a wide range of 10kpc<r<2Mpc. The model lensing profiles cover the intrinsic spread of halo mass and shape, and are projected over all inclinations when comparing with the data. We show that the model parameters are only weakly constrained and strongly degenerate mainly because of the lack of information along the line of sight. In particular, the limits on the concentration parameter become less restrictive with increasing triaxiality. Yet, by comparing the obtained constraints with expected probability distributions for the axis ratio and concentration parameters computed from numerical simulations, we find that ~6% of cluster-size halos in the CDM model can match the A1689 lensing observations at the 2-sigma level, corresponding to cases where the major-axis of the halo is aligned with the line of sight. Thus halo triaxiality could reduce the apparent discrepancy between theory and observation. This hypothesis needs to be further explored by a statistical lensing study for other clusters as well as by complementary three-dimensional information derived using X-ray, kinematics, and SZ effect observations.

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.

The K20 survey is a near infrared-selected, deep (Ks < 20) redshift survey targeting galaxies in two independent regions of the sky, the CDFS and the q0055-2659 field. The total Ks-selected sample includes 545 objects. Optical spectra for 525 of them have been obtained, providing 501 spectroscopic identifications (including 12 type-1 AGN and 45 stars). This paper describes the final K20 spectroscopic catalogue, along with the technique used to determine redshifts, measure the spectral features and characterize the spectra. The classification of the galaxy spectra has been performed according to a simple parametric recipe that uses the equivalent widths of the two main emission lines (O[II] and Halpha+[N II]) and two continuum indices (the 4000A break index, and a near-UV color index, C(28-39)). We defined three main spectroscopic classes: red early-type galaxies, blue emission-line galaxies and the intermediate galaxies, which show emission lines but a red continuum. More than 95% of the examined galaxies is included in one of these spectral types and a composite spectrum is built for each of the three galaxy classes. The blue emission-line and the early-type galaxies have been divided in redshift bins, and the corresponding composite spectra have been built, in order to investigate the evolution of the spectral properties of the K20 galaxies with redshift.

We study the dynamical origin of the structures observed in the scattered-light images of the resolved debris disk around HD 141569A. We explore the roles of radiation pressure from the central star, gas drag from the gas disk, and the tidal forces from two nearby stars in creating and maintaining these structures. We use a simple one-dimensional axisymmetric model to show that the presence of the gas helps confine the dust and that a broad ring of dust is produced if a central hole exists in the disk. This model also suggests that the disk is in a transient, excited dynamical state, as the observed dust creation rate applied over the age of the star is inconsistent with submillimeter mass measurements. We model in two dimensions the effects of a fly-by encounter between the disk and a binary star in a prograde, parabolic, coplanar orbit. We track the spatial distribution of the disk's gas, planetesimals, and dust. We conclude that the surface density distribution reflects the planetesimal distribution for a wide range of parameters. Our most viable model features a disk of initial radius 400 AU, a gas mass of 50 M_earth, and beta = 4 and suggests that the system is being observed within 4000 yr of the fly-by periastron. The model reproduces some features of HD 141569A's disk, such as a broad single ring and large spiral arms, but it does not reproduce the observed multiple spiral rings or disk asymmetries nor the observed clearing in the inner disk. For the latter, we consider the effect of a 5 M_Jup planet in an eccentric orbit on the planetesimal distribution of HD 141569A.

We present a weak lensing analysis of one of the most distant massive galaxy cluster known, RDCS 1252.9-2927 at z=1.24, using deep images from the Advanced Camera for Survey (ACS) on board the Hubble Space Telescope (HST). By taking advantage of the depth and of the angular resolution of the ACS images, we detect for the first time at z>1 a clear weak lensing signal in both the i (F775W) and z (F850LP) filters. We measure a 5-\sigma signal in the i band and a 3-\sigma signal in the shallower z band image. The two radial mass profiles are found to be in very good agreement with each other, and provide a measurement of the total mass of the cluster inside a 1Mpc radius of M(<1Mpc) = (8.0 +/- 1.3) x 10^14 M_\odot in the current cosmological concordance model h =0.70, \Omega_m=0.3, \Omega_\Lambda=0.7, assuming a redshift distribution of background galaxies as inferred from the Hubble Deep Fields surveys. A weak lensing signal is detected out to the boundary of our field (3' radius, corresponding to 1.5Mpc at the cluster redshift). We detect a small offset between the centroid of the weak lensing mass map and the brightest cluster galaxy, and we discuss the possible origin of this discrepancy. The cumulative weak lensing radial mass profile is found to be in good agreement with the X-ray mass estimate based on Chandr and XMM-Newton observations, at least out to R_500=0.5Mpc.

We present new measurements of the galaxy luminosity function (LF) and its dependence on local galaxy density, color, morphology, and clustocentric radius for the massive z=0.83 cluster MS1054-0321. Our analyses are based on imaging performed with the ACS onboard the HST in the F606W, F775W and F850LP passbands and extensive spectroscopic data obtained with the Keck LRIS. Our main results are based on a spectroscopically selected sample of 143 cluster members with morphological classifications derived from the ACS observations. Our three primary findings are (1) the faint-end slope of the LF is steepest in the bluest filter, (2) the LF in the inner part of the cluster (or highest density regions) has a flatter faint-end slope, and (3) the fraction of early-type galaxies is higher at the bright end of the LF, and gradually decreases toward fainter magnitudes. These characteristics are consistent with those in local galaxy clusters, indicating that, at least in massive clusters, the common characteristics of cluster LFs are established at z=0.83. We also find a 2sigma deficit of intrinsically faint, red galaxies (i-z>0.5, Mi>-19) in this cluster. This trend may suggest that faint, red galaxies (which are common in z<0.1 rich clusters) have not yet been created in this cluster at z=0.83. The giant-to-dwarf ratio in MS1054-0321 starts to increase inwards of the virial radius or when Sigma>30 Mpc^-2, coinciding with the environment where the galaxy star formation rate and the morphology-density relation start to appear. (abridged)

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.

We present the mass distribution in the central area of the cluster A1689 by fitting over 100 multiply lensed images with the non-parametric Strong Lensing Analysis Package (SLAP, Diego et al. 2004). The surface mass distribution is obtained in a robust way finding a total mass of 0.25E15 M_sun/h within a 70'' circle radius from the central peak. Our reconstructed density profile fits well an NFW profile with small perturbations due to substructure and is compatible with the more model dependent analysis of Broadhurst et al. (2004a) based on the same data. Our estimated mass does not rely on any prior information about the distribution of dark matter in the cluster. The peak of the mass distribution falls very close to the central cD and there is substructure near the center suggesting that the cluster is not fully relaxed. We also examine the effect on the recovered mass when we include the uncertainties in the redshift of the sources and in the original shape of the sources. Using simulations designed to mimic the data, we identify some biases in our reconstructed mass distribution. We find that the recovered mass is biased toward lower masses beyond 1 arcmin (150 kpc) from the central cD and that in the very center we may be affected by degeneracy problems. On the other hand, we confirm that the reconstructed mass between 25'' and 70'' is a robust, unbiased estimate of the true mass distribution and is compatible with an NFW profile.

We present the results of an extensive spectroscopic survey of RX J0152.7-1357, one of the most massive distant clusters of galaxies known. Multi-object spectroscopy, carried out with FORS1 and FORS2 on the ESO Very Large Telescope (VLT), has allowed us to measure more than 200 redshifts in the cluster field and to confirm 102 galaxies as cluster members. The mean redshift of the cluster is $z=0.837 \pm 0.001$ and we estimate the velocity dispersion of the overall cluster galaxy distribution to be $\sim 1600 \mathrm{km \ s^{-1}}$. The distribution of cluster members is clearly irregular, with two main clumps that follow the X-ray cluster emission mapped by Chandra. A third clump of galaxies to the east of the central structure and at the cluster redshift has also been identified. The two main clumps have velocity dispersions of $\sim919$ and $\sim737 \mathrm{km s^{-1}}$ respectively, and the peculiar velocity of the two clumps suggests that they will merge into a single more massive cluster. A segregation in the star formation activity of the member galaxies is observed. All star forming galaxies are located outside the high-density peaks, which are populated only by passive galaxies. A population of red galaxies (belonging to the cluster red sequence) with clear post-starburst spectral features and [OII] ($\lambda$3727) emission lines is observed in the outskirts of the cluster. Two AGNs, which were previously confused with the diffuse X-ray emission from the intracluster medium in ROSAT and BeppoSAX observations, are found to be cluster members.

(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\alpha 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\alpha 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\alpha line it is apparent that it is very difficult to explain the observed Ly\alpha 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 \alpha ^{-1}, where \alpha^{-1} is the stellar disk exponential scale-length. The cosmic web is the most likely origin for the detected Ly\alpha lines.

ACS observations of massive lensing clusters permit an order of magnitude increase in the numbers of multiply-lensed background galaxies identified behind a given cluster. We have developed a code to take the pixels belonging to any given image and generate counter-images with full resolution, so that multiple systems are convincingly and exhaustively identified. Over 130 images of 35 multiply lensed galaxies are found behind A1689, including many radial arcs and also tiny counter-images projected on the center of mass. The derived mass profile is found to flatten steadily towards the center, like an NFW profile, with a mean slope d(log(Sigma))/d(log(r)) ~ -0.55 +/- 0.1, over the range r<250 kpc/h, which is somewhat steeper than predicted for such a massive halo. We also clearly see the expected geometric increase of bend angles with redshift, however, given the low redshift of A1689, z=0.18, the dependence on cosmological parameters is weak, but using higher redshift clusters from our GTO program we may derive a more competitive constraint.

We analyse deep multi-colour Advanced Camera images of the largest known gravitational lens, A1689. Radial and tangential arcs delineate the critical curves in unprecedented detail and many small counter-images are found near the center of mass. We construct a flexible light deflection field to predict the appearance and positions of counter-images. The model is refined as new counter-images are identified and incorporated to improve the model, yielding a total of 106 images of 30 multiply lensed background galaxies, spanning a wide redshift range, 1.0$<$z$<$5.5. The resulting mass map is more circular in projection than the clumpy distribution of cluster galaxies and the light is more concentrated than the mass within $r<50kpc/h$. The projected mass profile flattens steadily towards the center with a shallow mean slope of $d\log\Sigma/d\log r \simeq -0.55\pm0.1$, over the observed range, r$<250kpc/h$, matching well an NFW profile, but with a relatively high concentration, $C_{vir}=8.2^{+2.1}_{-1.8}$. A softened isothermal profile ($r_{core}=20\pm2$\arcs) is not conclusively excluded, illustrating that lensing constrains only projected quantities. Regarding cosmology, we clearly detect the purely geometric increase of bend-angles with redshift. The dependence on the cosmological parameters is weak due to the proximity of A1689, $z=0.18$, constraining the locus, $\Omega_M+\Omega_{\Lambda} \leq 1.2$. This consistency with standard cosmology provides independent support for our model, because the redshift information is not required to derive an accurate mass map. Similarly, the relative fluxes of the multiple images are reproduced well by our best fitting lens model.

We measure the luminosity function of morphologically selected E/S0 galaxies from $z=0.5$ to $z=1.0$ using deep high resolution Advanced Camera for Surveys imaging data. Our analysis covers an area of $48\Box\arcmin$ (8$\times$ the area of the HDF-N) and extends 2 magnitudes deeper ($I\sim24$ mag) than was possible in the Deep Groth Strip Survey (DGSS). At $0.5<z<0.75$, we find $M_B^*-5\log h_{0.7}=-21.1\pm0.3$ and $\alpha=-0.53\pm0.2$, and at $0.75<z<1.0$, we find $M_B^*-5\log h_{0.7}=-21.4\pm0.2$. These luminosity functions are similar in both shape and number density to the luminosity function using morphological selection (e.g., DGSS), but are much steeper than the luminosity functions of samples selected using morphological proxies like the color or spectral energy distribution (e.g., CFRS, CADIS, or COMBO-17). The difference is due to the `blue', $(U-V)_0<1.7$, E/S0 galaxies, which make up to $\sim30%$ of the sample at all magnitudes and an increasing proportion of faint galaxies. We thereby demonstrate the need for {\it both morphological and structural information} to constrain the evolution of galaxies. We find that the `blue' E/S0 galaxies have the same average sizes and Sersic parameters as the `red', $(U-V)_0>1.7$, E/S0 galaxies at brighter luminosities ($M_B<-20.1$), but are increasingly different at fainter magnitudes where `blue' galaxies are both smaller and have lower Sersic parameters. Fits of the colors to stellar population models suggest that most E/S0 galaxies have short star-formation time scales ($\tau<1$ Gyr), and that galaxies have formed at an increasing rate from $z\sim8$ until $z\sim2$ after which there has been a gradual decline.

We use the exceptional depth of the Ultra Deep Field (UDF) and UDF-Parallel ACS fields to study the sizes of high redshift (z~2-6) galaxies and address long-standing questions about possible biases in the cosmic star formation rate due to surface brightness dimming. Contrasting B, V, and i-dropout samples culled from the deeper data with those obtained from the shallower GOODS fields, we demonstrate that the shallower data are essentially complete at bright magnitudes to z<5.5 and that the principal effect of depth is to add objects at the magnitude limit. This indicates that high redshift galaxies are compact in size (~0.1-0.3") and that large (>0.4", >3 kpc) low surface brightness galaxies are rare. A simple comparison of the half-light radii of the HDF-N + HDF-S U-dropouts with B, V, and i-dropouts from the UDF shows that the sizes follow a (1+z)^{-1.05+/-0.21} scaling towards high redshift. A more rigorous measurement compares different scalings of our U-dropout sample with the mean profiles for a set of intermediate magnitude (26.0<z_{850,AB}<27.5) i-dropouts from the UDF. The best-fit is found with a (1+z)^{-0.94_{-0.25} ^{+0.19}} size scaling (for fixed luminosity). This result is then verified by repeating this experiment with different size measures, low redshift samples, and magnitude ranges. Very similar scalings are found for all comparisons. A robust measurement of size evolution is thereby demonstrated for galaxies from z~6 to z~2.5 using data from the UDF.

We describe the feasibility of detecting the gravitational deflection of light emitted by stars moving under the influence of the massive object at the Galactic center. Light emitted by a star orbiting behind the central mass has a smaller impact parameter than the star itself, and suffers the effect of gravitational lensing, providing a closer probe of the central mass distribution and hence a stricter test of the black hole hypothesis. A mass of $4.3\times 10^{6} M_{\odot}$ causes a $0.1-2\rm mas$ deviation in the apparent position of orbiting stars projected within $10^{\circ}$ of the line of sight to the galactic center. In addtion, we may uniquely constrain the distance to the center of the galaxy because lensing deflections constrain the ratio $\rg/R_{0}$ of the Schwarzschild radius to the distance to the black hole, $R_{o}$, whereas the ratio $\rg/R_{o}^{3}$ is obtained by fitting the orbit.

The properties of Ultra Compact Dwarf (UCD) galaxy candidates in Abell 1689 (z=0.183) are investigated, based on deep high resolution ACS images. A UCD candidate has to be unresolved, have i<28 (M_V<-11.5) mag and satisfy color limits derived from Bayesian photometric redshifts. We find 160 UCD candidates with 22<i<28 mag. It is estimated that about 100 of these are cluster members, based on their spatial distribution and photometric redshifts. For i>26.8 mag, the radial and luminosity distribution of the UCD candidates can be explained well by Abell 1689's globular cluster (GC) system. For i<26.8 mag, there is an overpopulation of 15 +/- 5 UCD candidates with respect to the GC luminosity function. For i<26 mag, the radial distribution of UCD candidates is more consistent with the dwarf galaxy population than with the GC system of Abell 1689. The UCD candidates follow a color-magnitude trend with a slope similar to that of Abell 1689's genuine dwarf galaxy population, but shifted fainter by about 2-3 mag. Two of the three brightest UCD candidates (M_V ~ -17 mag) are slightly resolved. At the distance of Abell 1689, these two objects would have King-profile core radii of ~35 pc and r_eff ~300 pc, implying luminosities and sizes 2-3 times those of M32's bulge. Additional photometric redshifts obtained with late type stellar and elliptical galaxy templates support the assignment of these two resolved sources to Abell 1689. Our findings imply that in Abell 1689 there are at least 10 UCDs with M_V<-12.7 mag. Compared to the UCDs in the Fornax cluster they are brighter, larger and have colors closer to normal dwarf galaxies. This suggests that they may be in an intermediate stage of the stripping process. Spectroscopy is needed to definitely confirm the existence of UCDs in Abell 1689.

(Abridged) We study the internal color properties of a morphologically selected sample of spheroidal galaxies taken from HST/ACS ERO program of UGC 10214 (``The Tadpole''). By taking advantage of the unprecedented high resolution of the ACS in this very deep dataset we are able to characterize spheroids at sub-arcseconds scales. Using the V_606W and I_814W bands, we construct V-I color maps and extract color gradients for a sample of spheroids at I_814W < 24 mag. We investigate the existence of a population of morphologically classified spheroids which show extreme variation in their internal color properties similar to the ones reported in the HDFs. These are displayed as blue cores and inverse color gradients with respect to those accounted from metallicity variations. Following the same analysis we find a similar fraction of early-type systems (~30%-40%) that show non-homologous internal colors, suggestive of recent star formation activity. We present two statistics to quantify the internal color variation in galaxies and for tracing blue cores, from which we estimate the fraction of non-homogeneous to homogeneous internal colors as a function of redshift up to z<1.2. We find that it can be described as about constant as a function of redshift, with a small increase with redshift for the fraction of spheroids that present strong color dispersions. The implications of a constant fraction at all redshifts suggests the existence of a relatively permanent population of evolving spheroids up to z~1. We discuss the implications of this in the context of spheroidal formation.

We present a detailed analysis of the stellar mass content of galaxies up to z=2.5 in the K20 galaxy sample, that has a 92% spectroscopic completeness and a complete $UBVRIzJK_s$ multicolor coverage. We find that the M/L ratio decreases with redshift: in particular, the average M/L ratio of early type galaxies decreases with $z$, with a scatter that is indicative of a range of star--formation time-scales and redshift of formation. More important, the typical M/L of massive early type galaxies is larger than that of less massive ones, suggesting that their stellar population formed at higher z. The final K20 galaxy sample spans a range of stellar masses from M*=10^9Msun to M*=10^12Msun, with massive galaxies ($M*>10^11Msun) detected up to z~2. We compute the Galaxy Stellar Mass Function at various z, of which we observe only a mild evolution (i.e. by 20-30%) up to z~1. At z>1, the evolution of the GSMF appears to be much faster: at z~2, about 35% of the present day stellar mass in objects with M*~10^11Msun appear to have assembled. We also detect a change in the physical nature of the most massive galaxies, since at z>1 a population of massive star--forming galaxies progressively appears. We finally analyze our results in the framework of Lambda-CDM hierarchical models. First, we show that the large number of massive galaxies detected at high z does not violate any fundamental Lambda-CDM constraint based on the number of massive DM halos. Then, we compare our results with the predictions of renditions of both semianalytic and hydro-dynamical models, that range from severe underestimates to slight overestimates of the observed mass density at z<~2. We discuss how the differences among these models are due to the different implementation of the main physical processes. (Abridged)

We report on the i-dropouts detected in two exceptionally deep ACS fields (B_{435}, V_{606}, i_{775}, and z_{850} with 10 sigma limits of 28.8, 29.0, 28.5, and 27.8, respectively) taken in parallel with the UDF NICMOS observations. Using an i-z>1.4 cut, we find 30 i-dropouts over 21 arcmin^2 down to z_AB=28.1, or 1.4 i-dropouts arcmin^{-2}, with significant field-to-field variation (as expected from cosmic variance). This extends i-dropout searches some ~0.9^m further down the luminosity function than was possible in the GOODS field, netting a ~7x increase in surface density. An estimate of the size evolution for UV bright objects is obtained by comparing the composite radial flux profile of the bright i-dropouts (z<27.2) with scaled versions of the HDF-N + HDF-S U-dropouts. The best-fit is found with a (1+z)^{-1.57_{-0.53} ^{+0.50}} scaling in size (for fixed luminosity), extending lower redshift (1<z<5) trends to z~6. Adopting this scaling and the brighter i-dropouts from both GOODS fields, we make incompleteness estimates and construct a z~6 LF in the rest-frame continuum UV (~1350 A) over a 3.5 magnitude baseline, finding a shape consistent with that found at lower redshift. To evaluate the evolution in the LF from z~3.8, we make comparisons against different scalings of a lower redshift B-dropout sample. Though a strong degeneracy is found between luminosity and density evolution, our best-fit model scales as (1+z)^{-2.8} in number and (1+z)^0.1 in luminosity, suggesting a rest-frame continuum UV luminosity density at z~6 which is just 0.38_{-0.07} ^{+0.09}x that at z~3.8. Our inclusion of size evolution makes the present estimate lower than previous z~6 estimates.

We have discovered three globular clusters beyond the Holmberg radius in Hubble Space Telescope Advanced Camera for Surveys images of the gas-rich dark matter dominated blue compact dwarf galaxy NGC2915. The clusters, all of which start to resolve into stars, have M_{V606} = -8.9 to -9.8 mag, significantly brighter than the peak of the luminosity function of Milky Way globular clusters. Their colors suggest a metallicity [Fe/H] ~ -1.9 dex, typical of metal-poor Galactic globular clusters. The specific frequency of clusters is at a minimum normal, compared to spiral galaxies. However, since only a small portion of the system has been surveyed it is more likely that the luminosity and mass normalized cluster content is higher, like that seen in elliptical galaxies and galaxy clusters. This suggests that NGC2915 resembles a key phase in the early hierarchical assembly of galaxies - the epoch when much of the old stellar population has formed, but little of the stellar disk. Depending on the subsequent interaction history, such systems could go on to build-up larger elliptical galaxies, evolve into normal spirals, or in rare circumstances remain suspended in their development to become systems like NGC2915.

We present the analysis of the faint galaxy population in the Advanced Camera for Surveys (ACS) Early Release Observation fields VV 29 (UGC 10214) and NGC 4676. Here we attempt to thoroughly consider all aspects relevant for faint galaxy counting and photometry, developing methods which are based on public software and that are easily reproducible by other astronomers. Using simulations we determine the best SExtractor parameters for the detection of faint galaxies in deep HST observations, paying special attention to the issue of deblending, which significantly affects the normalization and shape of the number count distribution. We confirm, as claimed by Bernstein, Freedman and Madore (2002), that Kron-like magnitudes, such as the ones generated by SExtractor, can miss more than half of the light of faint galaxies, what dramatically affects the slope of the number counts. We present catalogs for the VV 29 and NGC 4676 fields with photometry in the g,V and I bands. We also show that combining the bayesian software BPZ with superb ACS data and new spectral templates enables us to estimate reliable photometric redshifts for a significant fraction of galaxies with as few as three filters. After correcting for selection effects, we measure slopes of 0.32+- 0.01 for 22 < g < 28, 0.34+-0.01 for 22< V <27.5 and 0.33+-0.01 for 22 < I < 27. The counts do not flatten (except perhaps in the g-band), up to the depth of our observations. We find that the faint counts m_{AB}> 25.5 can be well approximated in all our filters by a passive luminosity evolution model based on the COMBO-17 luminosity function (\alpha=-1.5), with a strong merging rate following the prescription of Glazebrook et al. (1994), \phi^*\propto (1+Qz), with Q=4.

We investigate the color-magnitude (CM) relation of galaxies in the distant X-ray selected cluster RDCS1252.9-2927 at z=1.24 using images obtained with the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope in the F775W and F850LP bandpasses. We select galaxies based on morphological classifications extending about 3.5 mag down the galaxy luminosity function, augmented by spectroscopic membership information. At the core of the cluster is an extensive early-type galaxy population surrounding a central pair of galaxies that show signs of dynamical interaction. The early-type population defines a tight sequence in the CM diagram, with an intrinsic scatter in observed (i-z) of $0.029\pm0.007$ mag based on 52 galaxies, or $0.024\pm0.008$ mag for $\sim30$ ellipticals. Simulations using the latest stellar population models indicate an age scatter for the ellipticals of about 34%, with a mean age $\taul\gta2.6$ Gyr (corresponding to $z_L\gta2.7$), and the last star formation occurring at $z_{end}\gta1.5$. Transforming to rest-frame $(U-B)$, we conclude that the slope and scatter in the CM relation for morphologically selected early-type galaxies show little or no evidence for evolution out to $z\approx1.2$. Thus, elliptical galaxies were already well established in X-ray luminous clusters when the universe was a third of its present age.

Spectroscopic redshifts have been measured for 9 K-band luminous galaxies at 1.7 < z < 2.3, selected with Ks < 20 in the "K20 survey" region of the Great Observatories Origins Deep Survey area. Star formation rates (SFRs) of ~100-500 Msun/yr are derived when dust extinction is taken into account. The fitting of their multi-color spectral energy distributions indicates stellar masses M ~ 10^11 Msun for most of the galaxies. Their rest-frame UV morphology is highly irregular, suggesting that merging-driven starbursts are going on in these galaxies. Morphologies tend to be more compact in the near-IR, a hint for the possible presence of older stellar populations. Such galaxies are strongly clustered, with 7 out of 9 belonging to redshift spikes, which indicates a correlation length r_0 ~ 9-17 h^-1 Mpc (1 sigma range). Current semianalytical models of galaxy formation appear to underpredict by a large factor (about 30) the number density of such a population of massive and powerful starburst galaxies at z ~ 2. The high masses and SFRs together with the strong clustering suggest that at z ~ 2 we may have started to explore the major formation epoch of massive early-type galaxies.

Using an i-z dropout criterion, we determine the space density of z~6 galaxies from two deep ACS GTO fields with deep optical-IR imaging. A total of 23 objects are found over 46 arcmin^2, or ~0.5 objects/arcmin^2 down to z~27.3 (6 sigma; all AB mag) (including one probable z~6 AGN). Combining deep ISAAC data for our RDCS1252-2927 field (J~25.7 and Ks~25.0 (5 sigma)) and NICMOS data for the HDF North (JH~27.3 (5 sigma)), we verify that these dropouts have flat spectral slopes. i-dropouts in our sample range in luminosity from ~1.5 L* (z~25.6) to ~0.3 L* (z~27.3) with the exception of one very bright candidate at z~24.2. The half-light radii vary from 0.09" to 0.29", or 0.5 kpc to 1.7 kpc. We derive the z~6 rest-frame UV luminosity density using three different procedures, each utilizing simulations based on a CDF South V dropout sample. First, we compare our findings with a no-evolution projection of this V-dropout sample. We find 23+/-25% more i-dropouts than we predict. Adopting previous results to z~5, this works out to a 20+/-29% drop in the luminosity density from z~3 to z~6. Second, we use these same V-dropout simulations to derive a selection function for our i-dropout sample and compute the UV-luminosity density (7.2+/-2.5 x 10^25 ergs/s/Hz/Mpc^3 down to z~27). We find a 39+/-21% drop over the same redshift range. This is our preferred value and suggests a star formation rate of 0.0090+/-0.0031 M_sol/yr/Mpc^3 to z~27, or ~0.036+/- 0.012 M_sol/yr/Mpc^3 extrapolating the LF to the faint limit. Third, we follow a very similar procedure, but assume no incompleteness, finding a luminosity density which is ~2-3X lower. This final estimate constitutes a lower limit. All three estimates are within the canonical range of luminosity densities necessary for reionization of the universe at this epoch. (abridged)

The evolution of high redshift galaxies in the two Hubble Deep Fields, HDF-N and HDF-S, is investigated using a cloning technique that replicates z~ 2-3 U dropouts to higher redshifts, allowing a comparison with the observed B and V dropouts at higher redshifts (z ~ 4-5). We treat each galaxy selected for replication as a set of pixels that are k-corrected to higher redshift, accounting for resampling, shot-noise, surface-brightness dimming, and the cosmological model. We find evidence for size evolution (a 1.7x increase) from z ~ 5 to z ~ 2.7 for flat geometries (Omega_M+Omega_LAMBDA=1.0). Simple scaling laws for this cosmology predict that size evolution goes as (1+z)^{-1}, consistent with our result. The UV luminosity density shows a similar increase (1.85x) from z ~ 5 to z ~ 2.7, with minimal evolution in the distribution of intrinsic colors for the dropout population. In general, these results indicate less evolution than was previously reported, and therefore a higher luminosity density at z ~ 4-5 (~ 50% higher) than other estimates. We argue the present technique is the preferred way to understand evolution across samples with differing selection functions, the most relevant differences here being the color cuts and surface brightness thresholds (e.g., due to the (1+z)^4 cosmic surface brightness dimming effect).

We present the galaxy rest-frame near-IR Luminosity Function (LF) and its cosmic evolution to z=1.5 based on a spectroscopic survey of a magnitude limited sample of galaxies with Ks<20 (the K20 survey, Cimatti et al. 2002b). The LFs have been derived in the rest-frame J and Ks bands using 3 z bins (z_mean= 0.5, 1, 1.5) and compared to the local near-IR LF. The faint-end of the LFs is consistent with the local estimates, with no evidence for a change either in the slope or normalization up to z<1.3. Viceversa, the density of luminous galaxies (M_Ks-5logh_70<-25.5) is higher than locally at all z and relatively constant or mildly increasing with z within our sample. The data are consistent with a mild luminosity evolution both in the J and Ks-band up to z=1.5 (DeltaM_J=-0.69+-0.12 and DeltaM_K=-0.54+-0.12 at z=1). Moreover, we find that red and early-type galaxies dominate the bright-end of the LF, and that their number density shows at most a small decrease (<30%) up to z=1, thus suggesting that massive elliptical galaxies were already in place at z=1 and they should have formed their stars and assembled their mass at higher z. There appears to be a correlation of the optical/near-IR colors with near-IR luminosities, the most luminous/massive galaxies being red/old. We find a slow evolution with z of the near-IR comoving luminosity density to z=1.5. Finally, we show that hierarchical models overpredict significantly the density of low luminosity galaxies at z<=1 and underpredict the density of luminous galaxies at z>=1, whereas PLE models are more consistent with the data up to z=1.5. The GIF model (Kaufmann et al. 1999) shows a clear deficiency of red luminous galaxies at z=1 compared to our observations and predicts a decrease of luminous galaxies with z not observed in our sample.

Multicolor coronagraphic images of the circumstellar disk around HD141569A have been obtained with the Hubble Space Telescope (HST) Advanced Camera for Surveys. B, V, and I images show that the disk's previously-described multiple-ring structure is actually a continuous distribution of dust with a tightly-wound spiral structure. Two, more open spiral arms extend from the disk, one of which appears to reach the nearby binary star HD141569BC. Diffuse dust is seen up to 1200AU from HD141569A. Although planets may exist in the inner region of the disk, tidal interaction with HD141569BC seems more likely to be the cause of these phenomena. The disk appears redder than the star (B-V = 0.21 and V-I = 0.25), and its color is spatially uniform. A scattering asymmetry factor of g = 0.25-0.35 is derived. The azimuthal density distribution is asymmetric, varying by a factor of ~3 at some radii.

We present observations of the first two supernovae discovered with the recently installed Advanced Camera for Surveys (ACS) on the Hubble Space Telescope. The supernovae were found in Wide Field Camera images of the Hubble Deep Field North taken with the F775W, F850LP, and G800L optical elements as part of the ACS guaranteed time observation program. Spectra extracted from the ACS G800L grism exposures confirm that the objects are Type Ia supernovae (SNe Ia) at redshifts z=0.47 and z=0.95. Follow-up HST observations have been conducted with ACS in F775W and F850LP and with NICMOS in the near-infrared F110W bandpass, yielding a total of 9 flux measurements in the 3 bandpasses over a period of 50 days in the observed frame. We discuss many of the important issues in doing accurate photometry with the ACS. We analyze the multi-band light curves using two different fitting methods to calibrate the supernovae luminosities and place them on the SNe Ia Hubble diagram. The resulting distances are consistent with the redshift-distance relation of the accelerating universe model, although evolving intergalactic grey dust remains as a less likely possibility. The relative ease with which these SNe Ia were found, confirmed, and monitored demonstrates the potential ACS holds for revolutionizing the field of high-redshift SNe Ia, and therefore of testing the accelerating universe cosmology and constraining the "epoch of deceleration".

We present the first Advanced Camera for Surveys (ACS) observations of young star clusters in the colliding/merging galaxy UGC 10214. The observations were made as part of the Early Release Observation (ERO) program for the newly installed ACS during service mission SM3B for the Hubble Space Telescope (HST). Many young star clusters can be identified in the tails of UGC 10214, with ages ranging from ~3 Myr to 10 Myr. The extreme blue V-I (F606W-F814W) colors of the star clusters found in the tail of UGC 10214 can only be explained if strong emission lines are included with a young stellar population. This has been confirmed by our Keck spectroscopy of some of these bright blue stellar knots. The most luminous and largest of these blue knots has an absolute magnitude of M_V = -14.45, with a half-light radius of 161 pc, and if it is a single star cluster, would qualify as a super star cluster (SSC). Alternatively, it could be a superposition of multiple scaled OB associations or clusters. With an estimated age of ~ 4-5 Myr, its derived mass is < 1.3 x 10^6 solar masses. Thus the young stellar knot is unbound and will not evolve into a normal globular cluster. The bright blue clusters and associations are much younger than the dynamical age of the tail, providing strong evidence that star formation occurs in the tail long after it was ejected. UGC 10214 provides a nearby example of processes that contributed to the formation of halos and intra-cluster media in the distant and younger Universe.

We present the redshift distribution of a complete sample of 480 galaxies with Ks<20 distributed over two independent fields covering a total area of 52 arcmin^2. The redshift completeness is 87% and 98% respectively with spectroscopic and high-quality and tested photometric redshifts. The redshift distribution of field galaxies has a median redshift z_{med}\sim 0.80, with \sim 32% and \sim 9% of galaxies at z>1 and z>1.5 respectively. A ``blind'' comparison is made with the predictions of a set of the most recent LambdaCDM hierarchical merging and pure luminosity evolution (PLE) models. The hierarchical merging models overpredict and underpredict the number of galaxies at low-z and high-z respectively, whereas the PLE models match the median redshift and the low-z distribution, still being able to follow the high-z tail of N(z). We briefly discuss the implications of this comparison and the possible origins of the observed discrepancies. We make the redshift distribution publicly available.

The K20 survey is an ESO VLT optical and near-infrared spectroscopic survey aimed at obtaining spectral information and redshifts of a complete sample of about 550 objects to K_s\leq20.0 over two independent fields with a total area of 52 arcmin^2. In this paper we discuss the scientific motivation of such a survey, we describe the photometric and spectroscopic properties of the sample, and we release the $K_s$-band photometric catalog. Extensive simulations showed that the sample is photometrically highly complete to K_s=20. The observed galaxy counts and the R-K_s color distribution are consistent with literature results. We observed spectroscopically 94% of the sample, reaching a spectroscopic redshift identification completeness of 92% to K_s\leq20.0 for the observed targets, and of 87% for the whole sample (i.e. counting also the unobserved targets). Deep spectroscopy was complemented with multi-band deep imaging in order to derive tested and reliable photometric redshifts for the galaxies lacking spectroscopic redshifts. The results show a very good agreement between the spectroscopic and the photometric redshifts with <z_{spe}-z_{phot}>=0.01 and with a dispersion of \sigma_{\Delta z}=0.09. Using both the spectroscopic and the photometric redshifts, we reached an overall redshift completeness of about 98%. The size of the sample, the redshift completeness, the availability of high quality photometric redshifts and multicolor spectral energy distributions make the K20 survey database one of the most complete samples available to date for constraining the currently competing scenarios of galaxy formation and for a variety of other galaxy evolution studies.

We compare the 3D clustering of old passively-evolving and dusty star-forming z~1 EROs from the K20 survey. With detailed simulations of clustering, the comoving correlation length of dusty star-forming EROs is constrained to be less than r_0~2.5 Mpc/h. In contrast, the old EROs are much more positively correlated, with 5.5<~r_0/(Mpc/h)<~16, consistent with previous claims for z~1 field early-type galaxies based on analyses of ERO angular clustering. The low level of clustering of dusty star-forming EROs does not support these to be major mergers building up an elliptical galaxy, or typical counterparts of SCUBA sources, but it is instead consistent with the weak clustering of high redshift blue galaxies and of luminous local IRAS galaxies. Current hierarchical merging models can explain the large r_0 for z~1 field early-type galaxies, but fail in matching their high number density and overall old ages.

We present discovery spectra of a sample of eight lensed galaxies at high redshift, 3.7<z<5.2, selected by their red colors in the fields of four massive clusters: A1689, A2219, A2390, and AC114. Metal absorption lines are detected and observed to be blueshifted by 300-800 km/s with respect to the centroid of Ly-alpha emission. A correlation is found between this blueshift and the equivalent width of the metal lines, which we interpret as a broadening of saturated absorption lines caused by a dispersion in the outflow velocity of interstellar gas. Local starburst galaxies show similar behavior, associated with obvious gas outflows. We also find a trend of increasing equivalent width of Ly-alpha emission with redshift, which may be a genuine evolutionary effect towards younger stellar populations at high redshift with less developed stellar continua. No obvious emission is detected below the Lyman limit in any of our spectra, nor in deep U or B-band images. The UV continua are reproduced well by early B-stars, although some dust absorption would allow a fit to hotter stars. After correcting for the lensing, we derive small physical sizes for our objects, ~0.5-5 kpc/h for a flat cosmology with Omega_m=0.3, Omega_Lambda=0.7. The lensed images are only marginally resolved in good seeing despite their close proximity to the critical curve, where large arcs are visible and hence high magnifications of up to ~20x are inferred. Two objects show a clear spatial extension of the Ly-alpha emission relative to the continuum starlight, indicating a ``breakout'' of the gas. The sizes of our galaxies together with their large gas motion suggests that outflows of gas are common at high redshift and associated with galaxy formation.

We present the results of VLT optical spectroscopy of a complete sample of 78 EROs with R-Ks\geq5 over a field of 52 arcmin^2. About 70% of the 45 EROs with Ks\leq19.2 have been spectroscopically identified with old passively evolving and dusty star-forming galaxies at 0.7<z<1.5. The two classes are about equally populated and for each of them we present and discuss the average spectrum. From the old ERO average spectrum and for Z=Z_{\odot} we derive a minimum age of \sim 3 Gyr, corresponding to a formation redshift of z_f \gtsima 2.4. PLE models with such formation redshifts well reproduce the density of old EROs (consistent with being passively evolving ellipticals), whereas the predictions of the current hierarchical merging models are lower than the observed densities by large factors (up to an order of magnitude). From the average spectrum of the star-forming EROs we estimate a substantial dust extinction with E(B-V) \gtsima 0.5. The star formation rates, corrected for the average reddening, suggest a significant contribution from EROs to the cosmic star-formation density at z \sim 1.

We examine the spatial clustering of extremely red objects (EROs) found in a relatively large survey of 700 arcmin^2, containing 400 galaxies with R-Ks>5 to Ks=19.2. A comoving correlation length r_0=12\pm3 /h Mpc is derived, under the assumption that the selection function is described by a passively evolving early-type galaxy population, with an effective redshift of z~1.2. This correlation length is very similar to that of local L* elliptical galaxies implying, at face value, no significant clustering evolution in comoving coordinates of early-type galaxies to the limiting depth of our sample, z~1.5. A rapidly evolving clustering bias can be designed to reproduce a null result; however, our data do not show the corresponding strong reduction in the average population density expected for consistency with underlying growth of the mass-function. We discuss our data in the context of recent ideas regarding bias evolution. The uncertainty we quote on r_0 accounts for the spikey redshift distribution expected along relatively narrow sightlines, which we quantify with detailed simulations. This is an improvement over the standard use of Limber's equation which, because of its implicit assumption of a smooth selection function, underestimates the true noise by a factor of ~ 3 for the parameters of our survey. We propose a general recipe for the analysis of angular clustering, suggesting that any measurement of the angular clustering amplitude, A, has an intrinsic additional uncertainty of \sigma_A/A=\sqrt{AC}, where AC is the appropriate integral constraint.

We compare the metallicity distribution of Galactic Halo stars with 3D realizations of hierarchical galaxy formation. Outflows from dwarf galaxies enrich the intergalactic medium inhomogeneously, at a rate depending on the local galaxy density. Consequently, the first stars created in small early-forming galaxies are less metal-rich that the first stars formed in more massive galaxies which typically form later. As most halo stars are likely to originate in accreted dwarfs, while disk stars formed out of outflow-enriched gas, this scenario naturally generates a ``metallicity floor'' for old disk stars, which we find to be roughly coincident with the higher end of our predicted metallicity distribution of halo stars, in agreement with observations. The broad and centrally peaked distribution of halo star metallicities is well reproduced in our models, with a natural dispersion depending on the exact accretion history. Our modeling includes the important ``baryonic stripping'' effect of early outflows, which brush away the tenuously held gas in neighboring pre-virialized density perturbations. This stripping process does not significantly modify the predicted shape of the halo star metal distribution but inhibits star-formation and hence the number of accreted stars, helping to reproduce the observed total Galactic halo luminosity and also the lack of low-luminosity local dwarf galaxies relative to N-body predictions.

We show that the winds identified with high-redshift low-mass galaxies may strongly affect the formation of stars in more massive galaxies that form later. With 3D realizations of a simple linear growth model we track gas shocking, metal enrichment, and cooling, together with dark halo formation. We show that outflows typically strip baryonic material out of collapsing intermediate mass halos, suppressing star formation. More massive halos can trap the heated gas but collapse late, leading to a broad bimodal redshift distribution, with a larger characteristic mass associated with the lower redshift peak. This scenario accounts for the observed bell-shaped luminosity function of early type galaxies, explains the small number of Milky Way satellite galaxies relative to Cold Dark Matter models predictions, and provides a possible explanation for the lack of metal poor G-dwarfs in the solar neighborhood and the more general lack of low-metallicity stars in massive galaxies relative to ``closed-box'' models of chemical enrichment. Intergalactic medium heating from outflows should produce spectral distortions in the cosmic microwave background that will be measurable with the next generation of experiments.

We show that the gas in growing density perturbations is vulnerable to the influence of winds outflowing from nearby collapsed galaxies that have already formed stars. This suggests that the formation of nearby galaxies with masses less than 10^9 solar masses is likely to be suppressed, irrespective of the details of galaxy formation. An impinging wind may shock heat the gas of a nearby perturbation to above the virial temperature, thereby mechanically evaporating the gas, or the baryons may be stripped from the perturbation entirely if they are accelerated to above the escape velocity. We show that baryonic stripping is the most effective of these two processes, because shock-heated clouds that are too large to be stripped are able to radiatively cool within a sound-crossing time, limiting evaporation. The intergalactic medium temperatures and star-formation rates required for outflows to have a significant influence on the formation of low-mass galaxies are consistent with current observations, but may soon be examined directly via associated distortions in the cosmic microwave background, and with near-infrared observations from the Next Generation Space Telescope, which may detect the supernovae from early-forming stars.

We examine the gravitational redshift of radiation emitted from within the potential of a cluster. Spectral lines from the intracluster medium (ICM) are redshifted in proportion to the emission-weighted mean potential along the line of sight, amounting to approximately 50 km/s at a radius of 100 kpc/h, for a cluster dispersion of 1200 km/s. We show that the relative redshifts of different ionization states of metals in the ICM provide a unique probe of the three-dimensional matter distribution. An examination of the reported peculiar velocities of cD galaxies in well studied Abell clusters reveals they are typically redshifted by an average of $\sim +200$ km/s. This can be achieved by gravity with the addition of a steep central potential associated with the cD galaxy. Note that in general gravitational redshifts cause a small overestimate of the recessional velocities of clusters by an average of $\sim$ 20 km/s.

HST images of a sample of 285 galaxies with measured z from the CFRS and Autofib-LDSS redshift surveys are analysed to derive the evolution of the merger fraction out to z~1. We have performed visual and machine-based merger identifications, as well as counts of bright pairs of galaxies with magnitude differences less than 1.5 mag. We find that the pair fraction increases with z, with up to ~20% of the galaxies being in physical pairs at z~0.75-1. We derive a merger fraction varying with z as (1+z)^{3.2 +/- 0.6}, after correction for line-of-sight contamination, in excellent agreement with the merger fraction derived from the visual classification of mergers for which m = 3.4 +/- 0.6. After correcting for seeing effects on the ground-based selection of survey galaxies, we conclude that the pair fraction evolves as (1+z)^{2.7 +/- 0.6}. This implies that an average L* galaxy will have undergone 0.8 to 1.8 merger events from z=1 to 0, with 0.5 to 1.2 merger events occuring in a 2 Gyr time span at z~0.9. This result is consistent with predictions from semi-analytical models of galaxy formation. From the simple co-addition of the observed luminosities of the galaxies in pairs, physical mergers are computed to lead to a brightening of 0.5 mag for each pair on average, and a boost in star formation rate of a factor of 2, as derived from the average [O II] equivalent widths. Mergers of galaxies are therefore contributing significantly to the evolution of both the luminosity function and luminosity density of the Universe out to z~1.

Surface photometry has been performed on a sample of 46 field elliptical galaxies. These galaxies are described well by a deVaucouleurs R^{1/4} profile. The sample was selected from the combined Canada-France and LDSS redshift surveys and spans the range 0.20 < z < 1.00. The relationship between galaxy half-light radius and luminosity evolves such that a galaxy of a given size is more luminous by Delta M_B=-0.97 \pm 0.14 mag at z=0.92 and the mean rest-frame color shifts blueward by Delta (U-V) =-0.68 \pm 0.11 at z=0.92 relative to the local cluster relations. Approximately 1/3 of these elliptical galaxies exhibit [OII] 3727 emission lines with equivalent widths > 15 angstroms indicating ongoing star formation. Estimated star-formation rates imply that \le 5% of the stellar mass in the elliptical galaxy population has been formed since z=1. We see no evidence for a decline in the space density of early-type galaxies with look-back time. The <V/V_{max}> statistics and a comparison with local luminosity functions are both consistent with the view that the population of massive early-type galaxies was largely in place by z~1. This implies that merging is not required since that time to produce the present-day space density of elliptical galaxies.

The large-scale maximum at k~0.05 identified in the power-spectrum of galaxy fluctuations provides a co-moving scale for measuring cosmological curvature. In shallow 3D surveys the peak is broad, but appears to be well resolved in 1D, at ~130 Mpc (k=0.048), comprising evenly spaced peaks and troughs. Surprisingly similar behaviour is evident at z=3 in the distribution of Lyman-break galaxies, for which we find a 5 sigma excess of pairs separated by Delta z=0.22pm0.02, equivalent to 85Mpc for Omega=1, increasing to 170 Mpc for Omega=0, with a number density contrast of 30% averaged over 5 independent fields. The combination, 3.2\Omega_m -\Omega_{\Lambda}=0.7, matches the local scale of 130 Mpc, i.e. Omega=0.2\pm0.1 or Omega_{m}=0.4\pm0.1 for the matter-dominated and flat models respectively, with an uncertainty given by the width of the excess correlation. The consistency here of the flat model with SNIa based claims is encouraging, but overshadowed by the high degree of coherence observed in 1D compared with conventional Gaussian models of structure formation. The appearance of this scale at high redshift and its local prominence in the distribution of Abell clusters lends support to claims that the high-z `spikes' represent young clusters. Finally we show that a spike in the primordial power spectrum of delta\rho/\rho=0.3 at k=0.05 has little effect on the CMB, except to exaggerate the first Doppler peak in flat matter-dominated models, consistent with recent observations. \\effect on the CMB, except to exaggerate the first Doppler peak in flat matter-dominated models, consistent with recent observations.

The optical-IR images of the Northern and Southern Hubble Deep Fields are used to measure the spectral and density evolution of early-type galaxies. The mean optical SED is found to evolve passively towards a mid F-star dominated spectrum by z ~ 2. We demonstrate with realistic simulations that hotter ellipticals would be readily visible if evolution progressed blueward and brightward at z > 2, following a standard IMF. The colour distributions are best fitted by a `red' IMF, deficient above ~2 M_solar and with a spread of formation in the range 1.5 < z_f < 2.5. Traditional age dating is spurious in this context, a distant elliptical can be young but appear red, with an apparent age >3 Gyrs independent of its formation redshift. Regarding density evolution, we demonstrate that the sharp decline in numbers claimed at z > 1 results from a selection bias against distant red galaxies in the optical, where the flux is too weak for morphological classification, but is remedied with relatively modest IR exposures revealing a roughly constant space density to z ~ 2. We point out that the lack of high mass star-formation inferred here and the requirement of metals implicates cooling-flows of pre-enriched gas in the creation of the stellar content of spheroidal galaxies. Deep-field X-ray images will be very helpful to examine this possibility.

We report on the spectroscopic identification and the long awaited redshift measurement of the heavily obscured, gravitationally lensed radio source PKS 1830-211, which was first observed as a radio Einstein ring. The NE component of the doubly imaged core is identified, in our infrared spectrum covering the wavelength range 1.5-2.5 microns, as an impressively reddened quasar at z=2.507. Our redshift measurement, together with the recently measured time delay (Lovell et al.), means that we are a step closer to determining the Hubble constant from this lens. Converting the time delay into the Hubble constant by using existing models leads to high values for the Hubble constant. Since the lensing galaxy lies very close to the center of the lensed ring, improving the error bars on the Hubble constant will require not only a more precise time delay measurement, but also very precise astrometry of the whole system.

We present a spectroscopic redshift of z=1.675 for the well-known multiply lensed system of arcs seen in the z=0.39 cluster Cl0024+16. In contrast to earlier work, we find that the lensed images are accurately reproduced by a projected mass distribution which traces the locations of the brightest cluster ellipticals, suggesting that the most significant minima of the cluster potential are not fully erased. The averaged mass profile is shallow and consistent with predictions of recent numerical simulations. The source redshift enables us to determine an enclosed cluster mass of M(<100kpc/h)=1.11+/-0.03 x 10^{14}M_{\odot}/h and a mass-to-light ratio of M/L_B(<100kpc/h)=320h(M/L_B)_{\odot}, after correction for passive stellar evolution. The arc spectrum contains many ionized absorption lines and closely resembles that of the local Wolf-Rayet galaxy NGC4217. Our lens model predicts a high magnification (~20) for each image and identifies a new pair of multiple images at a predicted redshift of z=1.3.

A substantial population of high redshift early-type galaxies is detected in very deep UBVRIJHK images towards the HDF-South. Four elliptical profile galaxies are identified in the redshift range z=1-2, all with very red SEDs, implying ages of >2 Gyrs for standard passive evolution. We also find later type IR-luminous galaxies at similarly high redshift, (10 objects with z>1, H<25), whose weak UV emission implies single burst ages of >1 Gyr. The number and luminosity-densities of these galaxies are comparable with the local E/SO-Sbc populations for \Omega_m>0.2, and in the absence of a significant cosmological constant, we infer that the major fraction of luminous Hubble-sequence galaxies have evolved little since z~2. A highly complete photometric redshift distribution is constructed to H=25 (69 galaxies) showing a broad spread of redshift, peaking at z~1.5, in reasonable agreement with some analyses of the HDF. Five `dropout' galaxies are detected at z~3.8, which are compact in the IR, ~0.5 kpc/h at rest 3500\AA. No example of a blue IR luminous elliptical is found, restricting the star-formation epoch of ellipticals to z>10 for a standard IMF and modest extinction.

We detect a distant cluster of galaxies centered on the QSO lens and luminous X-ray source AXJ2019+112, a.k.a. ``The Dark Cluster'' (Hattori et al 1997). Using deep V,I Keck images and wide-field K_s imaging from the NTT, a tight red sequence of galaxies is identified within a radius of 0.2 h^{-1} Mpc of the known z=1.01 elliptical lensing galaxy. The sequence, which includes the central elliptical galaxy, has a slope in good agreement with the model predictions of Kodama et al (1998) for z~1. We estimate the integrated rest-frame luminosity of the cluster to be L_V > 3.2 x 10^{11}h^{-2}L_{\sun} (after accounting for significant extinction at the low latitude of this field), more than an order of magnitude higher than previous estimates. The central region of the cluster is deconvolved using the technique of Magain, Courbin & Sohy (1998), revealing a thick central arc coincident with an extended radio source. All the observed lensing features are readily explained by differential magnification of a radio loud AGN by a shallow elliptical potential. The QSO must lie just outside the diamond caustic, producing two images, and the arc is a highly magnified image formed from a region close to the center of the host galaxy, projecting inside the caustic. The mass--to--light ratio within an aperture of 0.4 h ^{-1} Mpc is M_x/L_V= 224^{+112}_{-78}h(M/L_V)_{\sun}, using the X-ray temperature. The strong lens model yields a compatible value, M/L_V= 372^{+94}_{-94}h(M/L_V)_{\sun}, whereas an independent weak lensing analysis sets an upper limit of M/L_V <520 h(M/L_V)_{\sun}, typical of massive clusters.

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 all 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 data, resulting in entirely empirical and parameter-free ``no-evolution'' deep fields of galaxies for direct comparison with the HDF. In addition, the overcounting rate and the level of incompleteness can be accurately quantified by this procedure. 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 interesting 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.

In a companion paper we outlined a methodology for generating parameter-free, model-independent ``no-evolution'' fields of faint galaxy images, demonstrating the need for significant evolution in the HDF at faint magnitudes. Here we incorporate evolution into our procedure, by transforming the input bright galaxy images with redshift, for comparison with the HDF at faint magnitudes. Pure luminosity evolution is explored assuming that galaxy surface brightness evolves uniformly, at a rate chosen to reproduce the I-band counts. This form of evolution exacerbates the size discrepancy identified by our no-evolution simulations, by increasing the area of a galaxy visible to a fixed isophote. Reasonable dwarf-augmented models are unable to generate the count excess invoking moderate rates of stellar evolution. A plausible fit to the counts and sizes is provided by `mass-conserving' density-evolution, consistent with small-scale hierarchical growth, where the product of disk area and space density is conserved with redshift. Here the increased surface brightness generated by stellar evolution is accomodated by the reduced average galaxy size, for a wide range of geometries. These models are useful for assessing the limitations of the HDF images, by calculating their rates of incompleteness and the degree of over-counting. Finally we demonstrate the potential for improvement in quantifying evolution at fainter magnitudes using the HST Advanced Camera, with its superior UV and optical performance.

We present preliminary results of a neutral hydrogen (HI) redshift survey to find low surface brightness (LSB) galaxies in the very nearby universe. Our sample consists of all galaxies in the APM catalog (Maddox et al 1990) with a mean surface brightness of mu > 24 mag/arcsec^2, down to a magnitude limit of b_j < 17. With the Parkes 64m Radiotelescope 35 objects were detected at v < 4300 km/s. The resulting luminosity function, HI mass function, and for the first time field mass function are presented. It is found that LSBs make a negligible contribution to the overall integrated luminosity, HI mass, and total mass contained in galaxies.

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.

This paper presents the results of simulated NGST observations. The goal of this work is to study the detectability of very high redshift galaxies and to try to characterize Wide Deep-Field observations. The main points made in this paper are 1) high-z galaxies observable by NGST would be similar to local starburst galaxies 2) even if they exist, large spiral galaxies similar to local ones would not be observed by a NGST-like telescope and 3) the density of high-z galaxies detected so far is apparently moderate.

New high quality Keck and ESO images of PKS 1830-211 are presented. Applying a powerful new deconvolution algorithm to these optical and infrared data, both images of the flat spectrum core of the radio source have been identified. An extended source is also detected in the optical images, consistent with the expected location of the lensing galaxy. The source counterparts are very red at I-K~7, suggesting strong Galactic absorption with additional absorption by the lensing galaxy at z=0.885, and consistent with the detection of high redshift molecules in the lens.

We present new Keck observations of giant arcs in the cluster Abell 2390. High resolution two-dimensional spectra of two arcs show metal lines at z=4.040 +/- 0.005 with Lyman-alpha emission spatially separated from the stellar continuum. In addition, spectroscopy along the notorious `straight' arc reveals two unrelated galaxies, at z=0.913 and z=1.033, with absorption lines of MgII and FeII at z=0.913 seen against the more distant object, indicating the presence of a large gaseous halo.

We discuss BIMA (Berkeley Illinois Maryland Association) data and present new high quality optical and near-IR Keck images of the bright radio ring PKS 1830-211. Applying a powerful new deconvolution algorithm we have been able to identify both images of the radio source. In addition we recover an extended source in the optical, consistent with the expected location of the lensing galaxy. The source counterparts are very red, I-K=7, suggesting strong Galactic absorption with additional absorption by the lensing galaxy at z=0.885, and consistent with the detection of high redshift molecules in the lens.

We present the first application of lens magnification to measure the absolute mass of a galaxy cluster; Abell 1689. The absolute mass of a galaxy cluster can be measured by the gravitational lens magnification of a background galaxy population by the cluster potential. The lensing signal is complicated by the variation in number counts due to galaxy clustering and shot-noise, and by additional uncertainties in relating magnification to mass in the strong lensing regime. Clustering and shot-noise can be dealt with using maximum likelihood methods. Local approximations can then be used to estimate the mass from magnification. Alternatively if the lens is axially symmetric we show that the amplification equation can be solved nonlocally for the surface mass density and the tangential shear. In this paper we present the first maps of the total mass distribution in Abell 1689, measured from the deficit of lensed red galaxies behind the cluster. Although noisier, these reproduce the main features of mass maps made using the shear distortion of background galaxies but have the correct normalisation, finally breaking the ``sheet-mass'' degeneracy that has plagued lensing methods based on shear. We derive the cluster mass profile in the inner 4' (0.48 Mpc/h). These show a profile with a near isothermal surface mass density \kappa = (0.5+/-0.1)(\theta/1')^{-1} out to a radius of 2.4' (0.28Mpc/h), followed by a sudden drop into noise. We find that the projected mass interior to 0.24 h^{-1}$Mpc is M(<0.24 Mpc/h)=(0.50+/- 0.09) \times 10^{15} Msol/h. We compare our results with masses estimated from X-ray temperatures and line-of-sight velocity dispersions, as well as weak shear and lensing arclets and find all are in fair agreement for Abell 1698.

We have recently obtained deep near-IR Keck imaging of a newly-discovered z=4.04 galaxy (Frye & Broadhurst 1998). This is lensed by the rich foreground cluster Abell 2390 (z=0.23) into highly-magnified arcs 3-5" in length. Our H- and K'-band Keck/NIRC imaging allows us to map the Balmer+4000Ang break amplitude. In combination with high-quality archival HST/WFPC2 data, we can spatially resolve stellar populations along the arcs. The WFPC2 images clearly reveal several bright knots, which may correspond to sites of active star formation. Indeed, in some spatial regions the Keck/LRIS discovery spectra are consistent with OB-star spectral energy distributions in the rest-ultraviolet. However, there are considerable portions of the arcs which appear redder with no Ly-alpha emission, consistent with being post-starburst regions.

We describe the properties of a system of red arcs discovered at z=4.04 around the cluster A2390 ($z=0.23$). These arcs are images of a single galaxy, where the lensing is compounded by an elliptical cluster member lying close to the critical curve of the cluster. The combined magnification is estimated to be $\sim 20$, using HST images, depending the gradient of the model potentials, implying an unlensed magnitude for the source of $I_{AB}=25.5$. Keck spectroscopy reveals a continuum well fitted by B-stars and attenuated by the Lyman-series forest with an opacity consistent with z$\sim$4 QSO spectra, making the arcs relatively red. Damped Ly$\alpha$ absorption is observed at the source redshift with a relatively high column, n$_{HI}=3\times 10^{21}$. Curiously, Ly$\alpha$ emission is spatially separated ($\sim 0.5$kpc) from the bright stellar continuum and lies $\sim$ 300Km/s redward of interstellar absorption lines at z=4.04. Similar redwardl of young galaxies. Finally, we briefly comment on the notorious ``straight arc'' in A2390, which is resolved into two unrelated galaxies at z=0.913 and z=1.033.

(Abridged) We analyse Hubble Space Telescope images of a complete sample of 341 galaxies drawn from both the Canada France and Autofib/Low Dispersion Survey Spectrograph ground-based redshift surveys. We discuss morphological classifications of these galaxies, and quantify possible biases that may arise from various redshift-dependent effects. We then discuss these biases in the context of automated classifications, and quantify the expected misclassification in our system. After allowing for such biases, the redshift distribution for normal spirals, together with their luminosity function derived as a function of redshift, indicates approximately 1 magnitude of luminosity evolution in B(AB) by z=1. The elliptical sample is too small for precise evolutionary constraints. However, we find a substantial increase in the proportion of galaxies with irregular morphology at large redshift. These galaxies also appear to be the dominant cause of the rapid rise with redshift in the blue luminosity density identified in the redshift surveys. Although galaxies with irregular morphology may well comprise a mixture of different physical systems and might not correspond to present day irregulars, it is clear that the apparently declining abundance and luminosities of our distant ``irregulars'' holds an important key to understanding recent evolution in the star formation history of normal galaxies.

(Abridged) Several aspects of the evolution of star-forming galaxies are studied using measures of the 2-dimensional surface brightness profile of a sample of 341 faint objects selected from the CFRS and LDSS redshift surveys that have been observed with the Hubble Space Telescope. The size function of disk scale lengths in disk-dominated galaxies is found to stay roughly constant to z=1. This result, which is strengthened by inclusion of the local de Jong et al (1996) size function, suggests that the scale lengths of typical disks can not have grown substantially with cosmic epoch since z=1, unless a corresponding number of large disks have been destroyed through merging. In addition to a roughly constant number density, the galaxies with large disks have, as a set, properties consistent with the idea that they are similar galaxies observed at different cosmic epochs. However, on average, they show higher B-band disk surface brightnesses, bluer overall (U-V) colors, higher [OII] 3727 equivalent widths and less regular morphologies at high redshift than at low redshift, suggesting an increase in the star-formation rate by a factor of about 3 to z=0.7. This is consistent with the expectations of recent models for the evolution of the disk of the Milky Way and similar galaxies. The evolution of the large disk galaxies is probably not sufficient to account for the evolution of the overall luminosity function of galaxies over the interval 0<z<1, especially if Omega=1. Analysis of the half-light radii of all the galaxies in the sample and construction of the bivariate size-luminosity function suggests that larger changes in the galaxy population are due to smaller galaxies, those with half-light radii around 5 h_{50}^{-1} kpc

We present a novel measurement of the pairwise peculiar velocity distribution function of galaxies on scales $r<3200$ \kmsec in the Las Campanas Redshift Survey. The distribution is well described by a scale-independent exponential with a width $\tau$ where $\sig =\sqrt{2}\tau=363$ \kmsec. The signal is very stable. Results from the Northern and Southern sets of slices agree within $\pm13$ \kmsec and the fluctuations among the six individual survey slices vary as $\pm44$ \kmsec. The distribution was determined using a Fourier-space deconvolution of the redshift-space distortions in the correlation function. This technique is insensitive to the effect of rich clusters in the survey and recovers the entire distribution function rather than just its second moment. Taken together with the large effective volume of the survey $6.0\times10^6 \Mpc^3$, we believe this to be a definitive result for $r$-band selected galaxies with absolute magnitudes $-18.5 <M_r<-22.5$ and $z<0.2$.

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 \phi^*). 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.

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.

Two colour photometry of the cluster A1689 reveals a `relative magnification-bias' between lensed blue and red background galaxies, arising from a dependence of the faint galaxy count-slope on colour. The colour distribution is skewed blueward of the far field, allowing us to measure the cluster magnification and to understand the notorious blueness of large arcs. We show that the magnification information can be combined with the usual image distortion measurements to isolate the local `convergence' component of lensing and hence derive the projected mass. This is achieved through a simple local relation between the convergence and the observables, which can be applied generally over the surface a cluster. In the weak lensing limit, the convergence reduces to a dependence on the magnification alone, so that in the outskirts of clusters the surface-density of matter is obtained directly from the surface-density of background galaxies. Hence, useful lensing work requires colour information but not necessarily good seeing. Interestingly, convergence varies slowly at high redshift, saturating at a level depending on the Horizon distance, allowing a useful model-independent measurement of the Global Geometry.

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.

We show that the high redshift IRAS source F10214 is highly magnified by the gravitational field of an intervening elliptical galaxy, accounting for its many anomalous properties. Detailed radio and near-IR images identify the IRAS source with a symmetric arc, centered on a red object, or lensing galaxy. To explain the observed structures, the center of the source must much more highly magnified than its outer regions. Lensing predicts a small counterimage to the arc, which we find adjacent to the lensing galaxy. A red component in the observed spectrum suggests a lens redshift of unity, and the lens model yields a mass estimate of $M(r<3kpc)\approx10^{11}{M_solar}$, consistent with an ordinary elliptical galaxy. We present new high-resolution optical images which show a thin arc of emission, implying an intrinsically small source ($<0.5kpc$) which is highly magnified ($\sim20\times$). Since the optical is strongly polarized with a Seyfert~II spectrum, we propose that the optical arc is magnified image of the inner region of an obscured AGN. The obscuring ``torus'' will be similarly magnified, naturally accounting for the large IR flux. We show that finding objects like F10214+4724 in redshift surveys is probable, given the level of magnification bias expected for compact luminous IRAS sources. Such cases represent the obscured AGN counterparts to the lensed QSO population and, because of their extended sizes, are useful in determining the mass distribution in the lensing galaxies.

Using the multislit LDSS-2 spectrograph on the {\it William Herschel Telescope} we have completed a redshift survey in the magnitude range $22.5<B< 24$ which has produced 73 redshifts representing a 73\% complete sample uniformly-selected from four deep fields at high Galactic latitude. The survey extends out to $z>1$ and includes the highest redshift galaxy ($z=1.108$) yet discovered in a field sample. The median redshift, $\zmed=0.46$, and form of the redshift distribution constitute compelling evidence against simple luminosity evolution as an explanation of the large excess of faint galaxies ($\simeq\times$2--4 no-evolution) seen in this magnitude range. Rather we identify the excess population as blue objects with $z\sim 0.4$ and $B$\, luminosities similar to local $L^*$ galaxies indicating a dramatic decrease in the density of such objects over the last Hubble time, confirming the trends found in brighter redshift surveys. We also find a marked absence of {\it very} low redshift galaxies ($z<$0.1) at faint limits, severely constraining any significant steepening of the local field galaxy luminosity function at low luminosities.

Measurement of the gravitational distortion of images of distant galaxies is rapidly becoming established as a powerful probe of the dark mass distribution in clusters of galaxies. With the advent of large mosaics of CCD's these methods should provide a composite total mass profile for galaxy haloes and also measure the power spectrum of mass fluctuations on supercluster scales. We describe how HST observations have been used to place the observational measurement of the shear on a quantitative footing. By artifically stretching and then degrading WFPC2 images to simulate ground based observing it is now possible to directly calibrate the effect of atmospheric seeing. Similar experiments show that one can remove the effect of artificial image anisotropy arising in the atmosphere or telescope. There have also been important advances in the interpretation of the shear: reconstruction techniques have been extended to encompass the strong distortion regime of giant arcs etc., progress has been made in removing a bias present in earlier reconstruction techniques, and we describe new techniques for `aperture densitometry'. We present some new results on clusters of galaxies, and discuss the intimate connections between weak lensing and deep spectroscopy.

We develop and test a method for measuring the gravitational lensing induced distortion of faint background galaxies. We first describe how we locate the galaxies and measure a 2-component `polarisation' or ellipticity statistic $e_\alpha$ whose expectation value should be proportional to the gravitational shear $\gamma_\alpha$. We then show that an anisotropic instrumental psf perturbs the polarisation by $\delta e_\alpha = P^s_{\alpha\beta} p_\beta$, where $p_\alpha$ is a measure of the psf anisotropy and $P^s_{\alpha\beta}$ is the `linearised smear polarisability tensor'. By estimating $P^s_{\alpha\beta}$ for each object we can determine $p_\alpha$ from the foreground stars and apply a correction $-P^s_{\alpha\beta}p_\beta$ to the galaxies. We test this procedure using deep high-resolution images from HST which are smeared with an anisotropic psf and then have noise added to simulate ground-based observations. We find that the procedure works very well. A similar analysis yields a linear shear polarisability tensor $P^\gamma_{\alpha\beta}$ which describes the response to a gravitational shear. This calibrates the polarisation-shear relation, but only for galaxies which are well resolved. To empirically calibrate the effect of seeing on the smaller galaxies we artificially stretch HST images to simulate lensing and then degrade them as before. These experiments provide a rigorous and exacting test of the method under realistic conditions. They show that it is possible to remove the effect of instrumental psf anisotropy, and that the method provides an efficient and quantitative measurement of the gravitational shear.

We present a new method for measuring the projected mass distributions of galaxy clusters. The gravitational amplification is measured by comparing the joint distribution in redshift and magnitude of galaxies behind the cluster with that of field galaxies. We show that the total amplification is directly related to the surface mass density in the weak field limit, and so it is possible to map the mass distribution of the cluster. The method is shown to be limited by discreteness noise and galaxy clustering behind the lens. Galaxy clustering sets a lower limit to the error along the redshift direction, but a clustering independent lensing signature may be obtained from the magnitude distribution at fixed redshift. Statistical techniques are developed for estimating the surface mass density of the cluster. We extend these methods to account for any obscuration by cluster halo dust, which may be mapped independently of the dark matter. We apply the method to a series of numerical simulations and show the feasibility of the approach. We consider approximate redshift information, and show how the mass estimates are degraded.

We have used HRCam on the CFHT to obtain subarcsecond images of 26 galaxies with z=0.1--0.7. The primary sample of 17 galaxies have enhanced star formation indicated by [OII] equivalent widths greater than 20\AA, while the comparison sample of 9 galaxies have equivalent widths less than 10\AA. By fitting disks or bulges to BVI images we derive scalelengths for the blue and red stellar populations and locate the star-formation in each galaxy. These moderate-redshift galaxies generally have straightforward low-redshift analogues, in that their colours, sizes and luminosities are consistent with those of various types of z$\approx$0 galaxies. The star-forming objects have structures consistent with the full range of present-day disk galaxies, and absolute magnitudes M$^\ast$$-$1 to M$^\ast$+5. Some have star-formation concentrated in their nuclei but in most it occurs across the entire disk. We find companions at projected distances closer than 10 kpc for 30% of the star-forming galaxies, whereas none of the comparison sample have such close companions. This provides the first direct evidence linking interactions or mergers to the increased fraction of star-forming field galaxies at moderate redshifts.