We present Hubble Space Telescope 1.4-1.6 micron images of the hosts of ten extremely red quasars (ERQs) and six type 2 quasar candidates at z=2-3. ERQs, whose bolometric luminosities range between 10^47 and 10^48 erg/sec, show spectroscopic signs of powerful ionized winds, whereas type 2 quasar candidates are less luminous and show only mild outflows. After performing careful subtraction of the quasar light, we clearly detect almost all host galaxies. The median rest-frame B-band luminosity of the ERQ hosts in our sample is 10^11.2 L_Sun, or 4 L* at this redshift. Two of the ten hosts of ERQs are in ongoing mergers. The hosts of the type 2 quasar candidates are 0.6 dex less luminous, with 2/6 in likely ongoing mergers. Intriguingly, despite some signs of interaction and presence of low-mass companions, our objects do not show nearly as much major merger activity as do high-redshift radio-loud galaxies and quasars. In the absence of an overt connection to major ongoing gas-rich merger activity, our observations are consistent with a model in which the near-Eddington accretion and strong feedback of ERQs are associated with relatively late stages of mergers resulting in early-type remnants. These results are in some tension with theoretical expectations of galaxy formation models, in which rapid black hole growth occurs within a short time of a major merger. Type 2 quasar candidates are less luminous, so they may instead be powered by internal galactic processes.

Extended narrow-line regions (NLRs) around active galactic nuclei (AGN) are shaped by the distribution of gas in the host galaxy and by the geometry of the circumnuclear obscuration, and thus they can be used to test the AGN unification model. In this work, we quantify the morphologies of the narrow-line regions in 308 nearby AGNs ($z=0-0.14$, \lbol $\sim 10^{42.4-44.1}$ \erg{}) from the MaNGA survey. Based on the narrow-line region maps, we find that a large fraction (81\%) of these AGN have bi-conical NLR morphology. The distribution of their measured opening angles suggests that the intrinsic opening angles of the ionization cones has a mean value of 85--98$^\circ$ with a finite spread of 39-44$^\circ$ (1-$\sigma$). Our inferred opening angle distribution implies a number ratio of type I to type II AGN of 1:1.6--2.3, consistent with other measurements of the type I / type II ratio at low AGN luminosities. Combining these measurements with the WISE photometry data, we find that redder mid-IR color (lower effective temperature of dust) corresponds to stronger and narrower photo-ionized bicones. This relation is in agreement with the unification model that suggests that the bi-conical narrow-line regions are shaped by a toroidal dusty structure within a few pc from the AGN. Furthermore, we find a significant alignment between the minor axis of host galaxy disks and AGN ionization cones. Together, these findings suggest that obscuration on both circumnuclear ($\sim $pc) and galactic ($\sim$ kpc) scales are important in shaping and orienting the AGN narrow-line regions.

We present the source separation framework SCARLET for multi-band images, which is based on a generalization of the Non-negative Matrix Factorization to alternative and several simultaneous constraints. Our approach describes the observed scene as a mixture of components with compact spatial support and uniform spectra over their support. We present the algorithm to perform the matrix factorization and introduce constraints that are useful for optical images of stars and distinct stellar populations in galaxies, in particular symmetry and monotonicity with respect to the source peak position. We also derive the treatment of correlated noise and convolutions with band-dependent point spread functions, rendering our approach applicable to coadded images observed under variable seeing conditions. SCARLET thus yields a PSF-matched photometry measurement with an optimally chosen weight function given by the mean morphology in all available bands. We demonstrate the performance of SCARLET for deblending crowded extragalactic scenes and on an AGN jet -- host galaxy separation problem in deep 5-band imaging from the Hyper Suprime-Cam Stategic Survey Program. Using simulations with prominent crowding we show that SCARLET yields superior results to the HSC-SDSS deblender for the recovery of total fluxes, colors, and morphologies. Due to its non-parametric nature, a conceptual limitation of SCARLET is its sensitivity to undetected sources or multiple stellar population within detected sources, but an iterative strategy that adds components at the location of significant residuals appears promising. The code is implemented in Python with C++ extensions and is available at https://github.com/fred3m/scarlet

Narrow-line regions excited by active galactic nuclei (AGN) are important for studying AGN photoionization and feedback. Their strong [O III] lines can be detected with broadband images, allowing morphological studies of these systems with large-area imaging surveys. We develop a new technique to reconstruct the [O III] images using the Subaru Hyper Suprime-Cam (HSC) Survey aided with spectra from the Sloan Digital Sky Survey (SDSS). The technique involves a careful subtraction of the galactic continuum to isolate emission from the [O III]$\lambda$5007 and [O III]$\lambda$4959 lines. Compared to traditional targeted observations, this technique is more efficient at covering larger samples with less dedicated observational resources. We apply this technique to an SDSS spectroscopically selected sample of 300 obscured AGN at redshifts 0.1 - 0.7, uncovering extended emission-line region candidates with sizes up to tens of kpc. With the largest sample of uniformly derived narrow-line region sizes, we revisit the narrow-line region size-luminosity relation. The area and radii of the [O III] emission-line regions are strongly correlated with the AGN luminosity inferred from the mid-infrared (15 $\mu$m rest-frame) with a power-law slope of $0.62^{+0.05}_{-0.06} \pm 0.10$ (statistical and systemic errors), consistent with previous spectroscopic findings. We discuss the implications for the physics of AGN emission-line region and future applications of this technique, which should be useful for current and next-generation imaging surveys to study AGN photoionization and feedback with large statistical samples.

We report the discovery of an infrared (IR)-bright dust-obscured galaxy (DOG) that shows a strong ionized-gas outflow but no significant molecular gas outflow. Based on detail analysis of their optical spectra, we found some peculiar IR-bright DOGs that show strong ionized-gas outflow ([OIII]$\lambda$5007) from the central active galactic nucleus (AGN). For one of these DOGs (WISE J102905.90+050132.4) at $z_{\rm spec} = 0.493$, we performed follow-up observations using ALMA to investigate their CO molecular gas properties. As a result, we successfully detected $^{12}$CO($J$=2-1) and $^{12}$CO($J$=4-3) lines, and the continuum of this DOG. The intensity-weighted velocity map of both lines shows a gradient, and the line profile of those CO lines is well-fitted by a single narrow Gaussian, meaning that this DOG has no sign of strong molecular gas outflow. The IR luminosity of this object is $\log\,(L_{\rm IR}/L_{\odot})$ = 12.40 that is classified as ultraluminous IR galaxy (ULIRG). We found that (i) the stellar mass and star-formation rate relation and (ii) the CO luminosity and far-IR luminosity relation are consistent with those of typical ULIRGs at similar redshifts. These results indicate that the molecular gas properties of this DOG are normal despite that its optical spectrum showing a powerful AGN outflow. We conclude that a powerful ionized-gas outflow caused by the AGN does not necessarily affect the cold interstellar medium in the host galaxy at least for this DOG.

We present the ionized gas properties of infrared (IR)-bright dust-obscured galaxies (DOGs) that show an extreme optical/IR color, $(i - [22])_{\rm AB} > 7.0$, selected with the Sloan Digital Sky Survey (SDSS) and Wide-field Infrared Survey Explorer (WISE). For 36 IR-bright DOGs that show [OIII]$\lambda$5007 emission in the SDSS spectra, we performed a detailed spectral analysis to investigate their ionized gas properties. In particular, we measured the velocity offset (the velocity with respect to the systemic velocity measured from the stellar absorption lines) and the velocity dispersion of the [OIII] line. We found that the derived velocity offset and dispersion of most IR-bright DOGs are larger than those of Seyfert 2 galaxies (Sy2s) at $z < 0.3$, meaning that the IR-bright DOGs show relatively strong outflows compared to Sy2s. This can be explained by the difference of IR luminosity contributed from active galactic nucleus, $L_{\rm IR}$ (AGN), because we found that (i) $L_{\rm IR}$ (AGN) correlates with the velocity offset and dispersion of [OIII] and (ii) our IR-bright DOGs sample has larger $L_{\rm IR}$ (AGN) than Sy2s. Nevertheless, the fact that about 75% IR-bright DOGs have a large ($>$ 300 km s$^{-1}$) velocity dispersion, which is a larger fraction compared to other AGN populations, suggests that IR-bright DOGs are good laboratories to investigate AGN feedback. The velocity offset and dispersion of [OIII] and [NeIII]$\lambda$3869 are larger than those of [OII]$\lambda$3727, which indicates that the highly ionized gas tends to show more stronger outflows.

We performed near-diffraction-limited (~0.4 FWHM) N-band imaging of one of the nearest Active Galactic Nucleus (AGN) in M51 with 8.2m Subaru telescope to study the nuclear structure and spectral energy distribution (SED) at 8-13 um. We found that the nucleus is composed of an unresolved core (at ~13 pc resolution, or intrinsic size corrected for the instrumental effect of <6 pc) and an extended halo (at a few tens pc scale), and each of their SEDs is almost flat. We examined the SED by comparing with the archival Spitzer IRS spectrum processed to mimic our chopping observation of the nucleus, and the published radiative-transfer model SEDs of the AGN clumpy dusty torus. The halo SED is likely due to circumnuclear star formation showing little Polycyclic Aromatic Hydrocarbon (PAH) emission due to the AGN. The core SED is likely dominated by the AGN because of the following two reasons. Firstly, the clumpy torus model SEDs can reproduce the red mid-infrared continuum with apparently moderate silicate 9.7 um absorption. Secondly, the core 12 um luminosity and the absorption-corrected X-ray luminosity at 2-10 keV in the literature follow the mid-infrared-X-ray luminosity correlation known for the nearby AGNs including the Compton-thick ones.

To study the impact of active galactic nuclei (AGN) feedback on the galactic ISM, we present Magellan long-slit spectroscopy of 12 luminous nearby type 2 AGN (L_bol~10^{45.0-46.5} erg/s, z~0.1). These objects are selected from a parent sample of spectroscopically identified AGN to have high [OIII]{\lambda}5007 and WISE mid-IR luminosities and extended emission in the SDSS r-band images, suggesting the presence of extended [OIII]{\lambda}5007 emission. We find spatially resolved [OIII] emission (2-35 kpc from the nucleus) in 8 out of 12 of these objects. Combined with samples of higher luminosity type 2 AGN, we confirm that the size of the narrow-line region (R_NLR) scales with the mid-IR luminosity until the relation flattens at ~10 kpc. Nine out of 12 objects in our sample have regions with broad [OIII] linewidths (w_80>600 km/s), indicating outflows. We define these regions as the kinematically-disturbed region (KDR). The size of the KDR (R_KDR) is typically smaller than R_NLR by few kpc but also correlates strongly with the AGN mid-IR luminosity. Given the unknown density in the gas, we derive a wide range in the energy efficiency {\eta}=dot{E}/L_bol=0.01%-30%. We find no evidence for an AGN luminosity threshold below which outflows are not launched. To explain the sizes, velocity profiles, and high occurrence rates of the outflows in the most luminous AGN, we propose a scenario in which energy-conserving outflows are driven by AGN episodes with ~10^8-year durations. Within each episode the AGN flickers on shorter timescales, with a cadence of ~10^6 year active phases separated by ~10^7 years.

We present Atacama Large Millimeter/Submillimeter Array (ALMA) CO (1-0) and CO (3-2) observations of SDSS J135646.10+102609.0, an obscured quasar and ultra-luminous infrared galaxy (ULIRG) with two merging nuclei and a known 20-kpc-scale ionized outflow. The total molecular gas mass is M_{mol} ~ 9^{+19}_{-6} x 10^8 Msun, mostly distributed in a compact rotating disk at the primary nucleus (M_{mol} ~ 3 x 10^8 Msun) and an extended tidal arm (M_{mol} ~ 5 x 10^8 Msun). The tidal arm is one of the most massive molecular tidal features known; we suggest that it is due to the lower chance of shock dissociation in this elliptical/disk galaxy merger. In the spatially resolved CO (3-2) data, we find a compact (r ~ 0.3 kpc) high velocity (v ~ 500 km/s) red-shifted feature in addition to the rotation at the N nucleus. We propose a molecular outflow as the most likely explanation for the high velocity gas. The outflowing mass of M_{mol} ~ 7 x 10^7 Msun and the short dynamical time of t_{dyn} ~ 0.6 Myr yield a very high outflow rate of \dot{M}_{mol} ~ 350 Msun/yr and can deplete the gas in a million years. We find a low star formation rate (< 16 Msun/yr from the molecular content and < 21 Msun/yr from the far-infrared spectral energy distribution decomposition) that is inadequate to supply the kinetic luminosity of the outflow (\dot{E} ~ 3 x 10^43 erg/s). Therefore, the active galactic nucleus, with a bolometric luminosity of 10^46 erg/s, likely powers the outflow. The momentum boost rate of the outflow (\dot{p}/(Lbol/c) ~ 3) is lower than typical molecular outflows associated with AGN, which may be related to its compactness. The molecular and ionized outflows are likely two distinct bursts induced by episodic AGN activity that varies on a time scale of 10^7 yr.

We present observations of extended, 20-kpc scale soft X-ray gas around a luminous obscured quasar hosted by an ultra-luminous infrared galaxy caught in the midst of a major merger. The extended X-ray emission is well fit as a thermal gas with a temperature of kT ~ 280 eV and a luminosity of L_X ~ 10^42 erg/s and is spatially coincident with a known ionized gas outflow. Based on the X-ray luminosity, a factor of ~10 fainter than the [OIII] emission, we conclude that the X-ray emission is either dominated by photoionization, or by shocked emission from cloud surfaces in a hot quasar-driven wind.

We present the Green Bank 100m Telescope (GBT) mapping observations of CS 1-0, and the Submillimeter Array (SMA) 157-pointings mosaic of the 0.86 mm dust continuum emission as well as several warm and dense gas tracers, in the central ~20 pc area in Galactic Center. The unprecedentedly large field-of-view and the high angular resolution of our SMA dust image allow the identification of abundant 0.1-0.2 pc scale dense gas clumps. We found that in the Galactic Center, the Class I methanol masers are excellently correlated with the dense gas clumps. However, on the ~0.1 pc scale, these dense gas clumps still have a extremely large linewidth (FWHM~10-20 km/s). Simple calculations suggest that the identified clumps can be possibly the pressurized gas reservoir feeding the formation of 1-10 solar-mass stars. These gas clumps may be the most promising candidates for ALMA to resolve the high-mass star-formation in the Galactic center.

A fast and practical method is described for calculating the microlensed flux variance of an arbitrary source by uncorrelated stars. The required inputs are the mean convergence and shear due to the smoothed potential of the lensing galaxy, the stellar mass function, and the absolute square of the Fourier transform of the surface brightness in the source plane. The mathematical approach follows previous authors but has been generalized, streamlined, and implemented in publicly available code. Examples of its application are given for Dexter and Agol's inhomogeneous-disk models as well as the usual gaussian sources. Since the quantity calculated is a second moment of the magnification, it is only logarithmically sensitive to the sizes of very compact sources. However, for the inferred sizes of actual QSOs, it has some discriminatory power and may lend itself to simple statistical tests. At the very least, it should be useful for testing the convergence of microlensing simulations.

Black hole - galaxy scaling relations provide information about the coevolution of supermassive black holes and their host galaxies. We compare the black hole mass - circular velocity (MBH - Vc) relation with the black hole mass - bulge stellar velocity dispersion (MBH - sigma) relation, to see whether the scaling relations can passively emerge from a large number of mergers, or require a physical mechanism, such as feedback from an active nucleus. We present VLA H I observations of five galaxies, including three water megamaser galaxies, to measure the circular velocity. Using twenty-two galaxies with dynamical MBH measurements and Vc measurements extending to large radius, our best-fit MBH - Vc relation, log MBH = alpha + beta log(Vc /200 km s^-1), yields alpha = 7.43+/-0.13, beta = 3.68+1.23/-1.20, and intrinsic scatter epsilon_int = 0.51+0.11/-0.09. The intrinsic scatter may well be higher than 0.51, as we take great care to ascribe conservatively large observational errors. We find comparable scatter in the MBH - sigma relations, epsilon_int = 0.48+0.10/-0.08, while pure merging scenarios would likely result in a tighter scaling with the dark halo (as traced by Vc) than baryonic (sigma) properties. Instead, feedback from the active nucleus may act on bulge scales to tighten the MBH - sigma relation with respect to the MBH - Vc relation, as observed.

We examine the alignment between H_2O megamaser disks on sub-pc scales with circumnuclear disks and bars on <500 pc scales observed with HST/WFC3. The HST imaging reveals young stars, indicating the presence of gas. The megamaser disks are not well aligned with the circumnuclear bars or disks as traced by stars in the HST images. We speculate on the implications of the observed misalignments for fueling supermassive black holes in gas-rich spiral galaxies. In contrast, we find a strong preference for the rotation axes of the megamaser disks to align with radio continuum jets observed on >50 pc scales, in those galaxies for which radio continuum detections are available. Sub-arcsecond observations of molecular gas with ALMA will enable a more complete understanding of the interplay between circumnuclear structures.

We present the Submillimeter Array (SMA) 157-pointing mosaic in 0.86 mm dust continuum emission with 5$".1\times4"$.2 angular resolution, and the NRAO Green Bank 100m Telescope (GBT) observations of the CS/C$^{34}$S/$^{13}$CS 1-0 and SiO 1-0 emission with $\le20"$$\times18"$ angular resolution. The dust continuum image marginally resolves at least several tens of 10-10$^{2}$ $M_{\odot}$ dense clumps in the 5$'$ field including the circumnuclear disk (CND) and the exterior gas streamers. There is very good agreement between the high resolution dust continuum map of the CND and all previous molecular line observations. As the dust emission is the most reliable optically thin tracer of the mass, free from most chemical and excitation effects, we demonstrate the reality of the abundant localized structures within the CND, and their connection to external gas structures. From the spectral line data, the velocity dispersions of the dense clumps and their parent molecular clouds are $\sim$10-20 times higher than their virial velocity dispersions. This supports the idea that the CND and its immediate environment may not be stationary or stable structures. Some of the dense gas clumps are associated with 22 GHz water masers and 36.2 GHz and 44.1 GHz CH$_{3}$OH masers. However, we do not find clumps which are bound by the gravity of the enclosed molecular gas. Hence, the CH$_{3}$OH or H$_{2}$O maser emission may be due to strong (proto)stellar feedback, which may be dispersing some of the gas clumps.

The supermassive black hole (SMBH), Sgr A*, at the Galactic Center is surrounded by a molecular circumnuclear disk (CND) lying between 1.5-4 pc radii. The irregular and clumpy structures of the CND, suggest dynamical evolution and episodic feeding of gas towards the central SMBH. New sensitive data from the SMA and GBT, reveal several >5-10 pc scale molecular arms, which either directly connect to the CND, or may penetrate inside the CND. The CND appears to be the convergence of the innermost parts of largescale gas streamers, which are responding to the central gravitational potential well. Rather than being a quasi-stationary structure, the CND may be dynamically evolving, incorporating inflow via streamers, and feeding gas towards the center.

The importance of the far-infrared (FIR) mapping is demonstrated for a face-on spiral galaxy, M81, by analyzing its imaging data at 65, 90, and 140 {\mu}m taken by AKARI. Basic products are the dust temperature map, the dust optical depth map, and the colour-colour diagram. The main features are as follows. (i) The dust temperature derived from the total fluxes at 90 {\mu}m and 140 {\mu}m reflects the relatively low temperatures seen in the interarm and spiral arms excluding the warm spots, rather than the high temperatures in warm spots and the centre. This indicates that the total FIR luminosity is dominated by the dust heated by the general interstellar radiation field. (ii) The galaxy is more extended at 140 {\mu}m than at the other shorter wavelengths, which reflects the radial dust temperature gradient. (iii) The dust optical depth derived from the FIR mapping is broadly consistent with that estimated from the FIR-to-ultraviolet luminosity ratio. (iv) The FIR colour-colour diagramis useful to identify a 'contamination' of warm dust. The existence of small-scale warm star-forming regions is supported in the bright spots along the spiral arms. This contamination also leads to an underestimate of dust optical depth (or dust column density).