A series of gravitational instabilities in a circumnuclear gas disk (CND) are required to trigger gas transport to a central supermassive black hole (SMBH) and ignite Active Galactic Nuclei (AGNs). A test of this scenario is to investigate whether an enhanced molecular gas mass surface density ($\Sigma_{\rm mol}$) is found in the CND-scale of quasars relative to a comparison sample of inactive galaxies. Here we performed sub-kpc resolution CO(2-1) observations with ALMA of four low-redshift ($z \sim 0.06$), luminous ($\sim 10^{45}$ erg s$^{-1}$) quasars with each matched to a different star-forming galaxy, having similar redshift, stellar mass, and star-formation rate. We detected CO(2-1) emission from all quasars, which show diverse morphologies. Contrary to expectations, $\Sigma_{\rm mol}$ of the quasar sample, computed from the CO(2-1) luminosity, tends to be smaller than the comparison sample at $r < 500$ pc; there is no systematic enhancement of $\Sigma_{\rm mol}$ in our quasars. We discuss four possible scenarios that would explain the lower molecular gas content (or CO(2-1) luminosity as an actual observable) at the CND-scale of quasars, i.e., AGN-driven outflows, gas-rich minor mergers, time-delay between the onsets of a starburst-phase and a quasar-phase, and X-ray-dominated region (XDR) effects on the gas chemical abundance and excitation. While not extensively discussed in the literature, XDR effects can have an impact on molecular mass measurements particularly in the vicinity of luminous quasar nuclei; therefore higher resolution molecular gas observations, which are now viable using ALMA, need to be considered.

We used the Atacama Large Millimeter/submillimeter Array (ALMA) to map $^{12}$CO($J$ = 1-0), $^{12}$CO($J$ = 2-1), $^{12}$CO($J$ = 3-2), $^{13}$CO($J$ = 2-1), and [CI]($^3P_1$-$^3P_0$) emission lines around the type 1 active galactic nucleus (AGN) of NGC 7469 ($z = 0.0164$) at $\sim 100$ pc resolutions. The CO lines are bright both in the circumnuclear disk (central $\sim 300$ pc) and the surrounding starburst (SB) ring ($\sim 1$ kpc diameter), with two bright peaks on either side of the AGN. By contrast, the [CI]($^3P_1$-$^3P_0$) line is strongly peaked on the AGN. Consequently, the brightness temperature ratio of [CI]($^3P_1$-$^3P_0$) to $^{13}$CO(2-1) is $\sim 20$ at the AGN, as compared to $\sim 2$ in the SB ring. Our local thermodynamic equilibrium (LTE) and non-LTE models indicate that the enhanced line ratios (or CI enhancement) are due to an elevated C$^0$/CO abundance ratio ($\sim 3-10$) and temperature ($\sim 100-500$ K) around the AGN as compared to the SB ring (abundance ratio $\sim 1$, temperature $\lesssim 100$ K), which accords with the picture of the X-ray-dominated Region (XDR). Based on dynamical modelings, we also provide CO(1-0)-to- and [CI]($^3P_1$-$^3P_0$)-to-molecular mass conversion factors at the central $\sim 100$ pc of this AGN as $\alpha_{\rm CO} = 4.1$ and $\alpha_{\rm CI} = 4.4~M_\odot$ (K km s$^{-1}$ pc$^2$)$^{-1}$, respectively. Our results suggest that the CI enhancement is potentially a good marker of AGNs that could be used in a new submillimeter diagnostic method toward dusty environments.

We report spatial distributions of the Fe-K$\alpha$ line at 6.4 keV and the CO($J$ = 2--1) line at 230.538 GHz in NGC 2110, which are respectively revealed by $Chandra$ and ALMA at $\approx$ 0.5 arcsec. A $Chandra$ 6.2--6.5 keV-to-3.0--6.0 keV image suggests that the Fe-K$\alpha$ emission extends preferentially in a northwest-to-southeast direction out to $\sim$ 3 arcsec, or 500 pc, on each side. Spatially-resolved spectral analyses support this by finding significant Fe-K$\alpha$ emission lines only in northwest and southeast regions. Moreover, their equivalent widths are found $\sim$ 1.5 keV, indicative for the fluorescence by nuclear X-ray irradiation as the physical origin. By contrast, CO($J$ = 2--1) emission is weak therein. For quantitative discussion, we derive ionization parameters by following an X-ray dominated region (XDR) model. We then find them high enough to interpret the weakness as the result of X-ray dissociation of CO and/or H$_2$. Another possibility also remains that CO molecules follow a super-thermal distribution, resulting in brighter emission in higher-$J$ lines. Further follow-up observations are encouraged to draw a conclusion on what predominantly changes the inter-stellar matter properties, and whether the X-ray irradiation eventually affects the surrounding star formation as an AGN feedback.

We report the results of the Australia Telescope Compact Array (ATCA) 15 mm observation of the Phoenix galaxy cluster possessing an extreme star-burst brightest cluster galaxy (BCG) at the cluster center. We spatially resolved radio emission around the BCG, and found diffuse bipolar and bar-shape structures extending from the active galactic nucleus (AGN) of the BCG. They are likely radio jets/lobes, whose sizes are ~10-20 kpc and locations are aligned with X-ray cavities. If we assume that the radio jets/lobes expand with the sound velocity, their ages are estimated to be ~10 Myr. We also found compact radio emissions near the center and suggest that they are another young bipolar jets with ~1 Myr of age. Moreover, we found extended radio emission surrounding the AGN and discussed the possibility that the component is a product of the cooling flow, by considering synchrotron radiation partially absorbed by molecular clumps, free-free emission from the warm ionized gas, and the spinning dust emission from dusty circum-galactic medium.

Gravitational lensing sometimes dominates the observed properties of apparently very bright objects. We present morphological properties in the high-resolution (FWHM $\sim0.''15$) Atacama Large Millimeter/submillimeter Array (ALMA) 1-mm map for an ultra-luminous quasar (QSO) at $z=6.30$, SDSS J010013.02+280225.8 (hereafter J0100+2802), whose black hole mass $M_{\rm BH}$ is the most massive ($\sim$ 1.2$\times10^{10}M_{\odot}$) at $z>6$ ever known. We find that the continuum emission of J0100+2802 is resolved into a quadruple system within a radius of $0.''2$, which can be interpreted as either multiple dusty star-forming regions in the host galaxy or multiple images due to strong gravitational lensing. The Mg$_{\rm II}$ absorption and the potential Ly$\alpha$ line features have been identified at $z=2.33$ in the near-infrared spectroscopy towards J0100+2802, and a simple mass model fitting well reproduces the positions and flux densities of the quadruple system, both of which are consistent with the latter interpretation. Although a high-resolution map taken in the Advanced Camera for Survey (ACS) on board Hubble Space Telescope (HST) shows a morphology with an apparently single component, in our fiducial lens mass model it can simply be explained by a $\sim50$ pc scale offset between the ALMA and HST emission regions. In this case, the magnification factor for the observed HST emission is obtained to $\sim450$, reducing the intrinsic $M_{\rm BH}$ estimate to even below $10^{9}\,M_{\odot}$. The confirmation or the rejection of the gravitational lensing scenario is important for our understanding of the super-massive black holes in the early Universe.

We present ALMA [CII] line and far-infrared (FIR) continuum observations of three $z > 6$ low-luminosity quasars ($M_{\rm 1450} > -25$) discovered by our Subaru Hyper Suprime-Cam (HSC) survey. The [CII] line was detected in all three targets with luminosities of $(2.4 - 9.5) \times 10^8~L_\odot$, about one order of magnitude smaller than optically luminous ($M_{\rm 1450} \lesssim -25$) quasars. The FIR continuum luminosities range from $< 9 \times 10^{10}~L_\odot$ (3$\sigma$ limit) to $\sim 2 \times 10^{12}~L_\odot$, indicating a wide range in star formation rates in these galaxies. Most of the HSC quasars studied thus far show [CII]/FIR luminosity ratios similar to local star-forming galaxies. Using the [CII]-based dynamical mass ($M_{\rm dyn}$) as a surrogate for bulge stellar mass ($M_{\rm bulge}$), we find that a significant fraction of low-luminosity quasars are located on or even below the local $M_{\rm BH} - M_{\rm bulge}$ relation, particularly at the massive end of the galaxy mass distribution. In contrast, previous studies of optically luminous quasars have found that black holes are overmassive relative to the local relation. Given the low luminosities of our targets, we are exploring the nature of the early co-evolution of supermassive black holes and their hosts in a less biased way. Almost all of the quasars presented in this work are growing their black hole mass at much higher pace at $z \sim 6$ than the parallel growth model, in which supermassive black holes and their hosts grow simultaneously to match the local $M_{\rm BH} - M_{\rm bulge}$ relation at all redshifts. As the low-luminosity quasars appear to realize the local co-evolutionary relation even at $z \sim 6$, they should have experienced vigorous starbursts prior to the currently observed quasar phase to catch up with the relation.

We present a study of the gas kinematics of star-forming galaxies associated with protocluster 4C 23.56 at $z=2.49$ using $0''.4$ resolution CO (4-3) data taken with ALMA. Eleven H$\alpha$ emitters (HAEs) are detected in CO (4-3), including six HAEs that were previously detected in CO (3-2) at a coarser angular resolution. The detections in both CO lines are broadly consistent in the line widths and the redshifts, confirming both detections. With an increase in the number of spectroscopic redshifts, we confirm that the protocluster is composed of two merging groups with a total halo mass of $\log{(M_{\rm cl}/M_{\odot})} =13.4-13.6$, suggesting that the protocluster would evolve into a Virgo-like cluster ($>10^{14} M_{\odot}$). We compare the CO line widths and the CO luminosities with galaxies in other (proto)clusters ($n_{\rm gal}=91$) and general fields ($n_{\rm gal}=80$) from other studies. The 4C23.56 protocluster galaxies have CO line widths and luminosities comparable to other protocluster galaxies on average. On the other hand, the CO line widths are on average broader by $\approx50\%$ compared to field galaxies, while the median CO luminosities are similar. The broader line widths can be attributed to both effects of unresolved gas-rich mergers and/or compact gas distribution, which is supported by our limited but decent angular resolution observations and the size estimate of three galaxies. Based on these results, we argue that gas-rich mergers may play a role in the retention of the specific angular momentum to a value similar to that of field populations during cluster assembly, though we need to verify this with a larger number of samples.

We report the discovery of 28 quasars and 7 luminous galaxies at 5.7 $\le$ z $\le$ 7.0. This is the tenth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The total number of spectroscopically identified objects in SHELLQs has now grown to 93 high-z quasars, 31 high-z luminous galaxies, 16 [O III] emitters at z ~ 0.8, and 65 Galactic cool dwarfs (low-mass stars and brown dwarfs). These objects were found over 900 deg2, surveyed by HSC between 2014 March and 2018 January. The full quasar sample includes 18 objects with very strong and narrow Ly alpha emission, whose stacked spectrum is clearly different from that of other quasars or galaxies. While the stacked spectrum shows N V 1240 emission and resembles that of lower-z narrow-line quasars, the small Ly alpha width may suggest a significant contribution from the host galaxies. Thus these objects may be composites of quasars and star-forming galaxies.

We present deep near-infrared spectroscopy of six quasars at 6.1<z<6.7 with VLT/X-Shooter and Gemini-N/GNIRS. Our objects, originally discovered through a wide-field optical survey with the Hyper Suprime-Cam (HSC) Subaru Strategic Program (HSC-SSP), have the lowest luminosities (-25.5< M1450<-23.1 mag) of the z>5.8 quasars with measured black hole masses. From single-epoch mass measurements based on MgII2798, we find a wide range in black hole masses, from M_BH=10^7.6 to 10^9.3 Msun. The Eddington ratios L_bol/L_Edd range from 0.16 to 1.1, but the majority of the HSC quasars are powered by M_BH=10^9 Msun supermassive black holes (SMBHs) accreting at sub-Eddington rates. The Eddington ratio distribution of the HSC quasars is inclined to lower accretion rates than those of Willott et al. (2010a), who measured the black hole masses for similarly faint z=6 quasars. This suggests that the global Eddington ratio distribution is wider than has previously been thought. The presence of M_BH=10^9 Msun SMBHs at z=6 cannot be explained with constant sub-Eddington accretion from stellar remnant seed black holes. Therefore, we may be witnessing the first buildup of the most massive black holes in the first billion years of the universe, the accretion activity of which is transforming from active growth to a quiescent phase. Measurements of a larger complete sample of z>6 low-luminosity quasars, as well as deeper observations with future facilities will enable us to better understand the early SMBH growth in the reionization epoch.

We explore X-ray evidence for the presence of active galactic nuclei (AGNs) in the two most actively star-forming Green Pea galaxies (GPs), SDSS J0749+3337 and SDSS J0822+2241, which have star-formation rates (SFRs) of $123~M_\odot$ yr$^{-1}$ and $78~M_\odot$ yr$^{-1}$, respectively. The GPs have red mid-infrared (MIR) spectral energy distributions and higher 22 $\mu$m luminosities than expected from a proxy of the SFR (H$\alpha$ luminosity), consistent with hosting AGNs with 2-10 keV luminosities of $\sim10^{44}$ erg s$^{-1}$. We thus obtain and analyze the first hard ($>$ 10 keV) X-ray data observed with NuSTAR and archival XMM-Newton data below 10 keV. From the NuSTAR $\approx$20 ksec data, however, we find no significant hard X-ray emission. By contrast, soft X-ray emission with 0.5--8 keV luminosities of $\approx10^{42}$ erg s$^{-1}$ is significantly detected in both targets, which can be explained only by star formation (SF). A possible reason for the lack of clear evidence is that a putative AGN torus absorbs most of the X-ray emission. Applying a smooth-density AGN torus model, we determine minimum hydrogen column densities along the equatorial plane ($N_{\rm H}^{\rm eq}$) consistent with the non-detection. The results indicate $N_{\rm H}^{\rm eq} \gtrsim 2\times10^{24}$ cm$^{-2}$ for SDSS J0749+3337 and $N_{\rm H}^{\rm eq} \gtrsim 5\times10^{24}$ cm$^{-2}$ for SDSS J0822+2241. Therefore, the GPs may host such heavily obscured AGNs. Otherwise, no AGN exists and the MIR emission is ascribed to SF. Active SF in low-mass galaxies is indeed suggested to reproduce red MIR colors. This would imply that diagnostics based on MIR photometry data alone may misidentify such galaxies as AGNs.

We present ALMA observations of the CO(2-1) and CO(3-2) molecular gas transitions and associated (sub)-mm continua of the nearby Seyfert 1.5 galaxy NGC3227 with angular resolutions 0.085-0.21" (7-15pc). On large scales the cold molecular gas shows circular motions as well as streaming motions on scales of a few hundred parsecs associated with a large scale bar. We fitted the nuclear ALMA 1.3mm emission with an unresolved component and an extended component. The 850$\mu$m emission shows at least two extended components, one along the major axis of the nuclear disk and the other along the axis of the ionization cone. The molecular gas in the central region (1" ~73pc) shows several CO clumps with complex kinematics which appears to be dominated by non-circular motions. While we cannot demonstrate conclusively the presence of a warped nuclear disk, we also detected non-circular motions along the kinematic minor axis. They reach line-of-sight velocities of v-vsys =150-200km/s. Assuming that the radial motions are in the plane of the galaxy, then we interpret them as a nuclear molecular outflow due to molecular gas in the host galaxy being entrained by the AGN wind. We derive molecular outflow rates of $5\,M_\odot\,{\rm yr}^{-1}$ and $0.6\,M_\odot\,{\rm yr}^{-1}$ at projected distances of up to 30pc to the northeast and southwest of the AGN, respectively. At the AGN location we estimate a mass in molecular gas of $5\times 10^{5}\,M_\odot$ and an average column density $N({\rm H}_2) = 2-3\times 10^{23}\,{\rm cm}^{-2}$ in the inner 15pc. The nuclear molecular gas and sub-mm continuum emission of NGC3227 do not resemble the classical compact torus. Rather, these emissions extend for several tens of parsecs and appear connected with the circumnuclear ring in the host galaxy disk, as found in other local AGN. (Abridged)

We report a study of X-ray-irradiated gas in the central ~100 pc of the Circinus galaxy, hosting a Compton-thick active galactic nucleus (AGN), at 10-pc resolution using Chandra and ALMA. Based on ~200 ksec Chandra/ACIS-S data, we created an image of the Fe Kalpha line at 6.4 keV, tracing X-ray-irradiated dense gas. The ALMA data in Bands 6 (~270 GHz) and 7 (~350 GHz) cover five molecular lines: CO(J=3--2), HCN(J=3--2), HCN(J=4--3), HCO^+(J=3--2), and HCO^+(J=4--3). The detailed spatial distribution of dense molecular gas was revealed, and compared to the iron line image. The molecular gas emission appeared faint in regions with bright iron emission. Motivated by this, we quantitatively discuss the possibility that the molecular gas is efficiently dissociated by AGN X-ray irradiation (i.e., creating an X-ray-dominated region). Based on a non-local thermodynamic equilibrium model, we constrained the molecular gas densities and found that they are as low as interpreted by X-ray dissociation. Furthermore, judging from inactive star formation (SF) reported in the literature, we suggest that the X-ray emission has potential to suppress SF, particularly in the proximity of the AGN.

We report on the relation between the mass of supermassive black holes (SMBHs; M_BH) and that of hosting dark matter halos (M_h) for 49 z ~ 6 quasi-stellar objects (QSOs) with [CII]158um velocity-width measurements. Here, we estimate M_h assuming that the rotation velocity from FWHM_CII is equal to the circular velocity of the halo; we have tested this procedure using z ~ 3 QSOs that also have clustering-based M_h estimates. We find that a vast majority of the z ~ 6 SMBHs are more massive than expected from the local M_BH - M_h relation, with one-third of the sample by factors >~ 10^2. The median mass ratio of the sample, M_BH/M_h = 6 x 10^{-4}, means that 0.4% of the baryons in halos are locked up in SMBHs. The mass growth rates of our SMBHs amount to ~ 10% of the SFRs, or ~ 1% of the mean baryon accretion rates, of the hosting galaxies. A large fraction of the hosting galaxies are consistent with average galaxies in terms of SFR and perhaps of stellar mass and size. Our study indicates that the growth of SMBHs (M_BH ~ 10^{8-10} Msun) in luminous z ~ 6 QSOs greatly precedes that of hosting halos owing to efficient gas accretion even under normal star formation activities, although we cannot rule out the possibility that undetected SMBHs have local M_BH/M_h ratios. This preceding growth is in contrast to much milder evolution of the stellar-to-halo mass ratio.

We present the results of our ALMA HCN J=3-2 and HCO+ J=3-2 line observations of a uniformly selected sample (>25) of nearby ultraluminous infrared galaxies (ULIRGs) at z < 0.15. The emission of these dense molecular gas tracers and continuum are spatially resolved in the majority of observed ULIRGs for the first time with achieved synthesized beam sizes of ~0.2 arcsec or ~500 pc. In most ULIRGs, the HCN-to-HCO+ J=3-2 flux ratios in the nuclear regions within the beam size are systematically higher than those in the spatially extended regions. The elevated nuclear HCN J=3-2 emission could be related to (a) luminous buried active galactic nuclei, (b) the high molecular gas density and temperature in ULIRG's nuclei, and/or (c) mechanical heating by spatially compact nuclear outflows. A small fraction of the observed ULIRGs display higher HCN-to-HCO+ J=3-2 flux ratios in localized off-nuclear regions than those of the nuclei, which may be due to mechanical heating by spatially extended outflows. The observed nearby ULIRGs are generally rich in dense (>10^5 cm^-3) molecular gas, with an estimated mass of >10^9 Msun within the nuclear (a few kpc) regions, and dense gas can dominate the total molecular mass there. We find a low detection rate (<20%) regarding the possible signature of a vibrationally excited (v2=1f) HCN J=3-2 emission line in the vicinity of the bright HCO+ J=3-2 line that may be due, in part, to the large molecular line widths of ULIRGs.

We report the discovery of a quasar at z = 7.07, which was selected from the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. This quasar, HSC J124353.93+010038.5, has an order of magnitude lower luminosity than do the other known quasars at z > 7. The rest-frame ultraviolet absolute magnitude is M1450 = -24.13 +/- 0.08 mag and the bolometric luminosity is Lbol = (1.4 +/- 0.1) x 10^{46} erg/s. Its spectrum in the optical to near-infrared shows strong emission lines, and shows evidence for a fast gas outflow, as the C IV line is blueshifted and there is indication of broad absorption lines. The Mg II-based black hole mass is Mbh = (3.3 +/- 2.0) x 10^8 Msun, thus indicating a moderate mass accretion rate with an Eddington ratio 0.34 +/- 0.20. It is the first z > 7 quasar with sub-Eddington accretion, besides being the third most distant quasar, known to date. The luminosity and black hole mass are comparable to, or even lower than, those measured for the majority of low-z quasars discovered by the Sloan Digital Sky Survey, and thus this quasar likely represents a z > 7 counterpart to quasars commonly observed in the low-z universe.

Understanding the nuclear growth and feedback processes in galaxies requires investigating their often obscured central regions. One way to do this is to use (sub)millimeter line emission from vibrationally excited HCN (HCN-vib). It has been suggested that the most intense HCN-vib emission from a galaxy is connected to a phase of nuclear growth that occurs before the nuclear feedback processes have been fully developed. We aim to investigate if there is a connection between the presence of strong HCN-vib emission and the development of feedback in (U)LIRGs. We collected literature and archival data to compare the luminosities of rotational lines of HCN-vib, normalized to the total infrared luminosity, to the median velocities of 119 {\mu}m OH absorption lines, potentially indicating outflows, in a total of 17 (U)LIRGs. The most HCN-vib luminous systems all lack signatures of significant molecular outflows in the far-infrared OH absorption lines. However, at least some of the systems with bright HCN-vib emission do have fast and collimated outflows that can be seen in spectral lines at longer wavelengths. We conclude that the galaxy nuclei with the highest L(HCN-vib)/L(IR) do not drive wide-angle outflows detectable using the median velocities of far-infrared OH absorption lines. It is possible that this is due to an orientation effect where sources which are oriented in such a way that their outflows are not along our line of sight also radiate a smaller proportion of their infrared luminosity in our direction. It could also be that massive wide-angle outflows destroy the deeply embedded regions responsible for bright HCN-vib emission, so that the two phenomena cannot coexist. This would strengthen the idea that vibrationally excited HCN traces a heavily obscured stage of evolution before nuclear feedback mechanisms are fully developed.

UltraFast Outflows (UFOs), seen as X-ray blueshifted absorption lines in active galactic nuclei (AGNs), are considered to be a key mechanism for AGN feedback. In this scenario, UFO kinetic energy is transferred into the cold and extended molecular outflow observed at the mm/sub-mm wavelength, which blows away the gas and suppresses star formation and accretion onto the central black hole (BH). However, the energy transfer between the inner UFO and the outer molecular outflow has not yet fully studied mainly due to the limited sample. In this paper, we performed comparison of their kinetic energy using the mm/sub-mm published data and the X-ray archival data. Among fourteen Seyfert galaxies whose molecular outflows are detected in the IRAM/PdBI data, eight targets are bright enough to perform spectral fitting in X-ray, and we have detected UFO absorption lines in six targets with 90% significance level, using XMM-Newton and Suzaku satellites. The time-averaged UFO kinetic energy was derived from the spectral fitting. As a result, we have found that the energy-transfer rate (kinetic energy ratio of the molecular outflow to the UFO) ranges from $\sim7\times10^{-3}$ to $\sim$1, and has a negative correlation with the BH mass, which shows that the AGN feedback is more efficient in the lower mass BHs. This tendency is consistent with the theoretical prediction that the cooling time scale of the outflowing gas becomes longer than the flow time scale when the BH mass is smaller.

We present new measurements of the quasar luminosity function (LF) at $z \sim 6$, over an unprecedentedly wide range of the rest-frame ultraviolet luminosity $M_{1450}$ from $-30$ to $-22$ mag. This is the fifth in a series of publications from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The LF was calculated with a complete sample of 110 quasars at $5.7 \le z \le 6.5$, which includes 48 SHELLQs quasars discovered over 650 deg$^2$, and 63 brighter quasars discovered by the Sloan Digital Sky Survey and the Canada-France-Hawaii Quasar Survey (including one overlapping object). This is the largest sample of $z \sim 6$ quasars with a well-defined selection function constructed to date, and has allowed us to detect significant flattening of the LF at its faint end. A double power-law function fit to the sample yields a faint-end slope $\alpha = -1.23^{+0.44}_{-0.34}$, a bright-end slope $\beta = -2.73^{+0.23}_{-0.31}$, a break magnitude $M_{1450}^* = -24.90^{+0.75}_{-0.90}$, and a characteristic space density $\Phi^* = 10.9^{+10.0}_{-6.8}$ Gpc$^{-3}$ mag$^{-1}$. Integrating this best-fit model over the range $-18 < M_{1450} < -30$ mag, quasars emit ionizing photons at the rate of $\dot{n}_{\rm ion} = 10^{48.8 \pm 0.1}$ s$^{-1}$ Mpc$^{-3}$ at $z = 6.0$. This is less than 10 % of the critical rate necessary to keep the intergalactic medium ionized, which indicates that quasars are not a major contributor to cosmic reionization.

We used the Atacama Large Millimeter/Submillimeter Array (ALMA) to map the CO(3-2) and [CI](1-0) lines, as well as their underlying continuum emission, from the central $\sim 200$ pc region of the Circinus galaxy that hosts the nearest type 2 Seyfert-class active galactic nucleus (AGN), with a spatial resolution of $\sim 6-15$ pc. The lines and continuum-emitting regions consist of a circumnuclear disk (CND; 74 pc $\times$ 34 pc) and spiral arms. The distribution of the continuum emission revealed a temperature-dependent dust geometry and possibly polar dust elongation in the torus region. The molecular mass of the CND is $M_{\rm H2} \sim 3 \times 10^6~M_\odot$ with a beam-averaged H$_2$ column density of $\sim 5 \times 10^{23}$ cm$^{-2}$ toward the AGN position, which contributes significantly to the nuclear obscuration. The [CI](1-0)/CO(3-2) ratio at the AGN position is unusually high, suggesting an X-ray dominated region-type chemistry. We decomposed the observed velocity fields into rotational and dispersion components, and revealed multi-phase dynamic nature in the $r \lesssim 10$ pc torus region, i.e., the diffuse atomic gas is more spatially extended along the vertical direction of the disk than the dense molecular gas. Through comparisons with our model predictions based on the radiation-driven fountain scheme, we indicate that atomic outflows are the driver of the geometrical thickness of the atomic disk. This supports the validity of the radiation-driven fountain scheme in the vicinity of this AGN, which would explain the long-lasting mystery, the physical origin of the AGN torus.

Submillimeter bright galaxies in the early Universe are vigorously forming stars at ~1000 times higher rate than the Milky Way. A large fraction of stars is formed in the central 1 kiloparsec region, that is comparable in size to massive, quiescent galaxies found at the peak of the cosmic star formation history, and eventually the core of giant elliptical galaxies in the present-day Universe. However, the physical and kinematic properties inside a compact starburst core are poorly understood because dissecting it requires angular resolution even higher than the Hubble Space Telescope can offer. Here we report 550 parsec-resolution observations of gas and dust in the brightest unlensed submillimeter galaxy at z=4.3. We map out for the first time the spatial and kinematic structure of molecular gas inside the heavily dust-obscured core. The gas distribution is clumpy while the underlying disk is rotation-supported. Exploiting the high-quality map of molecular gas mass surface density, we find a strong evidence that the starburst disk is gravitationally unstable, implying that the self-gravity of gas overcomes the differential rotation and the internal pressure by stellar radiation feedback. The observed molecular gas would be consumed by star formation in a timescale of 100 million years, that is comparable to those in merging starburst galaxies. Our results suggest that the most extreme starburst in the early Universe originates from efficient star formation due to a gravitational instability in the central 2 kpc region.

RXJ1347.5-1145 (z = 0.451) is one of the most luminous X-ray galaxy clusters, which hosts a prominent cool core and exhibits a signature of a major merger. We present the first direct observational evidence for sub-sonic nature of sloshing motion of the cool core. We find that a residual X-ray image from the Chandra X-ray Observatory after removing the global emission shows a clear dipolar pattern characteristic of gas sloshing, whereas we find no significant residual in the Sunyaev-Zel'dovich effect (SZE) image from the Atacama Large Millimeter/submillimeter Array (ALMA). We estimate the equation of state of perturbations in the gas from the X-ray and SZE residual images. The inferred velocity is 420 +310 -420 km s-1, which is much lower than the adiabatic sound speed of the intracluster medium in the core. We thus conclude that the perturbation is nearly isobaric, and gas sloshing motion is consistent with being in pressure equilibrium. Next, we report evidence for gas stripping of an infalling subcluster, which likely shock-heats gas to high temperature well in excess of 20 keV. Using mass distribution inferred from strong lensing images of the Hubble Space Telescope (HST), we find that the mass peak is located away from the peak position of stripped gas with statistical significance of > 5{\sigma}. Unlike for the gas sloshing, the velocity inferred from the equation of state of the excess hot gas is comparable to the adiabatic sound speed expected for the 20 keV intracluster medium. All of the results support that the southeast substructure is created by a merger. On the other hand, the positional offset between the mass and the gas limits the self-interaction cross section of dark matter to be less than 3.7 h-1 cm2 g-1 (95% CL).

We report discovery of 41 new high-z quasars and luminous galaxies, which were spectroscopically identified at 5.7 < z < 6.9. This is the fourth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, based on the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. We selected the photometric candidates by a Bayesian probabilistic algorithm, and then carried out follow-up spectroscopy with the Gran Telescopio Canarias and the Subaru Telescope. Combined with the sample presented in the previous papers, we have now spectroscopically identified 137 extremely-red HSC sources over about 650 deg2, which include 64 high-z quasars, 24 high-z luminous galaxies, 6 [O III] emitters at z ~ 0.8, and 43 Galactic cool dwarfs (low-mass stars and brown dwarfs). The new quasars span the luminosity range from M1450 ~ -26 to -22 mag, and continue to populate a few magnitude lower luminosities than have been probed by previous wide-field surveys. In a companion paper, we derive the quasar luminosity function at z ~ 6 over an unprecedentedly wide range of M1450 ~ -28 to -21 mag, exploiting the SHELLQs and other survey outcomes.

The Eddington ratio ($\lambda_{\rm Edd}$) of supermassive black holes (SMBHs) is a fundamental parameter that governs the cosmic growth of SMBHs. Although gas mass accretion onto SMBHs is sustained when they are surrounded by large amounts of gas, little is known about the molecular content of galaxies, particularly those hosting super-Eddington SMBHs ($\lambda_{\rm Edd} > 1$: the key phase of SMBH growth). Here, we compiled reported optical and $^{12}$CO(1--0) data of local quasars to characterize their hosts. We found that higher $\lambda_{\rm Edd}$ SMBHs tend to reside in gas rich (i.e., high gas mass to stellar mass fraction = $f_{\rm gas}$) galaxies. We used two methods to make this conclusion: one uses black hole mass as a surrogate for stellar mass by assuming a local co-evolutionary relationship, and the other directly uses stellar masses estimated from near-infrared observations. The $f_{\rm gas}$--$\lambda_{\rm Edd}$ correlation we found concurs with the cosmic decreasing trend in $\lambda_{\rm Edd}$, as cold molecular gas is primarily consumed by star formation. This correlation qualitatively matches predictions of recent semi-analytic models about the cosmic downsizing of SMBHs as well. As the gas mass surface density would eventually be a key parameter controlling mass accretion, we need high-resolution observations to identify further differences in the molecular properties around super-Eddington and sub-Eddington SMBHs.

We present the results of our ALMA observations of eleven (ultra)luminous infrared galaxies ((U)LIRGs) at J=4-3 of HCN, HCO+, HNC and J=3-2 of HNC. This is an extension of our previously published HCN and HCO+ J=3-2 observations to multiple rotational J-transitions of multiple molecules, to investigate how molecular emission line flux ratios vary at different J-transitions. We confirm that ULIRGs that contain or may contain luminous obscured AGNs tend to show higher HCN-to-HCO+ flux ratios than starburst galaxies, both at J=4-3 and J=3-2. For selected HCN-flux-enhanced AGN-important ULIRGs, our isotopologue H13CN, H13CO+, and HN13C J=3-2 line observations suggest a higher abundance of HCN than HCO+ and HNC, which is interpreted to be primarily responsible for the elevated HCN flux in AGN-important galaxies. For such sources, the intrinsic HCN-to-HCO+ flux ratios after line opacity correction will be higher than the observed ratios, making the separation between AGNs and starbursts even larger. The signature of the vibrationally excited (v2=1f) HCN J=4-3 emission line is seen in one ULIRG, IRAS 12112-0305 NE. P Cygni profiles are detected in the HCO+ J=4-3 and J=3-2 lines toward IRAS 15250+3609, with an estimated molecular outflow rate of ~250-750 Mo/year. The SiO J=6-5 line also exhibits a P Cygni profile in IRAS 12112+0305 NE, suggesting the presence of shocked outflow activity. Shock tracers are detected in many sources, suggesting ubiquitous shock activity in the nearby ULIRG population.

We present our ALMA Cycle 4 measurements of the [CII] emission line and the underlying far-infrared (FIR) continuum emission from four optically low-luminosity ($M_{\rm 1450} > -25$) quasars at $z \gtrsim 6$ discovered by the Subaru Hyper Suprime Cam (HSC) survey. The [CII] line and FIR continuum luminosities lie in the ranges $L_{\rm [CII]} = (3.8-10.2) \times 10^8~L_\odot$ and $L_{\rm FIR} = (1.2-2.0) \times 10^{11}~L_\odot$, which are at least one order of magnitude smaller than those of optically-luminous quasars at $z \gtrsim 6$. We estimate the star formation rates (SFR) of our targets as $\simeq 23-40~M_\odot ~{\rm yr}^{-1}$. Their line and continuum-emitting regions are marginally resolved, and found to be comparable in size to those of optically luminous quasars, indicating that their SFR or likely gas mass surface densities (key controlling parameter of mass accretion) are accordingly different. The $L_{\rm [CII]}/L_{\rm FIR}$ ratios of the hosts, $\simeq (2.2-8.7) \times 10^{-3}$, are fully consistent with local star-forming galaxies. Using the [CII] dynamics, we derived their dynamical masses within a radius of 1.5-2.5 kpc as $\simeq (1.4-8.2) \times 10^{10}~M_\odot$. By interpreting these masses as stellar ones, we suggest that these faint quasar hosts are on or even below the star-forming main sequence at $z \sim 6$, i.e., they appear to be transforming into quiescent galaxies. This is in contrast to the optically luminous quasars at those redshifts, which show starburst-like properties. Finally, we find that the ratios of black hole mass to host galaxy dynamical mass of the most of low-luminosity quasars including the HSC ones are consistent with the local value. The mass ratios of the HSC quasars can be reproduced by a semi-analytical model that assumes merger-induced black hole-host galaxy evolution.

We present the results of our ALMA Cycle 4 high-spatial-resolution (0.04-0.07") observations, at HCN J=3-2 and HCO+ J=3-2 lines, of the nucleus of NGC 1068, the nearby prototypical type 2 active galactic nucleus (AGN). Our previous ALMA observations identified the compact emission of these lines at the putative location of the torus around a mass-accreting supermassive black hole. We now report that we have detected the rotation of this compact emission, with the eastern and western sides being redshifted and blueshifted, respectively. Unlike the previously reported CO J=6-5 emission, both the morphological and dynamical alignments of the HCN J=3-2 and HCO+ J=3-2 emission are roughly aligned along the east-west direction (i.e., the expected torus direction), suggesting that these molecular lines are better probes of a rotating dense molecular gas component in the torus. The western part of the torus exhibits larger velocity dispersion and stronger emission in the HCN J=3-2 and HCO+ J=3-2 lines than the eastern part, revealing a highly inhomogeneous molecular torus. The dense molecular gas in the torus and that of the host galaxy at 0.5-2.0" from the AGN along the torus direction are found to be counter-rotating, suggesting an external process happened in the past at the NGC 1068 nucleus.

In this study, we investigate the line emissions from cold molecular gas based on our previous "radiation-driven fountain model" (Wada 2016), which reliably explains the spectral energy distribution of the nearest type 2 Seyfert galaxy, the Circinus galaxy. Using a snapshot of the best-fit radiation-hydrodynamic model for the central r < 16 pc, in which non-equilibrium X-ray-dominated region chemistry is solved, we conduct post-processed, non-local thermodynamic equilibrium radiation transfer simulations for the CO lines. We obtain a spectral line energy distribution with a peak around J=6, and its distribution suggests that the lines are not thermalized. However, for a given line-of-sight, the optical depth distribution is highly non-uniform between $\tau \ll 1$ and $\tau \gg 1$. The CO-to-H2 conversion factor (X_CO), which can be directly obtained from the results, is not a constant and depends strongly on the integrated intensity, and it differs from the fiducial value for local objects. X_CO exhibits a large dispersion of more than one order of magnitude, reflecting the non-uniform internal structure of a "torus." We also found that the physical conditions differ between grid cells on a scale of a few parsecs along the observed lines of sight; therefore, a specific observed line ratio does not necessarily represent a single physical state of the ISM.

We report the study of far-IR sizes of submillimeter galaxies (SMGs) in relation to their dust-obscured star formation rate (SFR) and active galactic nuclei (AGN) presence, determined using mid-IR photometry. We determined the millimeter-wave ($\lambda_{\rm obs}=1100 \mu$m) sizes of 69 ALMA-identified SMGs, selected with $\geq10$$\sigma$ confidence on ALMA images ($F_{\rm 1100 \mu m}=1.7$--7.4 mJy). We found that all the SMGs are located above an avoidance region in the millimeter size-flux plane, as expected by the Eddington limit for star formation. In order to understand what drives the different millimeter-wave sizes in SMGs, we investigated the relation between millimeter-wave size and AGN fraction for 25 of our SMGs at $z=1$--3. We found that the SMGs for which the mid-IR emission is dominated by star formation or AGN have extended millimeter-sizes, with respective median $R_{\rm c,e} = 1.6^{+0.34}_{-0.21}$ and 1.5$^{+0.93}_{-0.24}$ kpc. Instead, the SMGs for which the mid-IR emission corresponds to star-forming/AGN composites have more compact millimeter-wave sizes, with median $R_{\rm c,e}=1.0^{+0.20}_{-0.20}$ kpc. The relation between millimeter-wave size and AGN fraction suggests that this size may be related to the evolutionary stage of the SMG. The very compact sizes for composite star-forming/AGN systems could be explained by supermassive black holes growing rapidly during the SMG coalescing, star-formation phase.

We present an 8 pc x 5 pc resolution view of the central ~ 200 pc region of the nearby starburst galaxy NGC 253, based on ALMA Band 7 ({\lambda} ~ 0.85 mm or {\nu} ~ 350 GHz) observations covering 11 GHz. We resolve the nuclear starburst of NGC 253 into eight dusty star-forming clumps, 10 pc in scale, for the first time. These clumps, each of which contains (4-10) x 10^4 M_sun of dust (assuming that the dust temperature is 25 K) and up to 6 x 10^2 massive (O5V) stars, appear to be aligned in two parallel ridges, while they have been blended in previous studies. Despite the similarities in sizes and dust masses of these clumps, their line spectra vary drastically from clump to clump although they are separated by only ~ 10 pc. Specifically, one of the clumps, Clump 1, exhibits line confusion-limited spectra with at least 36 emission lines from 19 molecules (including CH3OH, HNCO, H2CO, CH3CCH, H2CS, and H3O+) and a hydrogen recombination line (H26{\alpha}), while much fewer kinds of molecular lines are detected in some other clumps where fragile species, such as complex organic molecules and HNCO, completely disappear from their spectra. We demonstrate the existence of hot molecular gas ($T_rot (SO2) = 90 {\pm} 11 K) in the former clump, which suggests that the hot and chemically rich environments are localized within a 10-pc scale star-forming clump.

We present our ALMA multi-transition molecular line observational results for the ultraluminous infrared galaxy, IRAS 20551-4250, which is known to contain a luminous buried AGN and shows detectable vibrationally excited (v2=1f) HCN and HNC emission lines. The rotational J=1-0, 4-3, and 8-7 of HCN, HCO+, and HNC emission lines were clearly detected at a vibrational ground level (v=0). Vibrationally excited (v2=1f) J=4-3 emission lines were detected for HCN and HNC, but not for HCO+. Their observed flux ratios further support our previously obtained suggestion, based on J=3-2 data, that (1) infrared radiative pumping plays a role in rotational excitation at v=0, at least for HCN and HNC, and (2) HCN abundance is higher than HCO+ and HNC. The flux measurements of the isotopologue H13CN, H13CO+, and HN13C J=3-2 emission lines support the higher HCN abundance scenario. Based on modeling with collisional excitation, we constrain the physical properties of these line-emitting molecular gas, but find that higher HNC rotational excitation than HCN and HCO+ is difficult to explain, due to the higher effective critical density of HNC. We consider the effects of infrared radiative pumping using the available 5-30 micron infrared spectrum and find that our observational results are well explained if the radiation source is located at 30-100 pc from the molecular gas. The simultaneously covered very bright CO J=3-2 emission line displays a broad emission wing, which we interpret as being due to molecular outflow activity with the estimated rate of ~150 Msun/yr.

We used the Atacama Large Millimeter/submillimeter Array (ALMA) to map the CO(3-2) and the underlying continuum emissions around the type 1 low-luminosity active galactic nucleus (LLAGN; bolometric luminosity $\lesssim 10^{42}$ erg~s$^{-1}$) of NGC 1097 at $\sim 10$ pc resolution. These observations revealed a detailed cold gas distribution within a $\sim 100$ pc of this LLAGN. In contrast to the luminous Seyfert galaxy NGC 1068, where a $\sim 7$ pc cold molecular torus was recently revealed, a distinctively dense and compact torus is missing in our CO(3-2) integrated intensity map of NGC 1097. Based on the CO(3-2) flux, the gas mass of the torus of NGC 1097 would be a factor of $\gtrsim 2-3$ less than that found for NGC 1068 by using the same CO-to-H$_2$ conversion factor, which implies less active nuclear star formation and/or inflows in NGC 1097. Our dynamical modeling of the CO(3-2) velocity field implies that the cold molecular gas is concentrated in a thin layer as compared to the hot gas traced by the 2.12 $\mu$m H$_2$ emission in and around the torus. Furthermore, we suggest that NGC 1097 hosts a geometrically thinner torus than NGC 1068. Although the physical origin of the torus thickness remains unclear, our observations support a theoretical prediction that geometrically thick tori with high opacity will become deficient as AGNs evolve from luminous Seyferts to LLAGNs.

We present the results of a blind millimeter line emitter search using ALMA Band 6 data with a single frequency tuning toward four gravitational lensing clusters (RXJ1347.5-1145, Abell S0592, MACS J0416.1-2403, and Abell 2744). We construct three-dimensional signal-to-noise ratio (S/N) cubes with 60 MHz and 100 MHz binning, and search for millimeter line emitters. We do not detect any line emitters with a peak S/N > 5, although we do find a line emitter candidate with a peak S/N ~ 4.5. These results provide upper limits to the CO(3-2), CO(4-3), CO(5-4), and [CII] luminosity functions at z ~ 0.3, 0.7, 1.2, and 6, respectively. Because of the magnification effect of gravitational lensing clusters, the new data provide the first constraints on the CO and [CII] luminosity functions at unprecedentedly low luminosity levels, i.e., down to $<~ 10^{-3}-10^{-1} $Mpc$^{-3}$ dex$^{-1}$ at $L'_{CO} ~ 10^8-10^{10}$ K km $s^{-1}$ pc$^2$ and $<~ 10^{-3}-10^{-2}$ Mpc$^{-3}$ dex$^{-1}$ at $L_[CII] ~ 10^8-10^{10} L_\odot$, respectively. Although the constraints to date are not stringent yet, we find that the evolution of the CO and [CII] luminosity functions are broadly consistent with the predictions of semi-analytical models. This study demonstrates that the wide observations with a single frequency tuning toward gravitational lensing clusters are promising for constraining the CO and [CII] luminosity functions.

We present the distribution and kinematics of the molecular gas in the circumnuclear disk (CND, 400 pc x 200 pc) of Centaurus A with resolutions of ~5 pc (0.3 arcsec) and shed light onto the mechanism feeding the Active Galactic Nucleus (AGN) using CO(3-2), HCO+(4-3), HCN(4-3), and CO(6-5) observations obtained with ALMA. Multiple filaments or streamers of tens to a hundred parsec scale exist within the CND, which form a ring-like structure with an unprojected diameter of 9 x 6 arcsec (162pc x 108pc) and a position angle PA = 155deg. Inside the nuclear ring, there are two leading and straight filamentary structures with lengths of about 30-60pc at PA = 120deg on opposite sides of the AGN, with a rotational symmetry of 180deg and steeper position-velocity diagrams, which are interpreted as nuclear shocks due to non-circular motions. Along the filaments, and unlike other nearby AGNs, several dense molecular clumps present low HCN/HCO+(4-3) ratios (~0.5). The filaments abruptly end in the probed transitions at r = 20pc from the AGN, but previous near-IR H2 (J=1-0) S(1) maps show that they continue in an even ~1000 K), winding up in the form of nuclear spirals, and forming an inner ring structure with another set of symmetric filaments along the N-S direction and within r = 10pc. The molecular gas is governed primarily by non-circular motions, being the successive shock fronts at different scales where loss of angular momentum occurs, a mechanism which may feed efficiently powerful radio galaxies down to parsec scales.

We present spectroscopic identification of 32 new quasars and luminous galaxies discovered at 5.7 < z < 6.8. This is the second in a series of papers presenting the results of the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The photometric candidates were selected by a Bayesian probabilistic algorithm, and then observed with spectrographs on the Gran Telescopio Canarias and the Subaru Telescope. Combined with the sample presented in the previous paper, we have now identified 64 HSC sources over about 430 deg2, which include 33 high-z quasars, 14 high-z luminous galaxies, 2 [O III] emitters at z ~ 0.8, and 15 Galactic brown dwarfs. The new quasars have considerably lower luminosity (M1450 ~ -25 to -22 mag) than most of the previously known high-z quasars. Several of these quasars have luminous (> 10^(43) erg/s) and narrow (< 500 km/s) Ly alpha lines, and also a possible mini broad absorption line system of N V 1240 in the composite spectrum, which clearly separate them from typical quasars. On the other hand, the high-z galaxies have extremely high luminosity (M1450 ~ -24 to -22 mag) compared to other galaxies found at similar redshift. With the discovery of these new classes of objects, we are opening up new parameter spaces in the high-z Universe. Further survey observations and follow-up studies of the identified objects, including the construction of the quasar luminosity function at z ~ 6, are ongoing.

We investigate gas contents of star-forming galaxies associated with protocluster 4C23.56 at z = 2.49 by using the redshifted CO(3-2) and 1.1 mm dust continuum with the Atacama Large Millimeter/submillimeter Array. The observations unveil seven CO detections out of 22 targeted H$\alpha$ emitters (HAEs) and four out of 19 in 1.1 mm dust continuum. They have high stellar mass ($M_{\star}>4\times 10^{10}$ $M_{\odot}$) and exhibit a specific star-formation rate typical of main-sequence star forming galaxies at $z\sim2.5$. Different gas mass estimators from CO(3-2) and 1.1 mm yield consistent values for simultaneous detections. The gas mass ($M_{\rm gas}$) and gas fraction ($f_{\rm gas}$) are comparable to those of field galaxies, with $M_{\rm gas}=[0.3, 1.8]\times10^{11} \times (\alpha_{\rm CO}/(4.36\times A(Z)$)) M$_{\odot}$, where $\alpha_{\rm CO}$ is the CO-to-H$_2$ conversion factor and $A(Z)$ the additional correction factor for the metallicity dependence of $\alpha_{\rm CO}$, and $\langle f_{\rm gas}\rangle = 0.53 \pm 0.07$ from CO(3-2). Our measurements place a constraint on the cosmic gas density of high-$z$ protoclusters, indicating the protocluster is characterized by a gas density higher than that of the general fields by an order of magnitude. We found $\rho (H_2)\sim 5 \times 10^9 \,M_{\odot}\,{\rm Mpc^{-3}}$ with the CO(3-2) detections. The five ALMA CO detections occur in the region of highest galaxy surface density, where the density positively correlates with global star-forming efficiency (SFE) and stellar mass. Such correlations imply a potentially critical role of environment on early galaxy evolution at high-z protoclusters, although future observations are necessary for confirmation.

We report the identification of four millimeter line emitting galaxies with the Atacama Large Milli/submillimeter Array (ALMA) in SSA22 Field (ADF22). We analyze the ALMA 1.1 mm survey data, with an effective survey area of 5 arcmin$^2$, a frequency range of 253.1--256.8 and 269.1--272.8 GHz, angular resolution of 0".7 and RMS noise of 0.8 mJy beam$^{-1}$ at 36 km s$^{-1}$ velocity resolution. We detect four line emitter candidates with significance levels above $6 \sigma$. We identify one of the four sources as a CO(9-8) emitter at $z = 3.1$ in a member of the proto-cluster known in this field. Another line emitter with an optical counterpart is likely a CO(4-3) emitter at $z = 0.7$. The other two sources without any millimeter continuum or optical/near-infrared counterpart are likely to be [CII] emitter candidates at $z = 6.0$ and $6.5$. The equivalent widths of the [CII] candidates are consistent with those of confirmed high-redshift [CII] emitters and candidates, and are a factor of 10 times larger than that of the CO(9-8) emitter detected in this search. The [CII] luminosity of the candidates are $4-7 \times 10^8~\rm L_\odot$. The star formation rates (SFRs) of these sources are estimated to be $10-20~\rm M_{\odot}~yr^{-1}$ if we adopt an empirical [CII] luminosity - SFR relation. One of them has a relatively low-S/N ratio, but shows features characteristic of emission lines. Assuming that at least one of the two candidates is a [CII] emitter, we derive a lower limit of [CII]-based star formation rate density (SFRD) at $z~\sim~6$. The resulting value of $> 10^{-2}$ $\rm M_\odot yr^{-1} Mpc^{-3}$ is consistent with the dust-uncorrected UV-based SFRD. Future millimeter/submillimeter surveys can be used to detect a number of high redshift line emitters, with which to study the star formation history in the early Universe.

We present the detailed characterization of two extremely red submillimeter galaxies (SMGs), ASXDF1100.053.1 and 231.1, with the Atacama Large Millimeter/submillimeter Array (ALMA) and the Jansky Very Large Array (VLA). These SMGs were selected originally using AzTEC at 1100 micron, and are observed by Herschel to be faint at 100--500 micron. Their (sub)millimeter colors are as red as -- or redder -- than known z>~5 SMGs; indeed, ASXDF1100.053.1 is redder than HFLS 3, which lies at z=6.3. They are also faint and red in the near-/mid-infrared: ~1 microJy at IRAC 4.5 micron and <0.2 microJy in the Ks filter. These SMGs are also faint in the radio waveband, where F_6GHz=4.5 microJy for ASXDF1100.053.1 and F_1.4GHz=28 microJy for ASXDF1100.231.1, suggestive of z=6.5^{+1.4}_{-1.1} and z=4.1^{+0.6}_{-0.7} for ASXDF1100.053.1 and 231.1, respectively. ASXDF1100.231.1 has a flux excess in the 3.6-micron filter, probably due to H$\alpha$ emission at z=4--5. Derived properties of ASXDF1100.053.1 for z=5.5--7.5 and 231.1 for z=3.5--5.5 are as follows: their infrared luminosities are [6.5-7.4]x10^{12} and [4.2-4.5]x10^{12} L_sun; their stellar masses are [0.9-2]x10^{11} and [0.4-3]x10^{10} M_sun; their circularized half-light radii in the ALMA maps are ~1 and <~0.2 kpc (~2--3 kpc for 90% of the total flux). Lastly, their surface infrared luminosity densities, Sigma_IR, are ~1x10^{12} and >~1.5x10^{13} L_sun kpc^{-2}, similar to values seen for local (U)LIRGs. These data suggest that ASXDF1100.053.1 and 231.1 are compact SMGs at z>~4 and can plausibly evolve into z>~3 compact quiescent galaxies.

We report the result from observations conducted with the Atacama Large Millimeter/submillimeter Array (ALMA) to detect [CII] 158 um fine structure line emission from galaxies embedded in one of the most spectacular Lyman-alpha blobs (LABs) at z=3.1, SSA22-LAB1. Of three dusty star-forming galaxies previously discovered by ALMA 860 um dust continuum survey toward SSA22-LAB1, we detected the [CII] line from one, LAB1-ALMA3 at z=3.0993+/-0.0004. No line emission was detected, associated with the other ALMA continuum sources or from three rest-frame UV/optical selected z_spec~3.1 galaxies within the field of view. For LAB1-ALMA3, we find relatively bright [CII] emission compared to the infrared luminosity (L_[CII]/L_[CII]) and an extremely high [CII] 158 um and [NII] 205 um emission line ratio (L_[CII]/L_[NII]>55). The relatively strong [CII] emission may be caused by abundant photodissociation regions and sub-solar metallicity, or by shock heating. The origin of the unusually strong [CII] emission could be causally related to the location within the giant LAB, although the relationship between extended Lyman-alpha emission and ISM conditions of associated galaxies is yet to be understand.

We present new ALMA Band 7 ($\sim340$ GHz) observations of the dense gas tracers HCN, HCO$^+$, and CS in the local, single-nucleus, ultraluminous infrared galaxy IRAS 13120-5453. We find centrally enhanced HCN (4-3) emission, relative to HCO$^+$ (4-3), but do not find evidence for radiative pumping of HCN. Considering the size of the starburst (0.5 kpc) and the estimated supernovae rate of $\sim1.2$ yr$^{-1}$, the high HCN/HCO$^+$ ratio can be explained by an enhanced HCN abundance as a result of mechanical heating by the supernovae, though the active galactic nucleus and winds may also contribute additional mechanical heating. The starburst size implies a high $\Sigma_{IR}$ of $4.7\times10^{12}$ $L_{\odot}$ kpc$^{-2}$, slightly below predictions of radiation-pressure limited starbursts. The HCN line profile has low-level wings, which we tentatively interpret as evidence for outflowing dense molecular gas. However, the dense molecular outflow seen in the HCN line wings is unlikely to escape the galaxy and is destined to return to the nucleus and fuel future star formation. We also present modeling of Herschel observations of the H$_2$O lines and find a nuclear dust temperature of $\sim40$ K. IRAS 13120-5453 has a lower dust temperature and $\Sigma_{IR}$ than is inferred for the systems termed "compact obscured nuclei" (such as Arp 220 and Mrk 231). If IRAS 13120-5453 has undergone a compact obscured nucleus phase, we are likely witnessing it at a time when the feedback has already inflated the nuclear ISM and diluted star formation in the starburst/AGN core.

We present $1''.4$ (98 pc) resolution ALMA observations of $^{13}$CO($J$=1-0), C$^{18}$O($J$=1-0), CS($J$=2-1) and CH$_3$OH($J_K$=$2_K$-$1_K$) molecular lines in the central $1'$ (4.2 kpc) region of NGC 1068 to study the physical and chemical properties of giant molecular clouds (GMCs) and to test whether these GMC-scale properties are linked to the larger-scale galactic environment. Using the derived $^{13}$CO cube, we have identified 187 high-significance ($> 8\sigma$) GMCs by employing the CLUMPFIND algorithm. The molecular gas masses of GMCs $M_{\rm ^{13}CO}$, derived from the $^{13}$CO data, range from $1.8 \times 10^4 ~ M_\odot$ to $4.2 \times 10^7 ~M_\odot$. A mass function of GMCs in NGC 1068 has been obtained for the first time at $\sim$100 pc resolution. We find the slope of the mass function $\gamma = -1.25 \pm 0.07$ for a mass range of $M_{\rm ^{13}CO} \geq 10^5M_\odot$. This is shallower than the GMCs in the disk regions of the Milky Way, M51 and NGC 300. Further, we find that the high mass cut-off of the GMC mass function occurs at $M_{\rm ^{13}CO} \sim 6 \times 10^7 M_\odot$, which is an order of magnitude larger than that in the nuclear bar region of M51, indicating that the more massive clouds dominate the mass budget in NGC 1068. The observed C$^{18}$O)/$^{13}$CO intensity ratios are found to be fairly uniform (0.27 $\pm$ 0.05) among the identified GMCs. In contrast, the CH$_3$OH/$^{13}$CO ratios exhibit striking spatial variation across the disk, with the smallest values around the bar-end ($<$ 0.03), and larger ratios along the spiral arms ($\sim$ 0.1-0.2). We find that GMCs with detectable methanol emission tend to have systematically larger velocity widths than those without methanol emission, suggesting that (relatively weak) shocks are responsible for the enhancement of the CH$_3$OH/$^{13}$CO ratios of GMCs in the disk of NGC 1068.

The central structure in three of the brightest unlensed z=3-4 submillimeter galaxies are investigated through 0.015" - 0.05" (120 -- 360~pc) 860 micron continuum images obtained using the Atacama Large Millimeter/submillimeter Array (ALMA). The distribution in the central kpc in AzTEC1 and AzTEC8 are extremely complex, and they are composed of multiple ~200 pc clumps. AzTEC4 consists of two sources that are separated by ~1.5 kpc, indicating a mid-stage merger. The peak star formation rate densities in the central clumps are ~300 - 3000 Msun/yr/kpc^2, suggesting regions with extreme star formation near the Eddington Limit. By comparing the flux obtained by ALMA and Submillimeter Array (SMA), we find that 68-90% of the emission is extended (> 1 kpc) in AzTEC 4 and 8. For AzTEC1, we identify at least 11 additional compact (~200 pc) clumps in the extended 3 - 4 kpc region. Overall, the data presented here suggest that the luminosity surface densities observed at < 150 pc scales are roughly similar to that observed in local ULIRGs, as in the eastern nucleus of Arp 220. Between 10 to 30% of the 860 micron continuum is concentrated in clumpy structures in the central kpc while the remaining flux is distributed over > 1 kpc regions, some of which could also be clumpy. These sources can be explained by a rapid inflow of gas such as a merger of gas-rich galaxies, surrounded by extended and clumpy starbursts. However, the cold mode accretion model is not ruled out.

We present the results of our ALMA observations of three AGN-dominated nuclei in optical Seyfert 1 galaxies (NGC 7469, I Zw 1, and IC 4329 A) and eleven luminous infrared galaxies (LIRGs) with various levels of infrared estimated energetic contributions by AGNs at the HCN and HCO+ J=3-2 emission lines. The HCN and HCO+ J=3-2 emission lines are clearly detected at the main nuclei of all sources, except for IC 4329 A. The vibrationally excited (v2=1f) HCN J=3-2 and HCO+ J=3-2 emission lines are simultaneously covered, and HCN v2=1f J=3-2 emission line signatures are seen in the main nuclei of two LIRGs, IRAS 12112+0305 and IRAS 22491-1808, neither of which show clear buried AGN signatures in the infrared. If the vibrational excitation is dominated by infrared radiative pumping, through the absorption of infrared 14 um photons, primarily originating from AGN-heated hot dust emission, then these two LIRGs may contain infrared-elusive, but (sub)millimeter-detectable, extremely deeply buried AGNs. These vibrationally excited emission lines are not detected in the three AGN-dominated optical Seyfert 1 nuclei. However, the observed HCN v2=1f to v=0 flux ratios in these optical Seyferts are still consistent with the intrinsic flux ratios in LIRGs with detectable HCN v2=1f emission lines. The observed HCN-to-HCO+ J=3-2 flux ratios tend to be higher in galactic nuclei with luminous AGN signatures compared with starburst-dominated regions, as previously seen at J=1-0 and J=4-3.

We present the first image of the thermal Sunyaev-Zel'dovich effect (SZE) obtained by the Atacama Large Millimeter/submillimeter Array (ALMA). Combining 7-m and 12-m arrays in Band 3, we create an SZE map toward a galaxy cluster RXJ1347.5-1145 with 5 arc-second resolution (corresponding to the physical size of 20 kpc/h), the highest angular and physical spatial resolutions achieved to date for imaging the SZE, while retaining extended signals out to 40 arc-seconds. The 1-sigma statistical sensitivity of the image is 0.017 mJy/beam or 0.12 mK_CMB at the 5 arc-second full width at half maximum. The SZE image shows a good agreement with an electron pressure map reconstructed independently from the X-ray data and offers a new probe of the small-scale structure of the intracluster medium. Our results demonstrate that ALMA is a powerful instrument for imaging the SZE in compact galaxy clusters with unprecedented angular resolution and sensitivity. As the first report on the detection of the SZE by ALMA, we present detailed analysis procedures including corrections for the missing flux, to provide guiding methods for analyzing and interpreting future SZE images by ALMA.

We present a multi-wavelength analysis of five submillimeter sources (S_1.1mm = 0.54-2.02 mJy) that were detected during our 1.1-mm-deep continuum survey in the SXDF-UDS-CANDELS field (2 arcmin^2, 1sigma = 0.055 mJy beam^-1) using the Atacama Large Millimeter/submillimeter Array (ALMA). The two brightest sources correspond to a known single-dish (AzTEC) selected bright submillimeter galaxy (SMG), whereas the remaining three are faint SMGs newly uncovered by ALMA. If we exclude the two brightest sources, the contribution of the ALMA-detected faint SMGs to the infrared extragalactic background light is estimated to be ~ 4.1^{+5.4}_{-3.0} Jy deg^{-2}, which corresponds to ~ 16^{+22}_{-12}% of the infrared extragalactic background light. This suggests that their contribution to the infrared extragalactic background light is as large as that of bright SMGs. We identified multi-wavelength counterparts of the five ALMA sources. One of the sources (SXDF-ALMA3) is extremely faint in the optical to near-infrared region despite its infrared luminosity (L_IR ~ 1e12 L_sun or SFR ~ 100 M_sun yr^{-1}). By fitting the spectral energy distributions (SEDs) at the optical-to-near-infrared wavelengths of the remaining four ALMA sources, we obtained the photometric redshifts (z_photo) and stellar masses (M_*): z_photo ~ 1.3-2.5, M_* ~ (3.5-9.5)e10 M_sun. We also derived their star formation rates (SFRs) and specific SFRs (sSFRs) as ~ 30-200 M_sun yr^{-1} and ~ 0.8-2 Gyr^{-1}, respectively. These values imply that they are main-sequence star-forming galaxies.

We present a positive correlation between the mass of dense molecular gas ($M_{\rm dense}$) of $\sim 100$ pc scale circumnuclear disks (CNDs) and the black hole mass accretion rate ($\dot{M}_{\rm BH}$) in total 10 Seyfert galaxies, based on data compiled from the literature and an archive (median aperture $\theta_{\rm med}$ = 220 pc). A typical $M_{\rm dense}$ of CNDs is 10$^{7-8}$ $M_\odot$, estimated from the luminosity of the dense gas tracer, the HCN($1-0$) emission line. Because dense molecular gas is the site of star formation, this correlation is virtually equivalent to the one between nuclear star formation rate and $\dot{M}_{\rm BH}$ revealed previously. Moreover, the $M_{\rm dense}-\dot{M}_{\rm BH}$ correlation was tighter for CND-scale gas than for the gas on kpc or larger scales. This indicates that CNDs likely play an important role in fueling black holes, whereas $>$kpc scale gas does not. To demonstrate a possible approach for studying the CND-scale accretion process with the Atacama Large Millimeter/submillimeter Array (ALMA), we used a mass accretion model where angular momentum loss due to supernova explosions is vital. Based on the model prediction, we suggest that only the partial fraction of the mass accreted from the CND ($\dot{M}_{\rm acc}$) is consumed as $\dot{M}_{\rm BH}$. However, $\dot{M}_{\rm acc}$ agrees well with the total nuclear mass flow rate (i.e., $\dot{M}_{\rm BH}$ + outflow rate). Although these results are still tentative with large uncertainties, they support the view that star formation in CNDs can drive mass accretion onto supermassive black holes in Seyfert galaxies.

We present the results of ALMA Cycle 2 observations of the ultraluminous infrared galaxy, IRAS 20551-4250, at HCN/HCO+/HNC J=3-2 lines at both vibrational-ground (v=0) and vibrationally excited (v2=1) levels. This galaxy contains a luminous buried active galactic nucleus (AGN), in addition to starburst activity, and our ALMA Cycle 0 data revealed a tentatively detected vibrationally excited HCN v2=1f J=4-3 emission line. In our ALMA Cycle 2 data, the HCN/HCO+/HNC J=3-2 emission lines at v=0 are clearly detected. The HCN and HNC v2=1f J=3-2 emission lines are also detected, but the HCO+ v2=1f J=3-2 emission line is not. Given the high-energy level of v2=1 and the resulting difficulty of collisional excitation, we compared these results with those of the calculation of infrared radiative pumping, using the available infrared 5-35 micron spectrum. We found that all of the observational results were reproduced, if the HCN abundance was significantly higher than that of HCO+ and HNC. The flux ratio and excitation temperature between v2=1f and v=0, after correction for possible line opacity, suggests that infrared radiative pumping affects rotational (J-level) excitation at v=0 at least for HCN and HNC. The HCN-to-HCO+ v=0 flux ratio is higher than those of starburst-dominated regions, and will increase even more when thederived high HCN opacity is corrected. The enhanced HCN-to-HCO+ flux ratio in this AGN-hosting galaxy can be explained by the high HCN-to-HCO+ abundance ratio and sufficient HCN excitation at up to J=4, rather than the significantly higher efficiency of infrared radiative pumping for HCN than HCO+.

We present the results of our ALMA Cycle 2 high angular resolution (0.1-0.2 arcsec) observations of the nuclear region of the nearby well-studied type-2 active galactic nucleus (AGN), NGC 1068, at HCN J=3-2 and HCO+ J=3-2 emission lines. For the first time, due to a higher angular resolution than previous studies, we clearly detected dense molecular gas emission at the putative AGN location, identified as a ~1.1 mm (~266 GHz) continuum emission peak, by separating this emission from brighter emission located at 0.5-2.0 arcsec on the eastern and western sides of the AGN. The estimated intrinsic molecular emission size and dense molecular mass, which are thought to be associated with the putative dusty molecular torus around an AGN, were ~10 pc and ~several x 10^5 Msun, respectively. HCN-to-HCO+ J=3-2 flux ratios substantially higher than unity were found throughout the nuclear region of NGC 1068. The continuum emission displayed an elongated morphology along the direction of the radio jet located at the northern side of the AGN, as well as a weak spatially resolved component at ~2.0 arcsec on the southwestern side of the AGN. The latter component most likely originated from star formation, with the estimated luminosity more than one order of magnitude lower than the luminosity of the central AGN. No vibrationally excited (v2=1f) J=3-2 emission lines were detected for HCN and HCO+ across the field of view.

Hydrogen recombination lines at the submillimetre band (submm-RLs) can serve as probes of ionized gas without dust extinction. One therefore expects to probe the broad line region (BLR) of an obscured (type 2) active galactic nucleus (AGN) with those lines. However, admitting the large uncertainty in the continuum level, here we report on the non-detection of both broad and narrow H26$\alpha$ emission line (rest frequency = 353.62 GHz) towards the prototypical type 2 AGN of NGC 1068 with the Atacama Large Millimeter/submillimeter Array (ALMA). We also investigate the nature of BLR clouds that can emit submm-RLs with model calculations. As a result, we suggest that clouds with an electron density ($N_e$) of $\sim$ 10$^9$ cm$^{-3}$ can mainly contribute to broad submm-RLs in terms of the line flux. On the other hand, line flux from other density clouds would be insignificant considering their too large or too small line optical depths. However, even for the case of $N_e \sim 10^9$ cm$^{-3}$ clouds, we also suggest that the expected line flux is extremely low, which is impractical to detect even with ALMA.

In the central regions of active galaxies, dense molecular medium are exposed to various types of radiation and energy injections, such as UV, X-ray, cosmic ray, and shock dissipation. With the rapid progress of chemical models and implementation of new-generation mm/submm interferometry, we are now able to use molecules as powerful diagnostics of the physical and chemical processes in galaxies. Here we give a brief overview on the recent ALMA results to demonstrate how molecules can reveal underlying physical and chemical processes in galaxies. First, new detections of Galactic molecular absorption systems with elevated HCO/H$^{13}$CO$^+$ column density ratios are reported, indicating that these molecular media are irradiated by intense UV fields. Second, we discuss the spatial distributions of various types of shock tracers including HNCO, CH$_3$OH and SiO in NGC 253 and NGC 1068. Lastly, we provide an overview of proposed diagnostic methods of nuclear energy sources using ALMA, with an emphasis on the synergy with sensitive mid-infrared spectroscopy, which will be implemented by JWST and SPICA to disentangle the complex nature of heavily obscured galaxies across the cosmic time.

Compiling data from literature and the ALMA archive, we show enhanced HCN(4-3)/HCO$^+$(4-3) and/or HCN(4-3)/CS(7-6) integrated intensity ratios in circumnuclear molecular gas around active galactic nuclei (AGNs) compared to those in starburst (SB) galaxies (submillimeter HCN-enhancement). The number of sample galaxies is significantly increased from our previous work. We expect this feature could potentially be an extinction-free energy diagnostic tool of nuclear regions of galaxies. Non-LTE radiative transfer modelings of the above molecular emission lines involving both collisional and radiative excitation, as well as a photon trapping effect were conducted to investigate the cause of the high line ratios in AGNs. As a result, we found that enhanced abundance ratios of HCN-to-HCO$^+$ and HCN-to-CS in AGNs as compared to SB galaxies by a factor of a few to even $>$ 10 is a plausible explanation for the submillimeter HCN-enhancement. However, a counter argument of a systematically higher gas density in AGNs than in SB galaxies can also be a plausible scenario. Although we could not fully discriminate these two scenarios at this moment due to insufficient amount of multi-transition, multi-species data, the former scenario equivalently claims for abnormal chemical composition in AGNs. Regarding the actual mechanism to realize the composition, we suggest it is difficult with conventional gas phase X-ray dominated region (XDR) ionization models to reproduce the observed high line ratios. We might have to take into account other mechanisms such as neutral-neutral reactions that are efficiently activated at high temperature environments and/or mechanically heated regions to further understand the high line ratios in AGNs.

We report the results of $1^{\prime}.5 \times3^{\prime}$ mapping at 1.1~mm with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the central region of the $z=3.09$ SSA22 protocluster. By combining our source catalog with archival spectroscopic redshifts, we find that eight submillimeter galaxies (SMGs) with flux densities, $S_{\rm 1.1~mm}=0.7-6.4$~mJy ($L_{\rm IR}\sim10^{12.1}-10^{13.1}L_\odot$) are at $z=3.08-3.10$. Not only are these SMGs members of the protocluster but they in fact reside within the node at the junction of the 50 Mpc-scale filamentary three-dimensional structure traced by Lyman-$\alpha$ emitters (LAEs) in this field. The eight SMGs account for a star formation rate density (SFRD) $\sim$10 $M_\odot$ yr$^{-1}$ Mpc$^{-3}$ in the node, which is two orders of magnitudes higher than the global SFRD at this redshift. We find that four of the eight SMGs host a X-ray luminous active galactic nuclei (AGN). Our results suggest that the vigorous star formation activity and the growth of super massive black holes (SMBHs) occurred simultaneously in the densest regions at $z\sim3$, which may correspond to the most active historical phase of the massive galaxy population found in the core of the clusters in the present universe. Two SMGs are associated with Lyman-$\alpha$ blobs (LABs), implying that the two populations coexist in high density environments for a few cases.

We report on Atacama Large Millimeter/submillimeter Array (ALMA) detections of molecular absorption lines in Bands 3, 6 and 7 toward four radio-loud quasars, which were observed as the bandpass and complex gain calibrators. The absorption systems, three of which are newly detected, are found to be Galactic origin. Moreover, HCO absorption lines toward two objects are detected, which almost doubles the number of HCO absorption samples in the Galactic diffuse medium. In addition, high HCO to H13CO+ column density ratios are found, suggesting that the interstellar media (ISM) observed toward the two calibrators are in photodissociation regions, which observationally illustrates the chemistry of diffuse ISM driven by ultraviolet (UV) radiation. These results demonstrate that calibrators in the ALMA Archive are potential sources for the quest for new absorption systems and for detailed investigation of the nature of the ISM.

We present ALMA Cycle 1 observations of the central kpc region of the luminous type-1 Seyfert galaxy NGC 7469 with unprecedented high resolution (0.5$"$ $\times$ 0.4$"$ = 165 pc $\times$ 132 pc) at submillimeter wavelengths. Utilizing the wide-bandwidth of ALMA, we simultaneously obtained HCN(4-3), HCO$^+$(4-3), CS(7-6), and partially CO(3-2) line maps, as well as the 860 $\mu$m continuum. The region consists of the central $\sim$ 1$"$ component and the surrounding starburst ring with a radius of $\sim$ 1.5$"$-2.5$"$. Several structures connect these components. Except for CO(3-2), these dense gas tracers are significantly concentrated towards the central $\sim$ 1$"$, suggesting their suitability to probe the nuclear regions of galaxies. Their spatial distribution resembles well those of centimeter and mid-infrared continuum emissions, but it is anti-correlated with the optical one, indicating the existence of dust obscured star formation. The integrated intensity ratios of HCN(4-3)/HCO$^+$(4-3) and HCN(4-3)/CS(7-6) are higher at the AGN position than at the starburst ring, which is consistent to our previous findings (submm-HCN enhancement). However, the HCN(4-3)/HCO$^+$(4-3) ratio at the AGN position of NGC 7469 (1.11$\pm$0.06) is almost half of the corresponding value of the low-luminosity type-1 Seyfert galaxy NGC 1097 (2.0$\pm$0.2), despite the more than two orders of magnitude higher X-ray luminosity of NGC 7469. But the ratio is comparable to that of the close vicinity of the AGN of NGC 1068 ($\sim$ 1.5). Based on these results, we speculate that some other heating mechanisms than X-ray (e.g., mechanical heating due to AGN jet) can contribute significantly for shaping the chemical composition in NGC 1097.

We present first results from the SXDF-ALMA 1.5 arcmin^2 deep survey at 1.1 mm using Atacama Large Millimeter Array (ALMA). The map reaches a 1sigma depth of 55 uJy/beam and covers 12 Halpha-selected star-forming galaxies at z = 2.19 or z=2.53. We have detected continuum emission from three of our Halpha-selected sample, including one compact star-forming galaxy with high stellar surface density, NB2315-07. They are all red in the rest-frame optical and have stellar masses of log (M*/Msun)>10.9 whereas the other blue, main-sequence galaxies with log(M*/Msun)=10.0-10.8 are exceedingly faint, <290 uJy (2sigma upper limit). We also find the 1.1 mm-brightest galaxy, NB2315-02, to be associated with a compact (R_e=0.7+-0.1 kpc), dusty star-forming component. Given high gas fraction (44^{+20}_{-8}% or 37^{+25}_{-3}%) and high star formation rate surface density (126^{+27}_{-30} Msun yr^{-1}kpc^{-2}), the concentrated starburst can within less than 50^{+12}_{-11} Myr build up a stellar surface density matching that of massive compact galaxies at z~2, provided at least 19+-3% of the total gas is converted into stars in the galaxy centre. On the other hand, NB2315-07, which already has such a high stellar surface density core, shows a gas fraction (23+-8%) and is located in the lower envelope of the star formation main-sequence. This compact less star-forming galaxy is likely to be in an intermediate phase between compact dusty star-forming and quiescent galaxies.

We report the source size distribution, as measured by ALMA millimetric continuum imaging, of a sample of 13 AzTEC-selected submillimeter galaxies (SMGs) at z_photo ~ 3-6. Their infrared luminosities and star-formation rates (SFR) are L_IR ~ 2-6 x 10^12 L_sun and ~ 200-600 M_sun yr-1, respectively. The size of z ~ 3-6 SMGs ranges from 0".10 to 0".38 with a median of 0".20+0".03-0".05 (FWHM), corresponding to a median circularized effective radius (Rc,e) of 0.67+0.13-0.14 kpc, comparable to the typical size of the stellar component measured in compact quiescent galaxies at z ~ 2 (cQGs) --- R ~ 1 kpc. The median surface SFR density of our z ~ 3-6 SMGs is 100+42-26 M_sun yr-1 kpc-2, comparable to that seen in local merger-driven (U)LIRGsrather than in extended disk galaxies at low and high redshifts. The discovery of compact starbursts in z >~ 3 SMGs strongly supports a massive galaxy formation scenario wherein z ~ 3-6 SMGs evolve into the compact stellar components of z ~ 2 cQGs. These cQGs are then thought to evolve into the most massive ellipticals in the local Universe, mostly via dry mergers. Our results thus suggest that z >~ 3 SMGs are the likely progenitors of massive local ellipticals, via cQGs, meaning that we can now trace the evolutionary path of the most massive galaxies over a period encompassing ~ 90% of the age of the Universe.

We present an unprecedented measurement of the disc stability and local instability scales in the luminous infrared Seyfert 1 host, NGC7469, based on ALMA observations of dense gas tracers and with a synthesized beam of 165 x 132 pc. While we confirm that non-circular motions are not significant in redistributing the dense interstellar gas in this galaxy, we find compelling evidence that the dense gas is a suitable tracer for studying the origin of its intensely high-mass star forming ring-like structure. Our derived disc stability parameter accounts for a thick disc structure and its value falls below unity at the radii in which intense star formation is found. Furthermore, we derive the characteristic instability scale and find a striking agreement between our measured scale of ~ 180 pc, and the typical sizes of individual complexes of young and massive star clusters seen in high-resolution images.

We present 0.8-mm band molecular images and spectra obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) toward one of the nearest galaxies with an active galactic nucleus (AGN), NGC 1068. Distributions of CO isotopic species ($^{13}$CO and C$^{18}$O) $\it{J}$ = 3--2, CN $\it{N}$ = 3--2 and CS $\it{J}$ = 7--6 are observed toward the circumnuclear disk (CND) and a part of the starburst ring with an angular resolution of $\sim$1.$^{\prime\prime}$3 $\times$ 1.$^{\prime\prime}$2. The physical properties of these molecules and shock-related molecules such as HNCO, CH$_{3}$CN, SO, and CH$_{3}$OH detected in the 3-mm band were estimated using rotation diagrams under the assumption of local thermodynamic equilibrium. The rotational temperatures of the CO isotopic species and the shock-related molecules in the CND are, respectively, 14--22 K and upper limits of 20--40 K. Although the column densities of the CO isotopic species in the CND are only from one-fifth to one-third of that in the starburst ring, those of the shock-related molecules are enhanced by a factor of 3--10 in the CND. We also discuss the chemistry of each species, and compare the fractional abundances in the CND and starburst ring with those of Galactic sources such as cold cores, hot cores, and shocked molecular clouds in order to study the overall characteristics. We find that the abundances of shock-related molecules are more similar to abundances in hot cores and/or shocked clouds than to cold cores. The CND hosts relatively complex molecules, which are often associated with shocked molecular clouds or hot cores. Because a high X-ray flux can dissociate these molecules, they must also reside in regions shielded from X-rays.

The nearby Sy 1 galaxy NGC 1097 represents an ideal laboratory to explore the molecular chemistry in the presence and surroundings of an active galactic nucleus. Exploring the distribution of different molecular species allows us to understand the physical processes affecting the ISM both in the AGN vicinity as well as in the outer star forming molecular ring. We carried out 3 mm ALMA observations of HCN, HCO+, CCH, CS, HNCO, SiO, HC3N, and SO as well as the 13C isotopologues. All species were imaged over the central 2 kpc (~30") of the galaxy at a resolution of ~2.2"x1.5 (150 pc x 100 pc). HCO+ and CS appear to be slightly enhanced in the star forming ring. CCH, showing the largest variations across NGC 1097, is suggested to be a good tracer of both obscured and early stage star formation. HNCO, SiO and HC3N are significantly enhanced in the inner circumnuclear disk surrounding the AGN. Differences in the molecular abundances are observed between the star forming ring and the inner circumnuclear disk. We conclude that the HCN/HCO+ and HCN/CS differences observed between AGN dominated and starburst galaxies are not due to a HCN enhancement due to X-rays, but rather this enhancement is produced by shocked material at distances of 200 pc from the AGN. Additionally we claim the lower HCN/CS to be a combination of a small under-abundance of CS in AGNs together with excitation effects, where a high dense gas component (~10^6 cm^-3) may be more prominent in SB galaxies. However the most promising are the differences found among the dense gas tracers which, at our modest spatial resolution, seem to outline the physical structure of the molecular disk around the AGN. In this picture, HNCO probes the well shielded gas in the disk, surrounding the dense material moderately exposed to X-ray radiation traced by HC3N. Finally SiO might be the innermost molecule in the disk structure.

Sensitive observations with ALMA allow astronomers to observe the detailed distributions of molecules with relatively weak intensity in nearby galaxies. In particular, we report distributions of several molecular transitions including shock and dust related species ($^{13}$CO $J$ = 1--0, C$^{18}$O $J$ = 1--0, $^{13}$CN $N$ = 1--0, CS $J$ = 2--1, SO $J_N$ = 3$_2$--2$_1$, HNCO $J_{Ka,Kc}$ = 5$_{0,5}$--4$_{0,4}$, HC$_3$N $J$ = 11--10, 12--11, CH$_3$OH $J_K$ = 2$_K$--1$_K$, and CH$_3$CN $J_K$ = 6$_K$--5$_K$) in the nearby Seyfert 2 galaxy NGC 1068 observed with the ALMA early science program. The central $\sim$1 arcmin ($\sim$4.3 kpc) of this galaxy was observed in the 100 GHz region covering $\sim$96--100 GHz and $\sim$108--111 GHz with an angular resolution of $\sim4"\times2"$ (290 pc$\times$140 pc) to study the effects of an active galactic nucleus and its surrounding starburst ring on molecular abundances. Here, we present images and report a classification of molecular distributions into three main categories: (1) Molecules concentrated in the circumnuclear disk (CND) (SO $J_N$ = 3$_2$--2$_1$, HC$_3$N $J$ = 11--10, 12--11, and CH$_3$CN $J_K$ = 6$_K$--5$_K$), (2) Molecules distributed both in the CND and the starburst ring (CS $J$ = 2--1 and CH$_3$OH $J_K$ = 2$_K$--1$_K$), (3) Molecules distributed mainly in the starburst ring ($^{13}$CO $J$ = 1--0 and C$^{18}$O $J$ = 1--0). Since most of the molecules such as HC$_3$N observed in the CND are easily dissociated by UV photons and X-rays, our results indicate that these molecules must be effectively shielded. In the starburst ring, the relative intensity of methanol at each clumpy region is not consistent with those of $^{13}$CO, C$^{18}$O, and CS. This difference is probably caused by the unique formation and destruction mechanisms of CH$_3$OH.

We present the first 100 pc scale view of the dense molecular gas in the central ~ 1.3 kpc region of the type-1 Seyfert NGC 1097 traced by HCN (J=4-3) and HCO+ (J=4-3) lines afforded with ALMA band 7. This galaxy shows significant HCN enhancement with respect to HCO+ and CO in the low-J transitions, which seems to be a common characteristic in AGN environments. Using the ALMA data, we study the characteristics of the dense gas around this AGN and search for the mechanism of HCN enhancement. We find a high HCN (J=4-3) to HCO+ (J=4-3) line ratio in the nucleus. The upper limit of the brightness temperature ratio of HCN (v2=1^{1f}, J=4-3) to HCN (J=4-3) is 0.08, which indicates that IR pumping does not significantly affect the pure rotational population in this nucleus. We also find a higher HCN (J=4-3) to CS (J=7-6) line ratio in NGC 1097 than in starburst galaxies, which is more than 12.7 on the brightness temperature scale. Combined from similar observations from other galaxies, we tentatively suggest that this ratio appears to be higher in AGN-host galaxies than in pure starburst ones similar to the widely used HCN to HCO+ ratio. LTE and non-LTE modeling of the observed HCN and HCO+ lines using J=4-3 and 1-0 data from ALMA, and J=3-2 data from SMA, reveals a high HCN to HCO+ abundance ratio (5 < [HCN]/[HCO+] < 20: non-LTE analysis) in the nucleus, and that the high-J lines (J=4-3 and 3-2) are emitted from dense (10^{4.5} < n_H2 [/cc] < 10^6), hot (70 < Tkin [K] < 550) regions. Finally we propose that the high temperature chemistry is more plausible to explain the observed enhanced HCN emission in NGC 1097 than the pure gas phase PDR/XDR chemistry.

High resolution (0.4 arcsec) Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 0 observations of HCO+(4-3) and HCN(4-3) toward a mid-stage infrared bright merger VV114 have revealed compact nuclear (<200 pc) and extended (3 - 4 kpc) dense gas distribution across the eastern part of the galaxy pair. We find a significant enhancement of HCN(4-3) emission in an unresolved compact and broad (290km/s) component found in the eastern nucleus of VV114, and we suggest dense gas associated with the surrounding material around an Active Galactic Nucleus (AGN), with a mass upper limit of < 4 x 10^8 Msun. The extended dense gas is distributed along a filamentary structure with resolved dense gas concentrations (230pc; 10^6 Msun) separated by a mean projected distance of 600 pc, many of which are generally consistent with the location of star formation traced in Pa alpha emission. Radiative transfer calculations suggest moderately dense (10^5 - 10^6 cm^-3) gas averaged over the entire emission region. These new ALMA observations demonstrate the strength of the dense gas tracers in identifying both the AGN and star formation activity in a galaxy merger, even in the most dust enshrouded environments in the local universe.

We present a kinematic analysis of the dense molecular gas in the central 200 parsecs of the nearby galaxy NGC1097, based on Cycle 0 observations with the Atacama Large Millimeter/sub-millimeter Array (ALMA). We use the HCN(4-3) line to trace the densest interstellar molecular gas, and quantify its kinematics, and estimate an inflow rate for the molecular gas. We find a striking similarity between the ALMA kinematic data and the analytic spiral inflow model that we have previously constructed based on ionized gas velocity fields on larger scales. We are able to follow dense gas streaming down to 40 pc distance from the supermassive black hole in this Seyfert 1 galaxy. In order to fulfill marginal stability, we deduce that the dense gas is confined to a very thin disc, and we derive a dense gas inflow rate of 0.09 Msun/yr at 40 pc radius. Combined with previous values from the Ha and CO gas, we calculate a combined molecular and ionized gas inflow rate of 0.2 Msun/yr at 40 pc distance from the central supermassive black hole of NGC1097.