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 sensitive [\ion{C}{1}]~$^3P_1$--$^3P_0$ and $^{12}$CO~$J$=4--3 observations of the nearby merging galaxy NGC 6052 using the Morita (Atacama Compact) Array of ALMA. We detect $^{12}$CO~$J$=4--3 toward the northern part of NGC 6052, but [\ion{C}{1}]~$^3P_1$--$^3P_0$ is not detected with a [\ion{C}{1}]~$^3P_1$--$^3P_0$ to $^{12}$CO~$J$=4--3 line luminosity ratio of$~\lesssim0.07$. According to models of photodissociation regions, the unusual weakness of [\ion{C}{1}]~$^3P_1$--$^3P_0$ relative to $^{12}$CO~$J$=4--3 can be explained if the interstellar medium has a hydrogen density larger than $10^5\,{\rm cm}^{-3}$, conditions that might arise naturally in the ongoing merging process in NGC 6052. Its [\ion{C}{1}]~$^3P_1$--$^3P_0$ emission is also weaker than expected given the molecular gas mass inferred from previous measurements of $^{12}$CO~$J$=1--0 and $^{12}$CO~$J$=2--1. This suggests that [\ion{C}{1}]~$^3P_1$--$^3P_0$ may not be a reliable tracer of molecular gas mass in this galaxy. NGC 6052 is a unique laboratory to investigate how the merger process impacts the molecular gas distribution.

We investigate low-density accretion flows onto massive black holes (BHs) with masses of $\gtrsim 10^5~M_\odot$ orbiting around in the outskirts of their host galaxies, performing three-dimensional hydrodynamical simulations. Those wandering BHs are populated via ejection from the galactic nuclei through multi-body BH interactions and gravitational wave recoils associated with galaxy and BH coalescences. We find that when a wandering BH is fed with hot and diffuse plasma with density fluctuations, the mass accretion rate is limited at $\sim 10-20\%$ of the canonical Bondi-Hoyle-Littleton rate owing to a wide distribution of inflowing angular momentum. We further calculate radiation spectra from radiatively inefficient accretion flows onto the wandering BH and find that the spectra have a peak at the millimeter band, where the Atacama Large Millimeter/submillimeter Array (ALMA) has the highest sensitivity and spatial resolution. Millimeter observations with ALMA and future facilities such as the next generation Very Large Array (ngVLA) will enable us to hunt wandering BH population and push the detectable mass limit down to $M_\bullet \simeq 2\times10^7~M_\odot$ for massive nearby ellipticals, e.g., M87, and $M_\bullet \simeq 10^5~M_\odot$ for Milky Way. This radiation spectral model combined with numerical simulations will be applied to give physical interpretations of off-nuclear BHs detected in dwarf galaxies, which may constrain BH seed formation scenarios.

A galaxy-galaxy merger and the subsequent triggering of starburst activity are fundamental processes linked to the morphological transformation of galaxies and the evolution of star formation across the history of the Universe. Both nuclear and disk-wide starbursts are assumed to occur during the merger process. However, quantifying both nuclear and disk-wide star formation activity is non-trivial because the nuclear starburst is dusty in the most active merging starburst galaxies. This paper presents a new approach to this problem: combining hydrogen recombination lines in optical, millimeter, and free-free emission. Using NGC~3256 as a case study, H$\beta$, H40$\alpha$, and free-free emissions are investigated using the Multi Unit Spectroscopic Explorer at the Very Large Telescope of the European Southern Observatory (MUSE/VLT) and the Atacama Large Millimeter/submillimeter Array (ALMA). The H$\beta$ image obtained by MUSE identifies star-forming regions outside the nuclear regions, suggesting a disk-wide starburst. In contrast, the H40$\alpha$ image obtained by ALMA identifies a nuclear starburst where optical lines are undetected due to dust extinction ($A_{\rm V}\sim25$). Combining both MUSE and ALMA observations, we conclude that the total SFR is $49\pm2~M_{\odot}$~yr$^{-1}$ and the contributions from nuclear and disk-wide starbursts are $\sim34~\%$ and $\sim66~\%$, respectively. This suggests the dominance of disk-wide star formation in NGC~3256. In addition, pixel-by-pixel analyses for disk-wide star-forming regions suggest that shock gas tracers (e.g., CH$_3$OH) are enhanced where gas depletion time ($\tau_{\rm gas}$=$M_{\rm gas}/SFR$) is long. This possibly means that merger-induced shocks regulate disk-wide star formation activities.

We report the detection of a non-corotating gas component in a bright unlensed submillimeter galaxy at z=4.3, COSMOS-AzTEC-1, hosting a compact starburst. ALMA 0.17 and 0.09 arcsec resolution observations of [CII] emission clearly demonstrate that the gas kinematics is characterized by an ordered rotation. After subtracting the best-fit model of a rotating disk, we kinematically identify two residual components in the channel maps. Both observing simulations and analysis of dirty images confirm that these two subcomponents are not artificially created by noise fluctuations and beam deconvolution. One of the two has a velocity offset of 200 km/s and a physical separation of 2 kpc from the primary disk and is located along the kinematic minor axis of disk rotation. We conclude that this gas component is falling into the galaxy from a direction perpendicular to the disk rotation. The accretion of such small non-corotating gas components could stimulate violent disk instability, driving radial gas inflows into the center of galaxies and leading to formation of in-situ clumps such as identified in dust continuum and CO. We require more theoretical studies on high gas fraction mergers with mass ratio of 1:>10 to verify this process.

We present the results from ALMA observations of [NII]205 micron, [CII]158 micron, and [OII]88 micron lines in an unlensed submillimeter galaxy at z=4.3, COSMOS-AzTEC-1, hosting a compact starburst core with an effective radius of 1 kpc. The [CII] and [NII] emission are spatially-resolved in 0.3 arcsec-resolution (1 kpc in radius). The kinematic properties of the [NII] emission are consistent with those of the CO(4-3) and [CII] emission, suggesting that the ionized gas feels the same gravitational potential as the associated molecular gas and photodissociation regions (PDRs). On the other hand, the spatial extent is different among the lines and dust continuum: the [CII] emitting gas is the most extended and the dust is the most compact, leading to a difference of the physical conditions in the interstellar medium. We derive the incident far-ultraviolet flux and the hydrogen gas density through PDR modeling by properly subtracting the contribution of ionized gas to the total [CII] emission. The observed [CII] emission is likely produced by dense PDRs with nH(PDR)=10^(5.5-5.75) cm^-3 and G0=10^(3.5-3.75) in the central 1 kpc region and nH(PDR)=10^(5.0-5.25) cm^-3 and G0=10^(3.25-3.5) in the central 3 kpc region. We have also successfully measured the line ratio of [OIII/[NII] in the central 3 kpc region of COSMOS-AzTEC-1 at z=4.3, which is the highest redshift where both nitrogen and oxygen lines are detected. Under the most likely physical conditions, the measured luminosity ratio of L([OIII])/L([NII])=6.4+-2.2 indicates a near solar metallicity with Zgas=0.7-1.0 Zsol, suggesting a chemically evolved system at z=4.3.

We report ~2" resolution Atacama Large Millimeter/submillimeter Array observations of the HCN(1-0), HCO+(1-0), CO(1-0), CO(2-1), and CO(3-2) lines towards the nearby merging double-nucleus galaxy NGC 3256. We find that the high density gas outflow traced in HCN(1-0) and HCO+(1-0) emission is co-located with the diffuse molecular outflow emanating from the southern nucleus, where a low-luminosity active galactic nucleus (AGN) is believed to be the dominant source of the far-infrared luminosity. On the other hand, the same lines were undetected in the outflow region associated with the northern nucleus, whose primary heating source is likely related to starburst activity without obvious signs of AGN. Both HCO+(1-0)/CO(1-0) line ratio (i.e. dense gas fraction) and the CO(3-2)/CO(1-0) line ratio are larger in the southern outflow (0.20$\pm$0.04 and 1.3$\pm$0.2, respectively) than in the southern nucleus (0.08$\pm$0.01, 0.7$\pm$0.1, respectively). By investigating these line ratios for each velocity component in the southern outflow, we find that the dense gas fraction increases and the CO(3-2)/CO(1-0) line ratio decreases towards the largest velocity offset. This suggests the existence of a two-phase (diffuse and clumpy) outflow. One possible scenario to produce such a two-phase outflow is an interaction between the jet and the interstellar medium, which possibly triggers shocks and/or star formation associated with the outflow.

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.

We present CO J=4-3 line and 3 mm dust continuum observations of a 100 kpc-scale filamentary Ly{\alpha} nebula (SSA22 LAB18) at z=3.1 using the Atacama Large Millimeter/submillimeter Array (ALMA). We detected the CO J=4-3 line at a systemic z(CO)=3.093 {\pm} 0.001 at 11 {\sigma} from one of the ALMA continuum sources associated with the Ly{\alpha} filament. We estimated the CO J=4-3 luminosity of L'CO(4-3)=(2.3\pm0.2)x10^9 K km s^{-1} pc^2 for this CO source, which is one order of magnitude smaller than those of typical z>1 dusty star-forming galaxies (DSFGs) of similar far-infrared luminosity L(IR)~10^{12} Lsun. We derived a molecular gas mass of Mgas=(4.4^{+0.9}_{-0.6})x10^9 Msun and a star-formation rate of SFR=270\pm160 Msun yr^{-1}. We also estimated a gas depletion time of {\tau}(dep)=17\pm10 Myr, being shorter than those of typical DSFGs. It is suggested that this source is in a transition phase from DSFG to a gas-poor, early-type galaxy. From ALMA to Herschel multi-band dust continuum observations, we measured a dust emissivity index {\beta}=2.3\pm0.2, which is similar to those of local gas-poor, early-type galaxies. Such a high {\beta} can be reproduced by specific chemical compositions for interstellar dust at the submillimeter wavelengths from recent laboratory experiments. ALMA CO and multi-band dust continuum observations can constrain the evolutionary stage of high-redshift galaxies through {\tau}(dep) and {\beta}, and thus we can investigate dust chemical compositions even in the early Universe.

We present high-resolution observations (0".2-1".5) of multiple dense gas tracers, HCN and HCO$^+$ ($J$ = 1-0, 3-2, and 4-3), HNC ($J$ = 1-0), and CS ($J$ = 7-6) lines, toward the nearby luminous infrared galaxy VV 114 with the Atacama Large Millimeter/submillimeter Array. All lines are robustly detected at the central gaseous filamentary structure including the eastern nucleus and the Overlap region, the collision interface of the progenitors. We found that there is no correlation between star formation efficiency and dense gas fraction, indicating that the amount of dense gas does not simply control star formation in VV 114. We predict the presence of more turbulent and diffuse molecular gas clouds around the Overlap region compared to those at the nuclear region assuming a turbulence-regulated star formation model. The intracloud turbulence at the Overlap region might be excited by galaxy-merger-induced shocks, which also explains the enhancement of gas-phase CH$_3$OH abundance previously found there. We also present spatially resolved spectral line energy distributions of HCN and HCO$^+$ for the first time, and derive excitation parameters by assuming optically-thin and local thermodynamic equilibrium (LTE) conditions. The LTE model revealed that warmer, HCO$^+$-poorer molecular gas medium is dominated around the eastern nucleus, harboring an AGN. The HCN abundance is remarkably flat ($\sim$3.5 $\times$ 10$^{-9}$) independently of the various environments within the filament of VV 114 (i.e., AGN, star formation, and shock).

We present 0".97 $\times$ 0".53 (470 pc $\times$ 250 pc) resolution CO ($J$ = 2-1) observations toward the nearby luminous merging galaxy NGC 6240 with the Atacama Large Millimeter/submillimeter Array. We confirmed a strong CO concentration within the central 700 pc, which peaks between the double nuclei, surrounded by extended CO features along the optical dust lanes ($\sim$11 kpc). We found that the CO emission around the central a few kpc has extremely broad velocity wings with full width at zero intensity $\sim$ 2000 km s$^{-1}$, suggesting a possible signature of molecular outflow(s). In order to extract and visualize the high-velocity components in NGC 6240, we performed a multiple Gaussian fit to the CO datacube. The distribution of the broad CO components show four extremely large linewidth regions ($\sim$1000 km s$^{-1}$) located 1-2 kpc away from both nuclei. Spatial coincidence of the large linewidth regions with H$\alpha$, near-IR H$_2$, and X-ray suggests that the broad CO (2-1) components are associated with nuclear outflows launched from the double nuclei.

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 new single dish CO (3-2) emission data obtained toward 19 early stage and 7 late stage nearby merging galaxies using the Atacama Submillimeter Telescope Experiment (ASTE). Combining with the single dish and interferometric data of galaxies observed in previous studies, we investigate the relation between the CO (3-2) luminosity (L'CO(3-2)) and the far Infrared luminosity (LFIR) in a sample of 29 early stage and 31 late stage merging galaxies, and 28 nearby isolated spiral galaxies. We find that normal isolated spiral galaxies and merging galaxies have different slopes (alpha) in the log L'CO(3-2) - log LFIR plane (alpha ~ 0.79 for spirals and ~ 1.12 for mergers). The large slope (alpha > 1) for merging galaxies can be interpreted as an evidence for increasing Star Formation Efficiency (SFE=LFIR/L'CO(3-2)) as a function of LFIR. Comparing our results with sub-kpc scale local star formation and global star-burst activity in the high-z Universe, we find deviations from the linear relationship in the log L'CO(3-2) - log LFIR plane for the late stage mergers and high-z star forming galaxies. Finally, we find that the average SFE gradually increases from isolated galaxies, merging galaxies, and to high-z submillimeter galaxies / quasi-stellar objects (SMGs/QSOs). By comparing our findings with the results from numerical simulations, we suggest; (1) inefficient star-bursts triggered by disk-wide dense clumps occur in the early stage of interaction and (2) efficient star-bursts triggered by central concentration of gas occur in the final stage. A systematic high spatial resolution survey of diffuse and dense gas tracers is a key to confirm this scenario.