Sean T. Linden, Thomas Lai, Aaron S. Evans, Lee Armus, Kirsten L. Larson, Jeffrey A. Rich, Vivian U, George C. Privon, Hanae Inami, Yiqing Song, Marina Bianchin, Thomas Bohn, Victorine A. Buiten, Maria Sanchez-Garcia, Justin Kader, Laura Lenkic, Anne M. Medling, Torsten Boeker, Tanio Diaz-Santos, Vassilis Charmandaris, et al (11) We present the results of a James Webb Space Telescope (JWST) NIRCam and NIRSpec investigation into the young massive star cluster (YMC) population of NGC 3256, the most cluster-rich luminous infrared galaxy (LIRG) in the Great Observatories All Sky LIRG Survey. We detect 3061 compact YMC candidates with a $S/N \geq 3$ at F150W, F200W, and F335M. Based on yggdrasil stellar population models, we identify 116/3061 sources with F150W - F200W $> 0.47$ and F200W - F355M $> -1.37$ colors suggesting they are young (t $\leq 5$ Myr), dusty ($A_{V} = 5 - 15$), and massive ($M_{\odot} > 10^{5}$). This increases the sample of dust-enshrouded YMCs detected in this system by an order of magnitude relative to previous HST studies. With NIRSpec IFU pointings centered on the northern and southern nucleus, we extract the Pa$\alpha$ and 3.3$\mu$m PAH equivalent widths for 8 bright and isolated YMCs. Variations in both the F200W - F335M color and 3.3$\mu$m PAH emission with the Pa$\alpha$ line strength suggest a rapid dust clearing ($< 3 - 4$ Myr) for the emerging YMCs in the nuclei of NGC 3256. Finally, with both the age and dust emission accurately measured we use yggdrasil to derive the color excess (E(B - V)) for all 8 YMCs. We demonstrate that YMCs with strong 3.3$\mu$m PAH emission (F200W - F335M $> 0$) correspond to sources with E(B - V) $> 3$, which are typically missed in UV-optical studies. This underscores the importance of deep near-infrared imaging for finding and characterizing these very young and dust-embedded sources.
Katherine Alatalo, Andreea O. Petric, Lauranne Lanz, Kate Rowlands, Vivian U, Kirsten L. Larson, Lee Armus, Loreto Barcos-Muñoz, Aaron S. Evans, Jin Koda, Yuanze Luo, Anne M. Medling, Kristina E. Nyland, Justin A. Otter, Pallavi Patil, Fernando Peñaloza, Diane Salim, David B. Sanders, Elizaveta Sazonova, Maya Skarbinski, et al (3) We present the CO(1-0) maps of 28 infrared-bright galaxies from the Great Observatories All-Sky Luminous Infrared Galaxy Survey (GOALS) taken with the Combined Array for Research in Millimeter Astronomy (CARMA). We detect 100GHz continuum in 16 of 28 galaxies, which trace both active galactic nuclei (AGNs) and compact star-forming cores. The GOALS galaxies show a variety of molecular gas morphologies, though in the majority of cases, the average velocity fields show a gradient consistent with rotation. We fit the full continuum SEDs of each of the source using either MAGPHYS or SED3FIT (if there are signs of an AGN) to derive the total stellar mass, dust mass, and star formation rates of each object. We adopt a value determined from luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs) of $\alpha_{\rm CO}=1.5^{+1.3}_{-0.8}~M_\odot$ (K km s$^{-1}$ pc$^2)^{-1}$, which leads to more physical values for $f_{\rm mol}$ and the gas-to-dust ratio. Mergers tend to have the highest gas-to-dust ratios. We assume the cospatiality of the molecular gas and star formation, and plot the sample on the Schmidt-Kennicutt relation, we find that they preferentially lie above the line set by normal star-forming galaxies. This hyper-efficiency is likely due to the increased turbulence in these systems, which decreases the freefall time compared to star-forming galaxies, leading to "enhanced" star formation efficiency. Line wings are present in a non-negligible subsample (11/28) of the CARMA GOALS sources and are likely due to outflows driven by AGNs or star formation, gas inflows, or additional decoupled gas components.
Measurements of oxygen abundance throughout galaxies provide insight to the formation histories and ongoing processes. Here we present a study of the gas phase oxygen abundance in the HII regions and diffuse gas of the nearby starburst dwarf galaxy, IC 10. Using the Keck Cosmic Web Imager (KCWI) at W.M. Keck Observatory, we map the central region of IC 10 from 3500-5500A. The auroral [OIII]4363A line is detected with high signal-to-noise in 12 of 46 HII regions observed, allowing for direct measurement of the oxygen abundance, yielding a median and standard deviation of $\rm12+log(O/H)=8.37\pm0.25$. We investigate trends between these directly measured oxygen abundances and other HII region properties, finding weak negative correlations with the radius, velocity dispersion, and luminosity. We also find weak negative correlations between oxygen abundance and the derived quantities of turbulent pressure and ionized gas mass, and a moderate correlation with the derived dynamical mass. Strong line, $\rm R_{23}$ abundance estimates are used in the remainder of the HII regions and on a resolved spaxel-by-spaxel basis. There is a large offset between the abundances measured with $\rm R_{23}$ and the auroral line method. We find that the $\rm R_{23}$ method is unable to capture the large range of abundances observed via the auroral line measurements. The extent of this variation in measured abundances further indicates a poorly mixed interstellar medium (ISM) in IC 10, which is not typical of dwarf galaxies and may be partly due to the ongoing starburst, accretion of pristine gas, or a late stage merger.
Leonid Sajkov, Anna Sajina, Alexandra Pope, Stacey Alberts, Lee Armus, Duncan Farrah, Jamie Lin, Danilo Marchesini, Jed McKinney, Sylvain Veilleux, Lin Yan, Jason Young We present an analysis of 8 JWST Mid-Infrared Instrument 5.6 micron images with 5sigma depths of ~0.1uJy. We detect 2854 sources within our combined area of 18.4 sq.arcmin -- a >4x increase in source density over earlier IRAC channel 3 data. We compute the MIRI 5.6um number counts including an analysis of the field-to-field variation. Relative to earlier published MIRI 5.6micron counts, our counts have a more pronounced knee, at roughly 2\,$\mu$Jy. The location and amplitude of the counts at the knee are consistent with the Cowley et al. (2018) model predictions, although these models tend to overpredict the counts below the knee. In areas of overlap, 84% of the MIRI sources have a counterpart in the COSMOS2020 catalog. These MIRI sources have redshifts that are mostly in the z~0.5-2, with a tail out to z~5. They are predominantly moderate to low stellar masses $10^8-10^{10}$M$_{\odot}$) main sequence star-forming galaxies suggesting that with $\approx$2hr exposures, MIRI can reach well below $M^*$ at cosmic noon and reach higher mass systems out to z~5. Nearly 70% of the COSMOS2020 sources in areas of overlap now have a data point at 5.6micron (rest-frame near-IR at cosmic noon) which allows for more accurate stellar population parameter estimates. Finally, we discover 31 MIRI-bright sources not in COSMOS2020. A cross-match with IRAC channel 1 suggests that 10-20% of these are likely lower mass (M$_*\approx10^9$M$_{\odot}$), $z\sim1$ dusty galaxies. The rest (80--90%) are consistent with more massive, but still very dusty galaxies at z>3.
K. M. Dasyra, G. F. Paraschos, F. Combes, P. Patapis, G. Helou, M. Papachristou, J. A. Fernandez-Ontiveros, T. G. Bisbas, L. Spinoglio, L. Armus, M. Malkan We present James Webb Space Telescope MIRI data of the inner ~3x2 kpc^2 of the galaxy IC5063, in which the jets of a supermassive black hole interact with the gaseous disk they are crossing. Jet-driven outflows were known to be initiated along or near the jet path, and the stability conditions of clouds were known to vary because of these outflows. The MIRI data, of unprecedented resolution and sensitivity in the infrared, now reveal that there are more than ten discrete regions with outflows, nearly doubling the number of such known regions. Outflows exist near the radio lobes, at the nucleus, in a biconical structure perpendicular to the jet, and in a bubble moving against the disk. In some of them, velocities above escape velocity are observed. Stratification is also observed, with higher ionization or excitation gas attaining higher velocities. More outflows and bow shocks, found further away from the nucleus than the radio lobes, in regions without significant radio emission, reveal the existence of past or weak radio jets that interacted with the interstellar medium. The coincidence of the bow shocks with the optical extended emission line region (EELR) suggests that the jets also contributed to the gas ionization. Maps of the H2 gas excitation temperature, T_ex, indicate that the molecular gas is most excited in regions with radio emission. There, T_ex is more than 100K higher than in the EELR interior. We argue that a combination of jet-related shocks and cosmic rays is likely responsible for this excess molecular gas excitation.
Cory M. Whitcomb, J.-D. T. Smith, Karin Sandstrom, Carl A. Starkey, Grant P. Donnelly, Bruce T. Draine, Evan D. Skillman, Daniel A. Dale, Lee Armus, Brandon S. Hensley, Thomas S.-Y. Lai, Robert C. Kennicutt We use deep Spitzer mid-infrared spectroscopic maps of radial strips across three nearby galaxies with well-studied metallicity gradients (M101, NGC 628, and NGC 2403) to explore the physical origins of the observed deficit of polycyclic aromatic hydrocarbons (PAHs) at sub-solar metallicity (i.e. the PAH-metallicity relation or PZR). These maps allow us to trace the evolution of all PAH features from 5-18 $\mu$m as metallicity decreases continuously from solar ($Z_\odot$) to 0.2 $Z_\odot$. The total PAH to dust luminosity ratio remains relatively constant until reaching a threshold of $\sim$$\frac{2}{3}$$Z_\odot$, below which it declines smoothly but rapidly. The PZR has been attributed to preferential destruction of the smallest grains in the hard radiation environments found at low metallicity. In this scenario, a decrease in emission from the shortest wavelength PAH features is expected. In contrast, we find a steep decline in long wavelength power below $Z_\odot$, especially in the 17 $\mu$m feature, with the shorter wavelength PAH bands carrying an increasingly large fraction of power at low metallicity. We use newly developed grain models to reproduce the observed PZR trends, including these variations in fractional PAH feature strengths. The model that best reproduces the data employs an evolving grain size distribution that shifts to smaller sizes as metallicity declines. We interpret this as a result of inhibited grain growth at low metallicity, suggesting continuous replenishment in the interstellar medium is the dominant process shaping the PAH grain population in galaxies.
L. Bisigello, C. Gruppioni, A. Bolatto, L. Ciesla, A. Pope, L. Armus, L., J.D. Smith, R. Somerville, L.Y.A. Yung, R.J. Wright, C.M. Bradford, J. Glenn, A. Feltre The most active phases of star formation and black hole accretion are strongly affected by dust extinction, making far-infrared (far-IR) observations the best way to disentangle and study the co-evolution of galaxies and super massive black holes. The plethora of fine structure lines and emission features from dust, ionised and neutral atomic and warm molecular gas in the rest-frame mid- and far-IR provide unmatched diagnostic power to determine the properties of gas and dust, measure gas-phase metallicities and map cold galactic outflows in even the most obscured galaxies. By combining multi-band photometric surveys with low and high-resolution far-IR spectroscopy, the PRobe far-Infrared Mission for Astrophysics (PRIMA), a concept for a far-IR, 1.8m-diameter, cryogenically cooled observatory, will revolutionise the field of galaxy evolution by taking advantage of this IR toolkit to find and study dusty galaxies across galactic time. In this work, we make use of the phenomenological simulation SPRITZ and the Santa Cruz semi-analytical model to describe how a moderately deep multi-band PRIMA photometric survey can easily reach beyond previous IR missions to detect and study galaxies down to $10^{11}\,L_{\odot}$ beyond cosmic noon and at least up to z=4, even in the absence of gravitational lensing. By decomposing the spectral energy distribution (SED) of these photometrically selected galaxies, we show that PRIMA can be used to accurately measure the relative AGN power, the mass fraction contributed by polycyclic aromatic hydrocarbon (PAH) and the total IR luminosity. At the same time, spectroscopic follow up with PRIMA will allow to trace both the star formation and black hole accretion rates (SFR, BHAR), the gas phase metallicities and the mass outflow rates of cold gas in hundreds to thousands of individual galaxies to z=2.
Thomas Bohn, Hanae Inami, Aditya Togi, Lee Armus, Thomas S.-Y. Lai, Loreto Barcos-Munoz, Yiqing Song, Sean T. Linden, Jason Surace, Marina Bianchin, Vivian U, Aaron S. Evans, Torsten Böker, Matthew A. Malkan, Kirsten L. Larson, Sabrina Stierwalt, Victorine A. Buiten, Vassilis Charmandaris, Tanio Diaz-Santos, Justin H. Howell, et al (9) We present James Webb Space Telescope (JWST) Integral Field Spectrograph observations of NGC 3256, a local infrared-luminous late-stage merging system with two nuclei about 1 kpc apart, both of which have evidence of cold molecular outflows. Using JWST NIRSpec and MIRI datasets, we investigate this morphologically complex system on spatial scales of $<$100 pc, where we focus on the warm molecular H$_2$ gas surrounding the nuclei. We detect collimated outflowing warm H$_2$ gas originating from the southern nucleus, though we do not find significant outflowing warm H$_2$ gas surrounding the northern nucleus. Within the observed region, the maximum intrinsic velocities of the outflow reach up to $\sim$1,000 km s$^{-1}$, and extend out to a distance of 0.7 kpc. Based on H$_2$ S(7)/S(1) line ratios, we find a larger fraction of warmer gas near the S nucleus, which decreases with increasing distance from the nucleus, signifying the S nucleus as a primary source of H$_2$ heating. The gas mass of the warm H$_2$ outflow component is estimated to be $M\rm{_{warm,out}}$ = 8.9$\times$10$^5\;M_{\odot}$, as much as 4$\%$ of the cold H$_2$ mass as estimated using ALMA CO data. The outflow time scale is about $7\times10^5$ yr, resulting in a mass outflow rate of $\dot{M}\rm{_{warm,out}}$ = 1.3 M$_{\odot}$ yr$^{-1}$ and kinetic power of $P\rm{_{warm,out}}\;\sim\;2\times10^{41}$ erg s$^{-1}$. Lastly, the regions where the outflowing gas reside show high [FeII]/Pa$\beta$ and H$_2$/Br$\gamma$ line ratios, indicating enhanced mechanical heating caused by the outflows. At the same time, the 3.3 $\mu$m and 6.2 $\mu$m Polycyclic Aromatic Hydrocarbon fluxes in these regions are not significantly suppressed compared to those outside the outflows, suggesting the outflows have no clear negative feedback effect on the local star formation.
We present a focused study of radially-resolved varying PAH emission in the low-luminosity AGN-host NGC 4138 using deep Spitzer/IRS spectral maps. Using new model PAH spectra, we investigate whether these variations could be associated with changes to the PAH grain size distribution due to photodestruction by the AGN. Separately, we model the effects of the varying radiation field within NGC 4138, and we use this model to predict the corresponding changes in the PAH emission spectrum. We find that PAH band ratios are strongly variable with radius in this galaxy with short-to-long wavelength band ratios peaking in the starburst ring. The changing mix of starlight appears to have a considerable effect on the trends in these band ratios, and our radiation model predicts the shapes of these trends. However, the amplitude of observed variation is ~2.5 times larger than predicted for some ratios. A cutoff of small grains in the PAH size distribution, as has been suggested for AGN, together with changes in PAH ionization fraction could explain the behavior of the shorter bands, but this model fails to reproduce longer band behaviors. Additionally, we find that short-to-long wavelength PAH band ratios increase slightly within ~270pc of the center, suggesting that the AGN may directly influence PAH emission there.
Alberto D. Bolatto, Rebecca C. Levy, Elizabeth Tarantino, Martha L. Boyer, Deanne B. Fisher, Adam K. Leroy, Serena A. Cronin, Ralf S. Klessen, J. D. Smith, Dannielle A. Berg, Torsten Boeker, Leindert A. Boogaard, Eve C. Ostriker, Todd A. Thompson, Juergen Ott, Laura Lenkic, Laura A. Lopez, Daniel A. Dale, Sylvain Veilleux, Paul P. van der Werf, et al (20) We present new observations of the central 1 kpc of the M 82 starburst obtained with the James Webb Space Telescope (JWST) near-infrared camera (NIRCam) instrument at a resolution ~0.05"-0.1" (~1-2 pc). The data comprises images in three mostly continuum filters (F140M, F250M, and F360M), and filters that contain [FeII] (F164N), H2 v=1-0 (F212N), and the 3.3 um PAH feature (F335M). We find prominent plumes of PAH emission extending outward from the central starburst region, together with a network of complex filamentary substructure and edge-brightened bubble-like features. The structure of the PAH emission closely resembles that of the ionized gas, as revealed in Paschen alpha and free-free radio emission. We discuss the origin of the structure, and suggest the PAHs are embedded in a combination of neutral, molecular, and photoionized gas.
We present results from the ``Quasar hosts Unveiled by high Angular Resolution Techniques" (QUART) survey studying the Circumgalactic Medium (CGM) by observing rest-frame UV emission lines Ly$\alpha$, C IV and He II around two radio-loud quasars, 3C 9 (z=2.02) and 4C 05.84 (z=2.32), using Keck Cosmic Web Imager (KCWI). We detect large-scale Ly$\alpha$ nebulae around both quasars with projected diameters $\sim$ 100 kpc, with spatially resolved, embedded 15-30 kpc He II and C IV nebulae around both quasars as well as kinematically distinct He II and C IV nebulae at a physical separation of $\sim$ 15 kpc from both quasars. Observations of H$\alpha$, H$\beta$, and [O III] emission using Keck MOSFIRE spectroscopically confirm that the Ly$\alpha$ nebulae extend to companion galaxies and that these quasars are in a protogroup/protocluster environment. We confirm that the He II and C IV emission is kinematically and spatially coincident with the companion galaxies. We estimate the virial masses of the companion galaxies, their metallicities, and star formation rates, and investigate the sources of ionization. We measure the dynamical mass of the host dark matter halos and estimate that the dark matter halos of these systems will grow to a mass of 2 $\times 10^{14}$ M$_{\odot}$ (3C 9) and 2 $\times 10^{13}$ M$_{\odot}$ (4C 05.84) by z=0. The combined CGM and companion galaxies observations indicate Ly$\alpha$ substructure can indicate the presence of companion galaxies in the CGM.
Victorine A. Buiten, Paul P. van der Werf, Serena Viti, Lee Armus, Andrew G. Barr, Loreto Barcos-Muñoz, Aaron S. Evans, Hanae Inami, Sean T. Linden, George C. Privon, Yiqing Song, Jeffrey A. Rich, Susanne Aalto, Philip N. Appleton, Torsten Böker, Vassilis Charmandaris, Tanio Diaz-Santos, Christopher C. Hayward, Thomas S.-Y. Lai, Anne M. Medling, et al (2) The enormous increase in mid-IR sensitivity and spatial and spectral resolution provided by the JWST spectrographs enables, for the first time, detailed extragalactic studies of molecular vibrational bands. This opens an entirely new window for the study of the molecular interstellar medium in luminous infrared galaxies (LIRGs). We present a detailed analysis of rovibrational bands of gas-phase CO, H$_2$O, C$_2$H$_2$ and HCN towards the heavily-obscured eastern nucleus of the LIRG VV 114, as observed by NIRSpec and MIRI MRS. Spectra extracted from apertures of 130 pc in radius show a clear dichotomy between the obscured AGN and two intense starburst regions. We detect the 2.3 $\mu$m CO bandheads, characteristic of cool stellar atmospheres, in the star-forming regions, but not towards the AGN. Surprisingly, at 4.7 $\mathrm{\mu}$m we find highly-excited CO ($T_\mathrm{ex} \approx 700-800$ K out to at least rotational level $J = 27$) towards the star-forming regions, but only cooler gas ($T_\mathrm{ex} \approx 200$ K) towards the AGN. We conclude that only mid-infrared pumping through the rovibrational lines can account for the equilibrium conditions found for CO and H$_2$O in the deeply-embedded starbursts. Here the CO bands probe regions with an intense local radiation field inside dusty young massive star clusters or near the most massive young stars. The lack of high-excitation molecular gas towards the AGN is attributed to geometric dilution of the intense radiation from the bright point source. An overview of the relevant excitation and radiative transfer physics is provided in an appendix.
Jason Young, Alexandra Pope, Anna Sajina, Lin Yan, Thiago S Goncalves, Miriam Eleazer, Stacey Alberts, Lee Armus, Matteo Bonato, Daniel A. Dale, Duncan Farrah, Carl Ferkinhoff, Christopher C. Hayward, Jed McKinney, Eric J. Murphy, Nicole Nesvadba, Patrick Ogle, Leonid Sajkov, Sylvain Veilleux We present JWST/MIRI/MRS observations of an infrared luminous disk galaxy, FLS1, at z=0.54. With a lookback time of 5 Gyr, FLS1 is chronologically at the midpoint between the peak epoch of star formation and the present day. The MRS data provide maps of the atomic fine structure lines [Ar II]6.99 micron, [Ar III]8.99 micron, [Ne II]12.81 micron, and [Ne III]15.55 micron, polycyclic aromatic hydrocarbon (PAH) features at 3.3 micron, 6.2 micron, and 11.3 micron, and the warm molecular gas indicators H2S(5) and H2S(3); all these emission features are spatially resolved. We find that the PAH emission is more extended along the Northern side of the galaxy when compared to the well-studied star-formation tracer [Ne II]. The H2 rotational lines, which are shock indicators, are strongest and most extended on the Southern side of the galaxy. [Ar II] is the second brightest fine structure line detected in FLS1 and we show that it is a useful kinematic probe which can be detected with JWST out to z=3. Velocity maps of [Ar II] show a rotating disk with signs of turbulence. Our results provide an example of how spatially resolved mid-infrared spectroscopy can allow us to better understand the star formation and ISM conditions in a galaxy halfway back to the peak epoch of galaxy evolution.
Rebecca C. Levy, Alberto D. Bolatto, Elizabeth Tarantino, Adam K. Leroy, Lee Armus, Kimberly L. Emig, Rodrigo Herrera-Camus, Daniel P. Marrone, Elisabeth Mills, Oliver Ricken, Juergen Stutzki, Sylvain Veilleux, Fabian Walter M82 is an archetypal starburst galaxy in the local Universe. The central burst of star formation, thought to be triggered by M82's interaction with other members in the M81 group, is driving a multiphase galaxy-scale wind away from the plane of the disk that has been studied across the electromagnetic spectrum. Here, we present new velocity-resolved observations of the [CII] 158$\mu$m line in the central disk and the southern outflow of M82 using the upGREAT instrument onboard SOFIA. We also report the first detections of velocity-resolved ($\Delta V = 10$ km s$^{-1}$) [CII] emission in the outflow of M82 at projected distances of $\approx1-2$ kpc south of the galaxy center. We compare the [CII] line profiles to observations of CO and HI and find that likely the majority ($>55$%) of the [CII] emission in the outflow is associated with the neutral atomic medium. We find that the fraction of [CII] actually outflowing from M82 is small compared to the bulk gas outside the midplane (which may be in a halo or tidal streamers), which has important implications for observations of [CII] outflows at higher redshift. Finally, by comparing the observed ratio of the [CII] and CO intensities to models of photodissociation regions, we estimate that the far-ultraviolet (FUV) radiation field in the disk is $\sim10^{3.5}~G_0$, in agreement with previous estimates. In the outflow, however, the FUV radiation field is 2-3 orders of magnitudes lower, which may explain the high fraction of [CII] arising from the neutral medium in the wind.
Marina Bianchin, Vivian U, Yiqing Song, Thomas S.-Y. Lai, Raymond P. Remigio, Loreto Barcos-Munoz, Tanio Diaz-Santos, Lee Armus, Hanae Inami, Kirsten L. Larson, Aaron S. Evans, Torsten Boker, Justin A. Kader, Sean T. Linden, Vassilis Charmandaris, Matthew A. Malkan, Jeff Rich, Thomas Bohn, Anne M. Medling, Sabrina Stierwalt, et al (19) We present new JWST-NIRSpec IFS data for the luminous infrared galaxy NGC7469: a nearby (70.6Mpc) active galaxy with a Sy 1.5 nucleus that drives a highly ionized gas outflow and a prominent nuclear star-forming ring. Using the superb sensitivity and high spatial resolution of the JWST instrument NIRSpec-IFS, we investigate the role of the Seyfert nucleus in the excitation and dynamics of the circumnuclear gas. Our analysis focuses on the [Fe ii], H2, and hydrogen recombination lines that trace the radiation/shocked-excited molecular and ionized ISM around the AGN. We investigate the gas excitation through H2/Br\gamma and [Fe ii]/Pae̱ta emission line ratios and find that photoionization by the AGN dominates within the central 300 pc of the galaxy and together with a small region show ing signatures of shock-heated gas; these shock-heated regions are likely associated with a compact radio jet. In addition, the velocity field and velocity dispersion maps reveal complex gas kinematics. Rotation is the dominant feature, but we also identify non-circular motions consistent with gas inflows as traced by the velocity residuals and the spiral pattern in the Pa\alpha velocity dispersion map. The inflow is consistent with the mass outflow rate and two orders of magnitude higher than the AGN accretion rate. The compact nuclear radio jet has enough power to drive the highly ionized outflow. This scenario suggests that the inflow and outflow are in a self-regulating feeding-feedback process, with a contribution from the radio jet helping to drive the outflow.
Thomas S.-Y. Lai, Lee Armus, Marina Bianchin, Tanio Diaz-Santos, Sean T. Linden, George C. Privon, Hanae Inami, Vivian U, Thomas Bohn, Aaron S. Evans, Kirsten L. Larson, Brandon S. Hensley, J.-D.T. Smith, Matthew A. Malkan, Yiqing Song, Sabrina Stierwalt, Paul P. van der Werf, Jed McKinney, Susanne Aalto, Victorine A. Buiten, et al (19) We present James Webb Space Telescope (JWST) Near Infrared Spectrograph (NIRSpec) integral-field spectroscopy of the nearby luminous infrared galaxy, NGC 7469. We take advantage of the high spatial/spectral resolution and wavelength coverage of JWST /NIRSpec to study the 3.3 um neutral polycyclic aromatic hydrocarbon (PAH) grain emission on ~60 pc scales. We find a clear change in the average grain properties between the star-forming ring and the central AGN. Regions in the vicinity of the AGN, with [NeIII]/[NeII]>0.25, tend to have larger grain sizes and lower aliphatic-to-aromatic (3.4/3.3) ratios indicating that smaller grains are preferentially removed by photo-destruction in the vicinity of the AGN. We find an overall suppression of the total PAH emission relative to the ionized gas in the central 1 kpc region of the AGN in NGC 7469 compared to what has been observed with Spitzer on 3 kpc scales. However, the fractional 3.3 um to total PAH power is enhanced in the starburst ring, possibly due to a variety of physical effects on sub-kpc scales, including recurrent fluorescence of small grains or multiple photon absorption by large grains. Finally, the IFU data show that while the 3.3 um PAH-derived star formation rate (SFR) in the ring is 8% higher than that inferred from the [NeII] and [NeIII] emission lines, the integrated SFR derived from the 3.3 um feature would be underestimated by a factor of two due to the deficit of PAHs around the AGN, as might occur if a composite system like NGC 7469 were to be observed at high-redshift.
Jed McKinney, Alexandra Pope, Allison Kirkpatrick, Lee Armus, Tanio Diaz-Santos, Carlos Gomez-Guijarro, Maximilien Franco, David Elbaz, Christopher C. Hayward, Hanae Inami, Gergo Popping, Mengyuan Xiao Surface densities of gas, dust and stars provide a window into the physics of star-formation that, until the advent of high-resolution far-infrared/sub-millimeter observations, has been historically difficult to assess amongst dusty galaxies. To study the link between infrared (IR) surface densities and dust properties, we leverage the Atacama Large Millimetre/Submillimetre Array (ALMA) archive to measure the extent of cold dust emission in 15 $z\sim2$ IR selected galaxies selected on the basis of having available mid-IR spectroscopy from Spitzer. We use the mid-IR spectra to constrain the relative balance between dust heating from star-formation and active galactic nuclei (AGN), and to measure emission from Polycylic Aromatic Hydrocarbons (PAHs) -- small dust grains that play a key role in the photoelectric heating of gas. In general, we find that dust-obscured star-formation at high IR surface densities exhibits similar properties at low- and high-redshift, namely: local luminous IR galaxies have comparable PAH luminosity to total dust mass ratios as high-$z$ galaxies, and star-formation at $z\sim0-2$ is more efficient at high IR surface densities despite the fact that our sample of high$-z$ galaxies are closer to the main-sequence than local luminous IR galaxies. High star-formation efficiencies are coincident with a decline in the PAH/IR luminosity ratio reminiscent of the deficit observed in far-infrared fine-structure lines. Changes in the gas and dust conditions arising from high star-formation surface densities might help drive the star-formation efficiency up. This could help explain high efficiencies needed to reconcile star-formation and gas volume densities in dusty galaxies at cosmic noon.
J. Rich, S. Aalto, A.S. Evans, V. Charmandaris, G. C. Privon, T. Lai, H. Inami, S. Linden, L. Armus, T. Diaz-Santos, P. Appleton, L. Barcos-Muñoz, T. Böker, K. L. Larson, D. R. Law, M. A. Malkan, A. M. Medling, Y. Song, V. U, P. van der Werf, et al (16) We present results from the James Webb Space Telescope (JWST) Director's Discretionary Time Early Release Science (ERS) program 1328 targeting the nearby, Luminous Infrared Galaxy (LIRG), VV 114. We use the MIRI and NIRSpec instruments to obtain integral-field spectroscopy of the heavily obscured Eastern nucleus (V114E) and surrounding regions. The spatially resolved, high-resolution, spectra reveal the physical conditions in the gas and dust over a projected area of 2-3 kpc that includes the two brightest IR sources, the NE and SW cores. Our observations show for the first time spectroscopic evidence that the SW core hosts an AGN as evidenced by its very low 6.2 \mum and 3.3 \mum PAH equivalent widths (0.12 and 0.017 \mum respectively) and mid and near-IR colors. Our observations of the NE core show signs of deeply embedded star formation including absorption features due to aliphatic hydrocarbons, large quantities of amorphous silicates, as well as HCN due to cool gas along the line of sight. We detect elevated [Fe II]/Pf\alpha consistent with extended shocks coincident with enhanced emission from warm H$_{2}$, far from the IR-bright cores and clumps. We also identify broadening and multiple kinematic components in both H$_{2}$ and fine structure lines caused by outflows and previously identified tidal features.
Jed McKinney, Luke Finnerty, Caitlin Casey, Maximilien Franco, Arianna Long, Seiji Fujimoto, Jorge Zavala, Olivia Cooper, Hollis Akins, Alexandra Pope, Lee Armus, B. T. Soifer, Kirsten Larson, Keith Matthews, Jason Melbourne, Michael Cushing Selecting the first galaxies at z>7-10 from JWST surveys is complicated by z<6 contaminants with degenerate photometry. For example, strong optical nebular emission lines at z<6 may mimic JWST/NIRCam photometry of z>7-10 Lyman Break Galaxies (LBGs). Dust-obscured 3<z<6 galaxies in particular are potentially important contaminants, and their faint rest-optical spectra have been historically difficult to observe. A lack of optical emission line and continuum measures for 3<z<6 dusty galaxies now makes it difficult to test their expected JWST/NIRCam photometry for degenerate solutions with NIRCam dropouts. Towards this end, we quantify the contribution by strong emission lines to NIRCam photometry in a physically motivated manner by stacking 21 Keck II/NIRES spectra of hot, dust-obscured, massive ($\log\mathrm{M_*/M_\odot}\gtrsim10-11$) and infrared (IR) luminous galaxies at z~1-4. We derive an average spectrum and measure strong narrow (broad) [OIII]5007 and H$\alpha$ features with equivalent widths of $130\pm20$ A ($150\pm50$ A) and $220\pm30$ A ($540\pm80$ A) respectively. These features can increase broadband NIRCam fluxes by factors of 1.2-1.7 (0.2-0.6 mag). Due to significant dust-attenuation ($A_V\sim6$), we find H$\alpha$+[NII] to be significantly brighter than [OIII]+H$\beta$, and therefore find that emission-line dominated contaminants of high-z galaxy searches can only reproduce moderately blue perceived UV continua of $S_\lambda\propto\lambda^\beta$ with $\beta>-1.5$ and z>4. While there are some redshifts (z~3.75) where our stack is more degenerate with the photometry of z>10 LBGs between $\lambda_{rest}\sim0.3-0.8\,\mu$m, redder filter coverage beyond $\lambda_{obs}>3.5\,\mu$m and far-IR/sub-mm follow-up may be useful for breaking the degeneracy and making a crucial separation between two fairly unconstrained populations, dust-obscured galaxies at z~3-6 and LBGs at z>10.
We present Keck Cosmic Web Imager (KCWI) integral field spectroscopy (IFS) observations of rest-frame UV emission lines $\rm Ly\alpha$, C IV $\lambda \lambda$ 1548 Å, 1550Å and He II 1640 Å observed in the circumgalactic medium (CGM) of two $z=2$ radio-loud quasar host galaxies. We detect extended emission on 80-90 kpc scale in $\rm Ly\alpha$ in both systems with C IV, and He II emission also detected out to 30-50 kpc. All emission lines show kinematics with a blue and redshifted gradient pattern consistent with velocities seen in massive dark matter halos and similar to kinematic patterns of inflowing gas seen in hydrodynamical simulations. Using the kinematics of both resolved $\rm Ly\alpha$ emission and absorption, we can confirm that both kinematic structures are associated with accretion. Combining the KCWI data with molecular gas observations with Atacama Large Millimeter/submillimeter Array (ALMA) and high spatial resolution of ionized gas with Keck OSIRIS, we find that both quasar host galaxies reside in proto-group environments at $z=2$. We estimate $1-6\times10^{10}$M$_\odot$ of warm-ionized gas within 30-50 kpc from the quasar that is likely accreting onto the galaxy group. We estimate inflow rates of 60-200 M$_\odot$yr$^{-1}$, within an order of magnitude of the outflow rates in these systems. In the 4C 09.17 system, we detect narrow gas streams associated with satellite galaxies, potentially reminiscent of ram-pressure stripping seen in local galaxy groups and clusters. We find that the quasar host galaxies reside in dynamically complex environments, with ongoing mergers, gas accretion, ISM stripping, and outflows likely playing an important role in shaping the assembly and evolution of massive galaxies at cosmic noon.
Sean T. Linden, Aaron S. Evans, Lee Armus, Jeffrey A. Rich, Kirsten L. Larson, Thomas Lai, George C. Privon, Vivian U, Hanae Inami, Thomas Bohn, Yiqing Song, Loreto Barcos-Muñoz, Vassilis Charmandaris, Anne M. Medling, Sabrina Stierwalt, Tanio Diaz-Santos, Torsten Böker, Paul van der Werf, Susanne Aalto, Philip Appleton, et al (13) We present the results of a \it James Webb Space Telescope NIRCam investigation into the young massive star cluster (YMC) population in the luminous infrared galaxy VV 114. We identify 374 compact YMC candidates with a $S/N \geq 3$, 5, and 5 at F150W, F200W, and F356W respectively. A direct comparison with our \it HST cluster catalog reveals that $\sim 20\%$ of these sources are undetected at optical wavelengths. Based on \it yggdrasil stellar population models, we identify 17 YMC candidates in our \it JWST imaging alone with F150W-F200W and F200W-F356W colors suggesting they are all very young, dusty ($A_{V} = 5 - 15$), and massive ($10^{5.8} < M_{\odot} < 10^{6.1}$). The discovery of these `hidden' sources, many of which are found in the `overlap' region between the two nuclei, quadruples the number of $t < 3$ Myr clusters, and nearly doubles the number of $t < 6$ Myr clusters detected in VV 114. Now extending the cluster age distribution ($dN/d\tau \propto \tau^{\gamma}$) to the youngest ages, we find a slope of $\gamma = -1.30 \pm 0.39$ for $10^{6} < \tau (\mathrm{yr}) < 10^{7}$, which is consistent with the previously determined value from $10^{7} < \tau (\mathrm{yr}) < 10^{8.5}$, and confirms that VV 114 has a steep age distribution slope for all massive star clusters across the entire range of cluster ages observed. Finally, the consistency between our \it JWST- and \it HST-derived age distribution slopes indicates that the balance between cluster formation and destruction has not been significantly altered in VV 114 over the last 0.5 Gyr.
L. Armus, T. Lai, V. U, K.L. Larson, T. Diaz-Santos, A.S. Evans, M.A. Malkan, J. Rich, A.M. Medling, D.R. law, H. Inami, F. Muller-Sanchez, V. Charmandaris, P. can der Werf, S. Stierwalt, S. Linden, G.C. Privon, L. Barcos-Munoz, C. Hayward, Y. Song, et al (15) We present mid-infrared spectroscopic observations of the nucleus of the nearby Seyfert galaxy NGC 7469 taken with the MIRI instrument on the James Webb Space Telescope (JWST) as part of Directors Discretionary Time Early Release Science (ERS) program 1328. The high resolution nuclear spectrum contains 19 emission lines covering a wide range of ionization. The high ionization lines show broad, blueshifted emission reaching velocities up to 1700 km s$^{-1}$ and FWHM ranging from $\sim500 - 1100$ km s$^{-1}$. The width of the broad emission and the broad to narrow line flux ratios correlate with ionization potential. The results suggest a decelerating, stratified, AGN driven outflow emerging from the nucleus. The estimated mass outflow rate is one to two orders of magnitude larger than the current black hole accretion rate needed to power the AGN. Eight pure rotational H$_{2}$ emission lines are detected with intrinsic widths ranging from FWHM $\sim 125-330$ km s$^{-1}$. We estimate a total mass of warm H$_{2}$ gas of $\sim1.2\times10^{7}$M$_{\odot}$ in the central 100 pc. The PAH features are extremely weak in the nuclear spectrum, but a $6.2\mu$m PAH feature with an equivalent width $\sim0.07\mu$m and a flux of $2.7\times10^{-17}$ W m$^{-2}$ is detected. The spectrum is steeply rising in the mid-infrared, with a silicate strength $\sim0.02$, significantly smaller than seen in most PG QSOs, but comparable to other Seyfert 1's. These early MIRI mid-infrared IFU data highlight the power of JWST to probe the multi-phase interstellar media surrounding actively accreting supermassive black holes.
Thomas S. -Y. Lai, Lee Armus, Vivian U, Tanio Diaz-Santos, Kirsten L. Larson, Aaron Evans, Matthew A. Malkan, Philip Appleton, Jeff Rich, Francisco Muller-Sanchez, Hanae Inami, Thomas Bohn, Jed McKinney, Luke Finnerty, David R. Law, Sean Linden, Anne M. Medling, George C. Privon, Yiqing Song, Sabrina Stierwalt, et al (17) We present James Webb Space Telescope (JWST) Mid-InfraRed Instrument (MIRI) integral-field spectroscopy of the nearby merging, luminous infrared galaxy, NGC 7469. This galaxy hosts a Seyfert type-1.5 nucleus, a highly ionized outflow, and a bright, circumnuclear star-forming ring, making it an ideal target to study AGN feedback in the local Universe. We take advantage of the high spatial/spectral resolution of JWST/MIRI to isolate the star-forming regions surrounding the central active nucleus and study the properties of the dust and warm molecular gas on ~100 pc scales. The starburst ring exhibits prominent Polycyclic Aromatic Hydrocarbon (PAH) emission, with grain sizes and ionization states varying by only ~30%, and a total star formation rate of $\rm 10 - 30 \ M_\odot$/yr derived from fine structure and recombination emission lines. Using pure rotational lines of H2, we detect 1.2$\times$10$^{7} \rm \ M_\odot$ of warm molecular gas at a temperature higher than 200 K in the ring. All PAH bands get significantly weaker towards the central source, where larger and possibly more ionized grains dominate the emission. However, the bulk of the dust and molecular gas in the ring appears unaffected by the ionizing radiation or the outflowing wind from the AGN. These observations highlight the power of JWST to probe the inner regions of dusty, rapidly evolving galaxies for signatures of feedback and inform models that seek to explain the co-evolution of supermassive black holes and their hosts.
Thomas Bohn, Hanae Inami, Tanio Diaz-Santos, Lee Armus, Sean T. Linden, Vivian U, Jason Surace, Kirsten L. Larson, Aaron S. Evans, Shunshi Hoshioka, Thomas Lai, Yiqing Song, Joseph M. Mazzarella, Loreto Barcos-Munoz, Vassilis Charmandaris, Justin H. Howell, Anne M. Medling, George C. Privon, Jeffrey A. Rich, Sabrina Stierwalt, et al (13) We present James Webb Space Telescope (JWST) imaging of NGC 7469 with the Near-Infrared Camera (NIRCam) and the Mid-InfraRed Instrument (MIRI). NGC 7469 is a nearby, $z=0.01627$, luminous infrared galaxy (LIRG) that hosts both a Seyfert Type-1.5 nucleus and a circumnuclear starburst ring with a radius of $\sim$0.5 kpc. The new near-infrared (NIR) JWST imaging reveals 66 star-forming regions, 37 of which were not detected by HST observations. Twenty-eight of the 37 sources have very red NIR colors that indicate obscurations up to A$_{\rm{v}}\sim7$ and a contribution of at least 25$\%$ from hot dust emission to the 4.4$\mu$m band. Their NIR colors are also consistent with young ($<$5 Myr) stellar populations and more than half of them are coincident with the MIR emission peaks. These younger, dusty star-forming regions account for $\sim$6$\%$ and $\sim$17$\%$ of the total 1.5$\mu$m and 4.4$\mu$m luminosity of the starburst ring, respectively. Thanks to JWST, we find a significant number of young dusty sources that were previously unseen due to dust extinction. The newly identified 28 young sources are a significant increase compared to the number of HST-detected young sources (4-5). This makes the total percentage of the young population rise from $\sim$15$\%$ to 48$\%$. These results illustrate the effectiveness of JWST in identifying and characterizing previously hidden star formation in the densest star-forming environments around AGN.
Y. Song, S. T. Linden, A. S. Evans, L. Barcos-Munoz, E. J. Murphy, E. Momjian, T. Diaz-Santos, K. L. Larson, G. C. Privon, X. Huang, L. Armus, J. M. Mazzarella, V. U, H. Inami, V. Charmandaris, C. Ricci, K. L. Emig, J. McKinney, I. Yoon, D. Kunneriath, et al (6) We present the analysis of $\sim 100$pc-scale compact radio continuum sources detected in 63 local (Ultra) Luminous Infrared Galaxies (U/LIRGs; $L_{\rm IR} \ge 10^{11} L_\odot$), using FWHM $\lesssim 0''.1 - 0''.2$ resolution 15 and 33 GHz observations with the Karl G. Jansky Very Large Array. We identify a total of 133 compact radio sources with effective radii of 8 - 170pc, which are classified into four main categories -- "AGN" (AGN), "AGN/SBnuc" (AGN-starburst composite nucleus), "SBnuc" (starburst nucleus) and "SF" (star-forming clumps) -- based on ancillary datasets and the literature. We find that "AGN" and "AGN/SBnuc" more frequently occur in late-stage mergers and have up to 3 dex higher 33 GHz luminosities and surface densities compared with "SBnuc" and "SF", which may be attributed to extreme nuclear starburst and/or AGN activity in the former. Star formation rates (SFRs) and surface densities ($\Sigma_{\rm SFR}$) are measured for "SF" and "SBnuc" using both the total 33 GHz continuum emission (SFR $\sim 0.14 - 13$ M$_\odot$ yr$^{-1}$, $\Sigma_{\rm SFR} \sim 13 - 1600$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$) and the thermal free-free emission from HII regions (median SFR$_{\rm th} \sim 0.4$ M$_\odot$ yr$^{-1}$, $\Sigma_{\rm SFR_{th}} \sim 44$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$). These values are 1 - 2 dex higher than those measured for similar-sized clumps in nearby normal (non-U/LIRGs). The latter also have much flatter median 15 - 33 GHz spectral index ($\sim -0.08$) compared with "SBnuc" and "SF" ($\sim -0.46$), which may reflect higher non-thermal contribution from supernovae and/or ISM densities in local U/LIRGs that directly result from and/or lead to their extreme star-forming activities on 100\u2009pc scales.
Vivian U, Thomas Lai, Marina Bianchin, Raymond P. Remigio, Lee Armus, Kirsten L. Larson, Tanio Diaz-Santos, Aaron Evans, Sabrina Stierwalt, David R. Law, Matthew A. Malkan, Sean Linden, Yiqing Song, Paul P. van der Werf, Tianmu Gao, George C. Privon, Anne M. Medling, Loreto Barcos-Muñoz, Christopher C. Hayward, Hanae Inami, et al (18) The nearby, luminous infrared galaxy (LIRG) NGC 7469 hosts a Seyfert nucleus with a circumnuclear star-forming ring and is thus the ideal local laboratory for investigating the starburst--AGN connection in detail. We present integral-field observations of the central 1.3 kpc region in NGC 7469 obtained with the JWST Mid-InfraRed Instrument. Molecular and ionized gas distributions and kinematics at a resolution of ∼100 pc over the 4.9 - 7.6\mum region are examined to study gas dynamics influenced by the central AGN. The low-ionization [Fe II] \lambda5.34\mum and [Ar II] \lambda6.99\mum lines are bright on the nucleus and in the starburst ring, as opposed to H2 S(5) \lambda6.91\mum which is strongly peaked at the center and surrounding ISM. The high-ionization [Mg V] line is resolved and shows a broad, blueshifted component associated with the outflow. It has a nearly face-on geometry that is strongly peaked on the nucleus, where it reaches a maximum velocity of -650 km/s, and extends about 400 pc to the East. Regions of enhanced velocity dispersion in H2 and [Fe II] ∼180 pc from the AGN that also show high L(H2)/L(PAH) and L([Fe II])/L(Pf\alpha) ratios to the W and N of the nucleus pinpoint regions where the ionized outflow is depositing energy, via shocks, into the dense interstellar medium between the nucleus and the starburst ring. These resolved mid-infrared observations of the nuclear gas dynamics demonstrate the power of JWST and its high-sensitivity integral-field spectroscopic capability to resolve feedback processes around supermassive black holes in the dusty cores of nearby LIRGs.
Aaron S. Evans, David Frayer, Vassilis Charmandaris, Lee Armus, Hanae Inami, Jason Surace, Sean Linden, Baruch Soifer, Tanio Diaz-Santos, Kirsten Larson, Jeffrey Rich, Yiqing Song, Loreto Barcos-Munoz, Joseph Mazzarella, George Privon, Vivian U, Anne Medling, Torsten Boeker, Susanne Aalto, Kazushi Iwasawa, et al (15) James Webb Space Telescope (JWST) Mid-InfraRed Instrument (MIRI) images of the luminous infrared (IR) galaxy VV 114 are presented. This redshift ~ 0.020 merger has a western component (VV 114W) rich in optical star clusters and an eastern component (VV 114E) hosting a luminous mid-IR nucleus hidden at UV and optical wavelengths by dust lanes. With MIRI, the VV 114E nucleus resolves primarily into bright NE and SW cores separated by 630 pc. This nucleus comprises 45% of the 15um light of VV 114, with the NE and SW cores having IR luminosities, L_ IR (8-1000um) ~ 8+/-0.8x10^10 L_sun and ~ 5+/-0.5x10^10 L_sun, respectively, and IR densities, Sigma_IR >~ 2+/-0.2x10^13 L_sun / kpc^2 and >~ 7+/-0.7x10^12 L_sun / kpc^2, respectively -- in the range of Sigma_IR for the Orion star-forming core and the nuclei of Arp 220. The NE core, previously speculated to have an Active Galactic Nucleus (AGN), has starburst-like mid-IR colors. In contrast, the VV 114E SW has AGN-like colors. Approximately 40 star-forming knots with L_IR ~ 0.02-5x10^10 L_sun are identified, 25% of which have no optical counterpart. Finally, diffuse emission accounts for 40-60% of the mid-IR emission. Mostly notably, filamentary Poly-cyclic Aromatic Hydrocarbon (PAH) emission stochastically excited by UV and optical photons accounts for half of the 7.7um light of VV 114. This study illustrates the ability of JWST to detect obscured compact activity and distributed PAH emission in the most extreme starburst galaxies in the local Universe.
Hanae Inami, Jason Surace, Lee Armus, Aaron S. Evans, Kirsten L. Larson, Loreto Barcos-Munoz, Sabrina Stierwalt, Joseph M. Mazzarella, George C. Privon, Yiqing Song, Sean Linden, Christopher C. Hayward, Torsten Boker, Vivian U, Thomas Bohn, Vassilis Charmandaris, Tanio Diaz-Santos, Justin H. Howell, Thomas Lai, Anne M. Medling, et al (13) We have used the Mid-InfraRed Instrument (MIRI) on the James Webb Space Telescope (JWST) to obtain the first spatially resolved, mid-infrared (mid-IR) images of IIZw096, a merging luminous infrared galaxy (LIRG) at $z = 0.036$. Previous observations with the Spitzer Space Telescope suggested that the vast majority of the total IR luminosity (LIR) of the system originated from a small region outside of the two merging nuclei. New observations with JWST/MIRI now allow an accurate measurement of the location and luminosity density of the source that is responsible for the bulk of the IR emission. We estimate that 40-70% of the IR bolometric luminosity, or $3-5 \times 10^{11}\,{\rm{L_{\odot}}}$, arises from a source no larger than 175pc in radius, suggesting a luminosity density of at least $3-5 \times 10^{12} \, {\rm{L_{\odot} \, kpc^{-2}}}$. In addition, we detect 11 other star forming sources, five of which were previously unknown. The MIRI F1500W/F560W colors of most of these sources, including the source responsible for the bulk of the far-IR emission, are much redder than the nuclei of local LIRGs. These observations reveal the power of JWST to disentangle the complex regions at the hearts of merging, dusty galaxies.
We present the results of a stacking analysis performed on Spitzer/Infrared Spectrograph high-resolution mid-infrared spectra of luminous infrared galaxies (LIRGs) in the Great Observatories All-Sky LIRG Survey (GOALS). By binning on mid-infrared active galactic nucleus (AGN) fraction and stacking spectra, we detect bright emission lines [Ne II] and [Ne III], which trace star formation, and fainter emission lines [Ne V] and [O IV], which trace AGN activity, throughout the sample. We find the [Ne II] luminosity is fairly constant across all AGN fraction bins, while the [O IV] and [Ne V] luminosities increase by over an order of magnitude. Our measured average line ratios, [Ne V]/[Ne II] and [O IV]/[Ne II], at low AGN fraction are similar to H II galaxies while the line ratios at high AGN fraction are similar to LINERs and Seyferts. We decompose the [O IV] luminosity into star-formation and AGN components by fitting the [O IV] luminosity as a function of the [Ne II] luminosity and the mid-infrared AGN fraction. The [O IV] luminosity in LIRGs is dominated by star formation for mid-infrared AGN fractions $\lesssim0.3$. With the corrected [O IV] luminosity, we calculate black hole accretion rates ranging from $10^{-5}$ M$_{\odot}$/yr at low AGN fractions to 0.2 M$_{\odot}$/yr at the highest AGN fractions. We find that using the [O IV] luminosity, without correcting for star formation, can lead to an overestimate of the BHAR by up to a factor of 30 in starburst dominated LIRGs. Finally, we show the BHAR/SFR ratio increases by more than three orders of magnitude as a function of mid-infrared AGN fraction in LIRGs.
G. Martin, A. E. Bazkiaei, M. Spavone, E. Iodice, J. C. Mihos, M. Montes, J. A. Benavides, S. Brough, J. L. Carlin, C. A. Collins, P. A. Duc, F. A. Gómez, G. Galaz, H. M. Hernández-Toledo, R. A. Jackson, S. Kaviraj, J. H. Knapen, C. Martínez-Lombilla, S. McGee, D. O'Ryan, et al (33) Tidal features in the outskirts of galaxies yield unique information about their past interactions and are a key prediction of the hierarchical structure formation paradigm. The Vera C. Rubin Observatory is poised to deliver deep observations for potentially of millions of objects with visible tidal features, but the inference of galaxy interaction histories from such features is not straightforward. Utilising automated techniques and human visual classification in conjunction with realistic mock images produced using the NEWHORIZON cosmological simulation, we investigate the nature, frequency and visibility of tidal features and debris across a range of environments and stellar masses. In our simulated sample, around 80 per cent of the flux in the tidal features around Milky Way or greater mass galaxies is detected at the 10-year depth of the Legacy Survey of Space and Time (30-31 mag / sq. arcsec), falling to 60 per cent assuming a shallower final depth of 29.5 mag / sq. arcsec. The fraction of total flux found in tidal features increases towards higher masses, rising to 10 per cent for the most massive objects in our sample (M*~10^11.5 Msun). When observed at sufficient depth, such objects frequently exhibit many distinct tidal features with complex shapes. The interpretation and characterisation of such features varies significantly with image depth and object orientation, introducing significant biases in their classification. Assuming the data reduction pipeline is properly optimised, we expect the Rubin Observatory to be capable of recovering much of the flux found in the outskirts of Milky Way mass galaxies, even at intermediate redshifts (z<0.2).
Maren Cosens, Shelley A. Wright, Norman Murray, Lee Armus, Karin Sandstrom, Tuan Do, Kirsten Larson, Gregory Martinez, Sanchit Sabhlok, Andrey Vayner, James Wiley We present a survey of the central region of the nearest starburst galaxy, IC 10, using the W. M. Keck Observatory Keck Cosmic Web Imager (KCWI) at high spectral and spatial resolution. We map the central starburst of IC 10 to sample the kinematic and ionization properties of the individual star-forming regions. Using the low spectral resolution mode of KCWI we map the oxygen abundance and with the high spectral resolution mode we identify 46 individual H II regions. These H II regions have an average radius of 4.0 pc, star formation rate $\sim1.3\times10^{-4}$ M$_\odot$ yr$^{-1}$, and velocity dispersion $\sim$16 km s$^{-1}$. None of the H II regions appear to be virialized ($\rm \alpha_{vir}>>1$), and, on average, they show evidence of ongoing expansion. IC 10's H II regions are offset from the star forming region size-luminosity scaling relationships, as well as Larson's Law that relates size and velocity dispersion. We investigate the balance of inward and outward pressure, $\rm P_{in}$ and $\rm P_{out}$, finding $\rm P_{out}>P_{in}$ in 89% of H II regions, indicating feedback driven expansion even in these low mass H II regions. We find warm gas pressure ($\rm P_{gas}$) provides the dominant contribution to the outward pressure ($\rm P_{out}$). This counteracts the inward pressure which is dominated by turbulence in the surrounding gas rather than self-gravity. Five H II regions show evidence of outflows which are most likely supported by either stellar winds (2 regions) or champagne flows (3 regions). These observations provide new insights into the state of the star-forming regions in IC 10 and negative feedback from low mass clusters.
Yun Wang, Lee Armus, Andrew Benson, Emanuele Daddi, Andreas Faisst, Anthony Gonzalez, Casey Papovich, Zoran Ninkov, Massimo Robberto, Randall J. Rose, Thomas Rose, Claudia Scarlata, S. A. Stanford, Todd Veach, Zhongxu Zhai, Bradford Benson, L. E. Bleem, Michael W. Davis, George Helou, Lynne Hillenbrand ISCEA (Infrared Satellite for Cosmic Evolution Astrophysics) is a small astrophysics mission whose Science Goal is to discover how galaxies evolved in the cosmic web of dark matter at cosmic noon. Its Science Objective is to determine the history of star formation and its quenching in galaxies as a function of local density and stellar mass when the Universe was 3-5 Gyrs old (1.2<z<2.1). ISCEA is designed to test the Science Hypothesis that during the period of cosmic noon, at 1.7 < z < 2.1, environmental quenching is the dominant quenching mechanism for typical galaxies not only in clusters and groups, but also in the extended cosmic web surrounding these structures. ISCEA meets its Science Objective by making a 10% shot noise measurement of star formation rate down to 6 solar masses per year using H-alpha out to a radius > 10 Mpc in each of 50 protocluster (cluster and cosmic web) fields at 1.2 < z < 2.1. ISCEA measures the star formation quenching factor in those fields, and galaxy kinematics with a precision < 50 km/s to deduce the 3D spatial distribution in each field. ISCEA will transform our understanding of galaxy evolution at cosmic noon. ISCEA is a small satellite observatory with a 30cm equivalent diameter aperture telescope with a FoV of 0.32 deg^2, and a multi-object spectrograph with a digital micro-mirror device (DMD) as its programmable slit mask. ISCEA will obtain spectra of 1000 galaxies simultaneously at an effective resolving power of R=1000, with 2.8"x2.8" slits, over the NIR wavelength range of 1.1 to 2.0 microns, a regime not accessible from the ground without large gaps in coverage. ISCEA will achieve a pointing accuracy of <= 2" FWHM over 200s. ISCEA will be launched into a Low Earth Orbit, with a prime mission of 2.5 years. ISCEA's space-qualification of DMDs opens a new window for spectroscopy from space, enabling revolutionary advances in astrophysics.
Anne M. Medling, Lisa J. Kewley, Daniela Calzetti, George C. Privon, Kirsten Larson, Jeffrey A. Rich, Lee Armus, Mark G. Allen, Geoffrey V. Bicknell, Tanio Díaz-Santos, Timothy M. Heckman, Claus Leitherer, Claire E. Max, David S. N. Rupke, Ezequiel Treister, Hugo Messias, Alexander Y. Wagner We study the ionization and excitation structure of the interstellar medium in the late-stage gas-rich galaxy merger NGC 6240 using a suite of emission line maps at $\sim$25 pc resolution from the Hubble Space Telescope, Keck NIRC2 with Adaptive Optics, and ALMA. NGC 6240 hosts a superwind driven by intense star formation and/or one or both of two active nuclei; the outflows produce bubbles and filaments seen in shock tracers from warm molecular gas (H$_2$ 2.12$\mu$m) to optical ionized gas ([O III], [N II], [S II], [O I]) and hot plasma (Fe XXV). In the most distinct bubble, we see a clear shock front traced by high [O III]/H$\beta$ and [O III]/[O I]. Cool molecular gas (CO(2-1)) is only present near the base of the bubble, towards the nuclei launching the outflow. We interpret the lack of molecular gas outside the bubble to mean that the shock front is not responsible for dissociating molecular gas, and conclude that the molecular clouds are partly shielded and either entrained briefly in the outflow, or left undisturbed while the hot wind flows around them. Elsewhere in the galaxy, shock-excited H$_2$ extends at least $\sim$4 kpc from the nuclei, tracing molecular gas even warmer than that between the nuclei, where the two galaxies' interstellar media are colliding. A ridgeline of high [O III]/H$\beta$ emission along the eastern arm aligns with the south nucleus' stellar disk minor axis; optical integral field spectroscopy from WiFeS suggests this highly ionized gas is centered at systemic velocity and likely photoionized by direct line-of-sight to the south AGN.
Sean Linden, Aaron Evans, Kirsten Larson, George Privon, Lee Armus, Jeff Rich, Tanio Diaz-Santos, Eric Murphy, Yiqing Song, Loreto Barcos-Munoz, Justin Howell, Vassilis Charmandaris, Hanae Inami, Vivian U, Jason Surace, Joseph Mazzarella, Daniela Calzetti We present the results of a Hubble Space Telescope WFC3 near-UV and ACS/WFC optical study into the star cluster populations of 10 luminous and ultra-luminous infrared galaxies (U/LIRGs) in the Great Observatories All-Sky LIRG Survey (GOALS). Through integrated broadband photometry we have derived ages, masses, and extinctions for a total of 1027 star clusters in galaxies with $d_{L} <$ 110 Mpc in order to avoid issues related to cluster blending. The measured cluster age distribution slope of $dN/d\tau \propto \tau^{-0.5 +/- 0.2}$ is steeper than what has been observed in lower-luminosity star-forming galaxies. Further, differences in the slope of the observed cluster age distribution between inner- ($dN/d\tau \propto \tau^{-1.07 +/- 0.12}$) and outer-disk ($dN/d\tau \propto \tau^{-0.37 +/- 0.09}$) star clusters provides evidence of mass-dependent cluster destruction in the central regions of LIRGs driven primarily by the combined effect of strong tidal shocks and encounters with massive GMCs. Excluding the nuclear ring surrounding the Seyfert 1 nucleus in NGC 7469, the derived cluster mass function (CMF: $dN/dM \propto M^{\alpha}$) has marginal evidence for a truncation in the power-law (PL) at $M_{t} \sim 2$x$10^{6} M_{\odot}$ for our three most cluster-rich galaxies, which are all classified as early-stage mergers. Finally, we find evidence of a flattening of the CMF slope of $dN/dM \propto M^{-1.42 \pm 0.1}$ for clusters in late-stage mergers relative to early-stage ($\alpha = -1.65 \pm 0.02$), which we attribute to an increase in the formation of massive clusters over the course of the interaction.
We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of six radio-loud quasar host galaxies at $z=1.4-2.3$. We combine the kpc-scale resolution ALMA observations with high spatial-resolution adaptive optics integral field spectrograph data of the ionized gas. We detect molecular gas emission in five quasar host galaxies and resolve the molecular interstellar medium using the CO (3-2) or CO (4-3) rotational transitions. Clumpy molecular outflows are detected in four quasar host galaxies and in a merger system 21 kpc away from one quasar. Between the ionized and cold-molecular gas phases, the majority of the outflowing mass is in a molecular phase, while for three out of four detected multi-phase gas outflows, the majority of the kinetic luminosity and momentum flux is in the ionized phase. Combining the energetics of the multi-phase outflows, we find that their driving mechanism is consistent with energy-conserving shocks produced by the impact of the quasar jets with the gas in the galaxy. By assessing the molecular gas mass to the dynamics of the outflows, we estimate a molecular gas depletion time scale of a few Myr. The gas outflow rates exceed the star formation rates, suggesting that quasar feedback is a major mechanism of gas depletion at the present time. The coupling efficiency between the kinetic luminosity of the outflows and the bolometric luminosity of the quasar of 0.1-1% is consistent with theoretical predictions. Studying multi-phase gas outflows at high redshift is important for quantifying the impact of negative feedback in shaping the evolution of massive galaxies.
Jason Glenn, Charles M. Bradford, Erik Rosolowsky, Rashied Amini, Katherine Alatalo, Lee Armus, Andrew J. Benson, Tzu-Ching Chang, Jeremy Darling, Peter K. Day, Jeanette Domber, Duncan Farrah, Brandon Hensley, Sarah Lipscy, Bradley Moore, Seb Oliver, Joanna Perido, David Redding, Michael Rodgers, Raphael Shirley, et al (4) The Galaxy Evolution Probe (GEP) is a concept for a mid- and far-infrared space observatory to measure key properties of large samples of galaxies with large and unbiased surveys. GEP will attempt to achieve zodiacal light and Galactic dust emission photon background-limited observations by utilizing a 6 Kelvin, 2.0 meter primary mirror and sensitive arrays of kinetic inductance detectors. It will have two instrument modules: a 10 - 400 micron hyperspectral imager with spectral resolution R = 8 (GEP-I) and a 24 - 193 micron, R = 200 grating spectrometer (GEP-S). GEP-I surveys will identify star-forming galaxies via their thermal dust emission and simultaneously measure redshifts using polycyclic aromatic hydrocarbon emission lines. Galaxy luminosities derived from star formation and nuclear supermassive black hole accretion will be measured for each source, enabling the cosmic star formation history to be measured to much greater precision than previously possible. Using optically thin far-infrared fine-structure lines, surveys with GEP-S will measure the growth of metallicity in the hearts of galaxies over cosmic time and extraplanar gas will be mapped in spiral galaxies in the local universe to investigate feedback processes. The science case and mission architecture designed to meet the science requirements are described, and the kinetic inductance detector and readout electronics state of the art and needed developments are described. This paper supersedes the GEP concept study report cited in it by providing new content, including: a summary of recent mid-infrared KID development, a discussion of microlens array fabrication for mid-infrared KIDs, and additional context for galaxy surveys. The reader interested in more technical details may want to consult the concept study report.
C. Ricci, G. C. Privon, R. W. Pfeifle, L. Armus, K. Iwasawa, N. Torres-Alba, S. Satyapal, F. E. Bauer, E. Treister, L. C. Ho, S. Aalto, P. Arevalo, L. Barcos-Munoz, V. Charmandaris, T. Diaz-Santos, A. S. Evans, T. Gao, H. Inami, M. J. Koss, G. Lansbury, et al (8) The merger of two or more galaxies can enhance the inflow of material from galactic scales into the close environments of Active Galactic Nuclei (AGN), obscuring and feeding the supermassive black hole (SMBH). Both recent simulations and observations of AGN in mergers have confirmed that mergers are related to strong nuclear obscuration. However, it is still unclear how AGN obscuration evolves in the last phases of the merger process. We study a sample of 60 Luminous and Ultra-luminous IR galaxies (U/LIRGs) from the GOALS sample observed by NuSTAR. We find that the fraction of AGN that are Compton-thick (CT; $N_{\rm H}\geq 10^{24}\rm\,cm^{-2}$) peaks at $74_{-19}^{+14}\%$ at a late merger stage, prior to coalescence, when the nuclei have projected separations of $d_{\rm sep}\sim 0.4-6$ kpc. A similar peak is also observed in the median $N_{\rm H}$ [$(1.6\pm0.5)\times10^{24}\rm\,cm^{-2}$]. The vast majority ($85^{+7}_{-9}\%$) of the AGN in the final merger stages ($d_{\rm sep}\lesssim 10$ kpc) are heavily obscured ($N_{\rm H}\geq 10^{23}\rm\,cm^{-2}$), and the median $N_{\rm H}$ of the accreting SMBHs in our sample is systematically higher than that of local hard X-ray selected AGN, regardless of the merger stage. This implies that these objects have very obscured nuclear environments, with the $N_{\rm H}\geq 10^{23}\rm\,cm^{-2}$ gas almost completely covering the AGN in late mergers. CT AGN tend to have systematically higher absorption-corrected X-ray luminosities than less obscured sources. This could either be due to an evolutionary effect, with more obscured sources accreting more rapidly because they have more gas available in their surroundings, or to a selection bias. The latter scenario would imply that we are still missing a large fraction of heavily obscured, lower luminosity ($L_{2-10}\lesssim 10^{43}\rm\,erg\,s^{-1}$) AGN in U/LIRGs.
Y. Song, S. T. Linden, A. S. Evans, L. Barcos-Muñoz, G. C. Privon, I. Yoon, E. J. Murphy, K. L. Larson, T. Díaz-Santos, L. Armus, Joseph M. Mazzarella, J. Howell, H. Inami, N. Torres-Albà, V. U, V. Charmandaris, E. Momjian, J. McKinney, D. Kunneriath Nuclear rings are excellent laboratories for studying intense star formation. We present results from a study of nuclear star-forming rings in five nearby normal galaxies from the Star Formation in Radio Survey (SFRS) and four local LIRGs from the Great Observatories All-sky LIRG Survey (GOALS) at sub-kpc resolutions using VLA high-frequency radio continuum observations. We find that nuclear ring star formation (NRSF) contributes 49 - 60\% of the total star formation of the LIRGs, compared to 7 - 40\% for the normal galaxies. We characterize a total of 58 individual star-forming regions in these rings, and find that with measured sizes of 10 - 200 pc, NRSF regions in the LIRGs have SFR and $\Sigma_\mathrm{SFR}$ up to 1.7 M$_\odot$yr$^{-1}$ and 402 M$_\odot$yr$^{-1}$kpc$^{-2}$, respectively, which are about 10 times higher than NRSF regions in the normal galaxies with similar sizes, and comparable to lensed high-$z$ star-forming regions. At $\sim 100 - 300$ pc scales, we estimate low contributions ($< 50\%$) of thermal free-free emission to total radio continuum emission at 33 GHz in the NRSF regions in the LIRGs, but large variations possibly exist at smaller physical scales. Finally, using archival sub-kpc resolution CO (J=1-0) data of nuclear rings in the normal galaxies and NGC 7469 (LIRG), we find a large scatter in gas depletion times at similar molecular gas surface densities, which tentatively points to a multi-modal star formation relation on sub-kpc scales.
Ori D. Fox, Harish Khandrika, David Rubin, Chadwick Casper, Gary Z. Li, Tamas Szalai, Lee Armus, Alexei V. Filippenko, Michael F. Skrutskie, Lou Strolger, Schuyler D. Van Dyk Supernova (SN) rates serve as an important probe of star-formation models and initial mass functions. Near-infrared seeing-limited ground-based surveys typically discover a factor of 3-10 fewer SNe than predicted from far-infrared (FIR) luminosities owing to sensitivity limitations arising from both a variable point-spread function (PSF) and high dust extinction in the nuclear regions of star-forming galaxies. This inconsistency has potential implications for our understanding of star-formation rates and massive-star evolution, particularly at higher redshifts, where star-forming galaxies are more common. To resolve this inconsistency, a successful SN survey in the local universe must be conducted at longer wavelengths and with a space-based telescope, which has a stable PSF to reduce the necessity for any subtraction algorithms and thus residuals. Here we report on a two-year Spitzer/IRAC 3.6 um survey for dust-extinguished SNe in the nuclear regions of forty luminous infrared galaxies (LIRGs) within 200 Mpc. The asymmetric Spitzer PSF results in worse than expected subtraction residuals when implementing standard template subtraction. Forward-modeling techniques improve our sensitivity by ~1.5 magnitudes. We report the detection of 9 SNe, five of which were not discovered by optical surveys. After adjusting our predicted rates to account for the sensitivity of our survey, we find that the number of detections is consistent with the models. While this search is nonetheless hampered by a difficult-to-model PSF and the relatively poor resolution of Spitzer, it will benefit from future missions, such as Roman Space Telescope and JWST, with higher resolution and more symmetric PSFs.
We present observations of ionized gas outflows in eleven z$ =1.39-2.59$ radio-loud quasar host galaxies. Data was taken with the integral field spectrograph (IFS) OSIRIS and the adaptive optics system at the W.M. Keck Observatory targeting nebular emission lines (H$\beta$, [OIII], H$\alpha$, [NII] and [SII]) redshifted into the near-infrared (1-2.4 \micron). Outflows with velocities of 500 - 1700 km\u2009s$^{-1}$ are detected in 10 systems on scales ranging from $<1$ kpc to 10 kpc with outflow rates from 8-2400 M$_\odot$yr$^{-1}$. For five sources, the outflow momentum rates are 4-80 times $L_{AGN}$/c, consistent with outflows being driven by an energy conserving shock. The five other outflows are either driven by radiation pressure or an isothermal shock. The outflows are the dominant source of gas depletion, and we find no evidence for star formation along the outflow paths. For eight objects, the outflow paths are consistent with the orientation of the jets. Yet, given the calculated pressures, we find no evidence of the jets currently doing work on these galactic-scale ionized outflows. We find that galactic-scale feedback occurs well before galaxies establish a substantial fraction of their stellar mass, as expected from local scaling relationships.
L. Armus, S. T. Megeath, L. Corrales, M. Marengo, A. Kirkpatrick, J. D. Smith, M. Meyer, S. Gezari, R. P. Kraft, S. McCandliss, S. Tuttle, M. Elvis, M. Bentz, B. Binder, F. Civano, D. Dragomir, C. Espaillat, S. Finkelstein, D. B. Fox, M. Greenhouse, et al (13) NASA's Great Observatories have opened up the electromagnetic spectrum from space, providing sustained access to wavelengths not accessible from the ground. Together, Hubble, Compton, Chandra, and Spitzer have provided the scientific community with an agile and powerful suite of telescopes with which to attack broad scientific questions, and react to a rapidly changing scientific landscape. As the existing Great Observatories age, or are decommissioned, community access to these wavelengths will diminish, with an accompanying loss of scientific capability. This report, commissioned by the NASA Cosmic Origins, Physics of the Cosmos and Exoplanet Exploration Program Analysis Groups (PAGs), analyzes the importance of multi-wavelength observations from space during the epoch of the Great Observatories, providing examples that span a broad range of astrophysical investigations.
I. Heywood, E. J. Murphy, E. F. Jiménez-Andrade, L. Armus, W. D. Cotton, C. DeCoursey, M. Dickinson, T. J. W. Lazio, E. Momjian, K. Penner, I. Smail, O. M. Smirnov The Frontier Fields project is an observational campaign targeting six galaxy clusters, with the intention of using the magnification provided by gravitational lensing to study galaxies that are extremely faint or distant. We used the Karl G. Jansky Very Large Array (VLA) at 3 and 6 GHz to observe three Frontier Fields: MACSJ0416.1$-$2403 ($z$ = 0.396), MACSJ0717.5+3745 ($z$ = 0.545), and MACSJ1149.5+2223 ($z$ = 0.543). The images reach noise levels of $\sim$1 $\mu$Jy beam$^{-1}$ with sub-arcsecond resolution ($\sim$2.5 kpc at $z$ = 3), providing a high-resolution view of high-$z$ star-forming galaxies that is unbiased by dust obscuration. We generate dual-frequency continuum images at two different resolutions per band, per cluster, and derive catalogs totalling 1966 compact radio sources. Components within the areas of Hubble Space Telescope and Subaru observations are cross-matched, providing host galaxy identifications for 1296 of them. We detect 13 moderately-lensed (2.1 $<$ $\mu$ $<$ 6.5) sources, one of which has a demagnified peak brightness of 0.9 $\mu$Jy beam$^{-1}$, making it a candidate for the faintest radio source ever detected. There are 66 radio sources exhibiting complex morphologies, and 58 of these have host galaxy identifications. We reveal that MACSJ1149.5+2223 is not a cluster with a double relic, as the western candidate relic is resolved as a double-lobed radio galaxy associated with a foreground elliptical at $z$ = 0.24. The VLA Frontier Fields project is a public legacy survey. The image and catalog products from this work are freely available.
To investigate the growth history of galaxies, we measure the rest-frame radio, ultraviolet (UV), and optical sizes of 98 radio-selected, star-forming galaxies (SFGs) distributed over $0.3 \lesssim z \lesssim 3$ and median stellar mass of $\log(M_\star/ \rm M_\odot)\approx10.4$. We compare the size of galaxy stellar disks, traced by rest-frame optical emission, relative to the overall extent of star formation activity that is traced by radio continuum emission. Galaxies in our sample are identified in three Hubble Frontier Fields: MACSJ0416.1$-$2403, MACSJ0717.5+3745, and MACSJ1149.5+2223. Radio continuum sizes are derived from 3 GHz and 6 GHz radio images ($\lesssim 0$''$.6$ resolution, $\approx0.9\, \rm \mu Jy\, beam^{-1}$ noise level) from the Karl G. Jansky Very Large Array. Rest-frame UV and optical sizes are derived using observations from the Hubble Space Telescope and the ACS and WFC3 instruments. We find no clear dependence between the 3 GHz radio size and stellar mass of SFGs, which contrasts with the positive correlation between the UV/optical size and stellar mass of galaxies. Focusing on SFGs with $\log(M_\star/\rm M_\odot)>10$, we find that the radio/UV/optical emission tends to be more compact in galaxies with high star-formation rates ($\rm SFR\gtrsim 100\,M_\odot\,yr^{-1}$), suggesting that a central, compact starburst (and/or an Active Galactic Nucleus) resides in the most luminous galaxies of our sample. We also find that the physical radio/UV/optical size of radio-selected SFGs with $\log(M_\star/\rm M_\odot)>10$ increases by a factor of $1.5-2$ from $z\approx 3$ to $z\approx0.3$, yet the radio emission remains two-to-three times more compact than that from the UV/optical. These findings indicate that these massive, radio-selected SFGs at $0.3 \lesssim z \lesssim 3$ tend to harbor centrally enhanced star formation activity relative to their outer-disks.
Luigi Spinoglio, Sabrina Mordini, Juan Antonio Fernandez-Ontiveros, Almudena Alonso-Herrero, Lee Armus, Laura Bisigello, Francesco Calura, Francisco J. Carrera, Asantha Cooray, Helmut Dannerbauer, Roberto Decarli, Eiichi Egami, David Elbaz, Alberto Franceschini, Eduardo Gonzalez Alfonso, Luca Graziani, Carlotta Gruppioni, Evanthia Hatziminaoglou, Hidehiro Kaneda, Kotaro Kohno, et al (14) We use the SPace Infrared telescope for Cosmology and Astrophysics (SPICA) project as a template to demonstrate how deep spectrophotometric surveys covering large cosmological volumes over extended fields (1-15 square degrees) with a mid-IR imaging spectrometer (17-36 micron) in conjunction with deep 70 micron photometry with a far-IR camera, at wavelengths which are not affected by dust extinction can answer the most crucial questions in current galaxy evolution studies. A SPICA-like mission will be able for the first time to provide an unobscured three dimensional (3-D, i.e. x, y and redshift z) view of galaxy evolution back to an age of the Universe of less than ~2 Gyrs, in the mid-IR rest-frame. This survey strategy will produce a full census of the Star formation Rate (SFR) in the Universe, using Polycyclic Aromatic Hydrocarbons (PAH) bands and fine-structure ionic lines, reaching the characteristic knee of the galaxy luminosity function, where the bulk of the population is distributed, at any redshift up to z ~3.5. Deep follow-up pointed spectroscopic observations with grating spectrometers onboard the satellite, across the full IR spectral range (17-210 micron), would simultaneously measure Black Hole Accretion Rate (BHAR), from high-ionization fine-structure lines, and SFR, from PAH and low- to mid-ionization lines in thousands of galaxies from solar to low metallicities, down to the knee of their luminosity functions. The analysis of the resulting atlas of IR spectra will reveal the physical processes at play in evolving galaxies across cosmic time, especially its heavily dust-embedded phase during the activity peak at the cosmic noon (z ~1-3), through IR emission lines and features that are insensitive to the dust obscuration.
We present detailed observations of photoionization conditions and galaxy kinematics in eleven z$=1.39-2.59$ radio-loud quasar host galaxies. Data was taken with OSIRIS integral field spectrograph (IFS) and the adaptive optics system at the W.M. Keck Observatory that targeted nebular emission lines (H$\beta$,[OIII],H$\alpha$,[NII]) redshifted into the near-infrared (1-2.4 \micron). We detect extended ionized emission on scales ranging from 1-30 kpc photoionized by stars, shocks, and active galactic nuclei (AGN). Spatially resolved emission-line ratios indicate that our systems reside off the star formation and AGN-mixing sequence on the Baldwin, Phillips $\&$ Terlevich (BPT) diagram at low redshift. The dominant cause of the difference between line ratios of low redshift galaxies and our sample is due to lower gas-phase metallicities, which are 2-5$\times$ less compared to galaxies with AGN in the nearby Universe. Using gas velocity dispersion as a proxy to stellar velocity dispersion and dynamical mass measurement through inclined disk modeling we find that the quasar host galaxies are under-massive relative to their central supermassive black hole (SMBH) mass, with all systems residing off the local scaling ($M_{\bullet}-\sigma~$,$M_{\bullet}-M_{*}~$) relationship. These quasar host galaxies require substantial growth, up to an order of magnitude in stellar mass, to grow into present-day massive elliptical galaxies. Combining these results with part I of our sample paper (Vayner et al. 2021) we find evidence for winds capable of causing feedback before the AGN host galaxies land on the local scaling relation between black hole and galaxy stellar mass, and before the enrichment of the ISM to a level observed in local galaxies with AGN.
Star formation in galaxies is regulated by the heating and cooling in the interstellar medium. In particular, the processing of molecular gas into stars will depend strongly on the ratio of gas heating to gas cooling in the neutral gas around sites of recent star-formation. In this work, we combine mid-infrared (mid-IR) observations of Polycyclic Aromatic Hydrocarbons (PAHs), the dominant heating mechanism of gas in the interstellar medium (ISM), with [C II], [O I], and [Si II] fine-structure emission, the strongest cooling channels in dense, neutral gas. The ratio of IR cooling line emission to PAH emission measures the photoelectric efficiency, a property of the ISM which dictates how much energy carried by ultraviolet photons gets transferred into the gas. We find that star-forming, IR luminous galaxies in the Great Observatories All-Sky LIRG Survey (GOALS) with high IR surface densities have low photoelectric efficiencies. These systems also have, on average, higher ratios of radiation field strength to gas densities, and larger average dust grain size distributions. The data support a scenario in which the most compact galaxies have more young star-forming regions per unit area, which exhibit less efficient gas heating. These conditions may be more common at high-z, and may help explain the higher star-formation rates at cosmic noon. We make predictions on how this can be investigated with JWST.
Luke Finnerty, Kirsten Larson, B. T. Soifer, Lee Armus, Keith Matthews, Hyunsung D. Jun, Dae-Sik Moon, Jason Melbourne, Percy Gomez, Chao-Wei Tsai, Tanio Diaz-Santos, Peter Eisenhardt, Michael Cushing We present rest-frame optical spectroscopic observations of 24 Hot Dust-Obscured Galaxies (Hot DOGs) at redshifts 1.7-4.6 with KECK/NIRES. Our targets are selected based on their extreme red colors to be the highest luminosity sources from the WISE infrared survey. In 20 sources with well-detected emission we fit the key [O III], H$\beta$, H$\alpha$, [N II], and [S II] diagnostic lines to constrain physical conditions. Of the 17 targets with a clear detection of the [O III]$\rm \lambda$5007A emission line, 15 display broad blueshifted and asymmetric line profiles, with widths ranging from 1000 to 8000 $\rm km\ s^{-1}$ and blueshifts up to 3000 $\rm km\ s^{-1}$. These kinematics provide strong evidence for the presence of massive ionized outflows of up to $8000\ \rm M_\odot\ yr^{-1}$, with a median of $150\ \rm M_\odot\ yr^{-1}$. As many as eight sources show optical emission line ratios consistent with vigorous star formation. Balmer line star-formation rates, uncorrected for reddening, range from 30--1300 $\rm M_\odot\ yr^{-1}$, with a median of $50\ \rm M_\odot\ yr^{-1}$. Estimates of the SFR from SED fitting of mid and far-infrared photometry suggest significantly higher values. We estimate the central black hole masses to be of order $10^{8-10}\rm\ M_\odot$, assuming the present-day $\rm M_{BH}-\sigma_*$ relation. The bolometric luminosities and the estimated masses of the central black holes of these galaxies suggest that many of the AGN-dominated Hot DOGs are accreting at or above their Eddington limit. The combination of ongoing star formation, massive outflows, and high Eddington ratios suggest Hot DOGs are a transitional phase in galaxy evolution.
Lena Murchikova, Eric J. Murphy, Dariusz C. Lis, Lee Armus, Selma de Mink, Kartik Sheth, Nadia Zakamska, Frank Tramper, Angela Bongiorno, Martin Elvis, Lisa Kewley, Hugues Sana The extreme ultraviolet (EUV) spectra of distant star-forming regions cannot be probed directly using either ground- or space-based telescopes due to the high cross-section for interaction of EUV photons with the interstellar medium. This makes EUV spectra poorly constrained. The mm/submm recombination lines of H and He, which can be observed from the ground, can serve as a reliable probe of the EUV. Here we present a study based on ALMA observations of three Galactic ultra-compact HII regions and the starburst region Sgr B2(M), in which we reconstruct the key parameters of the EUV spectra using mm recombination lines of HI, HeI and HeII. We find that in all cases the EUV spectra between 13.6 and 54.4 eV have similar frequency dependence: L_\nu~ \nu^-4.5 +/- 0.4. We compare the inferred values of the EUV spectral slopes with the values expected for a purely single stellar evolution model (Starburst99) and the Binary Population and Spectral Synthesis code (BPASS). We find that the observed spectral slope differs from the model predictions. This may imply that the fraction of interacting binaries in HII regions is substantially lower than assumed in BPASS. The technique demonstrated here allows one to deduce the EUV spectra of star forming regions providing critical insight into photon production rates at \lambda < 912 A and can serve as calibration to starburst synthesis models, improving our understanding of star formation in distant universe and the properties of ionizing flux during reionization.
G. C. Privon, C. Ricci, S. Aalto, S. Viti, L. Armus, T. Díaz-Santos, E. González-Alfonso, K. Iwasawa, D. L. Jeff, E. Treister, F. Bauer, A. S. Evans, P. Garg, R. Herrero-Illana, J. M. Mazzarella, K. Larson, L. Blecha, L. Barcos-Muñoz, V. Charmandaris, S. Stierwalt, et al (1) Enhanced emission from the dense gas tracer HCN (relative to HCO$^+ $) has been proposed as a signature of active galactic nuclei (AGN). In a previous single-dish millimeter line survey we identified galaxies with HCN/HCO$ ^+ $ (1-0) intensity ratios consistent with those of many AGN but whose mid-infrared spectral diagnostics are consistent with little to no ( $\lesssim15\% $) contribution of an AGN to the bolometric luminosity. To search for putative heavily obscured AGN, we present and analyze \nustar hard X-ray (3-79 keV) observations of four such galaxies from the Great Observatories All-sky LIRG Survey. We find no X-ray evidence for AGN in three of the systems and place strong upper limits on the energetic contribution of any heavily obscured ($N_{\rm H}>10^{24}$ cm$^{-2}$) AGN to their bolometric luminosity. The X-ray flux upper limits are presently an order of magnitude below what XDR-driven chemistry model predict are necessary to drive HCN enhancements. In a fourth system we find a hard X-ray excess consistent with the presence of an AGN, but contributing only $\sim3\%$ of the bolometric luminosity. It is also unclear if the AGN is spatially associated with the HCN enhancement. We further explore the relationship between HCN/HCO$^+$ (for several $\mathrm{J}_\mathrm{upper}$ levels) and $L_\mathrm{AGN}/L_\mathrm{IR}$ for a larger sample of systems in the literature. We find no evidence for correlations between the line ratios and the AGN fraction derived from X-rays, indicating that HCN/HCO$^+$ intensity ratios are not driven by the energetic dominance of AGN, nor are they reliable indicators of whether SMBH accretion is ongoing.
Star formation depends critically on cooling mechanisms in the interstellar medium (ISM); however, thermal properties of gas in galaxies at the peak epoch of star formation (z ~ 2) remain poorly understood. A limiting factor in understanding the multiphase ISM is the lack of multiple tracers detected in the same galaxies, such as Polycyclic Aromatic Hydrocarbon (PAH) emission, a tracer of a critical photoelectric heating mechanism in interstellar gas, and [C II] 158\mum fine-structure emission, a principal coolant. We present ALMA Band 9 observations targeting [C II] in six z ~ 2 star-forming galaxies with strong Spitzer IRS detections of PAH emission. All six galaxies are detected in dust continuum and marginally resolved. We compare the properties of PAH and [C II] emission, and constrain their relationship as a function of total infrared luminosity (LIR) and IR surface density. [C II] emission is detected in one galaxy at high signal-to-noise (34\sigma), and we place a secure upper limit on a second source. The rest of our sample are not detected in [C II] likely due to redshift uncertainties and narrow ALMA bandpass windows. Our results are consistent with the deficit in [C II]/LIR and PAH/LIR observed in the literature. However, the ratio of [C II] to PAH emission at z ~ 2 is possibly much lower than what is observed in nearby dusty star-forming galaxies. This could be the result of enhanced cooling via [O I] at high-z, hotter gas and dust temperatures, and/or a reduction in the photoelectric efficiency, in which the coupling between interstellar radiation and gas heating is diminished.
Alison F. Crocker, Eric Pellegrini, J.-D. T. Smith, Bruce T. Draine, Christine D. Wilson, Mark Wolfire, Lee Armus, Elias Brinks, Daniel A. Dale, Brent Groves, Rodrigo Herrera-Camus, Leslie K. Hunt, Robert C. Kennicutt, Eric J. Murphy, Karin Sandstrom, Eva Schinnerer, Dimitra Rigopoulou, Erik Rosolowsky, Paul van der Werf We present resolved [CI] line intensities of 18 nearby galaxies observed with the SPIRE FTS spectrometer on the Herschel Space Observatory. We use these data along with resolved CO line intensities from $J_\mathrm{up} = 1$ to 7 to interpret what phase of the interstellar medium the [CI] lines trace within typical local galaxies. A tight, linear relation is found between the intensities of the CO(4-3) and [CI](2-1) lines; we hypothesize this is due to the similar upper level temperature of these two lines. We modeled the [CI] and CO line emission using large velocity gradient models combined with an empirical template. According to this modeling, the [CI](1-0) line is clearly dominated by the low-excitation component. We determine [CI] to molecular mass conversion factors for both the [CI](1-0) and [CI](2-1) lines, with mean values of $\alpha_{\mathrm{[CI](1-0)}} = 7.3$ M$_{\mathrm{sun}}$ K$^{-1}$ km$^{-1}$ s pc$^{-2}$ and $\alpha_{\mathrm{[CI](2-1)}} = 34 $ M$_{\mathrm{sun}}$ K$^{-1}$ km$^{-1}$ s pc$^{-2}$ with logarithmic root-mean-square spreads of 0.20 and 0.32 dex, respectively. The similar spread of $\alpha_{\mathrm{[CI](1-0)}}$ to $\alpha_{\mathrm{CO}}$ (derived using the CO(2-1) line) suggests that [CI](1-0) may be just as good a tracer of cold molecular gas as CO(2-1) in galaxies of this type. On the other hand, the wider spread of $\alpha_{\mathrm{[CI](2-1)}}$ and the tight relation found between [CI](2-1) and CO(4-3) suggest that much of the [CI](2-1) emission may originate in warmer molecular gas.
E. Treister, H. Messias, G. C. Privon, N. Nagar, A. M. Medling, V. U., F. E. Bauer, C. Cicone, L. Barcos Munoz, A. S. Evans, F. Muller-Sanchez, J. M. Comerford, L. Armus, C. Chang, M. Koss, G. Venturi, K. Schawinski, C. Casey, C. M. Urry, D. B. Sanders, et al (2) We present the highest resolution --- 15 pc (0.03'') --- ALMA $^{12}$CO(2-1) line emission and 1.3mm continuum maps, tracers of the molecular gas and dust, respectively, in the nearby merging galaxy system NGC 6240, that hosts two supermassive black holes growing simultaneously. These observations provide an excellent spatial match to existing Hubble optical and near-infrared observations of this system. A significant molecular gas mass, $\sim$9$\times$10$^9$M$_\odot$, is located in between the two nuclei, forming a clumpy stream kinematically dominated by turbulence, rather than a smooth rotating disk as previously assumed from lower resolution data. Evidence for rotation is seen in the gas surrounding the southern nucleus, but not in the northern one. Dynamical shells can be seen, likely associated with nuclear supernovae remnants. We further detect the presence of significant high velocity outflows, some of them reaching velocities $>$500 km/s, affecting a significant fraction, $\sim$11\% of the molecular gas in the nuclear region. Inside the spheres of influence of the northern and southern supermassive black holes we find molecular masses of 7.4$\times$10$^8$M$_\odot$ and 3.3$\times$10$^9$M$_\odot$, respectively. We are thus directly imaging the reservoir of gas that can accrete onto each supermassive black hole. These new ALMA maps highlight the critical need for high resolution observations of molecular gas in order to understand the feeding of supermassive black holes and its connection to galaxy evolution in the context of a major galaxy merger.
M. Meixner, A. Cooray, D. Leisawitz, J. Staguhn, L. Armus, C. Battersby, J. Bauer, E. Bergin, C.M. Bradford, K. Ennico-Smith, J. Fortney, T. Kataria, G. Melnick, S. Milam, D. Narayanan, D. Padgett, K. Pontoppidan, A. Pope, T. Roellig, K. Sandstrom, et al (49) The Origins Space Telescope (Origins) traces our cosmic history, from the formation of the first galaxies and the rise of metals to the development of habitable worlds and present-day life. Origins does this through exquisite sensitivity to infrared radiation from ions, atoms, molecules, dust, water vapor and ice, and observations of extra-solar planetary atmospheres, protoplanetary disks, and large-area extragalactic fields. Origins operates in the wavelength range 2.8 to 588 microns and is 1000 times more sensitive than its predecessors due to its large, cold (4.5 K) telescope and advanced instruments. Origins was one of four large missions studied by the community with support from NASA and industry in preparation for the 2020 Decadal Survey in Astrophysics. This is the final study report.
G. Aniano, B.T. Draine, L.K. Hunt, K. Sandstrom, D. Calzetti, R.C. Kennicutt, D.A. Dale, M. Galametz, K.D. Gordon, A.K. Leroy, J.-D.T. Smith, H. Roussel, M. Sauvage, F. Walter, L. Armus, A.D. Bolatto, M. Boquien, A. Crocker, I. De Looze, J. Donovan Meyer, et al (13) Dust and starlight are modeled for the KINGFISH project galaxies. With data from 3.6 micron to 500 micron, models are strongly constrained. For each pixel in each galaxy we estimate (1) dust surface density; (2) q_PAH, the dust mass fraction in PAHs; (3) distribution of starlight intensities heating the dust; (4) luminosity emitted by the dust; and (5) dust luminosity from regions with high starlight intensity. The models successfully reproduce both global and resolved spectral energy distributions. We provide well-resolved maps for the dust properties. As in previous studies, we find q_PAH to be an increasing function of metallicity, above a threshold Z/Z_sol approx 0.15. Dust masses are obtained by summing the dust mass over the map pixels; these "resolved" dust masses are consistent with the masses inferred from model fits to the global photometry. The global dust-to-gas ratios obtained from this study correlate with galaxy metallicities. Systems with Z/Z_sol > 0.5 have most of their refractory elements locked up in dust, whereas when Z/Z_sol < 0.3 most of these elements tend to remain in the gas phase. Within galaxies, we find that q_PAH is suppressed in regions with unusually warm dust with nu L_nu(70 um) > 0.4L_dust. With knowledge of one long-wavelength flux density ratio (e.g., f_160/f_500), the minimum starlight intensity heating the dust (U_min) can be estimated to within ~50%. For the adopted dust model, dust masses can be estimated to within ~0.07 dex accuracy using the 500 micron luminosity nu L_nu(500) alone. There are additional systematic errors arising from the choice of dust model, but these are hard to estimate. These calibrated prescriptions may be useful for studies of high-redshift galaxies.
K. L. Larson, T. Díaz-Santos, L. Armus, G. C. Privon, S. T. Linden, A. S. Evans, J. Howell, V. Charmandaris, V. U, D. B. Sanders, S. Stierwalt, L. Barcos-Muñoz, J. Rich, A. Medling, D. Cook, A. Oklopĉić, E. J. Murphy, P. Bonfini We present $HST$ narrow-band near-infrared imaging of Pa$\alpha$ and Pa$\beta$ emission of 48 local Luminous Infrared Galaxies (LIRGs) from the Great Observatories All-Sky LIRG Survey (GOALS). These data allow us to measure the properties of 810 spatially resolved star-forming regions (59 nuclei and 751 extra-nuclear clumps), and directly compare their properties to those found in both local and high-redshift star-forming galaxies. We find that in LIRGs, the star-forming clumps have radii ranging from $\sim90-900$ pc and star formation rates (SFRs) of $\sim1\times10^{-3}$ to 10 M$_\odot$yr$^{-1}$, with median values for extra-nuclear clumps of 170 pc and 0.03 M$_\odot$yr$^{-1}$. The detected star-forming clumps are young, with a median stellar age of $8.7$ Myrs, and a median stellar mass of $5\times10^{5}$ M$_\odot$. The SFRs span the range of those found in normal local star-forming galaxies to those found in high-redshift star-forming galaxies at $\rm{z}=1-3$. The luminosity function of the LIRG clumps has a flatter slope than found in lower-luminosity, star-forming galaxies, indicating a relative excess of luminous star-forming clumps. In order to predict the possible range of star-forming histories and gas fractions, we compare the star-forming clumps to those measured in the MassiveFIRE high-resolution cosmological simulation. The star-forming clumps in MassiveFIRE cover the same range of SFRs and sizes found in the local LIRGs and have total gas fractions that extend from 10 to 90%. If local LIRGs are similar to these simulated galaxies, we expect future observations with ALMA will find a large range of gas fractions, and corresponding star formation efficiencies, among the star-forming clumps in LIRGs.
Anne M. Medling, George C. Privon, Loreto Barcos-Muñoz, Ezequiel Treister, Claudia Cicone, Hugo Messias, David B. Sanders, Nick Scoville, Vivian U, Lee Armus, Franz E. Bauer, Chin-Shin Chang, Julia M. Comerford, Aaron S. Evans, Claire E. Max, Francisco Müller-Sánchez, Neil Nagar, Kartik Sheth Dynamical black hole mass measurements in some gas-rich galaxy mergers indicate that they are overmassive relative to their host galaxy properties. Overmassive black holes in these systems present a conflict with the standard progression of galaxy merger - quasar evolution; an alternative explanation is that a nuclear concentration of molecular gas driven inward by the merger is affecting these dynamical black hole mass estimates. We test for the presence of such gas near the two black holes in NGC 6240 using long-baseline ALMA Band 6 observations (beam size 0"06 $\times$ 0"03 or 30 pc$\times$15 pc). We find (4.2-9.8) $\times10^{7}$ M$_{\odot}$ and (1.2-7.7) $\times10^{8}$ M$_{\odot}$ of molecular gas within the resolution limit of the original black hole mass measurements for the north and south black holes, respectively. In the south nucleus, this measurement implies that 6-89% of the original black hole mass measurement actually comes from molecular gas, resolving the tension in the original black hole scaling relations. For the north, only 5% to 11% is coming from molecular gas, suggesting the north black hole is actually overmassive. Our analysis provides the first measurement of significant molecular gas masses contaminating dynamical black hole mass measurements. These high central molecular gas densities further present a challenge to theoretical black hole accretion prescriptions, which often assume accretion proceeds rapidly through the central 10 pc.
B. Mennesson, V. Bailey, J. Kasdin, J. Trauger, O. Absil, R. Akeson, L. Armus, J. L. Baudino, P. Baudoz, A. Bellini, D. Bennett, B. Berriman, A. Boccaletti, S. Calchi-Novati, K. Carpenter, C. Chen, W. Danchi, J. Debes, D. Defrere, S. Ertel, et al (43) The Wide Field Infrared Survey Telescope (WFIRST) Coronagraph Instrument (CGI) will be the first high-performance stellar coronagraph using active wavefront control for deep starlight suppression in space, providing unprecedented levels of contrast, spatial resolution, and sensitivity for astronomical observations in the optical. One science case enabled by the CGI will be taking images and(R~50)spectra of faint interplanetary dust structures present in the habitable zone of nearby sunlike stars (~10 pc) and within the snow-line of more distant ones(~20pc), down to dust density levels commensurate with that of the solar system zodiacal cloud. Reaching contrast levels below~10-7 for the first time, CGI will cross an important threshold in debris disks physics, accessing disks with low enough optical depths that their structure is dominated by transport phenomena than collisions. Hence, CGI results will be crucial for determining how exozodiacal dust grains are produced and transported in low-density disks around mature stars. Additionally, CGI will be able to measure the brightness level and constrain the degree of asymmetry of exozodiacal clouds around individual nearby sunlike stars in the optical, at the ~10x solar zodiacal emission level. This information will be extremely valuable for optimizing the observational strategy of possible future exo-Earth direct imaging missions, especially those planning to operate at optical wavelengths, such as Habitable Exoplanet Observatory (HabEx) and the Large Ultraviolet/Optical/Infrared Surveyor (LUVOIR).
Yun Wang, Mark Dickinson, Lynne Hillenbrand, Massimo Robberto, Lee Armus, Mario Ballardini, Robert Barkhouser, James Bartlett, Peter Behroozi, Robert A. Benjamin, Jarle Brinchmann, Ranga-Ram Chary, Chia-Hsun Chuang, Andrea Cimatti, Charlie Conroy, Robert Content, Emanuele Daddi, Megan Donahue, Olivier Dore, Peter Eisenhardt, et al (28) ATLAS (Astrophysics Telescope for Large Area Spectroscopy) is a concept for a NASA probe-class space mission. It is the spectroscopic follow-up mission to WFIRST, boosting its scientific return by obtaining deep NIR & MIR slit spectroscopy for most of the galaxies imaged by the WFIRST High Latitude Survey at z>0.5. ATLAS will measure accurate and precise redshifts for ~200M galaxies out to z=7 and beyond, and deliver spectra that enable a wide range of diagnostic studies of the physical properties of galaxies over most of cosmic history. ATLAS and WFIRST together will produce a definitive 3D map of the Universe over 2000 sq deg. ATLAS Science Goals are: (1) Discover how galaxies have evolved in the cosmic web of dark matter from cosmic dawn through the peak era of galaxy assembly. (2) Discover the nature of cosmic acceleration. (3) Probe the Milky Way's dust-enshrouded regions, reaching the far side of our Galaxy. (4) Discover the bulk compositional building blocks of planetesimals formed in the outer Solar System. These flow down to the ATLAS Scientific Objectives: (1A) Trace the relation between galaxies and dark matter with less than 10% shot noise on relevant scales at 1<z<7. (1B) Probe the physics of galaxy evolution at 1<z<7. (2) Obtain definitive measurements of dark energy and tests of General Relativity. (3) Measure the 3D structure and stellar content of the inner Milky Way to a distance of 25 kpc. (4) Detect and quantify the composition of 3,000 planetesimals in the outer Solar System. ATLAS is a 1.5m telescope with a FoV of 0.4 sq deg, and uses Digital Micro-mirror Devices (DMDs) as slit selectors. It has a spectroscopic resolution of R = 1000, and a wavelength range of 1-4 microns. ATLAS has an unprecedented spectroscopic capability based on DMDs, with a spectroscopic multiplex factor ~6,000. ATLAS is designed to fit within the NASA probe-class space mission cost envelope.
R. Herrero-Illana, G. C. Privon, A. S. Evans, T. Díaz-Santos, M. Á. Pérez-Torres, V. U, A. Alberdi, K. Iwasawa, L. Armus, S. Aalto, J. Mazzarella, J. Chu, D. B. Sanders, L. Barcos-Muñoz, V. Charmandaris, S. T. Linden, I. Yoon, D. T. Frayer, H. Inami, D.-C. Kim, et al (6) We present IRAM-30m Telescope $^{12}$CO and $^{13}$CO observations of a sample of 55 luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs) in the local universe. This sample is a subset of the Great Observatory All-Sky LIRG Survey (GOALS), for which we use ancillary multi-wavelength data to better understand their interstellar medium and star formation properties. Fifty-three (96%) of the galaxies are detected in $^{12}$CO, and 29 (52%) are also detected in $^{13}$CO above a 3$\sigma$ level. The median full width at zero intensity (FWZI) velocity of the CO line emission is 661km s$^{-1}$, and $\sim$54% of the galaxies show a multi-peak CO profile. Herschel photometric data is used to construct the far-IR spectral energy distribution of each galaxy, which are fit with a modified blackbody model that allows us to derive dust temperatures and masses, and infrared luminosities. We make the assumption that the gas-to-dust mass ratio of (U)LIRGs is comparable to local spiral galaxies with a similar stellar mass (i.e., gas/dust of mergers is comparable to their progenitors) to derive a CO-to-H$_2$ conversion factor of $\langle\alpha\rangle=1.8^{+1.3}_{-0.8}M_\odot$(K km s$^{-1}$pc$^{2}$)$^{-1}$; such a value is comparable to that derived for (U)LIRGs based on dynamical mass arguments. We derive gas depletion times of $400-600$Myr for the (U)LIRGs, compared to the 1.3Gyr for local spiral galaxies. Finally, we re-examine the relationship between the $^{12}$CO/$^{13}$CO ratio and dust temperature, confirming a transition to elevated ratios in warmer systems.
S. T. Linden, Y. Song, A. S. Evans, E. J. Murphy, L. Armus, L. Barcos-Muñoz, K. Larson, T. Díaz-Santos, G. C. Privon, J. Howell, J. A. Surace, V. Charmandaris, V. U, A. M. Medling, J. Chu, E. Momjian We present the first results of a high-resolution Karl G. Jansky Very Large Array (VLA) imaging survey of luminous and ultra-luminous infrared galaxies (U/LIRGs) in the Great Observatories All-Sky LIRG Survey (GOALS). From the full sample of 68 galaxies, we have selected 25 LIRGs that show resolved extended emission at sufficient sensitivity to image individual regions of star-formation activity beyond the nucleus.~With wideband radio continuum observations, which sample the frequency range from $3-33$ GHz, we have made extinction-free measurements of the luminosities and spectral indicies for a total of 48 individual star-forming regions identified as having de-projected galactocentric radii ($r_{G}$) that lie outside the 13.2$\mu$m core of the galaxy.~The median $3-33$ GHz spectral index and 33 GHz thermal fraction measured for these "extranuclear" regions is $-0.51 \pm 0.13$ and $65 \pm 11\%$ respectively.~These values are consistent with measurements made on matched spatial scales in normal star-forming galaxies, and suggests that these regions are more heavily-dominated by thermal free-free emission relative to the centers of local ULIRGs.~Further, we find that the median star-formation rate derived for these regions is $\sim 1 M_{\odot}$ yr$^{-1}$, and when we place them on the sub-galactic star-forming main sequence of galaxies (SFMS), we find they are offset from their host galaxies' globally-averaged specific star-formation rates (sSFRs).~We conclude that while nuclear starburst activity drives LIRGs above the SFMS, extranuclear star-formation still proceeds in a more extreme fashion relative to what is seen in local spiral galaxies.
Sarah Sadavoy, Mikako Matsuura, Lee Armus, Cara Battersby, Caitlin Casey, Christopher Clark, Asantha Cooray, Karine Demyk, Neal Evans, Karl Gordon, Frederic Galliano, Maryvonne Gerin, Benne Holwerda, Nia Imara, Doug Johnstone, Alvaro Labiano, David Leisawitz, Wanggi Lim, Leslie Looney, Margaret Meixner, et al (7) Dust offers a unique probe of the interstellar medium (ISM) across multiple size, density, and temperature scales. Dust is detected in outflows of evolved stars, star-forming molecular clouds, planet-forming disks, and even in galaxies at the dawn of the Universe. These grains also have a profound effect on various astrophysical phenomena from thermal balance and extinction in galaxies to the building blocks for planets, and changes in dust grain properties will affect all of these phenomena. A full understanding of dust in all of its forms and stages requires a multi-disciplinary investigation of the dust life cycle. Such an investigation can be achieved with a statistical study of dust properties across stellar evolution, star and planet formation, and redshift. Current and future instrumentation will enable this investigation through fast and sensitive observations in dust continuum, polarization, and spectroscopy from near-infrared to millimeter wavelengths.
Roberta Paladini, Matthew Povich, Lee Armus, Cara Battersby, Bruce Elmegreen, Adam Ginsburg, Doug Johnstone, David Leisawitz, Peregrine McGehee, Sarah Sadavoy, Marta Sewilo, Alessio Traficante, Martina Wiedner It is usually assumed that the stellar initial mass function (IMF) takes a universal form and that there exists a direct mapping between this and the distribution of natal core masses (the core mass function, CMF). The IMF and CMF have been best characterized in the Solar neighborhood. Beyond 500~pc from the Sun, in diverse environments where metallicity varies and massive star feedback may dominate, the IMF has been measured only incompletely and imprecisely, while the CMF has hardly been measured at all. In order to establish if the IMF and CMF are indeed universal and related to each other, it is necessary to: 1) perform multi-wavelength large-scale imaging and spectroscopic surveys of different environments across the Galaxy; 2) require an angular resolution of < 0.1'' in the optical/near-IR for stars and < 5'' in the far-IR for cores; 3) achieve far-IR sensitivities to probe 0.1~Msol cores at 2--3 kpc.
Alberto D. Bolatto, Lee Armus, Sylvain Veilleux, Adam K. Leroy, Fabian Walter, Richard Mushotzky, Karin M. Sandstrom, Paul Martini, Evan E. Schneider, Tony Wong, Roberto Decarli, Caitlin Casey, Dominik Riechers, David Meier, Desika Narayana There is wide consensus that galaxy outflows are one of the most important processes determining the evolution of galaxies through cosmic time, for example playing a key role in shaping the galaxy mass function. Our understanding of outflows and their drivers, however, is in its infancy --- this is particularly true for the cold (neutral atomic and molecular) phases of outflows, which present observational and modeling challenges. Here we outline several key open questions, briefly discussing the requirements of the observations necessary to make progress, and the relevance of several existing and planned facilities. It is clear that galaxy outflows, and particularly cold outflows, will remain a topic of active research for the next decade and beyond.
Philip Appleton, Lee Armus, Francois Boulanger, Charles M. Bradford, Jonathan Braine, Volker Bromm, Peter Capak, Michelle Cluver, Asantha Cooray, Tanio Diaz-Santos, Eiichi Egami, Bjorn Emonts, Pierre Guillard, George Helou, Lauranne Lanz, Susanne Madden, Anne Medling, Ewan O'Sullivan, Patrick Ogle, Alexandra Pope, et al (5) Galaxy formation depends on a complex interplay between gravitational collapse, gas accretion, merging, and feedback processes. Yet, after many decades of investigation, these concepts are poorly understood. This paper presents the argument that warm H$_2$ can be used as a tool to unlock some of these mysteries. Turbulence, shocks and outflows, driven by star formation, AGN activity or inflows, may prevent the rapid buildup of star formation in galaxies. Central to our understanding of how gas is converted into stars is the process by which gas can dissipate its mechanical energy through turbulence and shocks in order to cool. H$_2$ lines provide direct quantitative measurements of kinetic energy dissipation in molecular gas in galaxies throughout the Universe. Based on the detection of very powerful H$_2$ lines from z = 2 galaxies and proto-clusters at the detection limits of \it Spitzer, we are confident that future far-IR and UV H$_2$ observations will provide a wealth of new information and insight into galaxy evolution to high-z. Finally, at the very earliest epoch of star and galaxy formation, warm H$_2$ may also provide a unique glimpse of molecular gas collapse at 7 $<$ z $<$ 12 in massive dark matter (DM) halos on their way to forming the very first galaxies. Such measurements are beyond the reach of existing and planned observatories.
Alexandra Pope, Lee Armus, Eric Murphy, Susanne Aalto, David Alexander, Philip Appleton, Amy Barger, Matt Bradford, Peter Capak, Caitlin Casey, Vassilis Charmandaris, Ranga Chary, Asantha Cooray, Jim Condon, Tanio Diaz Santos, Mark Dickinson, Duncan Farrah, Carl Ferkinhoff, Norman Grogin, Ryan Hickox, et al (12) Galaxies grow their supermassive black holes in concert with their stars, although the relationship between these major galactic components is poorly understood. Observations of the cosmic growth of stars and black holes in galaxies suffer from disjoint samples and the strong effects of dust attenuation. The thermal infrared holds incredible potential for simultaneously measuring both the star formation and black hole accretion rates in large samples of galaxies covering a wide range of physical conditions. Spitzer demonstrated this potential at low redshift, and by observing some of the most luminous galaxies at z~2. JWST will apply these methods to normal galaxies at these epochs, but will not be able to generate large spectroscopic samples or access the thermal infrared at high-redshift. An order of magnitude gap in our wavelength coverage will persist between JWST and ALMA. A large, cold infrared telescope can fill this gap to determine when (in cosmic time), and where (within the cosmic web), stars and black holes co-evolve, by measuring these processes simultaneously in statistically complete and unbiased samples of galaxies to z>8. A next-generation radio interferometer will have the resolution and sensitivity to measure star-formation and nuclear accretion in even the dustiest galaxies. Together, the thermal infrared and radio can uniquely determine how stars and supermassive blackholes co-evolve in galaxies over cosmic time.
Caitlin M. Casey, Peter Capak, Johannes Staguhn, Lee Armus, Andrew Blain, Matthieu Bethermin, Jaclyn Champagne, Asantha Cooray, Kristen Coppin, Patrick Drew, Eli Dwek, Steven Finkelstein, Maximilien Franco, James Geach, Jacqueline Hodge, Maciej Koprowski, Claudia Lagos, Desika Narayanan, Alexandra Pope, David Sanders, et al (7) Two decades of effort have been poured into both single-dish and interferometric millimeter-wave surveys of the sky to infer the volume density of dusty star-forming galaxies (DSFGs, with SFR>100M$_\odot$ yr$^{-1}$) over cosmic time. Though obscured galaxies dominate cosmic star-formation near its peak at $z\sim2$, the contribution of such heavily obscured galaxies to cosmic star-formation is unknown beyond $z\sim2.5$ in contrast to the well-studied population of Lyman-break galaxies (LBGs) studied through deep, space- and ground-based pencil beam surveys in the near-infrared. Unlocking the volume density of DSFGs beyond $z>3$, particularly within the first 1 Gyr after the Big Bang is critical to resolving key open questions about early Universe galaxy formation: (1) What is the integrated star-formation rate density of the Universe in the first few Gyr and how is it distributed among low-mass galaxies (e.g. Lyman-break galaxies) and high-mass galaxies (e.g. DSFGs and quasar host galaxies)? (2) How and where do the first massive galaxies assemble? (3) What can the most extreme DSFGs teach us about the mechanisms of dust production (e.g. supernovae, AGB stars, grain growth in the ISM) <1 Gyr after the Big Bang? We summarize the types of observations needed in the next decade to address these questions.
Asantha Cooray, James Aguirre, Yacine Ali-Haimoud, Marcelo Alvarez, Phil Appleton, Lee Armus, George Becker, Jamie Bock, Rebecca Bowler, Judd Bowman, Matt Bradford, Patrick Breysse, Volker Bromm, Jack Burns, Karina Caputi, Marco Castellano, Tzu-Ching Chang, Ranga Chary, Hsin Chiang, Joanne Cohn, et al (54) The cosmic dawn and epoch of reionization mark the time period in the universe when stars, galaxies, and blackhole seeds first formed and the intergalactic medium changed from neutral to an ionized one. Despite substantial progress with multi-wavelength observations, astrophysical process during this time period remain some of the least understood with large uncertainties on our existing models of galaxy, blackhole, and structure formation. This white paper outlines the current state of knowledge and anticipated scientific outcomes with ground and space-based astronomical facilities in the 2020s. We then propose a number of scientific goals and objectives for new facilities in late 2020s to mid 2030s that will lead to definitive measurements of key astrophysical processes in the epoch of reionization and cosmic dawn.
Matteo Bonato, Gianfranco De Zotti, David Leisawitz, Mattia Negrello, Marcella Massardi, Ivano Baronchelli, Zhen-Yi Cai, Charles M. Bradford, Alexandra Pope, Eric J. Murphy, Lee Armus, Asantha Cooray We illustrate the extraordinary potential of the (far-IR) Origins Survey Spectrometer (OSS) on board the Origins Space Telescope (OST) to address a variety of open issues on the co-evolution of galaxies and AGNs. We present predictions for blind surveys, each of 1000 h, with different mapped areas (a shallow survey covering an area of 10 deg$^{2}$ and a deep survey of 1 deg$^{2}$) and two different concepts of the OST/OSS: with a 5.9 m telescope (Concept 2, our reference configuration) and with a 9.1 m telescope (Concept 1, previous configuration). In 1000 h, surveys with the reference concept will detect from $\sim 1.9 \times 10^{6}$ to $\sim 8.7 \times 10^{6}$ lines from $\sim 4.8 \times 10^{5}$-$2.7 \times 10^{6}$ star-forming galaxies and from $\sim 1.4 \times 10^{4}$ to $\sim 3.8 \times 10^{4}$ lines from $\sim 1.3 \times 10^{4}$-$3.5 \times 10^{4}$ AGNs. The shallow survey will detect substantially more sources than the deep one; the advantage of the latter in pushing detections to lower luminosities/higher redshifts turns out to be quite limited. The OST/OSS will reach, in the same observing time, line fluxes more than one order of magnitude fainter than the SPICA/SMI and will cover a much broader redshift range. In particular it will detect tens of thousands of galaxies at $z \geq 5$, beyond the reach of that instrument. The polycyclic aromatic hydrocarbons lines are potentially bright enough to allow the detection of hundreds of thousands of star-forming galaxies up to $z \sim 8.5$, i.e. all the way through the re-ionization epoch. The proposed surveys will allow us to explore the galaxy-AGN co-evolution up to $z\sim 5.5-6$ with very good statistics. OST Concept 1 does not offer significant advantages for the scientific goals presented here.
Rachel Akeson, Lee Armus, Etienne Bachelet, Vanessa Bailey, Lisa Bartusek, Andrea Bellini, Dominic Benford, David Bennett, Aparna Bhattacharya, Ralph Bohlin, Martha Boyer, Valerio Bozza, Geoffrey Bryden, Sebastiano Calchi Novati, Kenneth Carpenter, Stefano Casertano, Ami Choi, David Content, Pratika Dayal, Alan Dressler, et al (86) The Wide Field Infrared Survey Telescope (WFIRST) is a 2.4m space telescope with a 0.281 deg^2 field of view for near-IR imaging and slitless spectroscopy and a coronagraph designed for > 10^8 starlight suppresion. As background information for Astro2020 white papers, this article summarizes the current design and anticipated performance of WFIRST. While WFIRST does not have the UV imaging/spectroscopic capabilities of the Hubble Space Telescope, for wide field near-IR surveys WFIRST is hundreds of times more efficient. Some of the most ambitious multi-cycle HST Treasury programs could be executed as routine General Observer (GO) programs on WFIRST. The large area and time-domain surveys planned for the cosmology and exoplanet microlensing programs will produce extraordinarily rich data sets that enable an enormous range of Archival Research (AR) investigations. Requirements for the coronagraph are defined based on its status as a technology demonstration, but its expected performance will enable unprecedented observations of nearby giant exoplanets and circumstellar disks. WFIRST is currently in the Preliminary Design and Technology Completion phase (Phase B), on schedule for launch in 2025, with several of its critical components already in production.
Vanessa P. Bailey, Lee Armus, Bala Balasubramanian, Pierre Baudoz, Andrea Bellini, Dominic Benford, Bruce Berriman, Aparna Bhattacharya, Anthony Boccaletti, Eric Cady, Sebastiano Calchi Novati, Kenneth Carpenter, David Ciardi, Brendan Crill, William Danchi, John Debes, Richard Demers, Kjetil Dohlen, Robert Effinger, Marc Ferrari, et al (67) The Wide Field Infrared Survey Telescope (WFIRST) Coronagraph Instrument (CGI) is a high-contrast imager and integral field spectrograph that will enable the study of exoplanets and circumstellar disks at visible wavelengths. Ground-based high-contrast instrumentation has fundamentally limited performance at small working angles, even under optimistic assumptions for 30m-class telescopes. There is a strong scientific driver for better performance, particularly at visible wavelengths. Future flagship mission concepts aim to image Earth analogues with visible light flux ratios of more than 10^10. CGI is a critical intermediate step toward that goal, with a predicted 10^8-9 flux ratio capability in the visible. CGI achieves this through improvements over current ground and space systems in several areas: (i) Hardware: space-qualified (TRL9) deformable mirrors, detectors, and coronagraphs, (ii) Algorithms: wavefront sensing and control; post-processing of integral field spectrograph, polarimetric, and extended object data, and (iii) Validation of telescope and instrument models at high accuracy and precision. This white paper, submitted to the 2018 NAS Exoplanet Science Strategy call, describes the status of key CGI technologies and presents ways in which performance is likely to evolve as the CGI design matures.