Search

The growing volume of data produced by large astronomical surveys necessitates the development of efficient analysis techniques capable of effectively managing high-dimensional datasets. This study addresses this need by demonstrating some applications of manifold learning and dimensionality reduction techniques, specifically the Self-Organizing Map (SOM), on the optical+NIR SED space of galaxies, with a focus on sample comparison, selection biases, and predictive power using a small subset. To this end, we utilize a large photometric sample from the five CANDELS fields and a subset with spectroscopic measurements from the KECK MOSDEF survey in two redshift bins at $z\sim1.5$ and $z\sim2.2$. We trained SOM with the photometric data and mapped the spectroscopic data onto it as our study case. We found that MOSDEF targets do not cover all SED shapes existing in the SOM. Our findings reveal that Active Galactic Nuclei (AGN) within the MOSDEF sample are mapped onto the more massive regions of the SOM, confirming previous studies and known selection biases towards higher-mass, less dusty galaxies. Furthermore, SOM were utilized to map measured spectroscopic features, examining the relationship between metallicity variations and galaxy mass. Our analysis confirmed that more massive galaxies exhibit lower [OIII]/H$\beta$ and [OIII]/[OII] ratios and higher H$\alpha$/H$\beta$ ratios, consistent with the known mass-metallicity relation. These findings highlight the effectiveness of SOM in analyzing and visualizing complex, multi-dimensional datasets, emphasizing their potential in data-driven astronomical studies.

We examine star-formation and dust properties for a sample of 660 galaxies at $1.37\leq z\leq 2.61$ in the MOSDEF survey by dividing them into groups with similarly-shaped spectral energy distributions (SEDs). For each group, we combine the galaxy photometry into a finely-sampled composite SED, and stack their spectra. This method enables the study of more complete galaxy samples, including galaxies with very faint emission lines. We fit these composite SEDs with Prospector to measure the stellar attenuation and SED-based star-formation rates (SFRs). We also derive emission-line properties from the spectral stacks, including Balmer decrements, dust-corrected SFRs, and metallicities. We find that stellar attenuation correlates most strongly with mass, while nebular attenuation correlates strongly with both mass and SFR. Furthermore, the excess of nebular compared to stellar attenuation correlates most strongly with SFR. The highest SFR group has 2 mag of excess nebular attenuation. Our results are consistent with a model in which star-forming regions become more dusty as galaxy mass increases. To explain the increasing excess nebular attenuation, we require a progressively larger fraction of star formation to occur in highly-obscured regions with increasing SFR. This highly-obscured star formation could occur in dusty clumps or central starbursts. Additionally, as each galaxy group represents a different evolutionary stage, we study their locations on the UVJ and SFR-mass diagrams. As mass increases, metallicity and dust attenuation increase, while sSFR decreases. However, the most massive group moves towards the quiescent region of the UVJ diagram, while showing less obscuration, potentially indicating removal of dust., Comment: 21 pages, 7 figures, accepted for publication in ApJ

High-velocity outflows are ubiquitous in compact, massive (M$_* \sim$ 10$^{11}$ M$_{\odot}$), z $\sim$ 0.5 galaxies with extreme star formation surface densities ($\Sigma_{SFR} \sim$ 2000 M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$). We have previously detected and characterized these outflows using MgII absorption lines. To probe their full extent, we present Keck/KCWI integral field spectroscopy of the [OII] and MgII emission nebulae surrounding all of the 12 galaxies in this study. We find that [OII] is more effective than MgII in tracing low surface brightness, extended emission in these galaxies. The [OII] nebulae are spatially extended beyond the stars, with radial extent R$_{90}$ between 10 and 40 kpc. The nebulae exhibit non-gravitational motions, indicating galactic outflows with maximum blueshifted velocities ranging from -335 to -1920 km s$^{-1}$. The outflow kinematics correlate with the bursty star formation histories of these galaxies. Galaxies with the most recent bursts of star formation (within the last $<$ 3 Myr) exhibit the highest central velocity dispersions ($\sigma >$ 400 km s$^{-1}$), while the oldest bursts have the lowest-velocity outflows. Many galaxies exhibit both high-velocity cores and more extended, slower-moving gas indicative of multiple outflow episodes. The slower, larger outflows occurred earlier and have decelerated as they propagate into the CGM and mix on timescales $>$ 50 Myr., Comment: Accepted for publication in The Astrophysical Journal

Supermassive Black Holes (BHs) are known to efficiently grow through gas accretion, but even sustained and intense mass build-up through this mechanism struggles to explain the assembly of the most massive BHs observed in the local Universe. Using the Chandra Deep-Wide Field Survey (CDFWS) in the Bo\"otes field, we measure BH-galaxy assembly in massive galaxies ($M_\star\gtrsim10^{10}\, \rm M_\odot$) through the AGN fraction and specific Black Hole accretion rate (sBHAR) distribution as a function of redshift and stellar mass. We determine stellar masses and star formation rates for a parent sample of optically selected galaxies as well as those with X-ray detections indicating the presence of an AGN through Spectral Energy Distribution (SED) fitting. We derive a redshift-dependent mass completeness limit and extract X-ray information for every galaxy as to identify AGN in various types of galaxies at low and high redshift providing a comprehensive picture of the AGN population in massive galaxies. While X-ray AGN samples are dominated by moderately massive host galaxies of $M_{\star} \geqslant 10^{10}\, \rm M_{\odot}$, we do not find a strong stellar mass dependence in AGN fraction (to limits in sBHAR) or sBHAR distribution above this mass threshold, indicating a bias towards massive galaxies in the observed samples. A mild increase in AGN fraction is seen towards high stellar mass for low sBHAR AGN, as well as a suppression of high sBHAR events in the most massive galaxies. We derive BH-galaxy growth tracks over time, which reveal that while most BH mass has been accumulated since $z=4$ for lower mass BHs, the assembly of the most massive BHs is more complex, with little to no relative mass gain since $z=4$, implying that rapid and intense growth episodes prior to $z=4$ were necessary to form these massive BHs., Comment: 20 pages, 11 figures, submitted to MNRAS

We present Keck Cosmic Web Imager IFU observations around extended Ly$\alpha$ halos of 27 typical star-forming galaxies with redshifts $2.0 < z < 3.2$ drawn from the MOSFIRE Deep Evolution Field survey. We examine the average Ly$\alpha$ surface-brightness profiles in bins of star-formation rate (SFR), stellar mass ($M_*$), age, stellar continuum reddening, SFR surface density ($\rm \Sigma_{SFR}$), and $\rm \Sigma_{SFR}$ normalized by stellar mass ($\rm \Sigma_{sSFR}$). The scale lengths of the halos correlate with stellar mass, age, and stellar continuum reddening; and anti-correlate with star-formation rate, $\rm \Sigma_{SFR}$, and $\rm \Sigma_{sSFR}$. These results are consistent with a scenario in which the down-the-barrel fraction of Ly$\alpha$ emission is modulated by the low-column-density channels in the ISM, and that the neutral gas covering fraction is related to the physical properties of the galaxies. Specifically, we find that this covering fraction increases with stellar mass, age, and $E(B-V)$; and decreases with SFR, $\rm \Sigma_{SFR}$ and $\rm \Sigma_{sSFR}$. We also find that the resonantly scattered Ly$\alpha$ emission suffers greater attenuation than the (non-resonant) stellar continuum emission, and that the difference in attenuation increases with stellar mass, age, and stellar continuum reddening, and decreases with $\rm \Sigma_{sSFR}$. These results imply that more reddened galaxies have more dust in their CGM., Comment: 12 pages, 7 figures, 2 tables, submitted to ApJ

We use the large spectroscopic data set of the MOSFIRE Deep Evolution Field survey to investigate the kinematics and energetics of ionised gas outflows. Using a sample of 598 star-forming galaxies at redshift 1.4 < $z$ < 3.8, we decompose $\rm{H}\alpha$ and [OIII] emission lines into narrow and broad components, finding significant detections of broad components in 10% of the sample. The ionised outflow velocity from individual galaxies appears independent of galaxy properties, such as stellar mass, star-formation rate (SFR), and star-formation-rate surface density ($\Sigma_{\rm SFR}$). Adopting a simple outflow model, we estimate the mass-, energy- and momentum-loading factors of the ionised outflows, finding modest values with averages of 0.33, 0.04, and 0.22, respectively. The larger momentum- than energy-loading factors, for the adopted physical parameters, imply that these ionised outflows are primarily momentum-driven. We further find a marginal correlation (2.5$\sigma$) between the mass-loading factor and stellar mass in agreement with predictions by simulations, scaling as $\eta_{m}$ $\propto M_{\star}^{-0.45}$. This shallow scaling relation is consistent with these ionised outflows being driven by a combination of mechanical energy generated by supernovae explosions and radiation pressure acting on dusty material. In a majority of galaxies, the outflowing material does not appear to have sufficient velocity to escape the gravitational potential of their host, likely recycling back at later times. Together, these results suggest that the ionised outflows traced by nebular emission lines are negligible, with the bulk of mass and energy carried out in other gaseous phases., Comment: 17 pages, 9 figures, Submitted to Monthly Notices of the Royal Astronomical Society

[OII] emission maps obtained with the Keck Cosmic Web Imager (KCWI) are presented for four galaxies centered in cooling X-ray cluster atmospheres. Nebular emission extending tens of kpc is found in systems covering a broad range of atmospheric cooling rates, cluster masses, and dynamical states. Abell 262's central galaxy hosts a kpc-scale disk. The nebular gas in RXJ0820.9+0752 is offset and redshifted with respect to the central galaxy by $10-20$ kpc and 150 km s$^{-1}$, respectively. The nebular gases in PKS 0745-191 and Abell 1835 are being churned to higher velocity dispersion by X-ray bubbles and jets. The churned gas is enveloped by larger scale, lower velocity dispersion (quiescent) nebular emission. The mean line-of-sight speeds of the churned gas, quiescent gas, and the central galaxy each differ by up to $\sim 150$ km s$^{-1}$; nebular speeds upward of $800$ km s$^{-1}$ are found. Gases with outwardly-rising speeds upward of several hundred km s$^{-1}$ are consistent with being advected behind and being lifted by the rising bubbles. The peculiar motion between the galaxy, nebular gas, and perhaps the hot atmosphere from which it presumably condensed is affecting the bubble dynamics, and may strongly affect thermally unstable cooling, the dispersal of jet energy, and the angular momentum of gas accreting onto the galaxies and their nuclear black holes., Comment: 36 pages, 26 figures, Accepted for publication in ApJ

A new class of extragalactic astronomical sources discovered in 2021, named Odd Radio Circles (ORCs, Norris et al. 2021), are large rings of faint, diffuse radio continuum emission spanning ~1 arcminute on the sky. Galaxies at the centers of several ORCs have photometric redshifts of z~0.3-0.6, implying physical scales of several 100 kiloparsecs in diameter for the radio emission, the origin of which is unknown. Here we report spectroscopic data on an ORC including strong [OII] emission tracing ionized gas in the central galaxy of ORC4 at z=0.4512. The physical extent of the [OII] emission is ~40 kpc in diameter, larger than expected for a typical early-type galaxy (Pandya et al, 2017) but an order of magnitude smaller than the large-scale radio continuum emission. We detect a ~200 km/s velocity gradient across the [OII] nebula, as well as a high velocity dispersion of ~180 km/s. The [OII] equivalent width (EW, ~50 Ang) is extremely high for a quiescent galaxy. The morphology, kinematics, and strength of the [OII] emission are consistent with the infall of shock ionized gas near the galaxy, following a larger-scale, outward moving shock driven by a galactic wind. Both the extended optical and radio emission, while observed on very different scales, may therefore result from the same dramatic event., Comment: 7 figures, accepted to Nature

We report on the discovery of cool gas inflows towards three star-forming galaxies at $\left\sim$ 2.30. Analysis of Keck Low-Resolution Imaging Spectrometer spectroscopy reveals redshifted low-ionisation interstellar (LIS) metal absorption lines with centroid velocities of 60 - 130 km $\rm{s}^{-1}$. These inflows represent some of the most robust detections of inflowing gas into isolated, star-forming galaxies at high redshift. Our analysis suggests that the inflows are due to recycling metal-enriched gas from previous ejections. Comparisons between the galaxies with inflows and a larger parent sample of 131 objects indicate that galaxies with detected inflows may have higher specific star-formation rates (sSFR) and star-formation-rate surface densities. However, when additional galaxies without robustly detected inflows based on centroid velocity but whose LIS absorption line profiles indicate large red-wing velocities are considered, galaxies with inflows do not show unique properties relative to those lacking inflows. Additionally, we calculate the covering fraction of cool inflowing gas as a function of red-wing inflow velocity, finding an enhancement in high sSFR binned galaxies, likely due to an increase in the amount of recycling gas. Together, these results suggest that the low detection rate of galaxies with cool inflows is primarily related to the viewing angle rather than the physical properties of the galaxies., Comment: 10 pages, 7 figures, accepted for publication in MNRAS

We investigate dust attenuation and its dependence on viewing angle for 308 star-forming galaxies at $1.3\leq z\leq2.6$ from the MOSFIRE Deep Evolution Field (MOSDEF) survey. We divide galaxies with a detected H$\alpha$ emission line and coverage of H$\beta$ into eight groups by stellar mass, star formation rate (SFR), and inclination (i.e., axis ratio), then stack their spectra. From each stack, we measure Balmer decrement and gas-phase metallicity, then we compute median \AV and UV continuum spectral slope ($\beta$). First, we find that none of the dust properties (Balmer decrement, \AV, $\beta$) vary with axis ratio. Second, both stellar and nebular attenuation increase with increasing galaxy mass, showing little residual dependence on SFR or metallicity. Third, nebular emission is more attenuated than stellar emission, and this difference grows even larger at higher galaxy masses and SFRs. Based on these results, we propose a three-component dust model where attenuation predominantly occurs in star-forming regions and large, dusty star-forming clumps, with minimal attenuation in the diffuse ISM. In this model, nebular attenuation primarily originates in clumps, while stellar attenuation is dominated by star-forming regions. Clumps become larger and more common with increasing galaxy mass, creating the above mass trends. Finally, we argue that a fixed metal yield naturally leads to mass regulating dust attenuation. Infall of low-metallicity gas increases SFR and lowers metallicity, but leaves dust column density mostly unchanged. We quantify this idea using the Kennicutt-Schmidt and fundamental metallicity relations, showing that galaxy mass is indeed the primary driver of dust attenuation., Comment: 19 pages, 13 figures, accepted for publication in ApJ

This is the fifth paper in a series of investigations of the clustering properties of luminous, broad-emission-line active galactic nuclei (AGN) identified in the ROSAT All-Sky Survey (RASS) and Sloan Digital Sky Survey (SDSS). In this work we measure the cross-correlation function (CCF) between RASS/SDSS DR14 AGN with the SDSS CMASS galaxy sample at $0.44

We present results on the properties of extreme gas outflows in massive ($\rm M_* \sim$10$^{11} \ \rm M_{\odot}$), compact, starburst ($\rm SFR \sim$$200 \, \rm M_{\odot} \ yr^{-1}$) galaxies at z = $0.4-0.7$ with very high star formation surface densities ($\rm \Sigma_{SFR} \sim$$2000 \,\rm M_{\odot} \ yr^{-1} \ kpc^{-2}$). Using optical Keck/HIRES spectroscopy of 14 HizEA starburst galaxies we identify outflows with maximum velocities of $820 - 2860$ \kmps. High-resolution spectroscopy allows us to measure precise column densities and covering fractions as a function of outflow velocity and characterize the kinematics and structure of the cool gas outflow phase (T $\sim$10$^4$ K). We find substantial variation in the absorption profiles, which likely reflects the complex morphology of inhomogeneously-distributed, clumpy gas and the intricacy of the turbulent mixing layers between the cold and hot outflow phases. There is not a straightforward correlation between the bursts in the galaxies' star formation histories and their wind absorption line profiles, as might naively be expected for starburst-driven winds. The lack of strong \mgii \ absorption at the systemic velocity is likely an orientation effect, where the observations are down the axis of a blowout. We infer high mass outflow rates of $\rm \sim$50 $-$ 2200 $\rm M_{\odot} \, yr^{-1}$, assuming a fiducial outflow size of 5 kpc, and mass loading factors of $\eta\sim$5 for most of the sample. %with $\eta\sim$20 for two galaxies. While these values have high uncertainties, they suggest that starburst galaxies are capable of ejecting very large amounts of cool gas that will substantially impact their future evolution., Comment: Accepted for publication in The Astrophysical Journal

The Makani galaxy hosts the poster child of a galactic wind on scales of the circumgalactic medium. It consists of a two-episode wind in which the slow, outer wind originated 400 Myr ago (Episode I; R_I = 20-50 kpc) and the fast, inner wind is 7 Myr old (Episode II; R_II = 0-20 kpc). While this wind contains ionized, neutral, and molecular gas, the physical state and mass of the most extended phase--the warm, ionized gas--is unknown. Here we present Keck optical spectra of the Makani outflow. These allow us to detect hydrogen lines out to r = 30-40 kpc and thus constrain the mass, momentum, and energy in the wind. Many collisionally-excited lines are detected throughout the wind, and their line ratios are consistent with 200-400 km/s shocks that power the ionized gas, with v_shock = $\sigma$_wind. Combining shock models, density-sensitive line ratios, and mass and velocity measurements, we estimate that the ionized mass and outflow rate in the Episode II wind could be as high as that of the molecular gas: M_II(HII) ~ M_II(H_2) = (1-2)x10^9 Msun and dM/dt_II(HII) ~ dM/dt_II(H_2) = 170-250 Msun/yr. The outer wind has slowed, so that dM/dt_I(HII) ~ 10 Msun/yr, but it contains more ionized gas: M_I(HII) = 5x10^9 Msun. The momentum and energy in the recent Episode II wind imply a momentum-driven flow (p ``boost" ~ 7) driven by the hot ejecta and radiation pressure from the Eddington-limited, compact starburst. Much of the energy and momentum in the older Episode I wind may reside in a hotter phase, or lie further into the CGM., Comment: ApJ, in press

We use the large spectroscopic dataset of the MOSFIRE Deep Evolution Field (MOSDEF) survey to investigate some of the key factors responsible for the elevated ionization parameters (U) inferred for high-redshift galaxies, focusing in particular on the role of star-formation-rate surface density (Sigma_SFR). Using a sample of 317 galaxies with spectroscopic redshifts z~1.9-3.7, we construct composite rest-frame optical spectra in bins of Sigma_SFR and infer electron densities, n_e, using the ratio of the [OII] 3727, 3730 doublet. Our analysis suggests a significant (~3 sigma) correlation between n_e and Sigma_SFR. We further find significant correlations between U and Sigma_SFR for composite spectra of a subsample of 113 galaxies, and for a smaller sample of 25 individual galaxies with inferences of U. The increase in n_e -- and possibly also the volume filling factor of dense clumps in HII regions -- with Sigma_SFR appear to be important factors in explaining the relationship between U and Sigma_SFR. Further, the increase in n_e and SFR with redshift at a fixed stellar mass can account for most of the redshift evolution of U. These results suggest that the gas density, which sets n_e and the overall level of star-formation activity, may play a more important role than metallicity evolution in explaining the elevated ionization parameters of high-redshift galaxies., Comment: 11 pages, 9 figures, 1 table, submitted to ApJ

We present a measurement of the intrinsic space density of intermediate redshift ($z\sim0.5$), massive ($M_{*} \sim 10^{11} \ \text{M}_{\odot}$), compact ($R_{e} \sim 100$ pc) starburst ($\Sigma_{SFR} \sim 1000 \ \text{M}_{\odot} \ \text{yr}^{-1} \text{kpc}^{-1}$) galaxies with tidal features indicative of them having undergone recent major mergers. A subset of them host kiloparsec scale, $>1000 \ \text{km}\ \text{s}^{-1}$ outflows and have little indication of AGN activity, suggesting that extreme star formation can be a primary driver of large-scale feedback. The aim for this paper is to calculate their space density so we can place them in a better cosmological context. We do this by empirically modeling the stellar populations of massive, compact starburst galaxies. We determine the average timescale for which galaxies that have recently undergone an extreme nuclear starburst would be targeted and included in our spectroscopically selected sample. We find that massive, compact starburst galaxies targeted by our criteria would be selectable for $\sim 148 ^{+27}_{-24}$ Myr and have an intrinsic space density $n_{\text{CS}} \sim (1.1^{+0.5}_{-0.3}) \times 10^{-6} \ \ \text{Mpc}^{-3}$. This space density is broadly consistent with our $z\sim0.5$ compact starbursts being the most extremely compact and star forming low redshift analogs of the compact star forming galaxies in the early Universe as well as them being the progenitors to a fraction of intermediate redshift post starburst and compact quiescent galaxies., Comment: 27 pages, 13 figures, accepted for publication in the Astrophysical Journal

We define a new morphology metric called "patchiness" ($P$) that is sensitive to deviations from the average of a resolved distribution, does not require the galaxy center to be defined, and can be used on the spatially-resolved distribution of any galaxy property. While the patchiness metric has a broad range of applications, we demonstrate its utility by investigating the distribution of dust in the interstellar medium of 310 star-forming galaxies at spectroscopic redshifts $1.36

We analyze the completeness of the MOSDEF survey, in which z ~ 2 galaxies were selected for rest-optical spectroscopy from well-studied HST extragalactic legacy fields down to a fixed rest-optical magnitude limit (H_AB = 24.5). The subset of z ~ 2 MOSDEF galaxies with high signal-to-noise (S/N) emission-line detections analyzed in previous work represents a small minority (<10%) of possible z ~ 2 MOSDEF targets. It is therefore crucial to understand how representative this high S/N subsample is, while also more fully exploiting the MOSDEF spectroscopic sample. Using spectral-energy-distribution (SED) models and rest-optical spectral stacking, we compare the MOSDEF z ~ 2 high S/N subsample with the full MOSDEF sample of z ~ 2 star-forming galaxies with redshifts, the latter representing an increase in sample size of more than a factor of three. We find that both samples have similar emission-line properties, in particular in terms of the magnitude of the offset from the local star-forming sequence on the [N II] BPT diagram. There are small differences in median host galaxy properties, including the stellar mass (M_*), star-formation rate (SFR) and specific SFR (sSFR), and UVJ colors; however, these offsets are minor considering the wide spread of the distributions. Using SED modeling, we also demonstrate that the sample of z ~ 2 star-forming galaxies observed by the MOSDEF survey is representative of the parent catalog of available such targets. We conclude that previous MOSDEF results on the evolution of star-forming galaxy emission-line properties were unbiased relative to the parent z ~ 2 galaxy population., Comment: 18 pages, 10 figures, 3 tables

We present a detailed study of a galaxy merger taking place at $z=1.89$ in the GOODS-S field. Here we analyze Keck/MOSFIRE spectroscopic observations from the MOSFIRE Deep Evolution Field (MOSDEF) survey along with multi-wavelength photometry assembled by the 3D-HST survey. The combined dataset is modeled to infer the past star-formation histories (SFHs) of both merging galaxies. They are found to be massive, with log$_{10}(M_{\ast}/M_{\odot}) > 11$, with a close mass ratio satisfying the typical major-merger definition. Additionally, in the context of delayed-$\tau$ models, GOODS-S 43114 and GOODS-S 43683 have similar SFHs and low star-formation rates (log$_{10}$(SFR(SED)/$M_{\odot}/\rm{yr}^{-1}$) $<$ 1.0) compared to their past averages. The best-fit model SEDs show elevated H$\delta_{\rm{A}}$ values for both galaxies, indicating that their stellar spectra are dominated by A-type stars, and that star formation peaked $\sim0.5-1$ Gyr ago and has recently declined. Additionally, based on SED fitting both merging galaxies turned on and shut off star formation within a few hundred Myr of each other, suggesting that their bursts of star formation may be linked. Combining the SFHs and H$\delta_{\rm{A}}$ results with recent galaxy merger simulations, we infer that these galaxies have recently completed their first pericentric passage and are moving apart. Finally, the relatively low second velocity moment of GOODS-S 43114 given its stellar mass, suggests a disk-like structure. However, including the geometry of the galaxy in the modeling does not completely resolve the discrepancy between the dynamical and stellar masses. Future work is needed to resolve this inconsistency in mass., Comment: 12 pages, 7 figures, 1 table, Accepted for publication in MNRAS

We aim to quantify the relation between the dust-to-gas mass ratio (DTG) and gas-phase metallicity of $z=$2.1-2.5 luminous galaxies and contrast this high-redshift relation against analogous constraints at z$=$0. We present a sample of ten star-forming main-sequence galaxies in the redshift range $2.1

We present observations of CO(3-2) in 13 main-sequence $z=2.0-2.5$ star-forming galaxies at $\log(M_*/M_{\odot})=10.2-10.6$ that span a wide range in metallicity (O/H) based on rest-optical spectroscopy. We find that CO(3-2)/SFR decreases with decreasing metallicity, implying that the CO luminosity per unit gas mass is lower in low-metallicity galaxies at $z\sim2$. We constrain the CO-to-H$_2$ conversion factor ($\alpha_{\text{CO}}$) and find that $\alpha_{\text{CO}}$ inversely correlates with metallicity at $z\sim2$. We derive molecular gas masses ($M_{\text{mol}}$) and characterize the relations among $M_*$, SFR, $M_{\text{mol}}$, and metallicity. At $z\sim2$, $M_{\text{mol}}$ increases and molecular gas fraction ($M_{\text{mol}}$/$M_*$) decrease with increasing $M_*$, with a significant secondary dependence on SFR. Galaxies at $z\sim2$ lie on a near-linear molecular KS law that is well-described by a constant depletion time of 700 Myr. We find that the scatter about the mean SFR-$M_*$, O/H-$M_*$, and $M_{\text{mol}}$-$M_*$ relations is correlated such that, at fixed $M_*$, $z\sim2$ galaxies with larger $M_{\text{mol}}$ have higher SFR and lower O/H. We thus confirm the existence of a fundamental metallicity relation at $z\sim2$ where O/H is inversely correlated with both SFR and $M_{\text{mol}}$ at fixed $M_*$. These results suggest that the scatter of the $z\sim2$ star-forming main sequence, mass-metallicity relation, and $M_{\text{mol}}$-$M_*$ relation are primarily driven by stochastic variations in gas inflow rates. We place constraints on the mass loading of galactic outflows and perform a metal budget analysis, finding that massive $z\sim2$ star-forming galaxies retain only 30% of metals produced, implying that a large mass of metals resides in the circumgalactic medium., Comment: 35 pages, 17 figures, submitted to ApJ

The extent of black hole growth during different galaxy evolution phases and the connection between galaxy compactness and AGN activity remain poorly understood. We use Hubble Space Telescope imaging of the CANDELS fields to identify star-forming and quiescent galaxies at z=0.5-3 in both compact and extended phases and use Chandra X-ray imaging to measure the distribution of AGN accretion rates and track black hole growth within these galaxies. Accounting for the impact of AGN light changes ~20% of the X-ray sources from compact to extended galaxy classifications. We find that ~10-25% of compact star-forming galaxies host an AGN, a mild enhancement (by a factor ~2) compared to extended star-forming galaxies or compact quiescent galaxies of equivalent stellar mass and redshift. However, AGN are not ubiquitous in compact star-forming galaxies and this is not the evolutionary phase, given its relatively short timescale, where the bulk of black hole mass growth takes place. Conversely, we measure the highest AGN fractions (~10-30%) within the relatively rare population of extended quiescent galaxies. For massive galaxies that quench at early cosmic epochs, substantial black hole growth in this extended phase is crucial to produce the elevated black hole mass-to-galaxy stellar mass scaling relation observed for quiescent galaxies at z~0. We also show that AGN fraction increases with compactness in star-forming galaxies and decreases in quiescent galaxies within both the compact and extended sub-populations, demonstrating that AGN activity depends closely on the structural properties of galaxies., Comment: 30 pages, 18 figures, accepted for publication in MNRAS. Primary results are shown in Fig 7 and summarised by Fig 11. See Fig 15 and 16 for key interpretation/conclusions

The combination of the MOSDEF and KBSS-MOSFIRE surveys represents the largest joint investment of Keck/MOSFIRE time to date, with ~3000 galaxies at 1.4<=z<=3.8, roughly half of which are at z~2. MOSDEF is photometric- and spectroscopic-redshift selected with a rest-optical magnitude limit, while KBSS-MOSFIRE is primarily selected based on rest-UV colors and a rest-UV magnitude limit. Analyzing both surveys in a uniform manner with consistent spectral-energy-distribution (SED) models, we find that the MOSDEF z~2 targeted sample has a higher median M_* and redder rest U-V color than the KBSS-MOSFIRE z~2 targeted sample, and a smaller median SED-based SFR and sSFR (SFR(SED) and sSFR(SED)). Specifically, MOSDEF targeted a larger population of red galaxies with U-V and V-J >=1.25, while KBSS-MOSFIRE contains more young galaxies with intense star formation. Despite these differences in the z~2 targeted samples, the subsets of the surveys with multiple emission lines detected and analyzed in previously published work are much more similar. All median host-galaxy properties with the exception of stellar population age -- i.e., M_*, SFR(SED), sSFR(SED), A_V, and UVJ colors -- agree within the uncertainties. Additionally, when uniform emission-line fitting and stellar Balmer absorption correction techniques are applied, there is no significant offset between the two samples in the [OIII]$\lambda$5008/H$\beta$ vs. [NII]$\lambda$6585/H$\alpha$ diagnostic diagram, in contrast to previously-reported discrepancies. We can now combine the MOSDEF and KBSS-MOSFIRE surveys to form the largest z~2 sample with moderate-resolution rest-optical spectra and construct the fundamental scaling relations of star-forming galaxies during this important epoch., Comment: 23 pages, 14 figures, published in MNRAS

We perform an aperture-matched analysis of dust-corrected H$\alpha$ and UV SFRs using 303 star-forming galaxies with spectroscopic redshifts $1.36

We perform joint modeling of the composite rest-frame far-UV (FUV) and optical spectra of redshift 1.85

We present results on the nature of extreme ejective feedback episodes and the physical conditions of a population of massive ($\rm M_* \sim 10^{11} M_{\odot}$), compact starburst galaxies at z = 0.4-0.7. We use data from Keck/NIRSPEC, SDSS, Gemini/GMOS, MMT, and Magellan/MagE to measure rest-frame optical and near-IR spectra of 14 starburst galaxies with extremely high star formation rate surface densities (mean $\rm \Sigma_{SFR} \sim 3000 \,M_{\odot} yr^{-1} kpc^{-2}$) and powerful galactic outflows (maximum speeds v$_{98} \sim$ 1000-3000 km s$^{-1}$). Our unique data set includes an ensemble of both emission [OII]$\lambda\lambda$3726,3729, H$\beta$, [OIII]$\lambda\lambda$4959,5007, H$\alpha$, [NII]$\lambda\lambda$6548,6583, and [SII]$\lambda\lambda$6716,6731) and absorption MgII$\lambda\lambda$2796,2803, and FeII$\lambda$2586) lines that allow us to investigate the kinematics of the cool gas phase (T$\sim$10$^4$ K) in the outflows. Employing a suite of line ratio diagnostic diagrams, we find that the central starbursts are characterized by high electron densities (median n$_e \sim$ 530 cm$^{-3}$), and high metallicity (solar or super-solar). We show that the outflows are most likely driven by stellar feedback emerging from the extreme central starburst, rather than by an AGN. We also present multiple intriguing observational signatures suggesting that these galaxies may have substantial Lyman continuum (LyC) photon leakage, including weak [SII] nebular emission lines. Our results imply that these galaxies may be captured in a short-lived phase of extreme star formation and feedback where much of their gas is violently blown out by powerful outflows that open up channels for LyC photons to escape.

The complex structure of gas, metals, and dust in the interstellar and circumgalactic medium (ISM and CGM, respectively) in star-forming galaxies can be probed by Ly$\alpha$ emission and absorption, low-ionization interstellar (LIS) metal absorption, and dust reddening E(B-V). We present a statistical analysis of the mutual correlations among Ly$\alpha$ equivalent width (EW$_{Ly\alpha}$), LIS equivalent width (EW$_{LIS}$), and E(B-V) in a sample of 157 star-forming galaxies at $z\sim2.3$. With measurements obtained from individual, deep rest-UV spectra and spectral-energy distribution (SED) modeling, we find that the tightest correlation exists between EW$_{LIS}$ and E(B-V), although correlations among all three parameters are statistically significant. These results signal a direct connection between dust and metal-enriched HI gas, and that they are likely co-spatial. By comparing our results with the predictions of different ISM/CGM models, we favor a dusty ISM/CGM model where dust resides in HI gas clumps and Ly$\alpha$ photons escape through the low HI covering fraction/column density intra-clump medium. Finally, we investigate the factors that potentially contribute to the intrinsic scatter in the correlations studied in this work, including metallicity, outflow kinematics, Ly$\alpha$ production efficiency, and slit loss. Specifically, we find evidence that scatter in the relationship between EW$_{Ly\alpha}$ and E(B-V) reflects the variation in metal-to-HI covering fraction ratio as a function of metallicity, and the effects of outflows on the porosity of the ISM/CGM. Future simulations incorporating star-formation feedback and the radiative transfer of Ly$\alpha$ photons will provide key constraints on the spatial distributions of neutral hydrogen gas and dust in the ISM/CGM structure., Comment: 19 pages, 6 figures, 3 tables, accepted for publication in ApJ

We analyze the rest-optical emission-line ratios of z~1.5 galaxies drawn from the MOSFIRE Deep Evolution Field (MOSDEF) survey. Using composite spectra we investigate the mass-metallicity relation (MZR) at z~1.5 and measure its evolution to z=0. When using gas-phase metallicities based on the N2 line ratio, we find that the MZR evolution from z~1.5 to z=0 depends on stellar mass, evolving by $\Delta\rm log(\rm O/H)\sim0.25$ dex at $M_*<10^{9.75}M_{\odot}$ down to $\Delta\rm log(\rm O/H)\sim0.05$ at $M_*>10^{10.5}M_{\odot}$. In contrast, the O3N2-based MZR shows a constant offset of $\Delta\rm log(\rm O/H)\sim0.30$ across all masses, consistent with previous MOSDEF results based on independent metallicity indicators, and suggesting that O3N2 provides a more robust metallicity calibration for our z~1.5 sample. We investigated the secondary dependence of the MZR on SFR by measuring correlated scatter about the mean $M_*$-specific SFR and $M_*-\log(\rm O3N2)$ relations. We find an anti-correlation between $\log(\rm O/H)$ and sSFR offsets, indicating the presence of a $M_*$-SFR-Z relation, though with limited significance. Additionally, we find that our z~1.5 stacks lie along the z=0 metallicity sequence at fixed $\mu=\log(M_*/M_{\odot})-0.6\times\log(\rm SFR / M_{\odot} yr^{-1})$ suggesting that the z~1.5 stacks can be described by the z=0 fundamental metallicity relation (FMR). However, using different calibrations can shift the calculated metallicities off of the local FMR, indicating that appropriate calibrations are essential for understanding metallicity evolution with redshift. Finally, understanding how [NII]/H$\alpha$ scales with galaxy properties is crucial to accurately describe the effects of blended [NII] and H$\alpha$ on redshift and H$\alpha$ flux measurements in future large surveys utilizing low-resolution spectra such as with Euclid and the Roman Space Telescope., Comment: 14 Pages, 9 figures, Accepted to MNRAS

We present multi-band Hubble Space Telescope imaging that spans rest-frame near-ultraviolet through near-infrared wavelengths (0.3-1.1 $\mu$m) for 12 compact starburst galaxies at z=0.4-0.8. These massive galaxies (M_stellar ~ 10^11 M_Sun) are driving very fast outflows ($v_{max}$=1000-3000 km/s), and their light profiles are dominated by an extremely compact starburst component (half-light radius ~ 100 pc). Our goal is to constrain the physical mechanisms responsible for launching these fast outflows by measuring the physical conditions within the central kiloparsec. Based on our stellar population analysis, the central component typically contributes $\approx$25% of the total stellar mass and the central escape velocities $v_{esc,central}\approx900$ km/s are a factor of two smaller than the observed outflow velocities. This requires physical mechanisms that can accelerate gas to speeds significantly beyond the central escape velocities, and it makes clear that these fast outflows are capable of traveling into the circumgalactic medium, and potentially beyond. We find central stellar densities comparable to theoretical estimates of the Eddington limit, and we estimate $\Sigma_1$ surface densities within the central kpc comparable to those of compact massive galaxies at $0.5

Using the near-IR spectroscopy of the MOSFIRE Deep Evolution Field (MOSDEF) survey, we investigate the role of local environment in the gas-phase metallicity of galaxies. The local environment measurements are derived from accurate and uniformly calculated photometric redshifts with well-calibrated probability distributions. Based on rest-frame optical emission lines, [NII]$\lambda6584$ and H$\alpha$, we measure gas-phase oxygen abundance of 167 galaxies at $1.37\leq z\leq1.7$ and 303 galaxies at $2.09\leq z\leq2.61$, located in diverse environments. We find that at $z\sim1.5$, the average metallicity of galaxies in overdensities with $M_*\sim10^{9.8}M_\odot, 10^{10.2}M_\odot$ and $10^{10.8}M_\odot$ is higher relative to their field counterparts by $0.094\pm0.051$, $0.068\pm0.028$ and $0.052\pm0.043$ dex, respectively. However, this metallicity enhancement does not exist at higher redshift, $z\sim2.3$, where, compared to the field galaxies, we find $0.056\pm0.043$, $0.056\pm0.028$ and $0.096\pm 0.034$ dex lower metallicity for galaxies in overdense environments with $M_*\sim10^{9.8}M_\odot, 10^{10.2}M_\odot$ and $10^{10.7}M_\odot$, respectively. Our results suggest that, at $1.37\leq z\leq2.61$, the variation of mass-metallicity relation with local environment is small ($<0.1$dex), and reverses at $z\sim2$. Our results support the hypothesis that, at the early stages of cluster formation, owing to efficient gas cooling, galaxies residing in overdensities host a higher fraction of pristine gas with prominent primordial gas accretion, which lowers their gas-phase metallicity compared to their coeval field galaxies. However, as the Universe evolves to lower redshifts ($z\lesssim2$), shock-heated gas in overdensities cannot cool down efficiently, and galaxies become metal-rich rapidly due to the suppression of pristine gas inflow and re-accretion of metal-enriched outflows in overdensities., Comment: 18 pages, 8 figures, 3 tables, accepted for publication in ApJ

Ultraluminous infrared galaxies (ULIRGs) represent a critical stage in the merger-driven evolution of galaxies when AGN activity is common and AGN feedback is expected. We present high sensitivity and large field of view intergral field spectroscopy of the ULIRG Mrk 273 using new data from the Keck Cosmic Web Imager (KWCI). The KCWI data captures the complex nuclear region and the two extended nebulae in the northeast (NE) and southwest (SW) to $\sim 20$ kpc scales. Kinematics in the nuclear region show a fast, extended, bipolar outflow in the direction of the previously reported nuclear superbubbles spanning $\sim 5$ kpc, two to three times greater than the previously reported size. The larger scale extended nebulae on $\sim 20$ kpc show fairly uniform kinematics with FWHM $\sim 300 ~\kmps$ in the SW nebula and FWHM $\sim 120 ~\kmps$ in the NE nebula. We detect for the first time high ionization [NeV]3426, [OIII]4363 and HeII4684 emission lines in the extended NE nebula. Emission line ratios in the nuclear region correlate with the kinematic structures, with the bipolar outflow and the less collimated "outflow regions" showing distinct line ratio trends. Line ratio diagnostics of high ionization emission lines reveal non-trivial contribution from shock plus precursor ionization in the NE nebula and the nuclear region, mixed with AGN photoionization. These data are highly constraining for models of cool ionized gas existing ~20 kpc from a galactic nucleus., Comment: 17 pages, 16 figures, accepted for publication in ApJ

It is widely reported, based on clustering measurements of observed active galactic nuclei (AGN) samples, that AGN reside in similar mass host dark matter halos across the bulk of cosmic time, with log $M/M_\odot$~12.5-13.0 to z~2.5. We show that this is due in part to the AGN fraction in galaxies rising with increasing stellar mass, combined with AGN observational selection effects that exacerbate this trend. Here, we use AGN specific accretion rate distribution functions determined as a function of stellar mass and redshift for star-forming and quiescent galaxies separately, combined with the latest galaxy-halo connection models, to determine the parent and sub-halo mass distribution function of AGN to various observational limits. We find that while the median (sub-)halo mass of AGN, $\approx10^{12}M_\odot$, is fairly constant with luminosity, specific accretion rate, and redshift, the full halo mass distribution function is broad, spanning several orders of magnitude. We show that widely used methods to infer a typical dark matter halo mass based on an observed AGN clustering amplitude can result in biased, systematically high host halo masses. While the AGN satellite fraction rises with increasing parent halo mass, we find that the central galaxy is often not an AGN. Our results elucidate the physical causes for the apparent uniformity of AGN host halos across cosmic time and underscore the importance of accounting for AGN selection biases when interpreting observational AGN clustering results. We further show that AGN clustering is most easily interpreted in terms of the relative bias to galaxy samples, not from absolute bias measurements alone., Comment: 19 pages, 20 figures, accepted for publication in MNRAS. Updated to show median halo masses following referee's helpful comments

We use a sample of 532 star-forming galaxies at redshifts $z\sim 1.4-2.6$ with deep rest-frame optical spectra from the MOSFIRE Deep Evolution Field (MOSDEF) survey to place the first constraints on the nebular attenuation curve at high redshift. Based on the first five low-order Balmer emission lines detected in the composite spectra of these galaxies (${\rm H\alpha}$ through ${\rm H\epsilon}$), we derive a nebular attenuation curve that is similar in shape to that of the Galactic extinction curve, suggesting that the dust covering fraction and absorption/scattering properties along the lines-of-sight to massive stars at high redshift are similar to those of the average Milky Way sightline. The curve derived here implies nebular reddening values that are on average systematically larger than those derived for the stellar continuum. In the context of stellar population synthesis models that include the effects of stellar multiplicity, the difference in reddening of the nebular lines and stellar continuum may imply molecular cloud crossing timescales that are a factor of $\gtrsim 3\times$ longer than those inferred for local molecular clouds, star-formation rates that are constant or increasing with time such that newly-formed and dustier OB associations always dominate the ionizing flux, and/or that the dust responsible for reddening the nebular emission may be associated with non-molecular (i.e., ionized and neutral) phases of the ISM. Our analysis points to a variety of investigations of the nebular attenuation curve that will be enabled with the next generation of ground- and space-based facilities., Comment: 10 pages, 2 tables, 5 figures; submitted 2020 July 04, accepted 2020 September 06 to the Astrophysical Journal

We investigate the evolution of galaxy gas-phase metallicity (O/H) over the range $z=0-3.3$ using samples of $\sim300$ galaxies at $z\sim2.3$ and $\sim150$ galaxies at $z\sim3.3$ from the MOSDEF survey. This analysis crucially utilizes different metallicity calibrations at $z\sim0$ and $z>1$ to account for evolving ISM conditions. We find significant correlations between O/H and stellar mass ($M_*$) at $z\sim2.3$ and $z\sim3.3$. The low-mass power law slope of the mass-metallicity relation is remarkably invariant over $z=0-3.3$, such that $\textrm{O/H}\propto M_*^{0.30}$ at all redshifts in this range. At fixed $M_*$, O/H decreases with increasing redshift as dlog(O/H)/d$z=-0.11\pm0.02$. We find no evidence that the fundamental metallicity relation between $M_*$, O/H, and star-formation rate (SFR) evolves out to $z\sim3.3$, with galaxies at $z\sim2.3-3.3$ having O/H within 0.04~dex of local galaxies matched in $M_*$ and SFR on average. We employ analytic chemical evolution models to place constraints on the mass and metal loading factors of galactic outflows. The efficiency of metal removal increases toward lower $M_*$ at fixed redshift, and toward higher redshift at fixed $M_*$. These models suggest that the slope of the mass-metallicity relation is set by the scaling of the metal loading factor of outflows with $M_*$, not by the change in gas fraction as a function of $M_*$. The evolution toward lower O/H at fixed $M_*$ with increasing redshift is driven by both higher gas fraction (leading to stronger dilution of ISM metals) and higher metal removal efficiency, with models suggesting that both effects contribute approximately equally to the observed evolution. These results suggest that the processes governing the smooth baryonic growth of galaxies via gas flows and star formation hold in the same form over at least the past 12~Gyr., Comment: 36 pages, 19 figures, accepted to ApJ. Main results in Figures 8, 9, 12, and 13

Feedback through energetic outflows has emerged as a key physical process responsible for transforming star-forming galaxies into the quiescent systems observed in the local universe. To explore this process, this paper focuses on a sample of massive and compact merger remnant galaxies hosting high-velocity gaseous outflows ($|v| \gtrsim 10^{3}$ km s$^{-1}$), found at intermediate redshift ($z \sim 0.6$). From their mid-infrared emission and compact morphologies, these galaxies are estimated to have exceptionally large star formation rate (SFR) surface densities ($\Sigma_{SFR} \sim 10^{3}$ $\mathrm{M_{\odot}}$ yr$^{-1}$ kpc$^{-2}$), approaching the Eddington limit for radiation pressure on dust grains. This suggests that star formation feedback may be driving the observed outflows. However, these SFR estimates suffer from significant uncertainties. We therefore sought an independent tracer of star formation to probe the compact starburst activity in these systems. In this paper, we present SFR estimates calculated using 1.5 GHz continuum Jansky Very Large Array observations for 19 of these galaxies. We also present updated infrared (IR) SFRs calculated from WISE survey data. We estimate SFRs from the IR to be larger than those from the radio for 16 out of 19 galaxies by a median factor of 2.5. We find that this deviation is maximized for the most compact galaxies hosting the youngest stellar populations, suggesting that compact starbursts deviate from the IR-radio correlation. We suggest that this deviation stems either from free-free absorption of synchrotron emission, a difference in the timescale over which each indicator traces star formation, or exceptionally hot IR-emitting dust in these ultra-dense galaxies., Comment: 15 pages, 9 figures, 2 tables; accepted for publication in ApJ

We present new measurements of the clustering of stellar mass-complete samples of $\sim40,000$ SDSS galaxies at $z\sim0.03$ as a joint function of stellar mass and specific star formation rate (sSFR). Our results confirm what Coil et al. (2017) find at $z\sim0.7$: galaxy clustering is a stronger function of sSFR at fixed stellar mass than of stellar mass at fixed sSFR. We also find that galaxies above the star-forming main sequence (SFMS) with higher sSFR are less clustered than galaxies below the SFMS with lower sSFR, at a given stellar mass. A similar trend is present for quiescent galaxies. This confirms that main sequence scatter, and scatter within the quiescent sequence, is physically connected to the large-scale cosmic density field. We compare the resulting galaxy bias versus sSFR, and relative bias versus sSFR ratio, for different galaxy samples across ${0

We present results on the properties of neon emission in $z\sim2$ star-forming galaxies drawn from the MOSFIRE Deep Evolution Field (MOSDEF) survey. Doubly-ionized neon ([NeIII]3869) is detected at $\geq3\sigma$ in 61 galaxies, representing $\sim$25% of the MOSDEF sample with H$\alpha$, H$\beta$, and [OIII]$5007$ detections at similar redshifts. We consider the neon emission-line properties of both individual galaxies with [NeIII]3869 detections and composite $z\sim2$ spectra binned by stellar mass. With no requirement of [NeIII]3869 detection, the latter provide a more representative picture of neon emission-line properties in the MOSDEF sample. The [NeIII]3869/[OII]3727 ratio (Ne3O2) is anti-correlated with stellar mass in $z\sim2$ galaxies, as expected based on the mass-metallicity relation. It is also positively correlated with the [OIII]$5007$/[OII]$3727$ ratio (O32), but $z\sim2$ line ratios are offset towards higher Ne3O2 at fixed O32, compared with both local star-forming galaxies and individual H~II regions. Despite the offset towards higher Ne3O2 at fixed O32 at $z\sim2$, biases in inferred Ne3O2-based metallicity are small. Accordingly, Ne3O2 may serve as an important metallicity indicator deep into the reionization epoch. Analyzing additional rest-optical line ratios including [NeIII]$3869$/[OIII]$5007$ (Ne3O3) and [OIII]$5007$/H$\beta$ (O3H$\beta$), we conclude that the nebular emission-line ratios of $z\sim2$ star-forming galaxies suggest a harder ionizing spectrum (lower stellar metallicity, i.e., Fe/H) at fixed gas-phase oxygen abundance, compared to systems at $z\sim0$. These new results based on neon lend support to the physical picture painted by oxygen, nitrogen, hydrogen, and sulfur emission, of an ionized ISM in high-redshift star-forming galaxies irradiated by chemically young, $\alpha$-enhanced massive stars., Comment: 7 pages, 5 figures, accepted to ApJ Letters

We analyze the rest-optical emission-line spectra of $z\sim2.3$ star-forming galaxies in the complete MOSFIRE Deep Evolution Field (MOSDEF) survey. In investigating the origin of the well-known offset between the sequences of high-redshift and local galaxies in the [O III]5008/H$\beta$ vs. [N II]6585/H$\alpha$ ("[N II] BPT") diagram, we define two populations of $z\sim2.3$ MOSDEF galaxies. These include the "high" population that is offset towards higher [O III]5008/H$\beta$ and/or [N II]6585/H$\alpha$ with respect to the local SDSS sequence and the "low" population that overlaps the SDSS sequence. These two groups are also segregated within the [O III]5008/H$\beta$ vs. [S II]6718,6733/H$\alpha$ and the [O III]4960,5008/[O II]3727,3730 (O$_{32}$) vs. ([O III]4960,5008+[O II]3727,3730)/H$\beta$ (R$_{23}$) diagram, which suggests qualitatively that star-forming regions in the more offset galaxies are characterized by harder ionizing spectra at fixed nebular oxygen abundance. We also investigate many galaxy properties of the split sample and find that the "high" sample is on average smaller in size and less massive, but has higher specific star-formation rate and star-formation-rate surface density values and is slightly younger compared to the "low" population. From Cloudy+BPASS photoionization models, we estimate that the "high" population has a lower stellar metallicity (i.e., harder ionizing spectrum) but slightly higher nebular metallicity and higher ionization parameter compared to the "low" population. While the "high" population is more $\alpha$-enhanced (i.e., higher $\alpha$/Fe) than the "low" population, both samples are significantly more $\alpha$-enhanced compared to local star-forming galaxies with similar rest-optical line ratios. These differences must be accounted for in all high-redshift star-forming galaxies -- not only those "offset" from local excitation sequences., Comment: 16 pages, 7 figures

The observed empirical relation between the star-formation rates (SFR) of low-redshift galaxies and their radio continuum luminosity offers a potential means of measuring SFR in high redshift galaxies that is unaffected by dust obscuration. In this study, we make the first test for redshift evolution in the SFR-radio continuum relation at high redshift using dust-corrected H$\alpha$ SFR. Our sample consists of 178 galaxies from the MOSFIRE Deep Evolution Field (MOSDEF) Survey at $1.4 < z < 2.6$ with rest-frame optical spectroscopy and deep 1.5 GHz radio continuum observations from the Karl G. Jansky Very Large Array (VLA) GOODS North field. Using a stacking analysis we compare the observed radio continuum luminosities with those predicted from the dust-corrected H$\alpha$ SFR assuming a range of $z\sim0$ relations. We find no evidence for a systematic evolution with redshift, when stacking the radio continuum as a function of dust-corrected H$\alpha$ SFR and when stacking both optical spectroscopy and radio continuum as a function of stellar mass. We conclude that locally calibrated relations between SFR and radio continuum luminosity remain valid out to $z\sim 2$., Comment: 11 pages, 5 figures. Submitted for publication in MNRAS - comments welcome

We study the properties of 55 morphologically-identified merging galaxy systems at z~2. These systems are flagged as mergers based on features such as tidal tails, double nuclei, and asymmetry. Our sample is drawn from the MOSFIRE Deep Evolution Field (MOSDEF) survey, along with a control sample of isolated galaxies at the same redshift. We consider the relationships between stellar mass, star formation rate (SFR), and gas-phase metallicity for both merging and non-merging systems. In the local universe, merging systems are characterized by an elevated SFR and depressed metallicity compared to isolated systems at a given mass. Our results indicate SFR enhancement and metallicity deficit for merging systems relative to non-merging systems for a fixed stellar mass at z~2, though larger samples are required to establish these preliminary results with higher statistical significance. In future work, it will be important to establish if the enhanced SFR and depressed metallicity in high-redshift mergers deviate from the "fundamental metallicity relation," as is observed in mergers in the local universe, and therefore shed light on gas flows during galaxy interactions., Comment: 9 pages, 5 figures, 5 figures, accepted to MNRAS

We present constraints on the massive star and ionized gas properties for a sample of 62 star-forming galaxies at $z\sim2.3$. Using BPASS stellar population models, we fit the rest-UV spectra of galaxies in our sample to estimate age and stellar metallicity which, in turn, determine the ionizing spectrum. In addition to the median properties of well-defined subsets of our sample, we derive the ages and stellar metallicities for 30 high-SNR individual galaxies -- the largest sample of individual galaxies at high redshift with such measurements. Most galaxies in this high-SNR subsample have stellar metallicities of $0.001

We present spatially-resolved spectroscopy from the Keck Cosmic Web Imager (KCWI) of a star-forming galaxy at z=0.6942, which shows emission from the Mg II 2796, 2803 Angstrom doublet in the circumgalactic medium (CGM) extending ~37 kpc at 3-sigma significance in individual spaxels (1-sigma detection limit 4.8 x 10^{-19} erg s^-1 cm^-2 arcsec^-2). After deconvolution with the seeing, we obtain 5-sigma detections extending for ~31 kpc measured in 7-spaxel (1.1 arcsec^2) apertures. Spaxels covering the galaxy stellar regions show clear P-Cygni-like emission/absorption profiles with the blueshifted absorption extending to relative velocities of v = -800 km/s; however, the P-Cygni profiles give way to pure emission at large radii from the central galaxy. We have performed three-dimensional radiative transfer modeling to infer the geometry and velocity and density profiles of the outflowing gas. Our observations are most consistent with an isotropic outflow rather than biconical wind models with half-opening angles phi <= 80 deg. Furthermore, our modeling suggests that a wind velocity profile that decreases with radius is necessary to reproduce the velocity widths and strengths of Mg II line emission profiles at large circumgalactic radii. The extent of the Mg II emission we measure directly is further corroborated by our modeling, where we rule out outflow models with extent <30 kpc., Comment: Submitted to ApJ; comments welcome!

We derive the UV-optical stellar dust attenuation curve of galaxies at z=1.4-2.6 as a function of gas-phase metallicity. We use a sample of 218 star-forming galaxies, excluding those with very young or heavily obscured star formation, from the MOSFIRE Deep Evolution Field (MOSDEF) survey with H$\alpha$, H$\beta$, and [NII]$\lambda 6585$ spectroscopic measurements. We constrain the shape of the attenuation curve by comparing the average flux densities of galaxies sorted into bins of dust obscuration using Balmer decrements, i.e., H$\alpha$-to-H$\beta$ luminosities. The average attenuation curve for the high-metallicity sample (12+log(O/H)>8.5, corresponding to $M_*\gtrsim10^{10.4}\,M_{\odot}$) has a shallow slope, identical to that of the Calzetti local starburst curve, and a significant UV 2175A extinction bump that is $\sim 0.5\times$ the strength of the Milky Way bump. On the other hand, the average attenuation curve of the low-metallicity sample (12+log(O/H) $\sim 8.2-8.5$) has a steeper slope similar to that of the SMC curve, only consistent with the Calzetti slope at the $3\sigma$ level. The UV bump is not detected in the low-metallicity curve, indicating the relative lack of the small dust grains causing the bump at low metallicities. Furthermore, we find that on average the nebular reddening (E(B-V)) is a factor of 2 times larger than that of the stellar continuum for galaxies with low metallicities, while the nebular and stellar reddening are similar for galaxies with higher metallicities. The latter is likely due to a high surface density of dusty clouds embedding the star forming regions but also reddening the continuum in the high-metallicity galaxies., Comment: 20 pages and 9 figures and 1 appendix, accepted for publication in ApJ

We present a joint analysis of rest-UV and rest-optical spectra obtained using Keck/LRIS and Keck/MOSFIRE for a sample of 62 star-forming galaxies at $z\sim2.3$. We divide our sample into 2 bins based on their location in the [OIII]/Hb vs. [NII]/Ha BPT diagram, and perform the first differential study of the rest-UV properties of massive ionizing stars as a function of rest-optical emission-line ratios. Fitting BPASS stellar population synthesis models, including nebular continuum emission, to our rest-UV spectra, we find that high-redshift galaxies offset towards higher [OIII]/Hb and [NII]/Ha have younger ages (log(Age/yr)=$7.20^{+0.57}_{-0.20}$) and lower stellar metallicities ($Z_*=0.0010^{+0.0011}_{-0.0003}$) resulting in a harder ionizing spectrum, compared to the galaxies in our sample that lie on the local BPT star-forming sequence (log(Age/yr)=$8.57^{+0.88}_{-0.84}$, $Z_*=0.0019\pm0.0006$). Additionally, we find that the offset galaxies have an ionization parameter of log(U)=$-3.04^{+0.06}_{-0.11}$ and nebular metallicity of 12+log(O/H)=$8.40^{+0.06}_{-0.07}$, and the non-offset galaxies have log(U)=$-3.11\pm0.08$ and 12+log(O/H)=$8.30^{+0.05}_{-0.06}$. The stellar and nebular metallicities derived for our sample imply that the galaxies offset from the local BPT relation are more alpha-enhanced ($7.28^{+2.52}_{-2.82}$ O/Fe$_\odot$) compared to those consistent with the local sequence ($3.04^{+0.95}_{-0.54}$ O/Fe$_\odot$). However, even galaxies that are entirely consistent with the local nebular excitation sequence appear to be alpha-enhanced -- in contrast with typical local systems. Such differences must be considered when estimating gas-phase O/H at high redshift based on strong emission-line ratios. Specifically, a similarity in the location of high-redshift and local galaxies in the BPT diagram may not be indicative of a similarity in their physical properties., Comment: 16 pages, 12 figures. Accepted to MNRAS

We present measurements of [SIII]$\lambda\lambda$9069,9531 for a sample of $z\sim1.5$ star-forming galaxies, the first sample with measurements of these lines at z>0.1. We employ the line ratio S$_{32}$$\equiv$[SIII]$\lambda\lambda$9069,9531/[SII]$\lambda\lambda$6716,6731 as a novel probe of evolving ISM conditions. Since this ratio includes the low-ionization line [SII], it is crucial that the effects of diffuse ionized gas (DIG) on emission-line ratios be accounted for in $z\sim0$ integrated galaxy spectra, or else that comparisons be made to samples of local HII regions in which DIG emission is not present. We find that S$_{32}$ decreases with increasing stellar mass at both $z\sim1.5$ and $z\sim0$, but that the dependence is weak suggesting S$_{32}$ has a very shallow anticorrelation with metallicity, in contrast with O$_{32}$ that displays a strong metallicity dependence. As a result, S$_{32}$ only mildly evolves with redshift at fixed stellar mass. The $z\sim1.5$ sample is systematicallty offset towards lower S$_{32}$ and higher [SII]/H$\alpha$ at fixed [OIII]/H$\beta$ relative to $z=0$ HII regions. By comparing to photoionization model grids, we find that such trends can be explained by a scenario in which the ionizing spectrum is harder at fixed O/H with increasing redshift, but are inconsistent with an increase in ionization parameter at fixed O/H. This analysis demonstrates the advantages of expanding beyond the strongest rest-optical lines for evolutionary studies, and the particular utility of [SIII] for characterizing evolving ISM conditions and stellar compositions. These measurements provide a basis for estimating [SIII] line strengths for high-redshift galaxies, a line that the James Webb Space Telescope will measure out to z~5.5., Comment: 6 pages, 4 figures, accepted for publication in ApJ Letters

Population synthesis models predict that high-mass X-ray binary (HMXB) populations produced in low metallicity environments should be more X-ray luminous, a trend supported by studies of nearby galaxies. This trend may be responsible for the observed increase of the X-ray luminosity ($L_{\mathrm{X}}$) per star formation rate (SFR) with redshift due to the decrease of metallicity ($Z$) at fixed stellar mass as a function of redshift. To test this hypothesis, we use a sample of 79 $z\sim2$ star-forming galaxies with oxygen abundance measurements from the MOSDEF survey, which obtained rest-frame optical spectra for $\sim1500$ galaxies in the CANDELS fields at $1.37

We present results on the emission-line properties of 1.3<=z<=2.7 galaxies drawn from the complete MOSFIRE Deep Evolution Field (MOSDEF) survey. Specifically, we use observations of the emission-line diagnostic diagram of [OIII]5007/Hb vs. [SII]6717,6731/Ha, i.e., the "[SII] BPT diagram," to gain insight into the physical properties of high-redshift star-forming regions. High-redshift MOSDEF galaxies are offset towards lower [SII]6717,6731/Ha at fixed [OIII]5007/Hb, relative to local galaxies from the Sloan Digital Sky Survey (SDSS). Furthermore, at fixed [OIII]5007/Hb, local SDSS galaxies follow a trend of decreasing [SII]6717,6731/Ha as the surface density of star formation (Sigma_SFR) increases. We explain this trend in terms of the decreasing fractional contribution from diffuse ionized gas (f_DIG) as Sigma_SFR increases in galaxies, which causes galaxy-integrated line ratios to shift towards the locus of pure HII-region emission. The z~0 relationship between f_DIG and Sigma_SFR implies that high-redshift galaxies have lower f_DIG values than typical local systems, given their significantly higher typical Sigma_SFR. When an appropriate low-redshift benchmark with zero or minimal f_DIG is used, high-redshift MOSDEF galaxies appear offset towards higher [SII]6717,6731/Ha and/or [OIII]5007/Hb. The joint shifts of high-redshift galaxies in the [SII] and [NII] BPT diagrams are best explained in terms of the harder spectra ionizing their star-forming regions at fixed nebular oxygen abundance (expected for chemically-young galaxies), as opposed to large variations in N/O ratios or higher ionization parameters. The evolving mixture of HII regions and DIG is an essential ingredient to our description of the ISM over cosmic time., Comment: 6 pages, 4 figures, accepted to ApJ Letters

We present detections of [OIII]$\lambda$4363 and direct-method metallicities for star-forming galaxies at $z=1.7-3.6$. We combine new measurements from the MOSFIRE Deep Evolution Field (MOSDEF) survey with literature sources to construct a sample of 18 galaxies with direct-method metallicities at $z>1$, spanning $7.5<1$2+log(O/H$)<8.2$ and log(M$_*$/M$_{\odot})=7-10$. We find that strong-line calibrations based on local analogs of high-redshift galaxies reliably reproduce the metallicity of the $z>1$ sample on average. We construct the first mass-metallicity relation at $z>1$ based purely on direct-method O/H, finding a slope that is consistent with strong-line results. Direct-method O/H evolves by $\lesssim0.1$ dex at fixed M$_*$ and SFR from $z\sim0-2.2$. We employ photoionization models to constrain the ionization parameter and ionizing spectrum in the high-redshift sample. Stellar models with super-solar O/Fe and binary evolution of massive stars are required to reproduce the observed strong-line ratios. We find that the $z>1$ sample falls on the $z\sim0$ relation between ionization parameter and O/H, suggesting no evolution of this relation from $z\sim0$ to $z\sim2$. These results suggest that the offset of the strong-line ratios of this sample from local excitation sequences is driven primarily by a harder ionizing spectrum at fixed nebular metallicity compared to what is typical at $z\sim0$, naturally explained by super-solar O/Fe at high redshift caused by rapid formation timescales. Given the extreme nature of our $z>1$ sample, the implications for representative $z\sim2$ galaxy samples at $\sim10^{10}$ M$_{\odot}$ are unclear, but similarities to $z>6$ galaxies suggest that these conclusions can be extended to galaxies in the epoch of reionization., Comment: 27 pages + 4 page appendix, 18 figures, accepted for publication in MNRAS