For studies of galaxy formation and evolution, one of the major benefits of the James Webb Space Telescope is that space-based IFUs like those on its NIRSpec and MIRI instruments will enable spatially resolved spectroscopy of distant galaxies, including spectroscopy at the scale of individual star-forming regions in galaxies that have been gravitationally lensed. In the meantime, there is only a very small subset of lensed sources where work like this is possible even with the Hubble Space Telescope's Wide Field Camera 3 infrared channel grisms. We examine two of these sources, SDSS J1723+3411 and SDSS J2340+2947, using HST WFC3/IR grism data and supporting spatially-unresolved spectroscopy from several ground-based instruments to explore the size of spatial variations in observed strong emission line ratios like O32, R23, which are sensitive to ionization parameter and metallicity, and the Balmer decrement as an indicator of reddening. We find significant spatial variation in the reddening and the reddening-corrected O32 and R23 values which correspond to spreads of a few tenths of a dex in ionization parameter and metallicity. We also find clear evidence of a negative radial gradient in star formation in SDSS J2340+2947 and tentative evidence of one in SDSS J1723+3411, though its star formation is quite asymmetric. Finally, we find that reddening can vary enough spatially to make spatially-resolved reddening corrections necessary in order to characterize gradients in line ratios and the physical conditions inferred from them, necessitating the use of space-based IFUs for future work on larger, more statistically robust samples.

We study the highly magnified arc SGAS J122651.3+215220 caused by a star-forming galaxy at $z_s=2.93$ crossing the lensing caustic cast by the galaxy cluster SDSS J1226+2152 ($z_l=0.43$), using Hubble Space Telescope observations. We report in the arc several asymmetric surface brightness features whose angular separations are a fraction of an arcsecond from the lensing critical curve and appear to be highly but unequally magnified image pairs of underlying compact sources, with one brightest pair having clear asymmetry consistently across four filters. One explanation of unequal magnification is microlensing by intracluster stars, which induces independent flux variations in the images of individual or groups of source stars in the lensed galaxy. For a second possibility, intracluster dark matter subhalos invisible to telescopes effectively perturb lensing magnifications near the critical curve and give rise to persistently unequal image pairs. Our modeling suggests, at least for the most prominent identified image pair, that the microlensing hypothesis is in tension with the absence of notable asymmetry variation over a six-year baseline, while subhalos of $\sim 10^6$--$10^8\,M_\odot$ anticipated from structure formation with Cold Dark Matter typically produce stationary and sizable asymmetries. We judge that observations at additional times and more precise lens models are necessary to stringently constrain temporal variability and robustly distinguish between the two explanations. The arc under this study is a scheduled target of a Director's Discretionary Early Release Science program of the James Webb Space Telescope, which will provide deep images and a high-resolution view with integral field spectroscopy.

We present spatially-resolved echelle spectroscopy of an intervening MgII-FeII-MgI absorption-line system detected at $z_{\rm abs}=0.73379$ toward the giant gravitational arc PSZ1 G311.65-18.48. The absorbing gas is associated to an inclined disk-like star-forming galaxy, whose major axis is aligned with the two arc-segments reported here. We probe in absorption the galaxy's extended disk continuously, at $\approx 3$ kpc sampling, from its inner region out to $15\times$ the optical radius. We detect strong ($W_0^{2796}>0.3$ \r{A}) coherent absorption along $13$ independent positions at impact parameters $D=0$--$29$ kpc on one side of the galaxy, and no absorption at $D=28$--$57$ kpc on the opposite side (all de-lensed distances at $z_{\rm abs}$). We show that: (1) the gas distribution is anisotropic; (2) $W_0^{2796}$, $W_0^{2600}$, $W_0^{2852}$, and the ratio $W_0^{2600}\!/W_0^{2796}$, all anti-correlate with $D$; (3) the $W_0^{2796}$-$D$ relation is not cuspy and exhibits significantly less scatter than the quasar-absorber statistics; (4) the absorbing gas is co-rotating with the galaxy out to $D \lesssim 20$ kpc, resembling a `flat' rotation curve, but at $D\gtrsim 20$ kpc velocities decline below the expectations from a 3D disk-model extrapolated from the nebular [OII] emission. These signatures constitute unambiguous evidence for rotating extra-planar diffuse gas, possibly also undergoing enriched accretion at its edge. Arguably, we are witnessing some of the long-sought processes of the baryon cycle in a single distant galaxy expected to be representative of such phenomena.

During the epoch of reionisation, neutral gas in the early Universe was ionized by hard ultraviolet radiation emitted by young stars in the first galaxies. To do so, ionizing ultraviolet photons must escape from the host galaxy. We present Hubble Space Telescope observations of the gravitationally lensed galaxy PSZ1-ARC G311.6602-18.4624, revealing bright, multiply-imaged ionizing photon escape from a compact star-forming region through a narrow channel in an optically thick gas. The gravitational lensing magnification shows how ionizing photons escape this galaxy, contributing to the re-ionization of the Universe. The multiple sight lines to the source probe absorption by intergalactic neutral hydrogen on scales of no more than a few hundred, perhaps even less than ten, parsec.

In the past decade, our understanding of how stars and galaxies formed during the first 5 billion years after the Big Bang has been revolutionized by observations that leverage gravitational lensing by intervening masses, which act as natural cosmic telescopes to magnify background sources. Previous studies have harnessed this effect to probe the distant universe at ultraviolet, optical, infrared and millimeter wavelengths. However, strong lensing studies of young, star-forming galaxies have never extended into X-ray wavelengths, which uniquely trace high-energy phenomena. Here we report an X-ray detection of star formation in a highly magnified, strongly lensed galaxy. This lensed galaxy, seen during the first third of the history of the Universe, is a low--mass, low--metallicity starburst with elevated X-ray emission, and is a likely analog to the first generation of galaxies. Our measurements yield insight into the role that X-ray emission from stellar populations in the first generation of galaxies may play in re-ionizing the Universe. This observation paves the way for future strong lensing-assisted X-ray studies of distant galaxies reaching orders of magnitude below the detection limits of current deep fields, and previews the depths that will be attainable with future X-ray observatories.

We examine the diagnostic power of rest-frame ultraviolet (UV) nebular emission lines, and compare them to more commonly used rest-frame optical emission lines, using the test case of a single star-forming knot of the bright lensed galaxy RCSGA 032727-132609 at redshift z~1.7. This galaxy has complete coverage of all the major rest-frame UV and optical emission lines from Magellan/MagE and Keck/NIRSPEC. Using the full suite of diagnostic lines, we infer the physical properties: nebular electron temperature (T_e), electron density (n_e), oxygen abundance (log(O/H)), ionisation parameter (log(q)) and interstellar medium (ISM) pressure (log(P/k)). We examine the effectiveness of the different UV, optical and joint UV-optical spectra in constraining the physical conditions. Using UV lines alone we can reliably estimate log(q), but the same is difficult for log(O/H). UV lines yield a higher (~1.5 dex) log(P/k) than the optical lines, as the former probes a further inner nebular region than the latter. For this comparison, we extend the existing Bayesian inference code IZI, adding to it the capability to infer ISM pressure simultaneously with metallicity and ionisation parameter. This work anticipates future rest-frame UV spectral datasets from the James Webb Space Telescope (JWST) at high redshift and from the Extremely Large Telescope (ELT) at moderate redshift.

In response to the need for the Astro2020 Decadal Survey to explicitly engage early career astronomers, the National Academies of Sciences, Engineering, and Medicine hosted the Early Career Astronomer and Astrophysicist Focus Session (ECFS) on October 8-9, 2018 under the auspices of Committee of Astronomy and Astrophysics. The meeting was attended by fifty six pre-tenure faculty, research scientists, postdoctoral scholars, and senior graduate students, as well as eight former decadal survey committee members, who acted as facilitators. The event was designed to educate early career astronomers about the decadal survey process, to solicit their feedback on the role that early career astronomers should play in Astro2020, and to provide a forum for the discussion of a wide range of topics regarding the astrophysics career path. This white paper presents highlights and themes that emerged during two days of discussion. In Section 1, we discuss concerns that emerged regarding the coming decade and the astrophysics career path, as well as specific recommendations from participants regarding how to address them. We have organized these concerns and suggestions into five broad themes. These include (sequentially): (1) adequately training astronomers in the statistical and computational techniques necessary in an era of "big data", (2) responses to the growth of collaborations and telescopes, (3) concerns about the adequacy of graduate and postdoctoral training, (4) the need for improvements in equity and inclusion in astronomy, and (5) smoothing and facilitating transitions between early career stages. Section 2 is focused on ideas regarding the decadal survey itself, including: incorporating early career voices, ensuring diverse input from a variety of stakeholders, and successfully and broadly disseminating the results of the survey.

The ionizing photons produced by massive stars are key actors in galaxy evolution. Ionizing photon production and escape is poorly understood. Improved space-based, spatially-resolved, multiplexed spectroscopic capabilities covering observed wavelengths of 1000 to 3000 Angstrom, in concert with spectroscopy from the ELTs and JWST, would lead to definitive answers as to how ionizing photons are produced and leaked, what populations of galaxies are responsible for ionizing photon leakage, what determines whether escape is possible, and how ionizing galaxy populations evolve over cosmic time.

Cosmological cluster-scale strong gravitational lensing probes the mass distribution of the dense cores of massive dark matter halos and the structures along the line of sight from background sources to the observer. It is frequently assumed that the primary lens mass dominates the lensing, with the contribution of secondary masses along the line of sight being neglected. Secondary mass structures may, however, affect both the detectability of strong lensing in a given survey and modify the properties of the lensing that is detected. In this paper, we utilize a large cosmological N-body simulation and a multiple lens plane (and many source planes) ray-tracing technique to quantify the influence of line of sight halos on the detectability of cluster-scale strong lensing in a cluster sample with a mass limit that encompasses current cluster catalogs from the South Pole Telescope. We extract both primary and secondary halos from the "Outer Rim" simulation and consider two strong lensing realizations: one with only the primary halos included, and the other contains all secondary halos down to a mass limit. In both cases, we use the same source information extracted from the Hubble Ultra Deep Field, and create realistic lensed images consistent with moderately deep ground-based imaging. The results demonstrate that down to the mass limit considered the total number of lenses is boosted by about 13-21% when considering the complete multi-halo lightcone. The increment in strong lens counts peaks at lens redshifts of 0.6 approximately with no significant effect at z<0.3. The strongest trends are observed relative to the primary halo mass, with no significant impact in the most massive quintile of the halo sample, but increasingly boosting the observed lens counts toward small primary halo masses, with an enhancement greater than 50% in the least massive quintile of the halo masses considered.

We present new, deep observations of the Phoenix cluster from the Chandra X-ray Observatory, the Hubble Space Telescope, and the Karl Jansky Very Large Array. These data provide an order of magnitude improvement in depth and/or angular resolution at X-ray, optical, and radio wavelengths, yielding an unprecedented view of the core of the Phoenix cluster. We find that the one-dimensional temperature and entropy profiles are consistent with expectations for pure-cooling hydrodynamic simulations and analytic descriptions of homogeneous, steady-state cooling flow models. In the inner ~10 kpc, the cooling time is shorter by an order of magnitude than any other known cluster, while the ratio of the cooling time to freefall time approaches unity, signaling that the ICM is unable to resist multiphase condensation on kpc scales. When we consider the thermodynamic profiles in two dimensions, we find that the cooling is highly asymmetric. The bulk of the cooling in the inner ~20 kpc is confined to a low-entropy filament extending northward from the central galaxy. We detect a substantial reservoir of cool (10^4 K) gas (as traced by the [OII] doublet), which is coincident with the low-entropy filament. The bulk of this cool gas is draped around and behind a pair of X-ray cavities, presumably bubbles that have been inflated by radio jets, which are detected for the first time on kpc scales. These data support a picture in which AGN feedback is promoting the formation of a multiphase medium via a combination of ordered buoyant uplift and locally enhanced turbulence. These processes ought to counteract the tendency for buoyancy to suppress condensation, leading to rapid cooling along the jet axis. The recent mechanical outburst has sufficient energy to offset cooling, and appears to be coupling to the ICM via a cocoon shock, raising the entropy in the direction orthogonal to the radio jets.

We present strong gravitational lensing models for 37 galaxy clusters from the SDSS Giant Arcs Survey. We combine data from multi-band Hubble Space Telescope WFC3imaging, with ground-based imaging and spectroscopy from Magellan, Gemini, APO, and MMT, in order to detect and spectroscopically confirm new multiply-lensed background sources behind the clusters. We report spectroscopic or photometric redshifts of sources in these fields, including cluster galaxies and background sources. Based on all available lensing evidence, we construct and present strong lensing mass models for these galaxy clusters.

We present the first spatially resolved analysis of rest-frame optical and UV imaging and spectroscopy for a lensed galaxy at z = 2.39 hosting a Seyfert active galactic nucleus (AGN). Proximity to a natural guide star has enabled high signal-to-noise VLT SINFONI + adaptive optics observations of rest-frame optical diagnostic emission lines, which exhibit an underlying broad component with FWHM ~ 700 km/s in both the Balmer and forbidden lines. Measured line ratios place the outflow robustly in the region of the ionization diagnostic diagrams associated with AGN. This unique opportunity - combining gravitational lensing, AO guiding, redshift, and AGN activity - allows for a magnified view of two main tracers of the physical conditions and structure of the interstellar medium in a star-forming galaxy hosting a weak AGN at cosmic noon. By analyzing the spatial extent and morphology of the Ly-alpha and dust-corrected H-alpha emission, disentangling the effects of star formation and AGN ionization on each tracer, and comparing the AGN induced mass outflow rate to the host star formation rate, we find that the AGN does not significantly impact the star formation within its host galaxy.

We present $Hubble\ Space\ Telescope$ ($HST$) imaging and grism spectroscopy of a strongly lensed LIRG at $z=0.816$, SGAS 143845.1$+$145407, and use the magnification boost of gravitational lensing to study the distribution of star formation throughout this galaxy. Based on the $HST$ imaging data, we create a lens model for this system; we compute the mass distribution and magnification map of the $z=0.237$ foreground lens. We find that the magnification of the lensed galaxy ranges between $2$ and $10$, with a total magnification (measured over all the images of the source) of $\mu=11.8^{+4.6}_{-2.4}$. We find that the total projected mass density within $\sim34$ kpc of the brightest cluster galaxy is $6.0^{+0.3}_{-0.7}\times10^{12}\,M_{\odot}$. Using the lens model we create a source reconstruction for SGAS 143845.1$+$145407, which paired with a faint detection of H$\alpha$ in the grism spectroscopy, allows us to finally comment directly on the distribution of star formation in a $z\sim1$ LIRG. We find widespread star formation across this galaxy, in agreement with the current understanding of these objects. However, we note a deficit of H$\alpha$ emission in the nucleus of SGAS 143845.1$+$145407, likely due to dust extinction.

We report the detection of extended Lyman-$\alpha$ emission from the host galaxy of SDSS~J2222+2745, a strongly lensed quasar at $z = 2.8$. Spectroscopic follow-up clearly reveals extended Lyman-$\alpha$ in emission between two images of the central active galactic nucleus (AGN). We reconstruct the lensed quasar host galaxy in the source plane by applying a strong lens model to HST imaging, and resolve spatial scales as small as $\sim$200 parsecs. In the source plane we recover the host galaxy morphology to within a few hundred parsecs of the central AGN, and map the extended Lyman-$\alpha$ emission to its physical origin on one side of the host galaxy at radii $\sim$0.5-2 kpc from the central AGN. There are clear morphological differences between the Lyman-$\alpha$ and rest-frame ultraviolet stellar continuum emission from the quasar host galaxy. Furthermore, the relative velocity profiles of quasar Lyman-$\alpha$, host galaxy Lyman-$\alpha$, and metal lines in outflowing gas reveal differences in the absorbing material affecting the AGN and host galaxy. These data indicate the presence of patchy local intervening gas in front of the central quasar and its host galaxy. This interpretation is consistent with the central luminous quasar being obscured across a substantial fraction of its surrounding solid angle, resulting in strong anisotropy in the exposure of the host galaxy to ionizing radiation from the AGN. This work demonstrates the power of strong lensing-assisted studies to probe spatial scales that are currently inaccessible by other means.

For lensed galaxy SGAS J111020.0+645950.8 at redshift z=2.481, which is magnified by a factor of 28 +- 8, we analyze the morphology of star formation as traced by rest-frame ultraviolet emission, in both the highly-magnified source plane, and in simulations of how this galaxy would appear without lensing magnification. Were this galaxy not lensed but drawn from an HST deep field, we would conclude that almost all its star formation arises from an exponential disk (S\'ersic index of 1.0 +- 0.4) with an effective radius of r_e = 2.7 +- 0.3 kpc measured from two-dimensional fitting to F606W using Galfit, and r_e=1.9 +- 0.1 kpc measured by fitting a radial profile to F606W elliptical isophotes. At the normal spatial resolution of the deep fields, there is no sign of clumpy star formation within SGAS J111020.0+645950.8 . However, the enhanced spatial resolution enabled by gravitational lensing tells a very different story: much of the star formation arises in two dozen clumps with sizes of r=30--50 pc spread across the 7 kpc length of the galaxy. The color and spatial distribution of the diffuse component suggests that still smaller clumps are unresolved. Despite this clumpy, messy morphology, the radial profile is still well-characterized by an exponential profile. In this lensed galaxy, stars are forming in complexes with sizes well below 100 pc; such sizes are wholly unexplored by surveys of galaxy evolution at 1<z<3.

We present measurements of the surface density of star formation, the star-forming clump luminosity function, and the clump size distribution function, for the lensed galaxy SGAS J111020.0+645950.8 at a redshift of z=2.481. The physical size scales that we probe, radii r=30--50 pc, are considerably smaller scales than have yet been studied at these redshifts. The star formation surface density we find within these small clumps is consistent with surface densities measured previously for other lensed galaxies at similar redshift. Twenty-two percent of the rest-frame ultraviolet light in this lensed galaxy arises from small clumps, with r<100 pc. Within the range of overlap, the clump luminosity function measured for this lensed galaxy is remarkably similar to those of z~0 galaxies. In this galaxy, star-forming regions smaller than 100 pc---physical scales not usually resolved at these redshifts by current telescopes---are important locations of star formation in the distant universe. If this galaxy is representative, this may contradict the theoretical picture in which the critical size scale for star formation in the distant universe is of order 1 kiloparsec. Instead, our results suggest that current telescopes have not yet resolved the critical size scales of star-forming activity in galaxies over most of cosmic time.

Using the combined resolving power of the Hubble Space Telescope and gravitational lensing, we resolve star-forming structures in a z~2.5 galaxy on scales much smaller than the usual kiloparsec diffraction limit of HST. SGAS J111020.0+645950.8 is a clumpy, star forming galaxy lensed by the galaxy cluster SDSS J1110+6459 at z = 0.659, with a total magnification ~30x across the entire arc. We use a hybrid parametric/non-parametric strong lensing mass model to compute the deflection and magnification of this giant arc, reconstruct the light distribution of the lensed galaxy in the source plane, and resolve the star formation into two dozen clumps. We develop a forward-modeling technique to model each clump in the source plane. We ray trace the model to the image plane, convolve with the instrumental point spread function (PSF), and compare with the GALFIT model of the clumps in the image plane, which decomposes clump structure from more extended emission. This technique has the advantage, over ray tracing, by accounting for the asymmetric lensing shear of the galaxy in the image plane and the instrument PSF. At this resolution, we can begin to study star formation on a clump-by-clump basis, toward the goal of understanding feedback mechanisms and the buildup of exponential disks at high redshift.

Gravitational lensing has become one of the most powerful tools available for investigating the 'dark side' of the universe. Cosmological strong gravitational lensing, in particular, probes the properties of the dense cores of dark matter halos over decades in mass and offers the opportunity to study the distant universe at flux levels and spatial resolutions otherwise unavailable. Studies of strongly-lensed variable sources offer yet further scientific opportunities. One of the challenges in realizing the potential of strong lensing is to understand the statistical context of both the individual systems that receive extensive follow-up study, as well as that of the larger samples of strong lenses that are now emerging from survey efforts. Motivated by these challenges, we have developed an image-simulation pipeline, PICS (Pipeline for Images of Cosmological Strong lensing) to generate realistic strong gravitational lensing signals from group and cluster scale lenses. PICS uses a low-noise and unbiased density estimator based on (resampled) Delaunay Tessellations to calculate the density field; lensed images are produced by ray-tracing images of actual galaxies from deep Hubble Space Telescope observations. Other galaxies, similarly sampled, are added to fill in the light cone. The pipeline further adds cluster-member galaxies and foreground stars into the lensed images. The entire image ensemble is then observed using a realistic point spread function which includes appropriate detector artifacts for bright stars. Noise is further added, including such non-Gaussian elements as noise window-paning from mosaiced observations, residual bad pixels, and cosmic rays. The aim is to produced simulated images that appear identical - to the eye (expert or otherwise) - to real observations in various imaging surveys.

SDSS J2222+2745 is a galaxy cluster at z=0.49, strongly lensing a quasar at z=2.805 into six widely separated images. In recent HST imaging of the field, we identify additional multiply lensed galaxies, and confirm the sixth quasar image that was identified by Dahle et al. (2013). We used the Gemini North telescope to measure a spectroscopic redshift of z=4.56 of one of the secondary lensed galaxies. These data are used to refine the lens model of SDSS J2222+2745, compute the time delay and magnifications of the lensed quasar images, and reconstruct the source image of the quasar host and a second lensed galaxy at z=2.3. This second galaxy also appears in absorption in our Gemini spectra of the lensed quasar, at a projected distance of 34 kpc. Our model is in agreement with the recent time delay measurements of Dahle et al. (2015), who found tAB=47.7+/-6.0 days and tAC=-722+/-24 days. We use the observed time delays to further constrain the model, and find that the model-predicted time delays of the three faint images of the quasar are tAD=502+/-68 days, tAE=611+/-75 days, and tAF=415+/-72 days. We have initiated a follow-up campaign to measure these time delays with Gemini North. Finally, we present initial results from an X-ray monitoring program with Swift, indicating the presence of hard X-ray emission from the lensed quasar, as well as extended X-ray emission from the cluster itself, which is consistent with the lensing mass measurement and the cluster velocity dispersion.

Characterization of the morphology of strongly lensed galaxies is challenging because images of such galaxies are typically highly distorted. Lens modeling and source plane reconstruction is one approach that can provide reasonably undistorted images from which morphological measurements can be made, although at the expense of a highly spatially variable telescope PSF when mapped back to the source plane. Unfortunately, modeling the lensing mass is a time and resource intensive process, and in many cases there are too few constraints to precisely model the lensing mass. If, however, useful morphological measurements could be made in the image plane rather than the source plane, it would bypass this issue and obviate the need for a source reconstruction process for some applications. We examine the use of the Gini coefficient as one such measurement. Because it depends on the cumulative distribution of the light of a galaxy, but not the relative spatial positions, the fact that surface brightness is conserved by lensing means that the Gini coefficient may be well-preserved by strong gravitational lensing. Through simulations, we test the extent to which the Gini coefficient is conserved, including by effects due to PSF convolution and pixelization, to determine whether it is invariant enough under lensing to be used as a measurement of galaxy morphology that can be made in the image plane.

Using the Low Dispersion Survey Spectrograph 3 at the Magellan II Clay Telescope, we target {candidate absorption host galaxies} detected in deep optical imaging {(reaching limiting apparent magnitudes of 23.0-26.5 in $g, r, i,$ and $z$ filters) in the fields of three QSOs, each of which shows the presence of high metallicity, high $N_{\rm HI}$ absorption systems in their spectra (Q0826-2230: $z_{abs}$=0.9110, Q1323-0021: $z_{abs}=0.7160$, Q1436-0051: $z_{abs}=0.7377, 0.9281$). We confirm three host galaxies {at redshifts 0.7387, 0.7401, and 0.9286} for two of the Lyman-$\alpha$ absorption systems (one with two galaxies interacting). For these systems, we are able to determine the star formation rates (SFRs); impact parameters (from previous imaging detections); the velocity shift between the absorption and emission redshifts; and, for one system, also the emission metallicity.} Based on previous photometry, we find these galaxies have L$>$L$^{\ast}$. The [O II] SFRs for these galaxies are in the range $11-25$ M$_{\odot}$ yr$^{-1}$ {(uncorrected for dust)}, while the impact parameters lie in the range $35-54$ kpc. {Despite the fact that we have confirmed galaxies at 50 kpc from the QSO, no gradient in metallicity is indicated between the absorption metallicity along the QSO line of sight and the emission line metallicity in the galaxies.} We confirm the anti-correlation between impact parameter and $N_{\rm HI}$ from the literature. We also report the emission redshift of five other galaxies: three at $z_{em}>z_{QSO}$, and two (L$<$L$^{\ast}$) at $z_{em}<z_{QSO}$ not corresponding to any known absorption systems.

The sample of cosmological strong lensing systems has been steadily growing in recent years and with the advent of the next generation of space-based survey telescopes, the sample will reach into the thousands. The accuracy of strong lens models relies on robust identification of multiple image families of lensed galaxies. For the most massive lenses, often more than one background galaxy is magnified and multiply-imaged, and even in the cases of only a single lensed source, identification of counter images is not always robust. Recently, we have shown that the Gini coefficient in space-telescope-quality imaging is a measurement of galaxy morphology that is relatively well-preserved by strong gravitational lensing. Here, we investigate its usefulness as a diagnostic for the purposes of image family identification and show that it can remove some of the degeneracies encountered when using color as the sole diagnostic, and can do so without the need for additional observations since whenever a color is available, two Gini coefficients are as well.

We present the mass distribution at the core of SDSS J1531+3414, a strong-lensing cluster at z=0.335. We find that the mass distribution is well described by two cluster-scale halos with a contribution from cluster-member galaxies. New HST observations of SDSS J1531+3414 reveal a signature of ongoing star formation associated with the two central galaxies at the core of the cluster, in the form of a chain of star forming regions at the center of the cluster. Using the lens model presented here, we place upper limits on the contribution of a possible lensed image to the flux at the center region, and rule out that this emission is coming from a background source.

New Hubble Space Telescope ultraviolet and optical imaging of the strong-lensing galaxy cluster SDSS J1531+3414 (z=0.335) reveals two centrally dominant elliptical galaxies participating in an ongoing major merger. The interaction is at least somewhat rich in cool gas, as the merger is associated with a complex network of nineteen massive superclusters of young stars (or small tidal dwarf galaxies) separated by ~1 kpc in projection from one another, combining to an estimated total star formation rate of ~5 solar masses per year. The resolved young stellar superclusters are threaded by narrow H-alpha, [O II], and blue excess filaments arranged in a network spanning ~27 kpc across the two merging galaxies. This morphology is strongly reminiscent of the well-known "beads on a string" mode of star formation observed on kpc-scales in the arms of spiral galaxies, resonance rings, and in tidal tails between interacting galaxies. Nevertheless, the arrangement of this star formation relative to the nuclei of the two galaxies is difficult to interpret in a dynamical sense, as no known "beads on a string" systems associated with kpc-scale tidal interactions exhibit such lopsided morphology relative to the merger participants. In this Letter we present the images and follow-up spectroscopy, and discuss possible physical interpretations for the unique arrangement of the young stellar clusters. While we suggest that this morphology is likely to be dynamically short-lived, a more quantitative understanding awaits necessary multiwavelength follow-up, including optical integral field spectroscopy, ALMA sub-mm interferometry, and Chandra X-ray imaging.

We present optical and near-IR imaging and spectroscopy of SGAS J105039.6$+$001730, a strongly lensed galaxy at z $=$ 3.6252 magnified by $>$30$\times$, and derive its physical properties. We measure a stellar mass of log(M$_{*}$/M$_{\odot}$) $=$ 9.5 $\pm$ 0.35, star formation rates from [O II]$\lambda$$\lambda$3727 and H-$\beta$ of 55 $\pm$ 20 and 84 $\pm$ 17 M$_{\odot}$ yr$^{-1}$, respectively, an electron density of n$_{e} \leq$ 10$^{3}$ cm$^{-2}$, an electron temperature of T$_{e} \leq$ 14000 K, and a metallicity of 12+log(O/H) $=$ 8.3 $\pm$ 0.1. The strong C III]$\lambda$$\lambda$1907,1909 emission and abundance ratios of C, N, O and Si are consistent with well-studied starbursts at z $\sim$ 0 with similar metallicities. Strong P Cygni lines and He II$\lambda$1640 emission indicate a significant population of Wolf-Rayet stars, but synthetic spectra of individual populations of young, hot stars do not reproduce the observed integrated P Cygni absorption features. The rest-frame UV spectral features are indicative of a young starburst with high ionization, implying either 1) an ionization parameter significantly higher than suggest by rest-frame optical nebular lines, or 2) differences in one or both of the initial mass function and the properties of ionizing spectra of massive stars. We argue that the observed features are likely the result of a superposition of star forming regions with different physical properties. These results demonstrate the complexity of star formation on scales smaller than individual galaxies, and highlight the importance of systematic effects that result from smearing together the signatures of individual star forming regions within galaxies.