Search

Dual active galactic nuclei (AGN) offer a unique opportunity to probe the relationship between super massive black holes (SMBH) and their host galaxies as well as the role of major mergers in triggering AGN activity. The confirmed dual AGN Mrk 266 has been studied extensively with multi-wavelength imaging. Now, high spatial resolution IFU spectroscopy of Mrk 266 provides an opportunity to probe the kinematics of both the merger event and AGN feedback. We present for the first time high spatial resolution kinematic maps for both nuclei of Mrk 266 obtained with the Keck OSIRIS IFU spectrograph, utilizing adaptive optics to achieve a resolution of 0.31" and 0.20" for the NE and SW nuclei, respectively. Using the M-sigma relation for mergers, we infer a SMBH mass of approximately 7e7 solar masses for the southwestern nucleus. Additionally, we report that the molecular gas kinematics of the southwestern nucleus are dominated by rotation rather than large-scale chaotic motions. The southwest nucleus also contains both a circumnuclear ring of star formation from which an inflow of molecular gas is likely fueling the AGN and a compact, AGN-dominated outflow of highly ionized gas with a timescale of approximately 2 Myr, significantly shorter than the timescale of the merger. The northeastern nucleus, on the other hand, exhibits complex kinematics related to the merger, including molecular gas that appears to have decoupled from the rotation of the stars. Our results suggest that while the AGN activity in Mrk 266 was likely triggered during the merger, AGN feeding is currently the result of processes internal to each host galaxy, thus resulting in a strong asymmetry between the two nuclei., Comment: 18 pages, 18 figures, 3 tables. Accepted for publication in MNRAS

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., Comment: 18 pages, 7 figures, 1 table; Accepted for publication in ApJL

Our understanding of early-type galaxies (ETGs) has grown in the past decade with the advance of full-spectrum fitting techniques used to infer the properties of the stellar populations that make-up the galaxy. We present ages, central velocity dispersions, and abundance ratios relative to Fe of C, N, O, Mg, Si, Ca, Ti, Cr, Mn, Co, Ni, Cu, Sr, Ba, and Eu, derived using full-spectrum fitting techniques for three ETGs NGC 2865, NGC 3818, and NGC 4915. Each of these three galaxies were selected because they have optical, disturbed structures (fine structure) that are linked to major merger events that occurred 1, 7, and 6 Gyr ago, respectively. Two of the ETGs, NGC 3818 and NGC 4915, show chemical signatures similar to ETGs without fine structure, which is consistent with a gas-poor merger of elliptical galaxies in which substantial star formation is not expected. For NGC 2865, we find a statistically higher abundance of Ca (an $\alpha$-element) and Cr and Mn (Fe-peak elements). We show that for NGC 2865, a simple gas-rich merger scenario fails to explain the larger abundance ratios compared to ETGs without fine structure. These three early-type galaxies with fine structure exhibit a range of abundances, suggesting ETGs with fine structure can form via multiple pathways and types of galaxy mergers., Comment: 26 pages, 8 figures

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., Comment: 16 pages, 8 figures. Submitted to ApJ

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., Comment: 28 pages, 17 figures. Accepted for publication in ApJ. This version includes small revisions following the referee report

Low mass galaxy pair fractions are understudied, and it is unclear whether low mass pair fractions evolve in the same way as more massive systems over cosmic time. In the era of JWST, Roman, and Rubin, selecting galaxy pairs in a self-consistent way will be critical to connect observed pair fractions to cosmological merger rates across all mass scales and redshifts. Utilizing the Illustris TNG100 simulation, we create a sample of physically associated low mass ($\rm 10^81$. At $z=0$, the low mass major (1:4 mass ratio) pair fraction is 4$\times$ lower than high mass pairs, consistent with findings for cosmological merger rates. We show that separation limits that vary with the mass and redshift of the system, such as scaling by the virial radius of the host halo ($r_{\mathrm{sep}}< 1 R_{\rm vir}$), are critical for recovering pair fraction differences between low mass and high mass systems. Alternatively, static physical separation limits applied equivalently to all galaxy pairs do not recover the differences between low and high mass pair fractions, even up to separations of $300$ kpc. Finally, we place isolated mass-analogs of Local Group galaxy pairs, i.e., Milky Way (MW)--M31, MW--LMC, LMC--SMC, in a cosmological context, showing that isolated analogs of LMC--SMC-mass pairs and low-separation ($<50$ kpc) MW--LMC-mass pairs are $2-3\times$ more common at $z\gtrsim2-3$., Comment: 18 pages, 5 figures, accepted for publication in ApJ, updated Feb 2024 to reflect accepted version

We present a catalog of the millimeter-wave (mm-wave) continuum properties of 98 nearby ($z <$ 0.05) active galactic nuclei (AGNs) selected from the 70-month Swift/BAT hard X-ray catalog that have precisely determined X-ray spectral properties and subarcsec-resolution ALMA Band-6 (211--275 GHz) observations as of 2021 April. Due to the hard-X-ray ($>$ 10 keV) selection, the sample is nearly unbiased for obscured systems at least up to Compton-thick-level obscuration, and provides the largest number of AGNs with high physical resolution mm-wave data ($\lesssim$ 100--200 pc). Our catalog reports emission peak coordinates, spectral indices, and peak fluxes and luminosities at 1.3 mm (230 GHz). Additionally, high-resolution mm-wave images are provided. Using the images and creating radial surface brightness profiles of mm-wave emission, we identify emission extending from the central source and isolated blob-like emission. Flags indicating the presence of these emission features are tabulated. Among 90 AGNs with significant detections of nuclear emission, 37 AGNs ($\approx$ 41%) appear to have both or one of extended or blob-like components. We, in particular, investigate AGNs that show well-resolved mm-wave components and find that these seem to have a variety of origins (i.e., a jet, radio lobes, a secondary AGN, stellar clusters, a narrow line region, galaxy disk, active star-formation regions, and AGN-driven outflows), and some components have currently unclear origins., Comment: 49 pages, 7 figures, 3 tables, accepted for publication in ApJS

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]/Pa\b{eta} 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., Comment: 15 pages, 4 figures, 3 tables, accepted for publication in ApJ

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., Comment: 14 pages, 5 figures, 2 tables, Submitted to ApJL

In the current era of high-z galaxy discovery with JWST and ALMA, our ability to study the stellar populations and ISM conditions in a diverse range of galaxies at Cosmic Dawn has rapidly improved. At the same time, the need to understand the current limitations in modeling galaxy formation processes and physical properties in order to interpret these observations is critical. Here, we study the challenges in modeling galaxy dust temperatures, both in the context of forward modeling galaxy spectral properties from a hydrodynamical simulation and via backwards modeling galaxy physical properties from mock observations of far-infrared dust emission. Using the Simba model for galaxy formation combined with Powderday radiative transfer, we can accurately predict the evolution of dust at high redshift, though several aspects of the model are essentially free parameters (dust composition, sub-resolution dust in star-forming regions) that dull the predictive power of the model dust temperature distributions. We also highlight the uncertainties in the backwards modeling methods, where we find the commonly used models and assumptions to fit FIR SEDs and infer dust temperatures (e.g., single temperature, optically thin modified blackbody) largely fail to capture the complexity of high-z dusty galaxies. We caution that conclusions inferred from both simulations -- limited by resolution and post-processing techniques -- and observations -- limited by sparse data and simplistic model parameterizations -- are susceptible to unique and nuanced uncertainties that can limit the usefulness of current high-z dust measurements., Comment: 20 pages main text, 12 figures. Accepted to ApJ, newest version

Recent studies have proposed that the nuclear millimeter continuum emission observed in nearby active galactic nuclei (AGN) could be created by the same population of electrons that gives rise to the X-ray emission that is ubiquitously observed in accreting black holes. We present the results of a dedicated high spatial resolution ($\sim$60-100 milliarcsecond) ALMA campaign on a volume-limited ($<50$ Mpc) sample of 26 hard X-ray ($>10$ keV) selected radio-quiet AGN. We find an extremely high detection rate (25/26 or $94^{+3}_{-6}\%$), which shows that nuclear emission at mm-wavelengths is nearly ubiquitous in accreting SMBHs. Our high-resolution observations show a tight correlation between the nuclear (1-23 pc) 100GHz and the intrinsic X-ray emission (1$\sigma$ scatter of $0.22$ dex). The ratio between the 100GHz continuum and the X-ray emission does not show any correlation with column density, black hole mass, Eddington ratio or star formation rate, which suggests that the 100GHz emission can be used as a proxy of SMBH accretion over a very broad range of these parameters. The strong correlation between 100GHz and X-ray emission in radio-quiet AGN could be used to estimate the column density based on the ratio between the observed 2-10keV ($F^{\rm obs}_{2-10\rm\,keV}$) and 100GHz ($F_{100\rm\,GHz}$) fluxes. Specifically, a ratio $\log (F^{\rm obs}_{2-10\rm\,keV}/F_{100\rm\,GHz})\leq 3.5$ strongly suggests that a source is heavily obscured [$\log (N_{\rm H}/\rm cm^{-2})\gtrsim 23.8$]. Our work shows the potential of ALMA continuum observations to detect heavily obscured AGN (up to an optical depth of one at 100GHz, i.e. $N_{\rm H}\simeq 10^{27}\rm\,cm^{-2}$), and to identify binary SMBHs with separations $<100$ pc, which cannot be probed by current X-ray facilities., Comment: ApJL in press

We present a phenomenological framework for starburst-driven neutrino production via proton-proton collisions and apply it to (ultra)luminous infrared galaxies (U/LIRGs) in the Great Observatories All-Sky LIRG Survey (GOALS). The framework relates the infrared luminosity of a GOALS galaxy, derived from consistently available Herschel Space Observatory data, to the expected starburst-driven neutrino flux. The model parameters that define this relation can be estimated from multiwavelength data. We apply the framework in a case study to the LIRG NGC 3690 (Arp 299, Mrk 171) and compare the obtained neutrino fluxes to the current sensitivity of the IceCube Neutrino Observatory. Using our framework, we also conclude that the neutrino emission in the LIRG NGC 1068, recently presented as the first steady IceCube neutrino point source, cannot be explained by a starburst-driven scenario and is therefore likely dominated by the active galactic nucleus in this galaxy. In addition to the single-source investigations, we also estimate the diffuse starburst-driven neutrino flux from GOALS galaxies and the total LIRG population over cosmic history., Comment: 26 pages, 11 figures, published version in Phys. Rev. D

We present multi-wavelength high-spatial resolution (~0.1'', 70 pc) observations of UGC 4211 at z=0.03474, a late-stage major galaxy merger at the closest nuclear separation yet found in near-IR imaging (0.32'', ~230 pc projected separation). Using Hubble Space Telescope/STIS, VLT/MUSE+AO, Keck/OSIRIS+AO spectroscopy, and ALMA observations, we show that the spatial distribution, optical and NIR emission lines, and millimeter continuum emission are all consistent with both nuclei being powered by accreting supermassive black holes (SMBHs). Our data, combined with common black hole mass prescriptions, suggests that both SMBHs have similar masses, log MBH~8.1 (south) and log MBH~8.3 (north), respectively. The projected separation of 230 pc (~6X the black hole sphere of influence) represents the closest-separation dual AGN studied to date with multi-wavelength resolved spectroscopy and shows the potential of nuclear (<50 pc) continuum observations with ALMA to discover hidden growing SMBH pairs. While the exact occurrence rate of close-separation dual AGN is not yet known, it may be surprisingly high, given that UGC 4211 was found within a small, volume-limited sample of nearby hard X-ray detected AGN. Observations of dual SMBH binaries in the sub-kpc regime at the final stages of dynamical friction provide important constraints for future gravitational wave observatories., Comment: 19 pages, 7 figures, 2 tables, published in ApJL

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., Comment: 13 pages, 4 figures, 1 table, published in the Astrophysical Journal

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., Comment: 13 pages, 4 figures, 2 tables, Submitted to ApJL

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., Comment: 4 figures, 1 table, Accepted by ApJL

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 {\sim}100 pc over the 4.9 - 7.6{\mu}m region are examined to study gas dynamics influenced by the central AGN. The low-ionization [Fe II] {\lambda}5.34{\mu}m and [Ar II] {\lambda}6.99{\mu}m lines are bright on the nucleus and in the starburst ring, as opposed to H2 S(5) {\lambda}6.91{\mu}m 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] {\sim}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., Comment: 14 pages, 6 figures, 1 table; accepted for publication in ApJL

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., Comment: 8 pages, 4 figures, 1 table, Submitted to ApJL

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., Comment: 3 figures, 1 table, Accepted for publication in ApJL

To understand the origin of nuclear ($\lesssim$ 100 pc) millimeter-wave (mm-wave) continuum emission in active galactic nuclei (AGNs), we systematically analyzed sub-arcsec resolution Band-6 (211-275 GHz) ALMA data of 98 nearby AGNs ($z <$ 0.05) from the 70-month Swift/BAT catalog. The sample, almost unbiased for obscured systems, provides the largest number of AGNs to date with high mm-wave spatial resolution sampling ($\sim$ 1-200 pc), and spans broad ranges of 14-150 keV luminosity {$40 < \log[L_{\rm 14-150}/({\rm erg\,s^{-1}})] < 45$}, black hole mass [$5 < \log(M_{\rm BH}/M_\odot) < 10$], and Eddington ratio ($-4 < \log \lambda_{\rm Edd} < 2$). We find a significant correlation between 1.3 mm (230 GHz) and 14-150 keV luminosities. Its scatter is $\approx$ 0.36 dex, and the mm-wave emission may serve as a good proxy of the AGN luminosity, free of dust extinction up to $N_{\rm H} \sim 10^{26}$ cm$^{-2}$. While the mm-wave emission could be self-absorbed synchrotron radiation around the X-ray corona according to past works, we also discuss different possible origins of the mm-wave emission; AGN-related dust emission, outflow-driven shocks, and a small-scale ($<$ 200 pc) jet. The dust emission is unlikely to be dominant, as the mm-wave slope is generally flatter than expected. Also, due to no increase in the mm-wave luminosity with the Eddington ratio, a radiation-driven outflow model is possibly not the common mechanism. Furthermore, we find independence of the mm-wave luminosity on indicators of the inclination angle from the polar axis of the nuclear structure, which is inconsistent with a jet model whose luminosity depends only on the angle., Comment: 44 pages, 35 figures, 3 tables, accepted for publication in ApJ; modified format, added supplementary figure (Fig. 32)

We present new central stellar velocity dispersions for 484 Sy 1.9 and Sy 2 from the second data release of the Swift/BAT AGN Spectroscopic Survey (BASS DR2). This constitutes the largest study of velocity dispersion measurements in X-ray selected, obscured AGN with 956 independent measurements of the Ca H+K and Mg b region (3880-5550A) and the Ca triplet region (8350-8730A) from 642 spectra mainly from VLT/Xshooter or Palomar/DoubleSpec. Our sample spans velocity dispersions of 40-360 km/s, corresponding to 4-5 orders of magnitude in black holes mass (MBH=10^5.5-9.6 Msun), bolometric luminosity (LBol~10^{42-46 ergs/s), and Eddington ratio (L/Ledd~10^{-5}-2). For 281 AGN, our data provide the first published central velocity dispersions, including 6 AGN with low mass black holes (MBH=10^5.5-6.5 Msun), discovered thanks to our high spectral resolution observations (sigma~25 km/s). The survey represents a significant advance with a nearly complete census of hard-X-ray selected obscured AGN with measurements for 99% of nearby AGN (z<0.1) outside the Galactic plane. The BASS AGN have higher velocity dispersions than the more numerous optically selected narrow line AGN (i.e., ~150 vs. ~100 km/s), but are not biased towards the highest velocity dispersions of massive ellipticals (i.e., >250 km/s). Despite sufficient spectral resolution to resolve the velocity dispersions associated with the bulges of small black holes (~10^4-5 Msun), we do not find a significant population of super-Eddington AGN. Using estimates of the black hole sphere of influence, direct stellar and gas black hole mass measurements could be obtained with existing facilities for more than ~100 BASS AGN., Comment: 32 pages, 16 figures, 8 tables, published in ApJS as part of BASS DR2 special issue

We present the AGN catalog and optical spectroscopy for the second data release of the Swift BAT AGN Spectroscopic Survey (BASS DR2). With this DR2 release we provide 1425 optical spectra, of which 1181 are released for the first time, for the 858 hard X-ray selected AGN in the Swift BAT 70-month sample. The majority of the spectra (813/1425, 57%) are newly obtained from VLT/Xshooter or Palomar/Doublespec. Many of the spectra have both higher resolution (R>2500, N~450) and/or very wide wavelength coverage (3200-10000 A, N~600) that are important for a variety of AGN and host galaxy studies. We include newly revised AGN counterparts for the full sample and review important issues for population studies, with 44 AGN redshifts determined for the first time and 780 black hole mass and accretion rate estimates. This release is spectroscopically complete for all AGN (100%, 858/858) with 99.8% having redshift measurements (857/858) and 96% completion in black hole mass estimates of unbeamed AGN (outside the Galactic plane). This AGN sample represents a unique census of the brightest hard X-ray selected AGN in the sky, spanning many orders of magnitude in Eddington ratio (Ledd=10^-5-100), black hole mass (MBH=10^5-10^10 Msun), and AGN bolometric luminosity (Lbol=10^40-10^47 ergs/s)., Comment: 34 pages, 11 figures, 15 tables, published in ApJS as part of BASS DR2 special issue

The BAT AGN Spectroscopic Survey (BASS) is designed to provide a highly complete census of the key physical parameters of supermassive black holes (SMBHs) that power local active galactic nuclei (AGN) (z<0.3), including their bolometric luminosity, black hole mass, accretion rates, and line-of-sight gas obscuration, and the distinctive properties of their host galaxies (e.g., star formation rates, masses, and gas fractions). We present an overview of the BASS data release 2 (DR2), an unprecedented spectroscopic survey in spectral range, resolution, and sensitivity, including 1449 optical (3200-10000 A) and 233 NIR (1-2.5 um) spectra for the brightest 858 ultra-hard X-ray (14-195 keV) selected AGN across the entire sky and essentially all levels of obscuration. This release provides a highly complete set of key measurements (emission line measurements and central velocity dispersions), with 99.9% measured redshifts and 98% black hole masses estimated (for unbeamed AGN outside the Galactic plane). The BASS DR2 AGN sample represents a unique census of nearby powerful AGN, spanning over 5 orders of magnitude in AGN bolometric luminosity, black hole mass, Eddington ratio, and obscuration. The public BASS DR2 sample and measurements can thus be used to answer fundamental questions about SMBH growth and its links to host galaxy evolution and feedback in the local universe, as well as open questions concerning SMBH physics. Here we provide a brief overview of the survey strategy, the key BASS DR2 measurements, data sets and catalogs, and scientific highlights from a series of DR2-based works., Comment: 21 pages, 6 figures, 5 tables, published in ApJS as part of BASS DR2 special issue

We present measurements of broad emission lines and virial estimates of supermassive black hole masses ($M_{BH}$) for a large sample of ultra-hard X-ray selected active galactic nuclei (AGNs) as part of the second data release of the BAT AGN Spectroscopic Survey (BASS/DR2). Our catalog includes $M_{BH}$ estimates for a total 689 AGNs, determined from the H$\alpha$, H$\beta$, $MgII\lambda2798$, and/or $CIV\lambda1549$ broad emission lines. The core sample includes a total of 512 AGNs drawn from the 70-month Swift/BAT all-sky catalog. We also provide measurements for 177 additional AGNs that are drawn from deeper Swift/BAT survey data. We study the links between $M_{BH}$ estimates and line-of-sight obscuration measured from X-ray spectral analysis. We find that broad H$\alpha$ emission lines in obscured AGNs ($\log (N_{\rm H}/{\rm cm}^{-2})> 22.0$) are on average a factor of $8.0_{-2.4}^{+4.1}$ weaker, relative to ultra-hard X-ray emission, and about $35_{-12}^{~+7}$\% narrower than in unobscured sources (i.e., $\log (N_{\rm H}/{\rm cm}^{-2}) < 21.5$). This indicates that the innermost part of the broad-line region is preferentially absorbed. Consequently, current single-epoch $M_{BH}$ prescriptions result in severely underestimated ($>$1 dex) masses for Type 1.9 sources (AGNs with broad H$\alpha$ but no broad H$\beta$) and/or sources with $\log (N_{\rm H}/{\rm cm}^{-2}) > 22.0$. We provide simple multiplicative corrections for the observed luminosity and width of the broad H$\alpha$ component ($L[{\rm b}{\rm H}\alpha]$ and FWHM[bH$\alpha$]) in such sources to account for this effect, and to (partially) remedy $M_{BH}$ estimates for Type 1.9 objects. As key ingredient of BASS/DR2, our work provides the community with the data needed to further study powerful AGNs in the low-redshift Universe., Comment: published in ApJS

Ionised gas kinematics provide crucial evidence of the impact that active galactic nuclei (AGN) have in regulating star formation in their host galaxies. Although the presence of outflows in AGN host galaxies has been firmly established, the calculation of outflow properties such as mass outflow rates and kinetic energy remains challenging. We present the [OIII]5007 ionised gas outflow properties of 22 z$<$0.1 X-ray AGN, derived from the BAT AGN Spectroscopic Survey using MUSE/VLT. With an average spatial resolution of 1" (0.1-1.2 kpc), the observations resolve the ionised gas clouds down to sub-kiloparsec scales. Resolved maps show that the [OIII] velocity dispersion is, on average, higher in regions ionised by the AGN, compared to star formation. We calculate the instantaneous outflow rates in individual MUSE spaxels by constructing resolved mass outflow rate maps, incorporating variable outflow density and velocity. We compare the instantaneous values with time-averaged outflow rates by placing mock fibres and slits on the MUSE field-of-view, a method often used in the literature. The instantaneous outflow rates (0.2-275 $M_{\odot}$ yr$^{-1}$) tend to be 2 orders of magnitude higher than the time-averaged outflow rates (0.001-40 $M_{\odot}$ yr$^{-1}$). The outflow rates correlate with the AGN bolometric luminosity ($L_{\rm bol}\sim$ 10$^{42.71}$-10$^{45.62}$ erg/s) but we find no correlations with black hole mass (10$^{6.1}$-10$^{8.9}$ M$_{\odot}$), Eddington ratio (0.002-1.1) and radio luminosity (10$^{21}$-10$^{26}$ W/Hz). We find the median coupling between the kinetic energy and $L_{\rm bol}$ to be 1%, consistent with the theoretical predictions for an AGN-driven outflow., Comment: Main paper: 20 pages, 15 figures, 3 tables. Accepted for publication in MNRAS

The Baldwin, Philips, & Terlevich diagram of [O III]/H$\beta$ vs. [N II]/H$\alpha$ (hereafter N2-BPT) has long been used as a tool for classifying galaxies based on the dominant source of ionizing radiation. Recent observations have demonstrated that galaxies at $z\sim2$ reside offset from local galaxies in the N2-BPT space. In this paper, we conduct a series of controlled numerical experiments to understand the potential physical processes driving this offset. We model nebular line emission in a large sample of galaxies, taken from the SIMBA cosmological hydrodynamic galaxy formation simulation, using the CLOUDY photoionization code to compute the nebular line luminosities from H II regions. We find that the observed shift toward higher [O III]/H$\beta$ and [N II]/H$\alpha$ values at high redshift arises from sample selection: when we consider only the most massive galaxies $M_* \sim 10^{10-11} M_\odot$, the offset naturally appears, due to their high metallicities. We predict that deeper observations that probe lower-mass galaxies will reveal galaxies that lie on a locus comparable to $z\sim 0$ observations. Even when accounting for sample selection effects, we find that there is a subtle mismatch between simulations and observations. To resolve this discrepancy, we investigate the impact of varying ionization parameters, H II region densities, gas-phase abundance patterns, and increasing radiation field hardness on N2-BPT diagrams. We find that either decreasing the ionization parameter or increasing the N/O ratio of galaxies at fixed O/H can move galaxies along a self-similar arc in N2-BPT space that is occupied by high-redshift galaxies., Comment: 25 pages, 15 figures. Accepted for publication in ApJ

Observations of cold molecular gas reservoirs are critical for understanding the shutdown of star formation in massive galaxies. While dust continuum is an efficient and affordable tracer, this method relies upon the assumption of a "normal" molecular-gas to dust mass ratio, $\delta_{\mathrm{GDR}}$, typically of order one hundred. Recent null detections of quiescent galaxies in deep dust continuum observations support a picture where the cold gas and dust has been rapidly depleted or expelled. In this work, we present another viable explanation: a significant fraction of galaxies with low star formation per unit stellar mass are predicted to have extreme $\delta_{\mathrm{GDR}}$ ratios. We show that simulated massive quiescent galaxies at $0 < z < 3$ in the \textsc{simba} cosmological simulations have $\delta_{\mathrm{GDR}}$ values that extend $>$4 orders of magnitude. The dust in most simulated quiescent galaxies is destroyed significantly more rapidly than the molecular gas depletes, and cannot be replenished. The transition from star-forming to quiescent halts dust formation via star formation processes, with dust subsequently destroyed by supernova shocks and thermal sputtering of dust grains embedded in hot plasma. After this point, the dust growth rate in the models is not sufficient to overcome the loss of $>$3 orders of magnitude in dust mass to return to normal values of $\delta_{\mathrm{GDR}}$ despite having high metallicity. Our results indicate that it is not straight forward to use a single observational indicator to robustly pre-select exotic versus normal ratios. These simulations make strong predictions that can be tested with millimeter facilities., Comment: 9 pages, 5 figures, version accepted for publication in the Astrophysical Journal Letters

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., Comment: 27 pages, 18 figures; accepted for publication in ApJ

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., Comment: 70 pages, 21 figures, 5 tables, accepted for publication in ApJ

Post-starburst (PSB), or 'E+A', galaxies represent a rapid transitional phase between major, gas-rich mergers and gas-poor, quiescent early-type galaxies. Surprisingly, many PSBs have been shown to host a significant interstellar medium (ISM), despite theoretical predictions that the majority of star-forming gas should be expelled in AGN- or starburst-driven outflows. To-date, the resolved properties of this surviving ISM have remained unknown. We present high resolution ALMA continuum and CO(2$-$1) observations in six gas- and dust-rich PSBs, revealing for the first time the spatial and kinematic structure of their ISM on sub-kpc scales. We find extremely compact molecular reservoirs, with dust and gas surface densities rivaling those found in (ultra-)luminous infrared galaxies. We observe spatial and kinematic disturbances in all sources, with some also displaying disk-like kinematics. Estimates of the internal turbulent pressure in the gas exceed those of normal star-forming disks by at least 2 orders of magnitude, and rival the turbulent gas found in local interacting galaxies, such as the Antennae. Though the source of this high turbulent pressure remains uncertain, we suggest that the high incidence of tidal disruption events (TDEs) in PSBs could play a role. The star formation in these PSBs' turbulent central molecular reservoirs is suppressed, forming stars only 10% as efficiently as starburst galaxies with similar gas surface densities. "The fall" of star formation in these galaxies was not precipitated by complete gas expulsion or redistribution. Rather, this high-resolution view of PSBs' ISM indicates that star formation in their remaining compact gas reservoirs is suppressed by significant turbulent heating., Comment: 19 pages, 9 figures, 3 tables; Accepted for publication in ApJ

We perform a systematic X-ray spectroscopic analysis of 57 local ultra/luminous infrared galaxy systems (containing 84 individual galaxies) observed with Nuclear Spectroscopic Telescope Array and/or Swift/BAT. Combining soft X-ray data obtained with Chandra, XMM-Newton, Suzaku and/or Swift/XRT, we identify 40 hard ($>$10 keV) X-ray detected active galactic nuclei (AGNs) and constrain their torus parameters with the X-ray clumpy torus model XCLUMPY (Tanimoto et al. 2019). Among the AGNs at $z < 0.03$, for which sample biases are minimized, the fraction of Compton-thick ($N_{\rm H} \geq 10^{24}$ cm$^{-2}$) AGNs reaches 64$^{+14}_{-15}$% (6/9 sources) in late mergers, while 24$^{+12}_{-10}$% (3/14 sources) in early mergers, consistent with the tendency reported by Ricci et al. (2017). We find that the bolometric AGN luminosities derived from the infrared data increase, but the X-ray to bolometric luminosity ratios decrease, with merger stage. The X-ray weak AGNs in late mergers ubiquitously show massive outflows at sub-pc to kpc scales. Among them, the most luminous AGNs ($L_{\rm bol,AGN} \sim 10^{46}$ erg s$^{-1}$) have relatively small column densities of $\lesssim$10$^{23}$ cm$^{-2}$ and almost super-Eddington ratios ($\lambda_{\rm Edd} \sim$ 1.0). Their torus covering factors ($C_{\rm T}^{\rm (22)} \sim 0.6$) are larger than those of Swift/BAT selected AGNs with similarly high Eddington ratios. These results suggest a scenario that, in the final stage of mergers, multiphase strong outflows are produced due to chaotic quasi-spherical inflows and the AGN becomes extremely X-ray weak and deeply buried due to obscuration by inflowing and/or outflowing material., Comment: 85 pages, 88 figures, accepted in ApJS

The IceCube Neutrino Observatory has detected high-energy astrophysical neutrinos in the TeV-PeV range. These neutrinos have an isotropic distribution on the sky, and therefore, likely originate from extragalactic sources. Active Galactic Nuclei form a class of astronomical objects which are promising neutrino source candidates given their high electromagnetic luminosity and potential ability to accelerate cosmic rays up to energies greater than 10$^{16}$ eV. Interactions of these cosmic rays within the AGN environment are expected to produce both neutrinos and pionic gamma rays. Some hadronic models of AGN emission suggest that such gamma rays can in turn interact with the dense photon fields of AGN and cascade down to hard X-rays and MeV gamma rays. We present an update on the IceCube stacking analysis searching for high-energy neutrinos from hard X-ray sources sampled from the $\textit{Swift}$-BAT AGN Spectroscopic Survey., Comment: Presented at the 37th International Cosmic Ray Conference (ICRC 2021). See arXiv:2107.06966 for all IceCube contributions, 11 pages, 5 figures

We present Integral Field Spectroscopic (IFS) observations of the nearby ($z\sim0.03$) dual Active Galactic Nuclei (AGN) Mrk 739, whose projected nuclear separation is $\sim$3.4~kpc, obtained with the Multi Unit Spectroscopic Explorer (MUSE) at the Very Large Telescope (VLT). We find that the galaxy has an extended AGN-ionized emission-line region extending up to $\sim 20$ kpc away from the nuclei, while star-forming regions are more centrally concentrated within 2 - 3 kpc. We model the kinematics of the ionized gas surrounding the East nucleus using a circular disk profile, resulting in a peak velocity of $237^{+26}_{-28}$ km s$^{-1}$ at a distance of $\sim 1.2$ kpc. The enclosed dynamical mass within 1.2 kpc is $\log M(M_{\odot})=10.20\pm0.06$, $\sim$1,000 times larger than the estimated supermassive black hole (SMBH) mass of Mrk 739E. The morphology and dynamics of the system are consistent with an early stage of the collision, where the foreground galaxy (Mrk 739W) is a young star-forming galaxy in an ongoing first passage with its background companion (Mrk 739E). Since the SMBH in Mrk 739W does not show evidence of being rapidly accreting, we claim that the northern spiral arms of Mrk 739W are ionized by the nuclear activity of Mrk 739E., Comment: 26 pages, 15 figures. Accepted for publication in ApJ

Using parsec scale resolution hydrodynamical adaptive mesh refinement simulations we have studied the mass transport process throughout a galactic merger. The aim of such study is to connect both the peaks of mass accretion rate onto the BHs and star formation bursts with both gravitational and hydrodynamic torques acting on the galactic gaseous component. Our merger initial conditions were chosen to mimic a realistic system. The simulations include gas cooling, star formation, supernovae feedback, and AGN feedback. Gravitational and hydrodynamic torques near pericenter passes trigger gas funneling to the nuclei which is associated with bursts of star formation and black hole growth. Such episodes are intimately related with both kinds of torques acting on the galactic gas. Pericenters trigger both star formation and mass accretion rates of $\sim$ few $(1-10)$ $M_\odot$/yr. Such episodes last $\sim$ $(50-75)$ Myrs. Close passes also can produce black hole accretion that approaches and reaches the Eddington rate, lasting $\sim$ few Myrs. Our simulation shows that both gravitational and hydrodynamic torques are enhanced at pericenter passes with gravitational torques tending to have higher values than the hydrodynamic torques throughout the merger. We also find that in the closest encounters, hydrodynamic and gravitational torques can be comparable in their effect on the gas, the two helping in the redistribution of both angular momentum and mass in the galactic disc. Such phenomena allow inward mass transport onto the BH influence radius, fueling the compact object and lighting up the galactic nuclei., Comment: 13 pages, 9 figures

We present the host galaxy molecular gas properties of a sample of 213 nearby (0.0110^44 erg/s) increases by ~10-100 between a molecular gas mass of 10^8.7 Msun and 10^10.2 Msun. Higher Eddington ratio AGN galaxies tend to have higher molecular gas masses and gas fractions. Higher column density AGN galaxies (Log NH>23.4) are associated with lower depletion timescales and may prefer hosts with more gas centrally concentrated in the bulge that may be more prone to quenching than galaxy wide molecular gas. The significant average link of host galaxy molecular gas supply to SMBH growth may naturally lead to the general correlations found between SMBHs and their host galaxies, such as the correlations between SMBH mass and bulge properties and the redshift evolution of star formation and SMBH growth., Comment: 53 pages, 37 figures, accepted in ApJ

We present Powderday, a flexible, fast, open-source dust radiative transfer package designed to interface with galaxy formation simulations. Powderday builds on FSPS population synthesis models, Hyperion dust radiative transfer, and employs yt to interface between different software packages. We include our stellar population synthesis modeling on the fly, which allows for significant run-time flexibility in the assumed stellar physics. We include a model for nebular line emission that can employ either precomputed Cloudy lookup tables (for efficiency), or direct photoionization calculations for all young stars (for flexibility). The dust content follows either observationally-motivated prescriptions, direct modeling from galaxy formation simulations, or a novel approach that includes the dust content via learning-based algorithms from the SIMBA cosmological galaxy formation simulation. AGN can additionally be included via a range of prescriptions. The output of these models are broadband SEDs, as well as filter-convolved images. Powderday is designed to eliminate last-mile efforts by researchers that employ different hydrodynamic galaxy formation models, and seamlessly interfaces with GIZMO, AREPO, GASOLINE, CHANGA, and ENZO. We demonstrate the capabilities of the code via three applications: a model for the star formation rate (SFR) - infrared luminosity relation in galaxies (including the impact of AGN); the impact of circumstellar dust around AGB stars on the mid-infrared emission from galaxy SEDs; and the impact of galaxy inclination angle on dust attenuation laws., Comment: Submitted to ApJ; comments welcome

We present the Atacama Large Millimeter/submillimeter Array (ALMA) observation of the Sub-millimeter galaxy (SMG) ID 141 at z=4.24 in the [N II] 205 $\mu$m line (hereafter [N II]) and the underlying continuum at (rest-frame) 197.6 $\mu$m. Benefiting from lensing magnification by a galaxy pair at z=0.595, ID 141 is one of the brightest z$>4$ SMGs. At the angular resolutions of $\sim1.2"$ to $1.5"$ ($1" \sim6.9$ kpc), our observation clearly separates, and moderately resolves the two lensed images in both continuum and line emission at $\rm S/N>5$ . Our continuum-based lensing model implies an averaged amplification factor of $\sim5.8$ and reveals that the de-lensed continuum image has the S\'ersic index $\simeq 0.95$ and the S\'ersic radius of $\sim0.18" (\sim 1.24$ kpc). Furthermore, the reconstructed [N II] velocity field in the source plane is dominated by a rotation component with a maximum velocity of $\sim 300$ km/s at large radii, indicating a dark matter halo mass of $\sim 10^{12}M_{\odot}$. This, together with the reconstructed velocity dispersion field being smooth and modest in value ($<100$ km/s) over much of the outer parts of the galaxy, favours the interpretation of ID 141 being a disk galaxy dynamically supported by rotation. The observed [N II]/CO (7-6) and [N II]/[C II] 158 $\mu$m line luminosity ratios, which are consistent with the corresponding line ratio vs. far-infrared color correlation from local luminous infrared galaxies, imply a de-lensed star formation rate of ($1.8\pm 0.6)\times10^3M_\odot$/yr and provide an independent estimate on the size of the star-forming region $0.7^{+0.3}_{-0.3}$ kpc in radius., Comment: 13 pages, 6 figures, 2 tables. Lensing model code can be found here https://gitlab.com/cxylzlx/tiny_lens

Theory predicts that a supermassive black hole binary (SMBHB) could be observed as a luminous active galactic nucleus (AGN) that periodically varies on the order of its orbital timescale. In X-rays, periodic variations could be caused by mechanisms including relativistic Doppler boosting and shocks. Here we present the first systematic search for periodic AGNs using $941$ hard X-ray light curves (14-195 keV) from the first 105 months of the Swift Burst Alert Telescope (BAT) survey (2004-2013). We do not find evidence for periodic AGNs in Swift-BAT, including the previously reported SMBHB candidate MCG+11$-$11$-$032. We find that the null detection is consistent with the combination of the upper-limit binary population in AGNs in our adopted model, their expected periodic variability amplitudes, and the BAT survey characteristics. We have also investigated the detectability of SMBHBs against normal AGN X-ray variability in the context of the eROSITA survey. Under our assumptions of a binary population and the periodic signals they produce which have long periods of hundreds of days, up to $13$% true periodic binaries can be robustly distinguished from normal variable AGNs with the ideal uniform sampling. However, we demonstrate that realistic eROSITA sampling is likely to be insensitive to long-period binaries because longer observing gaps reduce their detectability. In contrast, large observing gaps do not diminish the prospect of detecting binaries of short, few-day periods, as 19% can be successfully recovered, the vast majority of which can be identified by the first half of the survey., Comment: 17 pages, including 8 figures and 4 tables. Accepted for publication in ApJ

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., Comment: accepted to ApJL

The role of feedback in triggering or quenching star formation and hence driving galaxy evolution can be directly studied with high resolution integral field observations. The manifestation of feedback in shocks is particularly important to examine in galaxy mergers, where violent interactions of gas takes place in the interstellar medium during the course of the galactic collision. As part of our effort to systematically study the local population of luminous infrared galaxies within the Great Observatories All-Sky LIRG Survey, we undertook the Keck OSIRIS AO LIRG Analysis observing campaign to study the gas dynamics in the inner kiloparsec regions of these systems at spatial scales of a few 10s of parsecs. With high-resolution near-infrared adaptive optics-assisted integral-field observations taken with OSIRIS on the Keck Telescopes, we employ near-infrared diagnostics such as Brg and the ro-vibrationally excited H2 lines to quantify the nuclear star formation rate and identify feedback associated with shocked molecular gas seen in 21 nearby luminous infrared galaxies. Shocked molecular gas is preferentially found in the ultraluminous infrared systems, but may also be triggered at a lower luminosity, earlier merging stage. On circumnuclear scales, AGN have a strong effect on heating the surrounding molecular gas, though their coupling is not simply driven by AGN strength but rather is complicated by orientation, dust shielding, density, and other factors. We find that the nuclear star formation correlates with merger class and diminishing projected nuclear separations. These trends are largely consistent with the picture of merger-induced starbursts within the center of galaxy mergers., Comment: Accepted to ApJ; 31 pages (with 5 tables and 10 embedded figures) + figure set

We quantify the frequency of companions of low redshift ($0.013 < z < 0.0252$), dwarf galaxies ($2 \times 10^8$ M$_\odot <$ M$_{*} < 5 \times 10^9$ M$_\odot$) that are isolated from more massive galaxies in SDSS and compare against cosmological expectations using mock observations of the Illustris simulation. Dwarf multiples are defined as 2 or more dwarfs that have angular separations > 55'', projected separations r$_p < 150$ kpc and relative line-of-sight velocities $\Delta V_{\rm LOS} < 150$ km/s. While the mock catalogs predict a factor of 2 more isolated dwarfs than observed in SDSS, the mean number of observed companions per dwarf is $N_c \sim 0.04$, in good agreement with Illustris when accounting for SDSS sensitivity limits. Removing these limits in the mock catalogs predicts $N_c\sim 0.06$ for future surveys (LSST, DESI), which will be complete to M$_* = 2\times 10^8$ M$_\odot$. The 3D separations of mock dwarf multiples reveal a contamination fraction of $\sim$40% in observations from projection effects. Most isolated multiples are pairs; triples are rare and it is cosmologically improbable that bound groups of dwarfs with more than 3 members exist within the parameter range probed in this study. We find that $<$1% of LMC-analogs in the field have an SMC-analog companion. The fraction of dwarf "Major Pairs'' (stellar mass ratio $>$1:4) steadily increases with decreasing Primary stellar mass, whereas the cosmological "Major Merger rate'' (per Gyr) has the opposite behaviour. We conclude that cosmological simulations can be reliably used to constrain the fraction of dwarf mergers across cosmic time., Comment: 21 pages, 16 Figures, submitted to MNRAS (Referee stage ; Comments are welcome)

We present ALMA and ACA [CI]$^{3}P_{1}-^{3}P_{0}$ ([CI](1-0)) observations of NGC6240, which we combine with ALMA CO(2-1) and IRAM PdBI CO(1-0) data to study the physical properties of the massive molecular (H$_2$) outflow. We discover that the receding and approaching sides of the H$_2$ outflow, aligned east-west, exceed 10 kpc in their total extent. High resolution ($0.24"$) [CI](1-0) line images surprisingly reveal that the outflow emission peaks between the two AGNs, rather than on either of the two, and that it dominates the velocity field in this nuclear region. We combine the [CI](1-0) and CO(1-0) data to constrain the CO-to-H$_2$ conversion factor ($\alpha_{\rm CO}$) in the outflow, which is on average $2.1\pm1.2~\rm M_{\odot} (K~km~s^{-1}~pc^2)^{-1}$. We estimate that $60\pm20$ % of the total H$_2$ gas reservoir of NGC6240 is entrained in the outflow, for a resulting mass-loss rate of $\dot{M}_{\rm out}=2500\pm1200~M_{\odot}~yr^{-1}$ $\equiv 50\pm30$ SFR. This energetics rules out a solely star formation-driven wind, but the puzzling morphology challenges a classic radiative-mode AGN feedback scenario. For the quiescent gas we compute $\langle\alpha_{\rm CO}\rangle = 3.2\pm1.8~\rm M_{\odot} (K~km~s^{-1}~pc^2)^{-1}$, which is at least twice the value commonly employed for (U)LIRGs. We observe a tentative trend of increasing $r_{21}\equiv L^{\prime}_{\rm CO(2-1)}/L^{\prime}_{\rm CO(1-0)}$ ratios with velocity dispersion and measure $r_{21}>1$ in the outflow, whereas $r_{21}\simeq1$ in the quiescent gas. We propose that molecular outflows are the location of the warmer, strongly unbound phase that partially reduces the opacity of the CO lines in (U)LIRGs, hence driving down their global $\alpha_{\rm CO}$ and increasing their $r_{21}$ values., Comment: Accepted for publication in ApJ

We present the results from our Atacama Large Millimeter/submillimeter Array (ALMA) imaging observations of the CO (7-6), [CI] 370 um (hereafter [CI]) and [NII] 205 um (hereafter [NI]I) lines and their underlying continuum emission of BRI 1335-0417, an infrared bright quasar at z = 4.407. At the achieved resolutions of 1.1" to 1.2" (or 7.5 to 8.2 kpc), the continuum at 205 and 372 um (rest-frame), the CO (7-6), and the [CI] emissions are at best barely resolved whereas the [NII] emission is well resolved with an ALMA beam de-convolved major axis of 1.3" (+/- 0.3") or 9 (+/-2) kpc. As a warm dense gas tracer, the CO (7-6) emission shows a more compact spatial distribution and a significantly higher peak velocity dispersion than the other two lines that probe lower density gas, a picture favoring a merger-triggered star formation (SF) scenario over an orderly rotating SF disk. The CO (7-6) data also indicate a possible QSO-driven gas outflow that reaches a maximum line-of-sight velocity of 500 to 600 km/s. The far-infrared (FIR) dust temperature (T_dust) of 41.5 K from a gray-body fit to the continuum agrees well with the average T_dust inferred from various line luminosity ratios. The resulting L_CO(7-6)/L_FIR luminosity ratio is consistent with that of local luminous infrared galaxies powered predominantly by SF. The CO(7-6) luminosity-inferred SF rate is 5.1 (+/-1.5) x 10^3 M_solar/yr . The system has an effective star-forming region of 1.7 (+1.7/-0.8) kpc in diameter and a molecular gas reservoir of ~5 x 10^{11} M_solar., Comment: 14 pages, 7 figures; accepted for publication in the Astrophysical Journal

Discoveries of low mass galaxy pairs and groups are increasing. Studies indicate that dwarf galaxy pairs are gas rich in the field and exhibit elevated star formation rates, suggestive of interactions. Lacking are dynamical models of observed dwarf galaxy pairs to disentangle the physical processes regulating their baryon cycles. We present new optical data and the first detailed theoretical model of an observed tidal encounter between two isolated low mass galaxies, NGC 4490 & NGC 4485. This system is an isolated analog of the Magellanic Clouds and is surrounded by a ~50 kpc extended HI envelope. We use hybrid $N$-body and test-particle simulations along with a visualization interface $Identikit$ to simultaneously reproduce the observed present-day morphology and kinematics. Our results demonstrate how repeated encounters between two dwarf galaxies can "park" baryons at very large distances, without the aid of environmental effects. Our best match to the data is an 8:1 mass ratio encounter where a one-armed spiral is induced in the NGC 4490-analog, which we postulate explains the nature of diffuse starlight presented in the new optical data. We predict that the pair will fully merge in ~370 Myr, but that the extended tidal features will continue to evolve and return to the merged remnant over ~5 Gyr. This pre-processing of baryons will affect the efficiency of gas stripping if such dwarf pairs are accreted by a massive host. In contrast, in isolated environments this study demonstrates how dwarf-dwarf interactions can create a long-lived supply of gas to the merger remnant., Comment: 30 pages, 11 figures, 2 tables, minor edits to reflect published version

Mid-infrared molecular hydrogen (H$_2$) emission is a powerful cooling agent in galaxy mergers and in radio galaxies; it is a potential key tracer of gas evolution and energy dissipation associated with mergers, star formation, and accretion onto supermassive black holes. We detect mid-IR H$_2$ line emission in at least one rotational transition in 91\% of the 214 Luminous Infrared Galaxies (LIRGs) observed with Spitzer as part of the Great Observatories All-sky LIRG Survey (GOALS). We use H$_2$ excitation diagrams to estimate the range of masses and temperatures of warm molecular gas in these galaxies. We find that LIRGs in which the IR emission originates mostly from the Active Galactic Nuclei (AGN) have about 100K higher H$_2$ mass-averaged excitation temperatures than LIRGs in which the IR emission originates mostly from star formation. Between 10 and 15\% of LIRGs have H$_2$ emission lines that are sufficiently broad to be resolved or partially resolved by the high resolution modules of Spitzer's Infrared Spectrograph (IRS). Those sources tend to be mergers and contain AGN. This suggests that a significant fraction of the H$_2$ line emission is powered by AGN activity through X-rays, cosmic rays, and turbulence. We find a statistically significant correlation between the kinetic energy in the H$_2$ gas and the H$_2$ to IR luminosity ratio. The sources with the largest warm gas kinetic energies are mergers. We speculate that mergers increase the production of bulk in-flows leading to observable broad H$_2$ profiles and possibly denser environments., Comment: 31 pages, 12 figures, 7 tables, accepted for publication in the Astrophysical Journal

We present the first spatially and spectrally resolved image of the molecular outflow in the western nucleus of Arp\,220. The outflow, seen in HCN~(1--0) by ALMA, is compact and collimated, with extension $\lesssim$ 120\,pc. Bipolar morphology emerges along the minor axis of the disk, with redshifted and blueshifted components reaching maximum inclination-corrected velocity of $\sim\,\pm$\,840\,km\,s$^{-1}$. The outflow is also seen in CO and continuum emission, the latter implying that it carries significant dust. We estimate a total mass in the outflow of $\geqslant$\,10$^{6}$\,M$_{\odot}$, a dynamical time of $\sim$\,10$^{5}$\,yr, and mass outflow rates of $\geqslant55$\,M$_{\odot}$\,yr$^{-1}$ and $\geqslant\,15$\,M$_{\odot}$\,yr$^{-1}$ for the northern and southern lobes, respectively. Possible driving mechanisms include supernovae energy and momentum transfer, radiation pressure feedback, and a central AGN. The latter could explain the collimated morphology of the HCN outflow, however we need more complex theoretical models, including contribution from supernovae and AGN, to pinpoint the driving mechanism of this outflow., Comment: Accepted for publication in ApJL

We present the results from Atacama Large Millimeter/submillimeter Array (ALMA) imaging in the [NII] 205 micron fine-structure line (hereafter [NII]) and the underlying continuum of BRI 1202-0725, an interacting galaxy system at $z =$ 4.7, consisting of an optical QSO, a sub-millimeter galaxy (SMG) and two Lyman-$\alpha$ emitters (LAEs), all within $\sim$25 kpc of the QSO. We detect the QSO and SMG in both [NII] and continuum. At the $\sim$$1"$ (or 6.6 kpc) resolution, both QSO and SMG are resolved in [NII], with the de-convolved major axes of $\sim$9 and $\sim$14 kpc, respectively. In contrast, their continuum emissions are much more compact and unresolved even at an enhanced resolution of $\sim$$0.7"$. The ratio of the [NII] flux to the existing CO (7$-$6) flux is used to constrain the dust temperature ($T_{\rm dust}$) for a more accurate determination of the FIR luminosity $L_{\rm FIR}$. Our best estimated $T_{\rm dust}$ equals $43 (\pm 2)$ K for both galaxies (assuming an emissivity index $\beta = 1.8$). The resulting $L_{\rm CO(7-6)}/L_{\rm FIR}$ ratios are statistically consistent with that of local luminous infrared galaxies, confirming that $L_{\rm CO(7-6)}$ traces the star formation (SF) rate (SFR) in these galaxies. We estimate that the on-going SF of the QSO (SMG) has a SFR of 5.1 $(6.9) \times 10^3 M_{\odot}$ yr$^{-1}$ ($\pm$ 30%) assuming Chabrier initial mass function, takes place within a diameter (at half maximum) of 1.3 (1.5) kpc, and shall consume the existing 5 $(5) \times 10^{11} M_{\odot}$ of molecular gas in 10 $(7) \times 10^7$ years., Comment: 4 pages, 1 table, 4 figures; accepted for publication in the ApJ Letters

We present ALMA $^{12}$CO (J=1-0, 3-2 and 6-5), $^{13}$CO (J=1-0) and C$^{18}$O (J=1-0) observations of the local Ultra Luminous Infrared Galaxy, IRAS 13120-5453 (dubbed "The Yo-yo"). The morphologies of the three isotopic species differ, where $^{13}$CO shows a hole in emission towards the center. We measure integrated brightness temperature line ratios of $^{12}$CO/$^{13}$CO $\geq$ 60 (exceeding 200) and $^{13}$CO/C$^{18}$O $\leq$ 1 in the central region. Assuming optical thin emission, C$^{18}$O is more abundant than $^{13}$CO in several regions. The abundances within the central 500 pc are consistent with enrichment of the ISM via a young starburst ($<$7Myr), a top-heavy initial mass function or a combination of both., Comment: 7 pages, 3 figures, accepted for publication in ApJL

We describe a Herschel Space Observatory 194-671 micron spectroscopic survey of a sample of 121 local luminous infrared galaxies and report the fluxes of the CO $J$ to $J$-1 rotational transitions for $4 \leqslant J \leqslant 13$, the [NII] 205 um line, the [CI] lines at 609 and 370 um, as well as additional and usually fainter lines. The CO spectral line energy distributions (SLEDs) presented here are consistent with our earlier work, which was based on a smaller sample, that calls for two distinct molecular gas components in general: (i) a cold component, which emits CO lines primarily at $J \lesssim 4$ and likely represents the same gas phase traced by CO (1-0), and (ii) a warm component, which dominates over the mid-$J$ regime ($4 < J < 10$) and is intimately related to current star formation. We present evidence that the CO line emission associated with an active galactic nucleus is significant only at $J > 10$. The flux ratios of the two [CI] lines imply modest excitation temperatures of 15 to 30 K; the [CI] 370 um line scales more linearly in flux with CO (4-3) than with CO (7-6). These findings suggest that the [CI] emission is predominately associated with the gas component defined in (i) above. Our analysis of the stacked spectra in different far-infrared (FIR) color bins reveals an evolution of the SLED of the rotational transitions of water vapor as a function of the FIR color in a direction consistent with infrared photon pumping., Comment: 244 pages, 23 figures, 7 tables, accepted for publication in the Astrophysical Journal Supplement Series