Multi-messenger astrophysics is a fast-growing, interdisciplinary field that combines data, which vary in volume and speed of data processing, from many different instruments that probe the Universe using different cosmic messengers: electromagnetic waves, cosmic rays, gravitational waves and neutrinos. In this Expert Recommendation, we review the key challenges of real-time observations of gravitational wave sources and their electromagnetic and astroparticle counterparts, and make a number of recommendations to maximize their potential for scientific discovery. These recommendations refer to the design of scalable and computationally efficient machine learning algorithms; the cyber-infrastructure to numerically simulate astrophysical sources, and to process and interpret multi-messenger astrophysics data; the management of gravitational wave detections to trigger real-time alerts for electromagnetic and astroparticle follow-ups; a vision to harness future developments of machine learning and cyber-infrastructure resources to cope with the big-data requirements; and the need to build a community of experts to realize the goals of multi-messenger astrophysics.

Statistical analyses from exoplanet surveys around low-mass stars indicate that super-Earth and Neptune-mass planets are more frequent than gas giants around such stars, in agreement with core accretion theory of planet formation. Using precise radial velocities derived from visual and near-infrared spectra, we report the discovery of a giant planet with a minimum mass of 0.46 Jupiter masses in an eccentric 204-day orbit around the very low-mass star GJ 3512. Dynamical models show that the high eccentricity of the orbit is most likely explained from planet-planet interactions. The reported planetary system challenges current formation theories and puts stringent constraints on the accretion and migration rates of planet formation and evolution models, indicating that disc instability may be more efficient in forming planets than previously thought.

Context. Teegarden's Star is the brightest and one of the nearest ultra-cool dwarfs in the solar neighbourhood. For its late spectral type (M7.0V), the star shows relatively little activity and is a prime target for near-infrared radial velocity surveys such as CARMENES. Aims. As part of the CARMENES search for exoplanets around M dwarfs, we obtained more than 200 radial-velocity measurements of Teegarden's Star and analysed them for planetary signals. Methods. We find periodic variability in the radial velocities of Teegarden's Star. We also studied photometric measurements to rule out stellar brightness variations mimicking planetary signals. Results. We find evidence for two planet candidates, each with $1.1M_\oplus$ minimum mass, orbiting at periods of 4.91 and 11.4 d, respectively. No evidence for planetary transits could be found in archival and follow-up photometry. Small photometric variability is suggestive of slow rotation and old age. Conclusions. The two planets are among the lowest-mass planets discovered so far, and they are the first Earth-mass planets around an ultra-cool dwarf for which the masses have been determined using radial velocities.

We outline a simple procedure designed for \emph{automatically} finding sets of multiple images in strong lensing (SL) clusters. We show that by combining (a) an arc-finding (or source extracting) program, (b) photometric redshift measurements, and (c) a preliminary light-traces-mass lens model, multiple-image systems can be identified in a fully automated (`blind') manner. The presented procedure yields an assessment of the likelihood of each arc to belong to one of the multiple-image systems, as well as the preferred redshift for the different systems. These could be then used to automatically constrain and refine the initial lens model for an accurate mass distribution. We apply this procedure to \emph{Cluster Lensing And Supernova with Hubble} observations of three galaxy clusters, MACS J0329.6-0211, MACS J1720.2+3536, and MACS J1931.8-2635, comparing the results to published SL analyses where multiple images were verified by eye on a particular basis. In the first cluster all originally identified systems are recovered by the automated procedure, and in the second and third clusters about half are recovered. Other known systems are not picked up, in part due to a crude choice of parameters, ambiguous photometric redshifts, or inaccuracy of the initial lens model. On top of real systems recovered, some false images are also mistakenly identified by the procedure, depending on the thresholds used. While further improvements to the procedure and a more thorough scrutinisation of its performance are warranted, the work constitutes another important step toward fully automatising SL analyses for studying mass distributions of large cluster samples.

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.

Supernovae (SNe) could be powerful probes of the properties of stars and galaxies at high redshifts in future surveys. Wide fields and longer exposure times are required to offset diminishing star formation rates and lower fluxes to detect useful numbers of events at high redshift. In principle, the Large Synoptic Survey Telescope (LSST) could discover large numbers of early SNe because of its wide fields but only at lower redshifts because of its AB mag limit of ~ 24. But gravitational lensing by galaxy clusters and massive galaxies could boost flux from ancient SNe and allow LSST to detect them at earlier times. Here, we calculate detection rates for lensed SNe at z ~ 5 - 7 for LSST. We find that the LSST Wide Deep Fast survey could detect up to 120 lensed Population (Pop) I and II SNe but no lensed Pop III SNe. Deep-drilling programs in a single 10 square degree FoV could detect Pop I and II core-collapse SNe and Pop III pair-instability SNe at AB magnitudes of 27 - 28 and 26, respectively. An alternative deep survey over 80 nights with a one-year cadence could find ~ 8 Pop III SNe.

[abridged] Characterizing the number counts of faint, dusty star-forming galaxies is currently a challenge even for deep, high-resolution observations in the FIR-to-mm regime. They are predicted to account for approximately half of the total extragalactic background light at those wavelengths. Searching for dusty star-forming galaxies behind massive galaxy clusters benefits from strong lensing, enhancing their measured emission while increasing spatial resolution. Derived number counts depend, however, on mass reconstruction models that properly constrain these clusters. We estimate the 1.1 mm number counts along the line of sight of three galaxy clusters, i.e. Abell 2744, MACSJ0416.1-2403 and MACSJ1149.5+2223, which are part of the ALMA Frontier Fields Survey. We perform detailed simulations to correct these counts for lensing effects. We use several publicly available lensing models for the galaxy clusters to derive the intrinsic flux densities of our sources. We perform Monte Carlo simulations of the number counts for a detailed treatment of the uncertainties in the magnifications and adopted source redshifts. We find an overall agreement among the number counts derived for the different lens models, despite their systematic variations regarding source magnifications and effective areas. Our number counts span ~2.5 dex in demagnified flux density, from several mJy down to tens of uJy. Our number counts are consistent with recent estimates from deep ALMA observations at a 3$\sigma$ level. Below $\approx$ 0.1 mJy, however, our cumulative counts are lower by $\approx$ 1 dex, suggesting a flattening in the number counts. In our deepest ALMA mosaic, we estimate number counts for intrinsic flux densities $\approx$ 4 times fainter than the rms level. This highlights the potential of probing the sub-10 uJy population in larger samples of galaxy cluster fields with deeper ALMA observations.

We present EasyCritics, an algorithm to detect strongly-lensing groups and clusters in wide-field surveys without relying on a direct recognition of arcs. EasyCritics assumes that light traces mass in order to predict the most likely locations of critical curves from the observed fluxes of luminous red early-type galaxies in the line of sight. The positions, redshifts and fluxes of these galaxies constrain the idealized gravitational lensing potential as a function of source redshift up to five free parameters, which are calibrated on few known lenses. From the lensing potential, EasyCritics derives the critical curves for a given, representative source redshift. The code is highly parallelized, uses fast Fourier methods and, optionally, GPU acceleration in order to process large datasets efficiently. The search of a $\smash{1 \, \mathrm{deg}^2}$ field of view requires less than 1 minute on a modern quad-core CPU, when using a pixel resolution of $0.25''/\mathrm{px}$. In this first part of a paper series on EasyCritics, we describe the main underlying concepts and present a first demonstration on data from the Canada-France-Hawaii-Telescope Lensing Survey. We show that EasyCritics is able to identify known group- and cluster-scale lenses, including a cluster with two giant arc candidates that were previously missed by automated arc detectors.

We report the results of $EasyCritics$, a fully automated algorithm for the efficient search of strong-lensing (SL) regions in wide-field surveys, applied to the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). By using only the photometric information of the brightest elliptical galaxies distributed over a wide redshift range ($\smash{0.2 \lesssim z \lesssim 0.9}$) and without requiring the identification of arcs, our algorithm produces lensing potential models and catalogs of critical curves of the entire survey area. We explore several parameter set configurations in order to test the efficiency of our approach. In a specific configuration, $EasyCritics$ generates only $\sim1200$ possibly super-critical regions in the CFHTLS area, drastically reducing the effective area for inspection from $154$ sq. deg to $\sim0.623$ sq. deg, $i.e.$ by more than two orders of magnitude. Among the pre-selected SL regions, we identify 32 of the 44 previously known lenses on the group and cluster scale, and discover 9 new promising lens candidates. The detection rate can be easily improved to $\sim82\%$ by a simple modification in the parameter set, but at the expense of increasing the total number of possible SL candidates. Note that $EasyCritics$ is fully complementary to other arc-finders since we characterize lenses instead of directly identifying arcs. Although future comparisons against numerical simulations are required for fully assessing the efficiency of $EasyCritics$, the algorithm seems very promising for upcoming surveys covering $\smash{10^{4}}$ sq. deg, such as the $Euclid$ mission and $LSST$, where the pre-selection of candidates for any kind of SL analysis will be indispensable due to the expected enormous data volume.

The CARMENES radial velocity (RV) survey is observing 324 M dwarfs to search for any orbiting planets. In this paper, we present the survey sample by publishing one CARMENES spectrum for each M dwarf. These spectra cover the wavelength range 520--1710nm at a resolution of at least $R > 80,000$, and we measure its RV, H$\alpha$ emission, and projected rotation velocity. We present an atlas of high-resolution M-dwarf spectra and compare the spectra to atmospheric models. To quantify the RV precision that can be achieved in low-mass stars over the CARMENES wavelength range, we analyze our empirical information on the RV precision from more than 6500 observations. We compare our high-resolution M-dwarf spectra to atmospheric models where we determine the spectroscopic RV information content, $Q$, and signal-to-noise ratio. We find that for all M-type dwarfs, the highest RV precision can be reached in the wavelength range 700--900nm. Observations at longer wavelengths are equally precise only at the very latest spectral types (M8 and M9). We demonstrate that in this spectroscopic range, the large amount of absorption features compensates for the intrinsic faintness of an M7 star. To reach an RV precision of 1ms$^{-1}$ in very low mass M dwarfs at longer wavelengths likely requires the use of a 10m class telescope. For spectral types M6 and earlier, the combination of a red visual and a near-infrared spectrograph is ideal to search for low-mass planets and to distinguish between planets and stellar variability. At a 4m class telescope, an instrument like CARMENES has the potential to push the RV precision well below the typical jitter level of 3-4ms$^{-1}$.

Context: The main goal of the CARMENES survey is to find Earth-mass planets around nearby M-dwarf stars. Seven M-dwarfs included in the CARMENES sample had been observed before with HIRES and HARPS and either were reported to have one short period planetary companion (GJ15A, GJ176, GJ436, GJ536 and GJ1148) or are multiple planetary systems (GJ581 and GJ876). Aims: We aim to report new precise optical radial velocity measurements for these planet hosts and test the overall capabilities of CARMENES. Methods: We combined our CARMENES precise Doppler measurements with those available from HIRES and HARPS and derived new orbital parameters for the systems. Bona-fide single planet systems are fitted with a Keplerian model. The multiple planet systems were analyzed using a self-consistent dynamical model and their best fit orbits were tested for long-term stability. Results: We confirm or provide supportive arguments for planets around all the investigated stars except for GJ15A, for which we find that the post-discovery HIRES data and our CARMENES data do not show a signal at 11.4 days. Although we cannot confirm the super-Earth planet GJ15Ab, we show evidence for a possible long-period ($P_{\rm c}$ = 7025$_{-629}^{+972}$ d) Saturn-mass ($m_{\rm c} \sin i$ = 51.8$_{-5.8}^{+5.5}M_\oplus$) planet around GJ15A. In addition, based on our CARMENES and HIRES data we discover a second planet around GJ1148, for which we estimate a period $P_{\rm c}$ = 532.6$_{-2.5}^{+4.1}$ d, eccentricity $e_{\rm c}$ = 0.34$_{-0.06}^{+0.05}$ and minimum mass $m_{\rm c} \sin i$ = 68.1$_{-2.2}^{+4.9}M_\oplus$. Conclusions: The CARMENES optical radial velocities have similar precision and overall scatter when compared to the Doppler measurements conducted with HARPS and HIRES. We conclude that CARMENES is an instrument that is up to the challenge of discovering rocky planets around low-mass stars.

We report on the first star discovered to host a planet detected by radial velocity (RV) observations obtained within the CARMENES survey for exoplanets around M dwarfs. HD 147379 ($V = 8.9$ mag, $M = 0.58 \pm 0.08$ M$_{\odot}$), a bright M0.0V star at a distance of 10.7 pc, is found to undergo periodic RV variations with a semi-amplitude of $K = 5.1\pm0.4$ m s$^{-1}$ and a period of $P = 86.54\pm0.06$ d. The RV signal is found in our CARMENES data, which were taken between 2016 and 2017, and is supported by HIRES/Keck observations that were obtained since 2000. The RV variations are interpreted as resulting from a planet of minimum mass $m_{\rm p}\sin{i} = 25 \pm 2$ M$_{\oplus}$, 1.5 times the mass of Neptune, with an orbital semi-major axis $a = 0.32$ au and low eccentricity ($e < 0.13$). HD 147379b is orbiting inside the temperate zone around the star, where water could exist in liquid form. The RV time-series and various spectroscopic indicators show additional hints of variations at an approximate period of 21.1d (and its first harmonic), which we attribute to the rotation period of the star.

Most sub-mm emission line studies of galaxies to date have targeted sources with known redshifts where the frequencies of the lines are well constrained. Recent blind line scans circumvent the spectroscopic redshift requirement, which could represent a selection bias. Our aim is to detect emission lines present in continuum oriented observations. The detection of such lines provides spectroscopic redshift and yields properties of the galaxies. We perform a search for emission lines in the ALMA observations of five Frontier Fields clusters and assess the reliability of our detection by associating line candidates with detected galaxies in deep near-infrared imaging. We find 26 significant emission lines candidates, with observed line fluxes between 0.2-4.6 Jy km s$^{-1}$ and velocity dispersions (FWHM) of 25-600 km s$^{-1}$. Nine of these candidates lie nearby to near-infrared sources, boosting their reliability; in six cases the observed line frequency and strength are consistent with expectations given the photometric redshift and properties of the galaxy counterparts. We present redshift identifications, magnifications and molecular gas estimates for the galaxies with identified lines. We show that two of these candidates likely originate from starburst galaxies, one of which is a jellyfish galaxy, while another two are consistent with being main sequence galaxies based in their depletion times. This work highlights the degree to which serendipitous emission lines can be discovered in large mosaic continuum observations when deep ancillary data are available. The low number of high-significance line detections, however, confirms that such surveys are not as optimal as blind line scans. We stress that Monte Carlo simulations should be used to assess the line detections significances, since using the negative noise suffers from stochasticity and incurs larger uncertainties.

[abridged] The Hubble and Spitzer Space Telescope surveys of the Frontier Fields (FF) provide extremely deep images around six massive, strong-lensing clusters of galaxies. The ALMA FF survey aims to cover the same fields at 1.1mm, with maps reaching (unlensed) sensitivities of $<$70$\mu$Jy, in order to explore the properties of background dusty star-forming galaxies. We report on the multi-wavelength photometric analysis of all 12 significantly detected ($>$5$\sigma$) sources in the first three FF clusters observed by ALMA, based on data from Hubble and Spitzer, the VLT and Herschel. We measure the total photometry in all available bands and determine the photometric redshifts and the physical properties of the counterparts via SED-fitting. In particular, we carefully estimate the FIR photometry using 1.1mm priors to limit the misidentification of blended FIR counterparts, which strongly affect some flux estimates in previous FIR catalogs. We identify robust near-infrared (NIR) counterparts for all 11 sources with K$_s$ detection, the majority of which are quite red, with eight having $F814W-K_s\gtrsim 4$ and five having $F160W-[4.5]\gtrsim3$. From the FIR point of view, all our objects have $z_{phot}$$\sim$1--3, whereas based on the optical SED one object prefers a high-$z$ solution ($z\geq\ $7). Five objects among our sample have spectroscopic redshifts from the GLASS survey for which we can reproduce their SEDs with existing templates. This verification confirms the validity of our photometric redshift methodology. The mean redshift of our sample is $z_{phot}$=1.99$\pm$0.27. All 1.1mm selected objects are massive (10.0$<\log[M_{\star}(M_{\odot})]<$ 11.5), with high star formation rates ($<\log[SFR(M_{\odot}/yr)]> \approx$1.6) and high dust contents (8.1 $<\log[M_{dust} (M_{\odot})]<$8.8), consistent with previous ALMA surveys.

We present an extensive spectroscopic follow-up campaign of 29 strong lensing (SL) selected galaxy clusters discovered primarily in the Second Red-Sequence Cluster Survey (RCS-2). Our spectroscopic analysis yields redshifts for 52 gravitational arcs present in the core of our galaxy clusters, which correspond to 35 distinct background sources that are clearly distorted by the gravitational potential of these clusters. These lensed galaxies span a wide redshift range of $0.8 \le z \le 2.9$, with a median redshift of $z_s = 1.8 \pm 0.1 $. We also measure reliable redshifts for 1004 cluster members, allowing us to obtain robust velocity dispersion measurements for 23 of these clusters, which we then use to determine their dynamical masses by using a simulation-based $\sigma_{DM} - M_{200}$ scaling relation. The redshift and mass ranges covered by our SL sample are $0.22 \le z \le 1.01$ and $5 \times10^{13} \le M_{200}/h^{-1}_{70}M_{\odot} \le 1.9\times10^{15}$, respectively. We analyze and quantify some possible effects that might bias our mass estimates, such as the presence of substructure, the region where cluster members are selected for spectroscopic follow-up, the final number of confirmed members, and line-of-sight effects. We find that 10 clusters of our sample with $N_{mem} \gtrsim 20$ show signs of dynamical substructure. However, the velocity data of only one system is inconsistent with a uni-modal distribution. We therefore assume that the substructures are only marginal and not of comparable size to the clusters themselves. Consequently, our velocity dispersion and mass estimates can be used as priors for SL mass reconstruction studies and also represent an important step toward a better understanding of the properties of the SL galaxy cluster population.

Dusty star-forming galaxies are among the most prodigious systems at high redshift (z>1), characterized by high star formation rates and huge dust reservoirs. The bright end of this population has been well characterized in recent years, but considerable uncertainties remain for fainter dusty star-forming galaxies, which are responsible for the bulk of star formation at high redshift and thus play a key role in galaxy growth and evolution. In this first paper of our series, we describe our methods for finding high redshift faint dusty galaxies using millimeter observations with ALMA. We obtained ALMA 1.1 mm mosaic images for three strong-lensing galaxy clusters from the Frontier Fields survey. The 2'x2' mosaics overlap with the deep HST WFC3/IR footprints and encompass the high magnification regions of each cluster. The combination of extremely high ALMA sensitivity and the magnification power of these clusters allows us to systematically probe the sub-mJy population of dusty star-forming galaxies over a large surveyed area. We present a description of the reduction and analysis of the ALMA continuum observations for the galaxy clusters Abell 2744 (z=0.308), MACSJ0416.1-2403 (z=0.396) and MACSJ1149.5+2223 (z=0.543), for which we reach observed rms sensitivities of 55, 59 and 71 $\mu$Jy beam$^{-1}$ respectively. We detect 12 dusty star-forming galaxies at S/N>5.0 across the three clusters, all of them presenting coincidence with NIR detected counterparts in the HST images. None of the sources fall close to the lensing caustics, thus they are not strongly lensed. The observed 1.1 mm flux densities for the total sample of galaxies range from 0.41 to 2.82 mJy, with observed effective radii spanning <0.05" to 0.37"$\pm$0.21". The lensing-corrected sizes of the detected sources appear to be in the same range as those measured in brighter samples, albeit with possibly larger dispersion.

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 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.

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 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.

The recent discovery of a large number of galaxy clusters using the Sunyaev-Zel'dovich (SZ) effect has opened a new era on the study of the most massive clusters in the Universe. Multi-wavelength analyses are required to understand the properties of these new sets of clusters, which are a sensitive probe of cosmology. We aim at a multi-wavelength characterization of PLCK G004.5$-$19.5, one of the most massive X-ray validated SZ effect-selected galaxy clusters discovered by the Planck satellite. We have observed PLCK G004.5$-$19.5 with GMOS on Gemini South for optical imaging and spectroscopy, and performed a strong lensing analysis. We also searched for associated radio emission in published catalogs. An analysis of the optical images confirms that this is a massive cluster, with a dominant central galaxy (the BCG) and an accompanying red sequence of galaxies, plus a $14''$-long strong lensing arc. Longslit pectroscopy of 6 cluster members shows that the cluster is at $z=0.516\pm0.002$. We also targeted the strongly lensed arc, and found $z_{\rm arc}=1.601$. We use LensTool to carry out a strong lensing analysis, from which we measure a median Einstein radius $\theta_E(z_s=1.6)\simeq30''$ and estimate an enclosed mass $M_E=2.45_{-0.47}^{+0.45}\times10^{14}\,M_\odot$. By extrapolating an NFW profile we find a total mass $M_{500}^{SL}=4.0_{-1.0}^{+2.1}\times10^{14}M_\odot$. Including a constraint on the mass from previous X-ray observations yields a slightly higher mass, $M_{500}^{SL+X}=6.7_{-1.3}^{+2.6}\times10^{14}M_\odot$, marginally consistent with the value from strong lensing alone. High-resolution radio images from TGSS at 150~MHz reveal that PLCK G004.5$-$19.5 hosts a powerful radio relic on scales $\lesssim500$ kpc. Emission at the same location is also detected in low resolution images at 843~MHz and 1.4~GHz. This is one of the higher redshift radio relics known to date.

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.

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 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.

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 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%.

We present the discovery of an extremely bright and extended lensed source from the second Red Sequence Cluster Survey (RCS2). RCSGA 032727-132609 is spectroscopically confirmed as a giant arc and counter-image of a background galaxy at $z=1.701$, strongly-lensed by the foreground galaxy cluster RCS2 032727-132623 at $z=0.564$. The giant arc extends over $\sim 38$\,\arcsec and has an integrated $g$-band magnitude of 19.15, making it $\sim 20$ times larger and $\sim 4$ times brighter than the prototypical lensed galaxy MS1512-cB58. This is the brightest distant lensed galaxy in the Universe known to date. Its location in the `redshift desert' provides unique opportunities to connect between the large samples of galaxies known at $z\sim3$ and $z\sim1$. We have collected photometry in 9 bands, ranging from $u$ to $K_s$, which densely sample the rest-frame UV and optical light, including the age-sensitive 4000\AA\ break. A lens model is constructed for the system, and results in a robust total magnification of $2.04 \pm 0.16$ for the counter-image; we estimate an average magnification of $17.2 \pm 1.4$ for the giant arc based on the relative physical scales of the arc and counter-image. Fits of single-component spectral energy distribution (SED) models to the photometry result in a moderately young age, $t = 115 \pm 65$\,Myr, small amounts of dust, $E(B-V) \le 0.035$, and an exponentially declining star formation history with \textit{e}-folding time $\tau = 10-100$\,Myr. After correcting for the lensing magnification, we find a stellar mass of $\log(\mathrm{M}/\mathrm{M}_\odot)=10.0 \pm 0.1$. Allowing for episodic star formation, an underlying old burst could contain up to twice the mass inferred from single-component modeling. This stellar mass estimate is consistent with the average stellar mass of a sample of `BM' galaxies ($1.4 < z < 2.0$) studied by Reddy et al. (2006).

CARMENES, Calar Alto high-Resolution search for M dwarfs with Exo-earths with a Near-infrared Echelle Spectrograph, is a study for a next-generation instrument for the 3.5m Calar Alto Telescope to be designed, built, integrated, and operated by a consortium of nine German and Spanish institutions. Our main objective is finding habitable exoplanets around M dwarfs, which will be achieved by radial velocity measurements on the m/s level in the near-infrared, where low-mass stars emit the bulk of their radiation.