Search

We present results from a Keck/DEIMOS survey to study satellite quenching in group environments at $z \sim 0.8$ within the Extended Groth Strip (EGS). We target $11$ groups in the EGS with extended X-ray emission. We obtain high-quality spectroscopic redshifts for group member candidates, extending to depths over an order of magnitude fainter than existing DEEP2/DEEP3 spectroscopy. This depth enables the first spectroscopic measurement of the satellite quiescent fraction down to stellar masses of $\sim 10^{9.5}~{\rm M}_{\odot}$ at this redshift. By combining an infall-based environmental quenching model, constrained by the observed quiescent fractions, with infall histories of simulated groups from the IllustrisTNG100-1-Dark simulation, we estimate environmental quenching timescales ($\tau_{\mathrm{quench}}$) for the observed group population. At high stellar masses (${M}_{\star}=10^{10.5}~{\rm M}_{\odot}$) we find that $\tau_{\mathrm{quench}} = 2.4\substack{+0.2 \\ -0.2}$ Gyr, which is consistent with previous estimates at this epoch. At lower stellar masses (${M}_{\star}=10^{9.5}~{\rm M}_{\odot}$), we find that $\tau_{\mathrm{quench}}=3.1\substack{+0.5 \\ -0.4}$ Gyr, which is shorter than prior estimates from photometry-based investigations. These timescales are consistent with satellite quenching via starvation, provided the hot gas envelope of infalling satellites is not stripped away. We find that the evolution in the quenching timescale between $0 \lt z \lt 1$ aligns with the evolution in the dynamical time of the host halo and the total cold gas depletion time. This suggests that the doubling of the quenching timescale in groups since $z\sim1$ could be related to the dynamical evolution of groups or a decrease in quenching efficiency via starvation with decreasing redshift., Comment: 10 figures; 2 tables; Accepted for publication in the Astrophysical Journal

Supermassive black holes (SMBHs) are observed in diverse galaxy populations across time yet a clear understanding of how they coevolve with their hosts has not been reached. Physically-motivated models of SMBH accretion and feedback vary widely between galaxy formation simulations due to the difficulty of modeling the range of scales important for galactic and SMBH processes. Here we use observational data to build an empirical model for SMBH growth. We apply observed specific accretion rate probability distributions as a function of galaxy star formation rate between $z = 0-2$ to the UniverseMachine galaxy formation model to determine SMBH accretion rates based on galaxy properties. We use observed $z = 0$ SMBH-stellar mass relations for the quiescent and star-forming populations to provide the local boundary conditions for SMBH growth histories. We then track the coevolutionary histories of galaxy stellar mass and their SMBHs backwards in time to $z = 2$. We find that the most massive SMBHs at $z = 0$ have grown very little of their total mass between $z = 0-2$, indicating early SMBH mass assembly for these systems. Conversely, lower mass SMBHs at $z = 0$ assembled their mass gradually across $z = 0-2$. This results in substantial evolution of the SMBH-stellar mass relation, shifting to higher normalization and shallower slope with increasing redshift. We find that the substantial scatter observed in the $z = 0$ SMBH-stellar mass relation results in the diversity of growth pathways found in our model, with some galaxies assembling their stellar mass before their SMBHs and others doing the opposite., Comment: 24 pages, 13 figures, submitted to ApJ

Active galactic nucleus (AGN) growth in disk-dominated, merger-free galaxies is poorly understood, largely due to the difficulty in disentangling the AGN emission from that of the host galaxy. By carefully separating this emission, we examine the differences between AGNs in galaxies hosting a (possibly) merger-grown, classical bulge, and AGNs in secularly grown, truly bulgeless disk galaxies. We use GALFIT to obtain robust, accurate morphologies of 100 disk-dominated galaxies imaged with the Hubble Space Telescope. Adopting an inclusive definition of classical bulges, we detect a classical bulge component in $53.3 \pm 0.5$ per cent of the galaxies. These bulges were not visible in Sloan Digital Sky Survey photometry, however these galaxies are still unambiguously disk-dominated, with an average bulge-to-total luminosity ratio of $0.1 \pm 0.1$. We find some correlation between bulge mass and black hole mass for disk-dominated galaxies, though this correlation is significantly weaker in comparison to the relation for bulge-dominated or elliptical galaxies. Furthermore, a significant fraction ($\gtrsim 90$ per cent) of our black holes are overly massive when compared to the relationship for elliptical galaxies. We find a weak correlation between total stellar mass and black hole mass for the disk-dominated galaxies, hinting that the stochasticity of black hole-galaxy co-evolution may be higher disk-dominated than bulge-dominated systems., Comment: 14 pages, 6 figures, 2 tables. Submitted to MNRAS

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 provide 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), indicating a bias towards massive galaxies in the observed samples. 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: 21 pages, 13 figures, accepted 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

The connection between active galactic nuclei (AGN) and their host dark matter halos provides powerful insights into how supermassive black holes (SMBHs) grow and coevolve with their host galaxies. Here we investigate the impact of observational AGN selection on the AGN halo occupation distribution (HOD) by forward-modeling AGN activity into cosmological N-body simulations. By assuming straightforward relationships between the SMBH mass, galaxy mass, and (sub)halo mass, as well as a uniform broken power law distribution of Eddington ratios, we find that luminosity-limited AGN samples result in biased HOD shapes. While AGN defined by an Eddington ratio threshold produce AGN fractions that are flat across halo mass (unbiased by definition), luminosity-limited AGN fractions peak around galaxy-group-sized halo masses and then decrease with increasing halo mass. With higher luminosities, the rise of the AGN fraction starts at higher halo masses, the peak is shifted towards higher halo masses, and the decline at higher halo masses is more rapid. These results are consistent with recent HOD constraints from AGN clustering measurements, which find (1) characteristic halo mass scales of $\log M_{Vir}\sim$ 12 - 13 [$h^{-1}M_{\odot}$] and (2) a shallower rise of the number of satellite AGN with increasing halo mass than for the overall galaxy population. Thus the observational biases due to AGN selection can naturally explain the constant, characteristic halo mass scale inferred from large-scale AGN clustering amplitudes over a range of redshifts, as well as the measured inconsistencies between AGN and galaxy HODs. We conclude that AGN selection biases can have significant impacts on the inferred AGN HOD, and can therefore lead to possible misinterpretations of how AGN populate dark matter halos and the AGN-host galaxy connection., Comment: Accepted for publication in A&A

A new class of radio source, the so-called Odd Radio Circles (ORCs), have been discovered by recent sensitive, large-area radio continuum surveys. The distances of these sources have so far relied on photometric redshifts of optical galaxies found at the centers of or near the ORCs. Here we present Gemini rest-frame optical spectroscopy of six galaxies at the centers of, or potentially associated with, the first five ORC discoveries. We supplement this with Legacy Survey imaging and Prospector fits to their griz+W1/W2 photometry. Of the three ORCs with central galaxies, all lie at distances (z = 0.27-0.55) that confirm the large intrinsic diameters of the radio circles (300-500 kpc). The central galaxies are massive ($M_*\sim10^{11}M_\odot$), red, unobscured ellipticals with old ($\gtrsim$1~Gyr) stellar populations. They have LINER spectral types that are shock- or AGN-powered. All three host low-luminosity, radio-quiet AGN. The similarity of their central galaxies are consistent with a common origin, perhaps as a blastwave from an ancient starburst. The other two ORCs are adjacent and have no prominent central galaxies. However, the z=0.25 disk galaxy that lies between them hosts a Type 2, moderate-luminosity AGN. They may instead be the lobes of a radio jet from this AGN., Comment: ApJ, in press. Changes to match accepted version; only significant change is addition of Fig. 3

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

Changing-look active galactic nuclei (CL-AGN) have been observed to change optical spectral type. Mrk 1018 is unique: first classified as a type 1.9 Seyfert galaxy, it transitioned to a type 1 before returning to its initial classification after approximately 30 years. We present a high-cadence monitoring programme that caught a major outburst in 2020. Due to sunblock, only the decline could be observed. We studied X-ray, UV, optical, and IR before and after the outburst to investigate the responses of the AGN structures. We derived a u'-band light curve of the AGN contribution alone. The flux increased by a factor of the order of 13. We confirmed this in other optical bands and determined the shape and speed of the decline in each waveband. The shapes of H beta and H alpha were analysed before and after the event. Two XMM-Newton observations from before and after the outburst were also exploited. The outburst is asymmetric, with a swifter rise than decline. The decline is best fit by a linear function, ruling out a tidal disruption event. The optical spectrum shows no change approximately 8 months before and 17 months after. The UV flux increased slightly after the outburst but the X-ray primary flux is unchanged. However, the 6.4 keV Iron line has doubled in strength. IR data taken 13 days after the observed optical peak show an increased emission level. Calculating the distance of the broad-line region and inner edge of the torus from the supermassive black hole can explain the multi-wavelength response to the outburst, in particular: i) the unchanged H beta and H alpha lines, ii) the unchanged primary X-ray spectral components, iii) the rapid and extended infrared response, as well as iv) the enhanced emission of the reflected 6.4 keV line. The outburst was due to a dramatic and short-lasting change in the intrinsic accretion rate. We discuss different models as potential causes., Comment: Accepted for publication in A&A

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

Noun compound interpretation is the task of expressing a noun compound (e.g. chocolate bunny) in a free-text paraphrase that makes the relationship between the constituent nouns explicit (e.g. bunny-shaped chocolate). We propose modifications to the data and evaluation setup of the standard task (Hendrickx et al., 2013), and show that GPT-3 solves it almost perfectly. We then investigate the task of noun compound conceptualization, i.e. paraphrasing a novel or rare noun compound. E.g., chocolate crocodile is a crocodile-shaped chocolate. This task requires creativity, commonsense, and the ability to generalize knowledge about similar concepts. While GPT-3's performance is not perfect, it is better than that of humans -- likely thanks to its access to vast amounts of knowledge, and because conceptual processing is effortful for people (Connell and Lynott, 2012). Finally, we estimate the extent to which GPT-3 is reasoning about the world vs. parroting its training data. We find that the outputs from GPT-3 often have significant overlap with a large web corpus, but that the parroting strategy is less beneficial for novel noun compounds., Comment: Findings of ACL 2023

Measurement of fusion excitation functions for stable nuclei has largely been restricted to nuclei with significant natural abundance. Typically, to investigate neighboring nuclei with low natural abundance has required obtaining isotopically enriched material. This restriction often limits the ability to perform such measurements. We report the measurement of a high quality fusion excitation function for a $^{17}$O beam produced from unenriched material with 0.038\% natural abundance. The measurement is enabled by using an active thick-target approach and the accuracy of the result is validated using its mirror nucleus $^{17}$F and resonances. The result provides important information about the average fusion cross-section for the oxygen isotopic chain as a function of neutron excess., Comment: 4 pages, 4 figures

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

The role of large-scale bars in the fuelling of active galactic nuclei (AGN) is still debated, even as evidence mounts that black hole growth in the absence of galaxy mergers cumulatively dominated and may substantially influence disc (i.e., merger-free) galaxy evolution. We investigate whether large-scale galactic bars are a good candidate for merger-free AGN fuelling. Specifically, we combine slit spectroscopy and Hubble Space Telescope imagery to characterise star formation rates (SFRs) and stellar masses of the unambiguously disc-dominated host galaxies of a sample of luminous, Type-1 AGN with 0.02 < z 0.024. After carefully correcting for AGN signal, we find no clear difference in SFR between AGN hosts and a stellar mass-matched sample of galaxies lacking an AGN (0.013 < z < 0.19), although this could be due to a small sample size (n_AGN = 34). We correct for SFR and stellar mass to minimise selection biases, and compare the bar fraction in the two samples. We find that AGN are marginally (1.7$\sigma$) more likely to host a bar than inactive galaxies, with AGN hosts having a bar fraction, fbar = 0.59^{+0.08}_{-0.09} and inactive galaxies having a bar fraction fbar = 0.44^{+0.08}_{-0.09}. However, we find no further differences between SFR- and mass-matched AGN and inactive samples. While bars could potentially trigger AGN activity, they appear to have no further, unique effect on a galaxy's stellar mass or SFR., Comment: 15 pages (9 figures). Accepted for publication in MNRAS

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

Recent observational and theoretical studies have suggested that supermassive black holes (SMBHs) grow mostly through non-merger (`secular') processes. Since galaxy mergers lead to dynamical bulge growth, the only way to observationally isolate non-merger growth is to study galaxies with low bulge-to-total mass ratio (e.g. B/T < 10%). However, bulge growth can also occur due to secular processes, such as disk instabilities, making disk-dominated selections a somewhat incomplete way to select merger-free systems. Here we use the Horizon-AGN simulation to select simulated galaxies which have not undergone a merger since z = 2, regardless of bulge mass, and investigate their location on typical black hole-galaxy scaling relations in comparison to galaxies with merger dominated histories. While the existence of these correlations has long been interpreted as co-evolution of galaxies and their SMBHs driven by galaxy mergers, we show here that they persist even in the absence of mergers. We find that the correlations between SMBH mass and both total mass and stellar velocity dispersion are independent of B/T ratio for both merger-free and merger-dominated galaxies. In addition, the bulge mass and SMBH mass correlation is still apparent for merger-free galaxies, the intercept for which is dependent on B/T. Galaxy mergers reduce the scatter around the scaling relations, with merger-free systems showing broader scatter. We show that for merger-free galaxies, the co-evolution is dominated by radio-mode feedback, and suggest that the long periods of time between galaxy mergers make an important contribution to the co-evolution between galaxies and SMBHs in all galaxies., Comment: RJS and RSB are joint first authors. 12 pages, 7 figures. Accepted 2023 February 13. Received 2023 February 2; in original form 2022 November 24

Here we use the Horizon-AGN simulation to test whether the spins of SMBHs in merger-free galaxies are higher. We select samples using an observationally motivated bulge-to-total mass ratio of < 0.1, along with two simulation motivated thresholds selecting galaxies which have not undergone a galaxy merger since z = 2, and those SMBHs with < 10% of their mass due to SMBH mergers. We find higher spins (> 5{\sigma} ) in all three samples compared to the rest of the population. In addition, we find that SMBHs with their growth dominated by BH mergers following galaxy mergers, are less likely to be aligned with their galaxy spin than those that have grown through accretion in the absence of galaxy mergers (3.4{\sigma} ). We discuss the implications this has for the impact of active galactic nuclei (AGN) feedback, finding that merger-free SMBHs spend on average 91% of their lifetimes since z = 2 in a radio mode of feedback (88% for merger-dominated galaxies). Given that previous observational and theoretical works have concluded that merger-free processes dominate SMBH-galaxy co-evolution, our results suggest that this co-evolution could be regulated by radio mode AGN feedback., Comment: RSB and RJS are joint first authors. 11 pages, 7 figures. Accepted 2023 June 9. Received 2023 April 25; in original form 2022 November 24

We report here on the recent upgrade of the TwinSol radioactive nuclear beam (RNB) facility at the University of Notre Dame. The new TriSol system includes a magnetic dipole to provide a second beamline and a third solenoid which acts to reduce the size of the radioactive beam on target., Comment: submitted to Nuclear Instruments and Methods A

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

Machine learning (ML), being now widely accessible to the research community at large, has fostered a proliferation of new and striking applications of these emergent mathematical techniques across a wide range of disciplines. In this paper, we will focus on a particular case study: the field of paleoanthropology, which seeks to understand the evolution of the human species based on biological and cultural evidence. As we will show, the easy availability of ML algorithms and lack of expertise on their proper use among the anthropological research community has led to foundational misapplications that have appeared throughout the literature. The resulting unreliable results not only undermine efforts to legitimately incorporate ML into anthropological research, but produce potentially faulty understandings about our human evolutionary and behavioral past. The aim of this paper is to provide a brief introduction to some of the ways in which ML has been applied within paleoanthropology; we also include a survey of some basic ML algorithms for those who are not fully conversant with the field, which remains under active development. We discuss a series of missteps, errors, and violations of correct protocols of ML methods that appear disconcertingly often within the accumulating body of anthropological literature. These mistakes include use of outdated algorithms and practices; inappropriate train/test splits, sample composition, and textual explanations; as well as an absence of transparency due to the lack of data/code sharing, and the subsequent limitations imposed on independent replication. We assert that expanding samples, sharing data and code, re-evaluating approaches to peer review, and, most importantly, developing interdisciplinary teams that include experts in ML are all necessary for progress in future research incorporating ML within anthropology.

The $\rm{H}\alpha$-to-UV luminosity ratio ($L(\rm H\alpha)/L(\rm UV)$) is often used to probe SFHs of star-forming galaxies and it is important to validate it against other proxies for burstiness. To address this issue, we present a statistical analysis of the resolved distribution of $\Sigma_{\rm{SFR}}$ as well as stellar age and their correlations with the globally measured $L(\rm H\alpha)/L(\rm UV)$ for a sample of 310 star-forming galaxies in two redshift bins of $1.37 < z < 1.70$ and $ 2.09 < z < 2.61$ observed by the MOSDEF survey. We use the multi-waveband CANDELS/3D-HST imaging of MOSDEF galaxies to construct $\Sigma_{\rm{SFR}}$ and stellar age maps. We analyze the composite rest-frame far-UV spectra of a subsample of MOSDEF targets obtained by the Keck/LRIS, which includes 124 star-forming galaxies (MOSDEF-LRIS) at redshifts $1.4 < z < 2.6$, to examine the average stellar population properties, and the strength of age-sensitive FUV spectral features in bins of $L(\rm H\alpha)/L(\rm UV)$. Our results show no significant evidence that individual galaxies with higher $L(\rm H\alpha)/L(\rm UV)$ are undergoing a burst of star formation based on the resolved distribution of $\Sigma_{\rm{SFR}}$ of individual star-forming galaxies. We segregate the sample into subsets with low and high $L(\rm H\alpha)/L(\rm UV)$. The high-$L(\rm H\alpha)/L(\rm UV)$ subset exhibits, on average, an age of $\log[\rm{Age/yr}]$ = 8.0, compared to $\log[\rm{Age/yr}]$ = 8.4 for the low-$L(\rm H\alpha)/L(\rm UV)$ galaxies, though the difference in age is significant at only the $2\sigma$ level. Furthermore, we find no variation in the strengths of Siiv$\lambda\lambda1393, 1402$ and Civ$\lambda\lambda1548, 1550$ P-Cygni features from massive stars between the two subsamples., Comment: 16 pages, 10 figures, published by the Monthly Notices of the Royal Astronomical Society

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

We investigate the conditions that facilitate galactic-scale outflows using a sample of 155 typical star-forming galaxies at $z$~2 drawn from the MOSFIRE Deep Evolution Field (MOSDEF) survey. The sample includes deep rest-frame UV spectroscopy from the Keck Low-Resolution Imaging Spectrometer (LRIS), which provides spectral coverage of several low-ionisation interstellar (LIS) metal absorption lines and Ly$\alpha$ emission. Outflow velocities are calculated from the centroids of the LIS absorption and/or Ly$\alpha$ emission, as well as the highest-velocity component of the outflow from the blue wings of the LIS absorption lines. Outflow velocities are found to be marginally correlated or independent of galaxy properties, such as star-formation rate (SFR) and star-formation rate surface density ($\Sigma_{\rm SFR}$). Outflow velocity scales with SFR as a power-law with index 0.24, which suggests that the outflows may be primarily driven by mechanical energy generated by supernovae explosions, as opposed to radiation pressure acting on dusty material. On the other hand, outflow velocity and $\Sigma_{\rm SFR}$ are not significantly correlated, which may be due to the limited dynamic range of $\Sigma_{\rm SFR}$ probed by our sample. The relationship between outflow velocity and $\Sigma_{\rm SFR}$ normalised by stellar mass ($\Sigma_{\rm sSFR}$), as a proxy for gravitational potential, suggests that strong outflows (e.g., > 200 km s$^{-1}$) appear ubiquitous above a threshold of log($\Sigma_{\rm sSFR}/\rm{yr}^{-1}\ \rm{kpc}^{-2}$) ~ -11.3, and that above this threshold, outflow velocity uncouples from $\Sigma_{\rm sSFR}$. These results highlight the need for higher resolution spectroscopic data and spatially resolved imaging to test the driving mechanisms of outflows predicted by theory., Comment: 15 pages, 5 figures, Submitted to Monthly Notices of the Royal Astronomical Society

We present a $z = 0.94$ quasar, SDSS J004846.45-004611.9, discovered in the SDSS-III BOSS survey. A visual analysis of this spectrum reveals highly broadened and blueshifted narrow emission lines, in particular [Ne~V]$\lambda3426$ and [O~III]$\lambda5007$, with outflow velocities of 4000 km s$^{-1}$, along with unusually large [Ne V]$\lambda3426$/[Ne III]$\lambda3869$ ratios. The gas shows higher ionization at higher outflow velocities, indicating a connection between the powerful outflow and the unusual strength of the high ionization lines. The SED and the $i - \text{W3}$ color of the source reveal that it is likely a "core" Extremely Red Quasar (core ERQ); a candidate population of young AGN that are violently "blowing out" gas and dust from their centers. The dominance of host galaxy light in its spectrum and its fortuitous position in the SDSS S82 region allows us to measure its star formation history and investigate for variability for the first time in an ERQ. Our analysis indicates that SDSS J004846.45-004611.9 underwent a short-lived starburst phase 400 Myr ago and was subsequently quenched, possibly indicating a time-lag between star formation quenching and the onset of AGN activity. We also find that the strong extinction can be uniquely attributed to the AGN and does not persist in the host galaxy, contradicting a scenario where the source has recently transitioned from being a dusty sub-mm galaxy. In our relatively shallow photometric data, the source does not appear to be variable at $0.24-2.4~\mu$m in the restframe, most likely due to the dominant contribution of host galaxy starlight at these wavelengths., Comment: 23 pages, 8 figures, accepted for publication in 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

We present results from ALMA 1.2mm continuum observations of a sample of 27 star-forming galaxies at z=2.1-2.5 from the MOSFIRE Deep Evolution Field (MOSDEF) survey. These galaxies have gas-phase metallicity and star-formation rate measurements from Hb, [OIII], Ha, and [NII]. Using stacks of Spitzer, Herschel, and ALMA photometry (rest-frame ~ 8-400$\mu$m), we examine the IR SED of high-redshift subsolar metallicity (~0.5 $Z_{\odot}$) LIRGs. We find that the data agree well with an average SED template of higher luminosity local low-metallicity dwarf galaxies (reduced $\chi^2$ of 1.8). When compared with the commonly used templates for solar-metallicity local galaxies or high-redshift LIRGs and ULIRGs, even in the most favorable case (with reduced $\chi^2$ of 2.8), the templates are rejected at >98% confidence level. The broader and hotter IR SED of both the local dwarfs and high-redshift subsolar metallicity galaxies may result from different grain properties, a clumpy dust geometry, or a harder/more intense ionizing radiation field that heats the dust to higher temperatures. The obscured SFR indicated by the FIR emission of the subsolar metallicity galaxies is only ~ 60% of the total SFR, which is considerably lower than that of the local LIRGs with ~ 96-97% obscured fractions. Due to the evolving IR SED shape, the local LIRG templates fit to mid-IR data can overestimate the Rayleigh-Jeans tail measurements at z~2 by a factor of 2-20, and these templates underestimate IR luminosities if fit to the observed ALMA fluxes by >0.4dex. At a given stellar mass or metallicity, dust masses at z~2.3 are an order of magnitude higher than those at z~0. Given the predicted molecular gas mass fractions, the observed z~2.3 dust-to-stellar mass ratios suggest lower dust-to-molecular gas masses than in local galaxies at the same metallicity., Comment: Accepted for publication in ApJ. 29 pages, 2 appendices, 10 figures, 3 tables

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