Identifying non-contaminated sample of high-redshift galaxies with escaping Lyman continuum (LyC) flux is important for understanding the sources and evolution of cosmic reionization. We present CLAUDS $u$-band photometry of the COSMOS field to probe LyC radiation from spectroscopically confirmed galaxies at $z\geq3.5$ and outside the standard Lyman-break galaxy colour selection expectations. Complementary to the CLAUDS data, we use Subaru multi-filter photometry, \textit {Hubble Space Telescope} ($HST$) multi-filter imaging, and the spectroscopic surveys D10K, VUDS and 3D-HST. We present a sample of Lyman continuum galaxy (LCG) candidates in the redshift range $3.5\lesssim z\lesssim5.1$. Here, we introduce 5 LCG candidates, where two are flagged quality 1 and three quality 2. The estimated $f_{\rm esc}^{\rm abs}$ for quality 1 candidates are in the range $\sim5\% - 73\%$ and $\sim30\% - 93\%$. These estimates are based on our derived parameters from individual galaxies as inputs to a range of BPASS models as well as mean intergalactic medium (IGM) and maximal intergalactic and circumgalactic media (IGM+CGM) transmission. We conclude that our search for LCGs is most likely biased to lines of sight with low HI densities or free from Lyman limit systems. Our two best LCG candidates have EW (Ly$\alpha)\leq50$\AA\ and we find no correlation or anti-correlation between EW (Ly$\alpha$), $f_{\rm esc}^{\rm abs}$, and $R_{\rm obs}$, the ratio of ionizing to non-ionizing observed flux in the measured passbands. Stacking candidates without solid LyC detections ($S/N<3$) results in an estimated $f_{\rm esc}^{\rm abs}$ from galaxies not greater than $1\%$.

Many recent integral integral field spectroscopy (IFS) survey teams have used stellar kinematic maps combined with imaging to statistically infer the underlying distributions of galaxy intrinsic shapes. With now several IFS samples at our disposal, the method, which was originally proposed by M. Franx and collaborators in 1991, is gaining in popularity, having been so far applied to ATLAS3D, SAMI, MANGA and MASSIVE. We present results showing that a commonly assumed relationship between dynamical and intrinsic shape alignment does not hold in Illustris, affecting our ability to recover accurate intrinsic shape distributions. A further implication is that so-called "prolate rotation", where the bulk of stars in prolate galaxies are thought to rotate around the projected major axis, is a misnomer.

The distribution of three dimensional intrinsic galaxy shapes has been a longstanding open question. The difficulty stems from projection effects meaning one must rely on statistical methods applied to galaxy samples to infer intrinsic shape distributions. Theoretical work using analytical galaxy potentials suggests a relationship between galaxy intrinsic shape (as defined by its "triaxiality", in practice a proxy for how prolate a galaxy is) and the intrinsic misalignment angle between kinematic and morphological axes ($\Psi_{\rm int}$). This relationship reduces the number of unknowns, providing more reliable inferred intrinsic shape distributions than methods using photometry alone. Here we explore the connection between intrinsic shape and stellar kinematics using cosmological hydrodynamical simulations from the Illustris project. The strongest relationship we find is that galaxy intrinsic flattening is correlated with specific angular momentum (j) with high j galaxies being flatter than galaxies with low specific angular momentum. Our analysis shows that, although the majority of kinematically misaligned galaxies exhibit prolate shapes, examples of kinematically aligned prolate galaxies are also present. Clearly a direct correspondence between prolate shape and minor-axis rotation (often referred to as "prolate rotation") is not present in Illustris. Thus, we demonstrate that the assumption of a simple relationship between $\Psi_{\rm int}$ and intrinsic shape commonly employed in shape recovery studies is not valid for Illustris galaxies. We suggest improvements on the method as well as some alternative methods for future work in this area.

We present the first results of our pilot study of 8 photometrically selected Lyman continuum (LyC) emitting galaxy candidates from the COSMOS field and focus on their optical emission line ratios. Observations were performed in the H and K bands using the Multi-Object Spectrometer for Infra-Red Exploration (MOSFIRE) instrument at the Keck Observatory, targeting the [OII], H$\beta$, and [OIII] emission lines. We find that photometrically selected LyC emitting galaxy candidates have high ionization parameters, based on their high [OIII]/[OII] ratios (O32), with an average ratio for our sample of 2.5$\pm$0.2. Preliminary results of our companion Low Resolution Imaging Spectrometer (LRIS) observations, targeting LyC and Ly$\alpha$, show that those galaxies with the largest O32 are typically found to also be Ly$\alpha$ emitters. High O32 galaxies are also found to have tentative non-zero LyC escape fractions ($f_{esc}(LyC)$) based on $u$ band photometric detections. These results are consistent with samples of highly ionized galaxies, including confirmed LyC emitting galaxies from the literature. We also perform a detailed comparison between the observed emission line ratios and simulated line ratios from density bounded H$_{\textrm{II}}$ regions modeled using the photoionization code MAPPINGS V. Estimates of $f_{esc}(LyC)$ for our sample fall in the range from 0.0-0.23 and suggest possible tension with published correlations between O32 and $f_{esc}(LyC)$, adding weight to dichotomy of arguments in the literature. We highlight the possible effects of clumpy geometry and mergers that may account for such tension.

Stars form in molecular complexes that are visible as giant clouds ($\sim 10^{5-6} \mathrm{M}_\odot$) in nearby galaxies and as giant clumps ($\sim 10^{8-9}\mathrm{M}_\odot$) in galaxies at redshifts $z\approx1$$-$$3$. Theoretical inferences on the origin and evolution of these complexes often require robust measurements of their characteristic size, which is hard to measure at limited resolution and often ill-defined due to overlap and quasi-fractal substructure. We show that maximum and luminosity-weighted sizes of clumps seen in star formation maps (e.g.\ H$\alpha$) can be recovered statistically using the two-point correlation function (2PCF), if an approximate stellar surface density map is taken as the normalizing random field. After clarifying the link between Gaussian clumps and the 2PCF analytically, we design a method for measuring the diameters of Gaussian clumps with realistic quasi-fractal substructure. This method is tested using mock images of clumpy disk galaxies at different spatial resolutions and perturbed by Gaussian white noise. We find that the 2PCF can recover the input clump scale at $\sim20\%$ accuracy, as long as this scale is larger than the spatial resolution. We apply this method to the local spiral galaxy NGC 5194, as well as to three clumpy turbulent galaxies from the DYNAMO-HST sample. In both cases, our statistical H$\alpha$-clump size measurements agree with previous measurements and with the estimated Jeans lengths. However, the new measurements are free from subjective choices when fitting individual clumps.

We present Keck/OSIRIS adaptive optics observations with 150-400 pc spatial sampling of 7 turbulent, clumpy disc galaxies from the DYNAMO sample ($0.07<z<0.2$). DYNAMO galaxies have previously been shown to be well matched in properties to main sequence galaxies at $z\sim1.5$. Integral field spectroscopy observations using adaptive optics are subject to a number of systematics including a variable PSF and spatial sampling, which we account for in our analysis. We present gas velocity dispersion maps corrected for these effects, and confirm that DYNAMO galaxies do have high gas velocity dispersion ($\sigma=40-80$\kms), even at high spatial sampling. We find statistically significant structure in 6 out of 7 galaxies. The most common distance between the peaks in velocity dispersion and emission line peaks is $\sim0.5$~kpc, we note this is very similar to the average size of a clump measured with HST H$\alpha$ maps. This could suggest that the peaks in velocity dispersion in clumpy galaxies likely arise due to some interaction between the clump and the surrounding ISM of the galaxy, though our observations cannot distinguish between outflows, inflows or velocity shear. Observations covering a wider area of the galaxies will be needed to confirm this result.

We examine the kinematic morphology of early-type galaxies (ETGs) in eight galaxy clusters in the Sydney-AAO Multi-object Integral field spectrograph (SAMI) Galaxy Survey. The clusters cover a mass range of 14.2<log(M_200/M_odot)<15.2 and we measure spatially-resolved stellar kinematics for 315 member galaxies with stellar masses 10.0<log(M_*/M_odot)<11.7 within 1R_200 of the cluster centers. We calculate the spin parameter, lambda_R and use that to classify the kinematic morphology of the galaxies as fast or slow rotators. The total fraction of slow rotators in the early-type galaxy population, F_SR=0.14+/-0.02 and does not depend on host cluster mass. Across the eight clusters, the fraction of slow rotators increases with increasing local overdensity. We also find that the slow-rotator fraction increases at small clustercentric radii (R_cl<0.3R_200), and note that there is also an increase in slow-rotator fraction at R_cl~0.6R_200. The slow rotators at these larger radii reside in cluster substructure. We find the strongest increase in slow-rotator fraction occurs with increasing stellar mass. After accounting for the strong correlation with stellar mass, we find no significant relationship between spin parameter and local overdensity in the cluster environment. We conclude that the primary driver for the kinematic morphology--density relationship in galaxy clusters is the changing distribution of galaxy stellar mass with local environment. The presence of slow rotators in substructure suggests that the cluster kinematic morphology--density relationship is a result of mass segregation of slow-rotating galaxies forming in groups that later merge with clusters and sink to the cluster center via dynamical friction.

The observation of counter rotation in galaxies (i.e. gas that rotates in the opposite direction to the stellar component or two co-spatial stellar populations with opposite rotation) is becoming more commonplace with modern integral field spectroscopic surveys. In this paper we explore the emergence of counter-rotation (both stellar and gaseous) in S0 galaxies from smoothed-particle hydrodynamics simulations of 1/10 mass ratio minor mergers between a 10^10.8 solar mass disk galaxy with a bulge-to-total ration of 0.17 and a gas rich companion (gas-to-stellar mass fraction of 5.0). These simulations include a self-consistent treatment of gas dynamics, star formation, the production/destruction of H2 and dust, and the time evolution of the interstellar radiation field. We explore the effect of retrograde versus prograde obits, gas and bulge mass fractions of the primary galaxy, and orbital parameters of the companion. The key requirement for producing counter rotation in stars or gas in a merger remnant is a retrograde primary, while the relative spin of the companion affects only the radial extent of the accreted gas. We also find that including a significant amount of gas in the primary can prevent the emergence of counter-rotating gas, although accreted stars retain counter-rotation. Bulge mass and orbit have a secondary effect, generally influencing the final distribution of accreted stars and gas within the framework outlined above. In addition to our primary focus of counter-rotating components in galaxies, we also make some predictions regarding the SFRs, H2 distributions, and dust in minor merger remnants.

Recently, large samples of visually classified early-type galaxies (ETGs) containing dust have been identified using space-based infrared observations with the Herschel Space Telescope. The presence of large quantities of dust in massive ETGs is peculiar as X-ray halos of these galaxies are expected to destroy dust in 10 Myr (or less). This has sparked a debate regarding the origin of the dust: is it internally produced by asymptotic giant branch (AGB) stars, or is it accreted externally through mergers? We examine the 2D stellar and ionised gas kinematics of dusty ETGs using IFS observations from the SAMI galaxy survey, and integrated star-formation rates, stellar masses, and dust masses from the GAMA survey. Only 8% (4/49) of visually-classified ETGs are kinematically consistent with being dispersion-supported systems. These "dispersion-dominated galaxies" exhibit discrepancies between stellar and ionised gas kinematics, either offsets in the kinematic position angle or large differences in the rotational velocity, and are outliers in star-formation rate at a fixed dust mass compared to normal star-forming galaxies. These properties are suggestive of recent merger activity. The remaining 90% of dusty ETGs have low velocity dispersions and/or large circular velocities, typical of "rotation-dominated galaxies". These results, along with the general evidence of published works on X-ray emission in ETGs, suggest that they are unlikely to host hot, X-ray gas consistent with their low stellar mass when compared to dispersion-dominated galaxies. This means dust will be long lived and thus these galaxies do not require external scenarios for the origin of their dust content.

Dust attenuation in galaxies has been extensively studied nearby, however, there are still many unknowns regarding attenuation in distant galaxies. We contribute to this effort using observations of star-forming galaxies in the redshift range z = 0.05-0.15 from the DYNAMO survey. Highly star-forming DYNAMO galaxies share many similar attributes to clumpy, star-forming galaxies at high redshift. Considering integrated Sloan Digital Sky Survey observations, trends between attenuation and other galaxy properties for DYNAMO galaxies are well matched to star-forming galaxies at high redshift. Integrated gas attenuations of DYNAMO galaxies are 0.2-2.0 mags in the V-band, and the ratio of stellar E(B-V) and gas E(B-V) is 0.78-0.08 (compared to 0.44 at low redshift). Four highly star-forming DYNAMO galaxies were observed at H-alpha using the Hubble Space Telescope and at Pa-alpha using integral field spectroscopy at Keck. The latter achieve similar resolution (~0.8-1 kpc) to our HST imaging using adaptive optics, providing resolved observations of gas attenuations of these galaxies on sub-kpc scales. We find < 1.0 mag of variation in attenuation (at H-alpha) from clump to clump, with no evidence of highly attenuated star formation. Attenuations are in the range 0.3-2.2 mags in the V band, consistent with attenuations of low redshift star-forming galaxies. The small spatial variation on attenuation suggests that a majority of the star-formation activity in these four galaxies occurs in relatively unobscured regions and, thus, star-formation is well characterised by our H-alpha observations.

We present $\sim$100 pc resolution Hubble Space Telescope H$\alpha$ images of 10 galaxies from the DYnamics of Newly-Assembled Massive Objects (DYNAMO) survey of low-$z$ turbulent disk galaxies, and use these to undertake the first detailed systematic study of the effects of resolution and clump clustering on observations of clumps in turbulent disks. In the DYNAMO-{\em HST} sample we measure clump diameters spanning the range $d_{clump} \sim 100-800$~pc, and individual clump star formation rates as high as $\sim5$~M$_{\odot}$~yr$^{-1}$. DYNAMO clumps have very high SFR surface densities, $\Sigma_{SFR}\sim 15$~M$_{\odot}$~yr$^{-1}$~kpc$^{-2}$, $\sim100\times$ higher than in H{\sc ii} regions of nearby spirals. Indeed, SFR surface density provides a simple dividing line between massive star forming clumps and local star forming regions, where massive star forming clumps have $\Sigma_{SFR}> 0.5$~M$_{\odot}$~yr$^{-1}$~kpc$^{-2}$. When degraded to match the observations of galaxies in $z\sim 1-3$ surveys, DYNAMO galaxies are similar in morphology and measured clump properties to clumpy galaxies observed in the high-$z$ Universe. Emission peaks in the simulated high-redshift maps typically correspond to multiple clumps in full resolution images. This clustering of clumps systematically increases the apparent size and SFR of clumps in 1~kpc resolution maps, and decreases the measured SFR surface density of clumps by as much as a factor of 20$\times$. From these results we can infer that clump clustering is likely to strongly effect the measured properties of clumps in high-$z$ galaxies, which commonly have kiloparsec scale resolution.

We present post-cryogenic Spitzer imaging at 3.6 and 4.5 micron with the Infrared Array Camera (IRAC) of the Spitzer/HETDEX Exploratory Large-Area (SHELA) survey. SHELA covers $\sim$deg$^2$ of the Sloan Digital Sky Survey "Stripe 82" region, and falls within the footprints of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) and the Dark Energy Survey. The HETDEX blind R $\sim$ 800 spectroscopy will produce $\sim$ 200,000 redshifts from the Lyman-$\alpha$ emission for galaxies in the range 1.9 < z < 3.5, and an additional $\sim$200,000 redshifts from the [OII] emission for galaxies at z < 0.5. When combined with deep ugriz images from the Dark Energy Camera, K-band images from NEWFIRM, and other ancillary data, the IRAC photometry from Spitzer will enable a broad range of scientific studies of the relationship between structure formation, galaxy stellar mass, halo mass, AGN, and environment over a co-moving volume of $\sim$0.5 Gpc$^3$ at 1.9 < z < 3.5. Here, we discuss the properties of the SHELA IRAC dataset, including the data acquisition, reduction, validation, and source catalogs. Our tests show the images and catalogs are 80% (50%) complete to limiting magnitudes of 22.0 (22.6) AB mag in the detection image, which is constructed from the weighted sum of the IRAC 3.6 and 4.5 micron images. The catalogs reach limiting sensitivities of 1.1 $\mu$Jy at both 3.6 and 4.5 micron (1$\sigma$, for R=2 arcsec circular apertures). As a demonstration of science, we present IRAC number counts, examples of highly temporally variable sources, and galaxy surface density profiles of rich galaxy clusters. In the spirit of Spitzer Exploratory programs we provide all images and catalogs as part of the publication.

We measure the stellar specific angular momentum jstar=Jstar/Mstar in four nearby (z~0.1) disk galaxies that have stellar masses Mstar near the break M* of the galaxy mass function, but look like typical star-forming disks at z~2 in terms of their low stability (Q~1), clumpiness, high ionized gas dispersion (40-50 km/s), high molecular gas fraction (20-30%) and rapid star formation (~20 Msun/yr). Combining high-resolution (Keck-OSIRIS) and large-radius (Gemini-GMOS) spectroscopic maps, only available at low z, we discover that these targets have about three times less stellar angular momentum than typical local spiral galaxies of equal stellar mass and bulge fraction. Theoretical considerations show that this deficiency in angular momentum is the main cause of their low stability, while the high gas fraction plays a complementary role. Interestingly, the low jstar values of our targets are similar to those expected in the M*-population at higher z from the approximate theoretical scaling jstar~(1+z)^(-1/2) at fixed Mstar. This suggests that a change in angular momentum, driven by cosmic expansion, is the main cause for the remarkable difference between clumpy M*-disks at high z (which likely evolve into early-type galaxies) and mass-matched local spirals.

We have undertaken an ambitious program to visually classify all galaxies in the five CANDELS fields down to H<24.5 involving the dedicated efforts of 65 individual classifiers. Once completed, we expect to have detailed morphological classifications for over 50,000 galaxies up to z<4 over all the fields. Here, we present our detailed visual classification scheme, which was designed to cover a wide range of CANDELS science goals. This scheme includes the basic Hubble sequence types, but also includes a detailed look at mergers and interactions, the clumpiness of galaxies, $k$-corrections, and a variety of other structural properties. In this paper, we focus on the first field to be completed -- GOODS-S. The wide area coverage spanning the full field includes 7634 galaxies that have been classified by at least three different people. In the deep area of the field, 2534 galaxies have been classified by at least five different people at three different depths. With this paper, we release to the public all of the visual classifications in GOODS-S along with the GUI that we developed to classify galaxies. We find that the level of agreement among classifiers is good and depends on both the galaxy magnitude and the galaxy type, with disks showing the highest level of agreement and irregulars the lowest. A comparison of our classifications with the Sersic index and rest-frame colors shows a clear separation between disk and spheroid populations. Finally, we explore morphological k-corrections between the V-band and H-band observations and find that a small fraction (84 galaxies in total) are classified as being very different between these two bands. These galaxies typically have very clumpy and extended morphology or are very faint in the V-band.

We study the relationship between the structure and star-formation rate (SFR) of X-ray selected low and moderate luminosity active galactic nuclei (AGNs) in the two Chandra Deep Fields, using Hubble Space Telescope imaging from the Cosmic Assembly Near Infrared Extragalactic Legacy Survey (CANDELS) and deep far-infrared maps from the PEP+GOODS-Herschel survey. We derive detailed distributions of structural parameters and FIR luminosities from carefully constructed control samples of galaxies, which we then compare to those of the AGNs. At z~1, AGNs show slightly diskier light profiles than massive inactive (non-AGN) galaxies, as well as modestly higher levels of gross galaxy disturbance (as measured by visual signatures of interactions and clumpy structure). In contrast, at z~2, AGNs show similar levels of galaxy disturbance as inactive galaxies, but display a red central light enhancement, which may arise due to a more pronounced bulge in AGN hosts or due to extinguished nuclear light. We undertake a number of tests of these alternatives, but our results do not strongly favour one interpretation over the other. The mean SFR and its distribution among AGNs and inactive galaxies are similar at z>1.5. At z<1, however, clear and significant enhancements are seen in the SFRs of AGNs with bulge-dominated light profiles. These trends suggest an evolution in the relation between nuclear activity and host properties with redshift, towards a minor role for mergers and interactions at z>1.5.

In this letter we report the discovery of CO fluxes, suggesting very high gas fractions in three disk galaxies seen in the nearby Universe (z ~ 0.1). These galaxies were investigated as part of the DYnamics of Newly Assembled Massive Objects (DYNAMO) survey. High-resolution Hubble Space Telescope imaging of these objects reveals the presence of large star forming clumps in the bodies of the galaxies, while spatially resolved spectroscopy of redshifted Halpha reveals the presence of high dispersion rotating disks. The internal dynamical state of these galaxies resembles that of disk systems seen at much higher redshifts (1 < z < 3). Using CO(1-0) observations made with the Plateau de Bure Interferometer, we find gas fractions of 20-30% and depletion times of tdep ~ 0.5 Gyr (assuming a Milky Way-like CO conversion factor). These properties are unlike those expected for low- redshift galaxies of comparable specific star formation rate, but they are normal for their high-z counterparts. DYNAMO galaxies break the degeneracy between gas fraction and redshift, and we show that the depletion time per specific star formation rate for galaxies is closely tied to gas fraction, independent of redshift. We also show that the gas dynamics of two of our local targets corresponds to those expected from unstable disks, again resembling the dynamics of high-z disks. These results provide evidence that DYNAMO galaxies are local analogues to the clumpy, turbulent disks, which are often found at high redshift.

We study the spatially resolved stellar kinematics of two star-forming galaxies at z = 0.1 from the larger DYnamics of Newly Assembled Massive Objects (DYNAMO) sample. These galaxies, which have been characterized by high levels of star formation and large ionized gas velocity dispersions, are considered possible analogs to high-redshift clumpy disks. They were observed using the GMOS instrument in integral field spectroscopy (IFS) mode at the Gemini Observatory with high spectral resolution (R=5400, equivalent to 24 km/s at the observed wavelengths) and 6 hour exposure times in order to measure the resolved stellar kinematics via absorption lines. We also obtain higher-quality emission line kinematics than previous observations. The spatial resolution (1.2 kpc) is sufficient to show that the ionized gas in these galaxies (as traced by H-beta emission) is morphologically irregular, forming multiple giant clumps while stellar continuum light is smooth and well described by an exponential profile. Clumpy gas morphologies observed in IFS data are confirmed by complementary narrow band H-alpha imaging from the Hubble Space Telescope. Morphological differences between the stars and ionized gas are not reflected dynamically as stellar kinematics are found the be closely coupled to the kinematics of the ionized gas: both components are smoothly rotating with large velocity dispersions (~40 km/s) suggesting that the high gas dispersions are not primarily driven by star-formation feedback. In addition, the stellar population ages of these galaxies are estimated to be quite young (60-500 Myr). The large velocity dispersions measured for these young stars suggest that we are seeing the formation of thick disks and/or stellar bulges in support of recent models which produce these from clumpy galaxies at high redshift.

We study the environmental dependence of color, stellar mass, and morphology by comparing galaxies in a forming cluster to those in the field at z = 1:6 with Hubble Space Telescope near-infrared imaging in the CANDELS/UDS field. We quantify the morphology of the galaxies using the effective radius, reff, and S\'ersic index, n. In both the cluster and field, approximately half of the bulge-dominated galaxies (n > 2) reside on the red sequence of the color-magnitude diagram, and most disk-dominated galaxies (n < 2) have colors expected for star-forming galaxies. There is weak evidence that cluster galaxies have redder rest-frame U - B colors and higher stellar masses compared to the field. Star-forming galaxies in both the cluster and field show no significant differences in their morphologies. In contrast, there is evidence that quiescent galaxies in the cluster have larger median effective radii and smaller S\'ersic indices compared to the field with a significance of 2?. These differences are most pronounced for galaxies at clustercentric distances 1 Mpc < Rproj < 1.5 Mpc, which have low S\'ersic indices and possibly larger effective radii, more consistent with star-forming galaxies at this epoch and in contrast to other quiescent galaxies. We argue that star-forming galaxies are processed under the influence of the cluster environment at distances greater than the cluster-halo virial radius. Our results are consistent with models where gas accretion onto these galaxies is suppressed from processes associated with the cluster environment.

We study the evolution of galaxy rest-frame ultraviolet (UV) colors in the epoch 4 < z < 8. We use new wide-field near-infrared data in GOODS-S from the CANDELS, HUDF09 and ERS programs to select galaxies via photometric redshift measurements. Our sample consists of 2812 candidate galaxies at z > 3.5, including 113 at z = 7 to 8. We fit the observed spectral energy distribution to a suite of synthetic stellar population models, and measure the value of the UV spectral slope (beta) from the best-fit model spectrum. The median value of beta evolves significantly from -1.82 (+0.00,-0.04) at z = 4, to -2.37 (+0.26,-0.06) at z = 7. Additionally, we find that faint galaxies at z = 7 have beta = -2.68 (+0.39,-0.24) (~ -2.4 after correcting for observational bias); this is redder than previous claims in the literature, and does not require "exotic" stellar populations to explain their colors. This evolution can be explained by an increase in dust extinction, with the timescale consistent with low-mass AGB stars forming the bulk of the dust. We find no significant (< 2-sigma) correlation between beta and M_UV when measuring M_UV at a consistent rest-frame wavelength of 1500 A. This is particularly true at bright magnitudes, though our results do show evidence for a weak correlation at faint magnitudes when galaxies in the HUDF are considered separately, hinting that dynamic range in sample luminosities may play a role. We do find a strong correlation between beta and the stellar mass at all redshifts, in that more massive galaxies exhibit redder colors. The most massive galaxies in our sample have red colors at each redshift, implying that dust can build up quickly in massive galaxies, and that feedback is likely removing dust from low-mass galaxies at z > 7. Thus the stellar-mass - metallicity relation, previously observed up to z ~ 3, may extend out to z = 7 - 8.

We discuss the structural and morphological properties of galaxies in a z=1.62 proto-cluster using near-IR imaging data from Hubble Space Telescope Wide Field Camera 3 data of the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). The cluster galaxies exhibit a clear color-morphology relation: galaxies with colors of quiescent stellar populations generally have morphologies consistent with spheroids, and galaxies with colors consistent with ongoing star formation have disk-like and irregular morphologies. The size distribution of the quiescent cluster galaxies shows a deficit of compact (< 1kpc), massive galaxies compared to CANDELS field galaxies at z=1.6. As a result the cluster quiescent galaxies have larger average effective sizes compared to field galaxies at fixed mass at greater than 90% significance. Combined with data from the literature, the size evolution of quiescent cluster galaxies is relatively slow from z~1.6 to the present, growing as (1+z)^(-0.6+/-0.1). If this result is generalizable, then it implies that physical processes associated with the denser cluster region seems to have caused accelerated size growth in quiescent galaxies prior to z=1.6 and slower subsequent growth at z<1.6 compared to galaxies in the lower density field. The quiescent cluster galaxies at z=1.6 have higher ellipticities compared to lower redshift samples at fixed mass, and their surface-brightness profiles suggest that they contain extended stellar disks. We argue the cluster galaxies require dissipationless (i.e., gas-poor or "dry") mergers to reorganize the disk material and to match the relations for ellipticity, stellar mass, size, and color of early-type galaxies in z<1 clusters.