Physical properties of the interstellar medium (ISM) at sub-galactic ($\sim$kpc) scales play an indispensable role in controlling the ability of gas to form stars. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, in this paper, we use the TNG50 cosmological simulation of the IllustrisTNG project to explore the physical parameter space comprised of 8 resolved ISM properties in star-forming regions across a wide range of host galaxy mass (M$_\star = 10^{7-11}$ M$_\odot$) and redshift ($0 \leq z \leq 3$). We deconstruct our simulated galaxies into regions of 1 kpc in size, and measure the gas/stellar surface densities, gas metallicity, vertical stellar velocity dispersion, epicyclic frequency and dark-matter volumetric density representative of each region, and study them in the context of their star formation activity and galactic environment (radial galactocentric location). By examining the star formation rate weighted distributions of these properties, we show that stars primarily form in two distinct regimes, which are brought about by an underlying bi-component radial star formation rate profile in galaxies. We examine how the relative prominence of these two regimes depends on host galaxy mass and cosmic time, and compare our findings with observations from the MaNGA IFU survey. Further, using principal component analysis, we characterise the aforementioned parameter space to reveal a high-degree of multicollinearity in relationships between ISM properties that drive the distribution of star formation at kiloparsec scales. Based on this, we find that a reduced three dimensional representation underpinned in essence by a multi-variate radius relationship is sufficient to capture most of the variance in the original 8D space.

Semi-analytic models (SAMs) are a promising means of tracking the physical processes associated with galaxy formation, but many of their approximations have not been rigorously tested. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, we compare predictions from the FIRE-2 hydrodynamical "zoom-in" simulations to those from the Santa Cruz SAM run on the same halo merger trees, with an emphasis on the global mass flow cycle. Our study includes 13 halos spanning low-mass dwarfs (M_vir~10^10 M_sun at z=0), intermediate-mass dwarfs (M_vir~10^11 M_sun) and Milky Way-mass galaxies (M_vir~10^12 M_sun). The SAM and FIRE-2 predictions agree relatively well with each other in terms of stellar and interstellar mass, but differ dramatically on circumgalactic mass (the SAM is lower than FIRE-2 by ~3 orders of magnitude for dwarfs). Strikingly, the SAM predicts higher gas accretion rates for dwarfs compared to FIRE-2 by factors of ~10-100, and this is compensated for with higher mass outflow rates in the SAM. We argue that the most severe model discrepancies are caused by the lack of preventative stellar feedback and the assumptions for halo gas cooling and recycling in the SAM. As a first step towards resolving these model tensions, we present a simple yet promising new preventative stellar feedback model in which the energy carried by supernova-driven winds is allowed to heat some fraction of gas outside of halos to at least the virial temperature such that accretion is suppressed.

Large scale outflows in star-forming galaxies are observed to be ubiquitous, and are a key aspect of theoretical modeling of galactic evolution in a cosmological context, the focus of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project. Gas blown out from galactic disks, similar to gas within galaxies, consists of multiple phases with large contrasts of density, temperature, and other properties. To study multiphase outflows as emergent phenomena, we run a suite of ~pc-resolution local galactic disk simulations using the TIGRESS framework. Explicit modeling of the interstellar medium (ISM), including star formation and self-consistent radiative heating plus supernova feedback, regulates ISM properties and drives the outflow. We investigate the scaling of outflow mass, momentum, energy, and metal loading factors with galactic disk properties, including star formation rate (SFR) surface density (\Sigma_SFR~10^{-4}-1 M_sun/kpc^2/yr), gas surface density (~1-100 M_sun/pc^2), and total midplane pressure (or weight) (~10^3-10^6 k_B cm^{-3} K). The main components of outflowing gas are mass-delivering cool gas (T~10^4 K) and energy/metal-delivering hot gas (T~10^6 K). Cool mass outflow rates measured at outflow launch points (one or two scale heights) are 1-100 times the SFR (decreasing with \Sigma_SFR), although in massive galaxies most mass falls back due to insufficient outflow velocity. The hot galactic outflow carries mass comparable to 10% of the SFR, together with 10-20% of the energy and 30-60% of the metal mass injected by SN feedback. The characteristic outflow velocities of both phases scale very weakly with SFR, as v_out \propto \Sigma_SFR^{0.1~0.2}, consistent with observations. Importantly, our analysis demonstrates that in any physically-motivated cosmological wind model, it is crucial to include at least two distinct thermal wind components.

Existing models of galaxy formation have not yet explained striking correlations between structure and star-formation activity in galaxies, notably the sloped and moving boundaries that divide star-forming from quenched galaxies in key structural diagrams. This paper uses these and other relations to ``reverse-engineer'' the quenching process for central galaxies. The basic idea is that star-forming galaxies with larger radii (at a given stellar mass) have lower black-hole masses due to lower central densities. Galaxies cross into the green valley when the cumulative effective energy radiated by their black hole equals $\sim4\times$ their halo-gas binding energy. Since larger-radii galaxies have smaller black holes, one finds they must evolve to higher stellar masses in order to meet this halo-energy criterion, which explains the sloping boundaries. A possible cause of radii differences among star-forming galaxies is halo concentration. The evolutionary tracks of star-forming galaxies are nearly parallel to the green-valley boundaries, and it is mainly the sideways motions of these boundaries with cosmic time that cause galaxies to quench. BH-scaling laws for star-forming, quenched, and green-valley galaxies are different, and most BH mass growth takes place in the green valley. Implications include: the radii of star-forming galaxies are an important second parameter in shaping their black holes; black holes are connected to their halos but in different ways for star-forming, quenched, and green-valley galaxies; and the same BH-halo quenching mechanism has been in place since $z \sim 3$. We conclude with a discussion of black hole-galaxy co-evolution, the origin and interpretation of BH scaling laws.

This paper studies pseudo-bulges (P-bulges) and classical bulges (C-bulges) in Sloan Digital Sky Survey central galaxies using the new bulge indicator $\Delta\Sigma_1$, which measures relative central stellar-mass surface density within 1 kpc. We compare $\Delta\Sigma_1$ to the established bulge-type indicator $\Delta\langle\mu_e\rangle$ from Gadotti (2009) and show that classifying by $\Delta\Sigma_1$ agrees well with $\Delta\langle\mu_e\rangle$. $\Delta\Sigma_1$ requires no bulge-disk decomposition and can be measured on SDSS images out to $z = 0.07$. Bulge types using it are mapped onto twenty different structural and stellar-population properties for 12,000 SDSS central galaxies with masses 10.0 < log $M_*$/$M_{\odot}$ < 10.4. New trends emerge from this large sample. Structural parameters show fairly linear log-log relations vs. $\Delta\Sigma_1$ and $\Delta\langle\mu_e\rangle$ with only moderate scatter, while stellar-population parameters show a highly non-linear "elbow" in which specific star-formation rate remains roughly flat with increasing central density and then falls rapidly at the elbow, where galaxies begin to quench. P-bulges occupy the low-density end of the horizontal arm of the elbow and are universally star-forming, while C-bulges occupy the elbow and the vertical branch and exhibit a wide range of star-formation rates at fixed density. The non-linear relation between central density and star-formation rate has been seen before, but this mapping onto bulge class is new. The wide range of star-formation rates in C-bulges helps to explain why bulge classifications using different parameters have sometimes disagreed in the past. The elbow-shaped relation between density and stellar indices suggests that central structure and stellar-populations evolve at different rates as galaxies begin to quench.

This paper is one in a series that explores the importance of radius as a second parameter in galaxy evolution. The topic investigated here is the relationship between star formation rate (SFR) and galaxy radius ($R_{\rm e}$) for main-sequence star-forming galaxies. The key observational result is that, over a wide range of stellar mass and redshift in both CANDELS and SDSS, there is little trend between SFR and $R_{\rm e}$ at fixed stellar mass. The Kennicutt-Schmidt law, or any similar density-related star formation law, then implies that smaller galaxies must have lower gas fractions than larger galaxies (at fixed $M_{\ast}$), and this is supported by observations of local star-forming galaxies. We investigate the implication by adopting the equilibrium "bathtub" model: the ISM gas mass is assumed to be constant over time and the net star formation rate is the difference between the accretion rate of gas onto the galaxy from the halo and the outflow rate due to winds. To match the observed null correlation between SFR and radius, the bathtub model requires that smaller galaxies at fixed mass have weaker galactic winds. Our hypothesis is that galaxies are a 2-dimensional family whose properties are set mainly by halo mass and concentration. Galaxy radius and accretion rate plausibly both depend on halo concentration, which predicts how wind strength should vary with $R_{\rm e}$ and SFR.

We present a WFC3 F160W ($H$-band) selected catalog in the CANDELS/GOODS-N field containing photometry from the ultraviolet (UV) to the far-infrared (IR), photometric redshifts and stellar parameters derived from the analysis of the multi-wavelength data. The catalog contains 35,445 sources over the 171 arcmin$^{2}$ of the CANDELS F160W mosaic. The 5$\sigma$ detection limits (within an aperture of radius 0\farcs17) of the mosaic range between $H=27.8$, 28.2 and 28.7 in the wide, intermediate and deep regions, that span approximately 50\%, 15\% and 35\% of the total area. The multi-wavelength photometry includes broad-band data from UV (U band from KPNO and LBC), optical (HST/ACS F435W, F606W, F775W, F814W, and F850LP), near-to-mid IR (HST/WFC3 F105W, F125W, F140W and F160W, Subaru/MOIRCS Ks, CFHT/Megacam K, and \spitzer/IRAC 3.6, 4.5, 5.8, 8.0 $\mu$m) and far IR (\spitzer/MIPS 24$\mu$m, HERSCHEL/PACS 100 and 160$\mu$m, SPIRE 250, 350 and 500$\mu$m) observations. In addition, the catalog also includes optical medium-band data (R$\sim50$) in 25 consecutive bands, $\lambda=500$ to 950~nm, from the SHARDS survey and WFC3 IR spectroscopic observations with the G102 and G141 grisms (R$\sim210$ and 130). The use of higher spectral resolution data to estimate photometric redshifts provides very high, and nearly uniform, precision from $z=0-2.5$. The comparison to 1,485 good quality spectroscopic redshifts up to $z\sim3$ yields $\Delta z$/(1+$z_{\rm spec}$)$=$0.0032 and an outlier fraction of $\eta=$4.3\%. In addition to the multi-band photometry, we release added-value catalogs with emission line fluxes, stellar masses, dust attenuations, UV- and IR- based star formation rates and rest-frame colors.

Hubble Space Telescope observations show that low-mass ($M_*=10^9-10^{10}M_{\odot}$) galaxies at high redshift ($z=1.0-2.5$) tend to be elongated (prolate) rather than disky (oblate) or spheroidal. This is explained in zoom-in cosmological hydrodynamical simulations by the fact that these galaxies are forming in cosmic web filaments where accretion happens preferentially along the direction of elongation. We ask whether the elongated morphology of these galaxies allows them to be used as effective tracers of cosmic web filaments at high redshift via their intrinsic alignments. Using mock lightcones and spectroscopically-confirmed galaxy pairs from the CANDELS survey, we test two types of alignments: (1) between the galaxy major axis and the direction to nearby galaxies of any mass, and (2) between the major axes of nearby pairs of low-mass, likely prolate, galaxies. The mock lightcones predict strong signals in 3D real space, 3D redshift space, and 2D projected redshift space for both types of alignments (assuming prolate galaxy orientations are the same as those of their host prolate halos), but we do not detect significant alignment signals in CANDELS observations. However, we show that spectroscopic redshifts have been obtained for only a small fraction of highly elongated galaxies, and accounting for spectroscopic incompleteness and redshift errors significantly degrades the 2D mock signal. This may partly explain the alignment discrepancy and highlights one of several avenues for future work.

In this paper we study the molecular gas content of a representative sample of 67 of the most massive early-type galaxies in the local universe, drawn uniformly from the MASSIVE survey. We present new IRAM-30m telescope observations of 30 of these galaxies, allowing us to probe the molecular gas content of the entire sample to a fixed molecular-to-stellar mass fraction of 0.1%. The total detection rate in this representative sample is 25$^{+5.9}_{-4.4}$%, and by combining the MASSIVE and ATLAS$^{\rm 3D}$ molecular gas surveys we find a joint detection rate of 22.4$^{+2.4}_{-2.1}$%. This detection rate seems to be independent of galaxy mass, size, position on the fundamental plane, and local environment. We show here for the first time that true slow rotators can host molecular gas reservoirs, but the rate at which they do so is significantly lower than for fast-rotators. Objects with a higher velocity dispersion at fixed mass (a higher kinematic bulge fraction) are less likely to have detectable molecular gas, and where gas does exist, have lower molecular gas fractions. In addition, satellite galaxies in dense environments have $\approx$0.6 dex lower molecular gas-to-stellar mass ratios than isolated objects. In order to interpret these results we created a toy model, which we use to constrain the origin of the gas in these systems. We are able to derive an independent estimate of the gas-rich merger rate in the low-redshift universe. These gas rich mergers appear to dominate the supply of gas to ETGs, but stellar mass loss, hot halo cooling and transformation of spiral galaxies also play a secondary role.

Ultra diffuse galaxies are a population of extended galaxies but with relatively low luminosities. The origin of these objects remains unclear, largely due to the observational challenges of the low surface brightness Universe. We present here a detailed stellar population analysis of a relatively isolated UDG, DGSAT I, based on spectroscopic data from the Keck Cosmic Web Imager integral field unit. The star formation history of DGSAT I seems to be extended, with a mean luminosity-weighted age of ~3 Gyr, in agreement with previous photometric studies. However, we find a very high [Mg/Fe] abundance ratio, which is extreme even in the context of the highly alpha-enhanced massive ellipticals and ultra-faint dwarfs. The [Mg/Fe]-enhancement of DGSAT I appears to be 10 times higher than the most magnesium-enhanced stellar systems discovered to date, and suggests that the chemical enrichment of this object was dominated by core-collapse supernovae. Intriguingly, this breaks the canonical relation between [Mg/Fe] and star formation time-scale. With a measured velocity dispersion of 56 +/- 10 km/s, DGSAT I also shows a high mass-to-light ratio, which indicates that it is highly dark matter-dominated. The metal-poor conditions of DGSAT I may have enhanced the formation of massive stars, while at the same time, additional mechanisms are needed to prevent iron-rich yields from being recycled into stars. These results suggest that some ultra-diffuse galaxies could have experienced chemical enrichment episodes similar to the first building blocks of galaxies.

A tightly correlated star formation rate-stellar mass relation of star forming galaxies, or star-forming sequence (SFS), is a key feature in galaxy property-space that is predicted by modern galaxy formation models. We present a flexible data-driven approach for identifying this SFS over a wide range of star formation rates and stellar masses using Gaussian mixture modeling (GMM). Using this method, we present a consistent comparison of the $z=0$ SFSs of central galaxies in the Illustris, EAGLE, and Mufasa hydrodynamic simulations and the Santa Cruz semi-analytic model (SC-SAM), alongside data from the Sloan Digital Sky Survey. We find, surprisingly, that the amplitude of the SFS varies by up to ${\sim} 0.7\,\mathrm{dex}$ (factor of ${\sim} 5$) among the simulations with power-law slopes range from $0.7$ to $1.2$. In addition to the SFS, our GMM method also identifies sub-components in the star formation rate-stellar mass relation corresponding to star-burst, transitioning, and quiescent sub-populations. The hydrodynamic simulations are similarly dominated by SFS and quiescent sub-populations unlike the SC-SAM, which predicts substantial fractions of transitioning and star-burst galaxies at stellar masses above and below $10^{10} M_\odot$, respectively. All of the simulations also produce an abundance of low-mass quiescent central galaxies in apparent tension with observations. These results illustrate that, even among models that well reproduce many observables of the galaxy population, the $z=0$ SFS and other sub-populations still show marked differences that can provide strong constraints on galaxy formation models.

In this second paper of the series, we study with new Keck/DEIMOS spectra the stellar populations of seven spectroscopically confirmed ultra-diffuse galaxies (UDGs) in the Coma cluster. We find intermediate to old ages (~ 7Gyr), low metallicities ([Z/H] ~ -0.7dex) and mostly super-solar abundance patterns ([Mg/Fe] ~ 0.13dex). These properties are similar to those of low-luminosity (dwarf) galaxies inhabiting the same area in the cluster and are most consistent with being the continuity of the stellar mass scaling relations of more massive galaxies. These UDGs' star formation histories imply a relatively recent infall into the Coma cluster, consistent with the theoretical predictions for a dwarf-like origin. However, considering the scatter in the resulting properties and including other UDGs in Coma, together with the results from the velocity phase-space study of the Paper I in this series, a mixed-bag of origins is needed to explain the nature of all UDGs. Our results thus reinforce a scenario in which many UDGs are field dwarfs that become quenched through their later infall onto cluster environments, whereas some UDGs could be genuine primordial galaxies that failed to develop due to an early quenching phase. The unknown proportion of dwarf-like to primordial-like UDGs leaves the enigma of the nature of UDGs still open.

We apply the Jeans Anisotropic MGE (JAM) dynamical modelling method to SAGES Legacy Unifying Globulars and GalaxieS (SLUGGS) survey data of early-type galaxies in the stellar mass range $10^{10}<M_*/{\rm M}_{\odot}<10^{11.6}$ that cover a large radial range of $0.1-4.0$ effective radii. We combine SLUGGS and ATLAS$^{\rm 3D}$ datasets to model the total-mass profiles of a sample of 21 fast-rotator galaxies, utilising a hyperparameter method to combine the two independent datasets. The total-mass density profile slope values derived for these galaxies are consistent with those measured in the inner regions of galaxies by other studies. Furthermore, the total-mass density slopes ($\gamma_{\rm tot}$) appear to be universal over this broad stellar mass range, with an average value of $\gamma_{\rm tot}=-2.12\,\pm\,0.05$, i.e. slightly steeper than isothermal. We compare our results to model galaxies from the Magneticum and EAGLE cosmological hydrodynamic simulations, in order to probe the mechanisms that are responsible for varying total-mass density profile slopes. The simulated-galaxy slopes are shallower than the observed values by $\sim0.1-0.3$, indicating that the physical processes shaping the mass distributions of galaxies in cosmological simulations are still incomplete. For galaxies with $M_*>10^{10.7}{\rm M}_{\odot}$ in the Magneticum simulations, we identify a significant anticorrelation between total-mass density profile slopes and the fraction of stellar mass formed ex situ (i.e. accreted), whereas this anticorrelation is weaker for lower stellar masses, implying that the measured total mass density slopes for low-mass galaxies are less likely to be determined by merger activity.

We present results from an optical study of two young Balmer-dominated remnants of SNIa in the Large Magellanic Cloud, 0509$-$67.5 and 0519$-$69.0, in an attempt to search for signatures of efficient cosmic-ray (CR) acceleration. We combine proper motion measurements from HST with corresponding optical spectroscopic measurements of the H$\alpha$ line at multiple rim positions from VLT/FORS2 and SALT/RSS and compare our results to published Balmer shock models. Analysis of the optical spectra result in broad H$\alpha$ widths in the range of 1800-4000 km s$^{-1}$ for twelve separate Balmer-dominated filaments that show no evidence for forbidden line emission, the corresponding shock speeds from proper motion measurements span a range of 1700-8500 km s$^{-1}$. Our measured values of shock speeds and broad H$\alpha$ widths in 0509$-$67.5 and 0519$-$69.0 are fit well with a Balmer shock model that does not include effects of efficient CR acceleration. We determine an upper limit of 7%/$\chi$ (95% confidence) on the CR acceleration efficiency for our ensemble of data points, where $\chi$ is the ionization fraction of the pre-shock gas. The upper limits on the individual remnants are 6%/$\chi$ (0509$-$67.5) and 11%/$\chi$ (0519$-$69.0). These upper limits are below the integrated CR acceleration efficiency in the Tycho supernova remnant, where the shocks predominantly show little H$\alpha$ emission, indicating that Balmer-dominated shocks are not efficient CR accelerators.

This is the first in a series of papers examining the demographics of star-forming galaxies at $0.2<z<2.5$ in CANDELS. We study 9,100 galaxies from GOODS-S and UDS having published values of redshifts, masses, star-formation rates (SFRs), and dust attenuation ($A_V$) derived from UV-optical SED fitting. In agreement with previous works, we find that the $UVJ$ colors of a galaxy are closely correlated with its specific star-formation rate (SSFR) and $A_V$. We define rotated $UVJ$ coordinate axes, termed $S_\mathrm{SED}$ and $C_\mathrm{SED}$, that are parallel and perpendicular to the star-forming sequence and derive a quantitative calibration that predicts SSFR from $C_\mathrm{SED}$ with an accuracy of ~0.2 dex. SFRs from UV-optical fitting and from UV+IR values based on Spitzer/MIPS 24 $\mu\mathrm{m}$ agree well overall, but systematic differences of order 0.2 dex exist at high and low redshifts. A novel plotting scheme conveys the evolution of multiple galaxy properties simultaneously, and dust growth, as well as star-formation decline and quenching, exhibit "mass-accelerated evolution" ("downsizing"). A population of transition galaxies below the star-forming main sequence is identified. These objects are located between star-forming and quiescent galaxies in $UVJ$ space and have lower $A_V$ and smaller radii than galaxies on the main sequence. Their properties are consistent with their being in transit between the two regions. The relative numbers of quenched, transition, and star-forming galaxies are given as a function of mass and redshift.

We present observational constraints on the stellar populations of two ultra-diffuse galaxies (UDGs) using optical through near-infrared (NIR) spectral energy distribution (SED) fitting. Our analysis is enabled by new $Spitzer$-IRAC 3.6 $\mu$m and 4.5 $\mu$m imaging, archival optical imaging, and the prospector fully Bayesian SED fitting framework. Our sample contains one field UDG (DGSAT I), one Virgo cluster UDG (VCC 1287), and one Virgo cluster dwarf elliptical for comparison (VCC 1122). We find that the optical--NIR colors of the three galaxies are significantly different from each other. We infer that VCC 1287 has an old ($\gtrsim7.7$ Gyr) and surprisingly metal-poor ($[Z/Z_{\odot}]\lesssim-1.0$) stellar population, even after marginalizing over uncertainties on diffuse interstellar dust. In contrast, the field UDG DGSAT I shows evidence of being younger than the Virgo UDG, with an extended star formation history and an age posterior extending down to $\sim3$ Gyr. The stellar metallicity of DGSAT I is sub-solar but higher than that of the Virgo UDG, with $[Z/Z_{\odot}]=-0.63^{+0.35}_{-0.62}$; in the case of exactly zero diffuse interstellar dust, DGSAT I may even have solar metallicity. With VCC 1287 and several Coma UDGs, a general picture is emerging where cluster UDGs may be "failed" galaxies, but the field UDG DGSAT I seems more consistent with a stellar feedback-induced expansion scenario. In the future, our approach can be applied to a large and diverse sample of UDGs down to faint surface brightness limits, with the goal of constraining their stellar ages, stellar metallicities, and circumstellar and diffuse interstellar dust content.

We use Keck/DEIMOS spectroscopy to confirm the cluster membership of 16 ultra-diffuse galaxies (UDGs) in the Coma cluster, bringing the total number of spectroscopically con- firmed UDGs to 24. We also identify a new cluster background UDG. In this pilot study of Coma UDGs in velocity phase-space, we find evidence that most present-day Coma UDGs have a recent infall epoch while a few may be ancient infalls. These recent infall UDGs have higher absolute relative line-of-sight velocities, bluer optical colors, and are smaller in size, unlike the ancient infalls. The kinematics of the spectroscopically confirmed Coma UDG sample is similar to that of the cluster late-type galaxy population. Our velocity phase-space analysis suggests that present-day cluster UDGs have a predominantly accretion origin from the field, acquire velocities corresponding to the mass of the cluster at accretion as they are accelerated towards the cluster center, and become redder and bigger as they experience the various physical processes at work within the cluster.

We estimate the cosmic number density of the recently identified class of HI-bearing ultra-diffuse sources (HUDs) based on the completeness limits of the ALFALFA survey. These objects fall in the range $8.5 < \log M_{\rm{HI}}/\rm{M_{\odot}} < 9.5$, have average $r$-band surface brightnesses fainter than 24 mag arsec$^{-2}$, half-light radii greater than 1.5 kpc, and are separated from neighbours by at least 350 kpc. We find that HUDs contribute at most 6% of the population of HI-bearing dwarfs, have a total cosmic number density of $(1.5 \pm 0.6) \times 10^{-3}$ $\rm{Mpc^{-3}}$, and an HI mass density of $(6.0 \pm 0.8) \times 10^{5}$ $\rm{M_{\odot}\,Mpc^{-3}}$. We estimate that this is similar to the total cosmic number density of UDGs in groups and clusters, and conclude that the relation between the number of UDGs hosted in a halo and the halo mass, must have a break below $M_{200} \sim 10^{12}$ $\rm{M_{\odot}}$ in order to account for the abundance of HUDs. The distribution of the velocity widths of HUDs rises steeply towards low values, indicating a preference for slow rotation rates. These objects have been absent from measurements of the galaxy stellar mass function owing to their low surface brightness. However, we estimate that due to their low number density, their inclusion would constitute a correction of less than 1%. Comparison with the Santa Cruz SAM shows that it produces HI-rich central UDGs that have similar colours to HUDs, but these are currently produced in much great a number. While previous results from this sample have favoured formation scenarios where HUDs form in high spin parameter halos, comparisons with the results of Rong et al. 2017, which invokes that formation mechanism, reveal that this model produces an order of magnitude more field UDGs than we observe in the HUD population.(Abridged)

We present high spatial resolution (FWHM$\sim$0.14'') observations of the CO($8-7$) line in GDS-14876, a compact star-forming galaxy at $z=2.3$ with total stellar mass of $\log(M_{\star}/M_{\odot})=10.9$. The spatially resolved velocity map of the inner $r\lesssim1$~kpc reveals a continous velocity gradient consistent with the kinematics of a rotating disk with $v_{\rm rot}(r=1\rm kpc)=163\pm5$ km s$^{-1}$ and $v_{\rm rot}/\sigma\sim2.5$. The gas-to-stellar ratios estimated from CO($8-7$) and the dust continuum emission span a broad range, $f^{\rm CO}_{\rm gas}=M_{\rm gas}/M_{\star}=13-45\%$ and $f^{\rm cont}_{\rm gas}=50-67\%$, but are nonetheless consistent given the uncertainties in the conversion factors. The dynamical modeling yields a dynamical mass of$\log(M_{\rm dyn}/M_{\odot})=10.58^{+0.5}_{-0.2}$ which is lower, but still consistent with the baryonic mass, $\log$(M$_{\rm bar}$= M$_{\star}$ + M$^{\rm CO}_{\rm gas}$/M$_{\odot}$)$=11.0$, if the smallest CO-based gas fraction is assumed. Despite a low, overall gas fraction, the small physical extent of the dense, star-forming gas probed by CO($8-7$), $\sim3\times$ smaller than the stellar size, implies a strong concentration that increases the gas fraction up to $f^{\rm CO, 1\rm kpc}_{\rm gas}\sim 85\%$ in the central 1 kpc. Such a gas-rich center, coupled with a high star-formation rate, SFR$\sim$ 500 M$_{\odot}$ yr$^{-1}$, suggests that GDS-14876 is quickly assembling a dense stellar component (bulge) in a strong nuclear starburst. Assuming its gas reservoir is depleted without replenishment, GDS-14876 will quickly ($t_{\rm depl}\sim27$ Myr) become a compact quiescent galaxy that could retain some fraction of the observed rotational support.

The relation between X-ray luminosity (L_X) and ambient gas temperature (T) among massive galactic systems is an important cornerstone of both observational cosmology and galaxy-evolution modeling. In the most massive galaxy clusters, the relation is determined primarily by cosmological structure formation. In less massive systems, it primarily reflects the feedback response to radiative cooling of circumgalactic gas. Here we present a simple but powerful model for the L_X-T relation as a function of physical aperture R within which those measurements are made. The model is based on the precipitation framework for AGN feedback and assumes that the circumgalactic medium is precipitation-regulated at small radii and limited by cosmological structure formation at large radii. We compare this model with many different data sets and show that it successfully reproduces the slope and upper envelope of the L_X-T-R relation over the temperature range from ~0.2 keV through >10 keV. Our findings strongly suggest that the feedback mechanisms responsible for regulating star formation in individual massive galaxies have much in common with the precipitation-triggered feedback that appears to regulate galaxy-cluster cores.

We explore observational and theoretical constraints on how galaxies might transition between the "star-forming main sequence" (SFMS) and varying "degrees of quiescence" out to $z=3$. Our analysis is focused on galaxies with stellar mass $M_*>10^{10}M_{\odot}$, and is enabled by GAMA and CANDELS observations, a semi-analytic model (SAM) of galaxy formation, and a cosmological hydrodynamical "zoom in" simulation with momentum-driven AGN feedback. In both the observations and the SAM, transition galaxies tend to have intermediate S\'ersic indices, half-light radii, and surface stellar mass densities compared to star-forming and quiescent galaxies out to $z=3$. We place an observational upper limit on the average population transition timescale as a function of redshift, finding that the average high-redshift galaxy is on a "fast track" for quenching whereas the average low-redshift galaxy is on a "slow track" for quenching. We qualitatively identify four physical origin scenarios for transition galaxies in the SAM: oscillations on the SFMS, slow quenching, fast quenching, and rejuvenation. Quenching timescales in both the SAM and the hydrodynamical simulation are not fast enough to reproduce the quiescent population that we observe at $z\sim3$. In the SAM, we do not find a clear-cut morphological dependence of quenching timescales, but we do predict that the mean stellar ages, cold gas fractions, SMBH masses, and halo masses of transition galaxies tend to be intermediate relative to those of star-forming and quiescent galaxies at $z<3$.

The bulk of the stellar growth over cosmic time is dominated by IR luminous galaxies at cosmic noon (z=1-2), many of which harbor a hidden active galactic nucleus (AGN). We use state of the art infrared color diagnostics, combining Spitzer and Herschel observations, to separate dust-obscured AGN from dusty star forming galaxies (SFGs) in the CANDELS and COSMOS surveys. We calculate 24 micron counts of SFGs, AGN/star forming "Composites", and AGN. AGN and Composites dominate the counts above 0.8 mJy at 24 micron, and Composites form at least 25% of an IR sample even to faint detection limits. We develop methods to use the Mid-Infrared Instrument (MIRI) on JWST to identify dust-obscured AGN and Composite galaxies from z~1-2. With the sensitivity and spacing of MIRI filters, we will detect >4 times as many AGN hosts than with Spitzer/IRAC criteria. Any star formation rates based on the 7.7 micron PAH feature (likely to be applied to MIRI photometry) must be corrected for the contribution of the AGN, or the SFR will be overestimated by ~35% for cases where the AGN provides half the IR luminosity and ~50% when the AGN accounts for 90% of the luminosity. Finally, we demonstrate that our MIRI color technique can select AGN with an Eddington ratio of $\lambda_{\rm Edd}\sim0.01$ and will identify AGN hosts with a higher sSFR than X-ray techniques alone. JWST/MIRI will enable critical steps forward in identifying and understanding dust-obscured AGN and the link to their host galaxies.

We explore empirical constraints on the statistical relationship between the radial size of galaxies and the radius of their host dark matter halos from $z\sim 0.1$--3 using the GAMA and CANDELS surveys. We map dark matter halo mass to galaxy stellar mass using relationships from abundance matching, applied to the Bolshoi-Planck dissipationless N-body simulation. We define SRHR$\equiv r_e/R_h$ as the ratio of galaxy radius to halo virial radius, and SRHR$\lambda \equiv r_e/(\lambda R_h)$ as the ratio of galaxy radius to halo spin parameter times halo radius. At $z\sim 0.1$, we find an average value of SRHR $\simeq 0.018$ and SRHR$\lambda \simeq 0.5$ with very little dependence on stellar mass. SRHR and SRHR$\lambda$ have a weak dependence on cosmic time since $z\sim 3$. SRHR shows a mild decrease over cosmic time for low mass galaxies, but increases slightly or does not evolve for more massive galaxies. We find hints that at high redshift ($z\sim 2$--3), SRHR$\lambda$ is lower for more massive galaxies, while it shows no significant dependence on stellar mass at $z\lesssim 0.5$. We find that for both the GAMA and CANDELS samples, at all redshifts from $z\sim 0.1$--3, the observed conditional size distribution in stellar mass bins is remarkably similar to the conditional distribution of $\lambda R_h$. We discuss the physical interpretation and implications of these results.

We present the first systematic investigation of the existence, spatial distribution, and kinematics of warm ionized gas as traced by the [O II] 3727AA emission line in 74 of the most massive galaxies in the local Universe. All of our galaxies have deep integral field spectroscopy from the volume- and magnitude-limited MASSIVE survey of early-type galaxies with stellar mass log(M_*/M_sun) > 11.5 (M_K < -25.3 mag) and distance D < 108 Mpc. Of the 74 galaxies in our sample, we detect warm ionized gas in 28, which yields a global detection fraction of 38\pm6% down to a typical [O II] equivalent width limit of 2AA. MASSIVE fast rotators are more likely to have gas than MASSIVE slow rotators with detection fractions of 80\pm10% and 28\pm6%, respectively. The spatial extents span a wide range of radii (0.6 - 18.2 kpc; 0.1 - 4R_e), and the gas morphologies are diverse, with 17/28 = 61\pm9% being centrally concentrated, 8/28 = 29\pm9% exhibiting clear rotation out to several kpc, and 3/28 = 11\pm6% being extended but patchy. Three out of four fast rotators show kinematic alignment between the stars and gas, whereas the two slow rotators with robust kinematic measurements available exhibit kinematic misalignment. Our inferred warm ionized gas masses are roughly ~10^5M_sun. The emission line ratios and radial equivalent width profiles are generally consistent with excitation of the gas by the old underlying stellar population. We explore different gas origin scenarios for MASSIVE galaxies and find that a variety of physical processes are likely at play, including internal gas recycling, cooling out of the hot gaseous halo, and gas acquired via mergers.

We study the correlation of galaxy structural properties with their location relative to the SFR-M* correlation, also known as the star formation "main sequence" (SFMS), in the CANDELS and GAMA surveys and in a semi-analytic model (SAM) of galaxy formation. We first study the distribution of median Sersic index, effective radius, star formation rate (SFR) density and stellar mass density in the SFR-M* plane. We then define a redshift dependent main sequence and examine the medians of these quantities as a function of distance from this main sequence, both above (higher SFRs) and below (lower SFRs). Finally, we examine the distributions of distance from the main sequence in bins of these quantities. We find strong correlations between all of these galaxy structural properties and the distance from the SFMS, such that as we move from galaxies above the SFMS to those below it, we see a nearly monotonic trend towards higher median Sersic index, smaller radius, lower SFR density, and higher stellar density. In the semi-analytic model, bulge growth is driven by mergers and disk instabilities, and is accompanied by the growth of a supermassive black hole which can regulate or quench star formation via Active Galactic Nucleus (AGN) feedback. We find that our model qualitatively reproduces the trends described above, supporting a picture in which black holes and bulges co-evolve, and AGN feedback plays a critical role in moving galaxies off of the SFMS.

We report evidence for excess blue light from the Type Ia supernova SN 2012cg at fifteen and sixteen days before maximum B-band brightness. The emission is consistent with predictions for the impact of the supernova on a non-degenerate binary companion. This is the first evidence for emission from a companion to a SN Ia. Sixteen days before maximum light, the B-V color of SN 2012cg is 0.2 mag bluer than for other normal SN~Ia. At later times, this supernova has a typical SN Ia light curve, with extinction-corrected M_B = -19.62 +/- 0.02 mag and Delta m_{15}(B) = 0.86 +/- 0.02. Our data set is extensive, with photometry in 7 filters from 5 independent sources. Early spectra also show the effects of blue light, and high-velocity features are observed at early times. Near maximum, the spectra are normal with a silicon velocity v_{Si} = -10,500$ km s^{-1}. Comparing the early data with models by Kasen (2010) favors a main-sequence companion of about 6 solar masses. It is possible that many other SN Ia have main-sequence companions that have eluded detection because the emission from the impact is fleeting and faint.

We present the first comprehensive archival study of the X-ray properties of ultracompact dwarf (UCD) galaxies, with the goal of identifying weakly-accreting central black holes in UCDs. Our study spans 578 UCDs distributed across thirteen different host systems, including clusters, groups, fossil groups, and isolated galaxies. Of the 336 spectroscopically-confirmed UCDs with usable archival Chandra imaging observations, 21 are X-ray-detected. Imposing a completeness limit of $L_X>2\times10^{38}$ erg s$^{-1}$, the global X-ray detection fraction for the UCD population is $\sim3\%$. Of the 21 X-ray-detected UCDs, seven show evidence of long-term X-ray time variability on the order of months to years. X-ray-detected UCDs tend to be more compact than non-X-ray-detected UCDs, and we find tentative evidence that the X-ray detection fraction increases with surface luminosity density and global stellar velocity dispersion. The X-ray emission of UCDs is fully consistent with arising from a population of low-mass X-ray binaries (LMXBs). In fact, there are fewer X-ray sources than expected using a naive extrapolation from globular clusters. Invoking the fundamental plane of black hole activity for SUCD1 near the Sombrero galaxy, for which archival Jansky Very Large Array imaging at 5 GHz is publicly available, we set an upper limit on the mass of a hypothetical central black hole in that UCD to be $\lesssim10^5M_{\odot}$. While the majority of our sources are likely LMXBs, we cannot rule out central black holes in some UCDs based on X-rays alone, and so we address the utility of follow-up radio observations to find weakly-accreting central black holes.

In this work we present CO(1-0) and CO(2-1) observations of a pilot sample of 15 early-type galaxies (ETGs) drawn from the MASSIVE galaxy survey, a volume-limited integral-field spectroscopic study of the most massive ETGs ($M_* >10^{11.5}M_\odot$) within 108 Mpc. These objects were selected because they showed signs of an interstellar medium and/or star formation. A large amount of gas ($>$2$\times$10$^8$ M$_{\odot}$) is present in 10 out of 15 objects, and these galaxies have gas fractions higher than expected based on extrapolation from lower mass samples. We tentatively interpret this as evidence that stellar mass loss and hot halo cooling may be starting to play a role in fuelling the most massive galaxies. These MASSIVE ETGs seem to have lower star-formation efficiencies (SFE=SFR/M$_{\rm H2}$) than spiral galaxies, but the SFEs derived are consistent with being drawn from the same distribution found in other lower mass ETG samples. This suggests that the SFE is not simply a function of stellar mass, but that local, internal processes are more important for regulating star formation. Finally we used the CO line profiles to investigate the high-mass end of the Tully-Fisher relation (TFR). We find that there is a break in the slope of the TFR for ETGs at high masses (consistent with previous studies). The strength of this break correlates with the stellar velocity dispersion of the host galaxies, suggesting it is caused by additional baryonic mass being present in the centre of massive ETGs. We speculate on the root cause of this change and its implications for galaxy formation theories.

The light curves of Type Ia supernovae (SNe Ia) are powered by the radioactive decay of $^{56}$Ni to $^{56}$Co at early times, and the decay of $^{56}$Co to $^{56}$Fe from ~60 days after explosion. We examine the evolution of the [Co III] 5892 A emission complex during the nebular phase for SNe Ia with multiple nebular spectra and show that the line flux follows the square of the mass of $^{56}$Co as a function of time. This result indicates both efficient local energy deposition from positrons produced in $^{56}$Co decay, and long-term stability of the ionization state of the nebula. We compile 77 nebular spectra of 25 SN Ia from the literature and present 17 new nebular spectra of 7 SNe Ia, including SN2014J. From these we measure the flux in the [Co III] 5892 A line and remove its well-behaved time dependence to infer the initial mass of $^{56}$Ni ($M_{Ni}$) produced in the explosion. We then examine $^{56}$Ni yields for different SN Ia ejected masses ($M_{ej}$ - calculated using the relation between light curve width and ejected mass) and find the $^{56}$Ni masses of SNe Ia fall into two regimes: for narrow light curves (low stretch s~0.7-0.9), $M_{Ni}$ is clustered near $M_{Ni}$ ~ 0.4$M_\odot$ and shows a shallow increase as $M_{ej}$ increases from ~1-1.4$M_\odot$; at high stretch, $M_{ej}$ clusters at the Chandrasekhar mass (1.4$M_\odot$) while $M_{Ni}$ spans a broad range from 0.6-1.2$M_\odot$. This could constitute evidence for two distinct SN Ia explosion mechanisms.

We examine the spheroid growth and star formation quenching experienced by galaxies from z~3 to the present by studying the evolution with redshift of the quiescent and spheroid-dominated fractions of galaxies from the CANDELS and GAMA surveys. We compare the observed fractions with predictions from a semi-analytic model which includes prescriptions for bulge growth and AGN feedback due to mergers and disk instabilities. We facilitate direct morphological comparison by converting our model bulge-to-total stellar mass ratios to Sersic indices. We then subdivide our population into the four quadrants of the sSFR-Sersic index plane and study the buildup of each of these subpopulations. We find that the fraction of star forming disks declines steadily, while the fraction of quiescent spheroids builds up over cosmic time. The fractions of star forming spheroids and quiescent disks are both non-negligible, and stay nearly constant over the period we have studied, at about 10% and 15-20% respectively. Our model is qualitatively successful at reproducing the evolution of the two "main" populations (star forming disk-dominated galaxies and quiescent spheroid-dominated galaxies), and approximately reproduces the relative fractions of all four types, but predicts a stronger decline in star forming spheroids, and increase in quiescent disks, than seen in the observations. A model with an additional channel for bulge growth via disk instabilities agrees better overall with the observations than a model in which bulges may grow only through mergers. We study evolutionary tracks of some individual galaxies as they experience morphological transformation and quenching, and examine the importance of different physical drivers of this transformation (major and minor mergers and disk instabilities). We find that complex histories with multiple transformative events are the norm.

We present UV through NIR broad-band photometry, and optical and NIR spectroscopy of Type Iax supernova 2012Z. The data set consists of both early and late-time observations, including the first late phase NIR spectrum obtained for a spectroscopically classified SN Iax. Simple model calculations of its bolometric light curve suggest SN 2012Z produced ~0.3 M_sun of (56)Ni, ejected about a Chandrasekhar mass of material, and had an explosion energy of ~10^51 erg, making it one of the brightest and most energetic SN Iax yet observed. The late phase NIR spectrum of SN 2012Z is found to broadly resemble similar epoch spectra of normal SNe Ia; however, like other SNe Iax, corresponding visual-wavelength spectra differ substantially compared to all supernova types. Constraints from the distribution of IMEs, e.g. silicon and magnesium, indicate that the outer ejecta did not experience significant mixing during or after burning, and the late phase NIR line profiles suggests most of the (56)Ni is produced during high density burning. The various observational properties of SN 2012Z are found to be consistent with the theoretical expectations of a Chandrasekhar mass white dwarf progenitor that experiences a pulsational delayed detonation, which produced several tenths of a solar mass of (56)Ni during the deflagration burning phase and little (or no) (56)Ni during the detonation phase. Within this scenario only a moderate amount of Rayleigh-Taylor mixing occurs both during the deflagration and fallback phase of the pulsation, and the layered structure of the IMEs is a product of the subsequent denotation phase. The fact that the SNe Iax population does not follow a tight brightness-decline relation similar to SNe Ia can then be understood in the framework of variable amounts of mixing during pulsational rebound and variable amounts of (56)Ni production during the early subsonic phase of expansion.