Resolving the interstellar medium in the nuclear region of two z=5.78 quasar host galaxies with ALMA

TL;DR

This study uses high-resolution ALMA observations to analyze the interstellar medium in the nuclear regions of two z=5.78 quasars, revealing complex gas dynamics, potential interactions, and insights into galaxy-black hole co-evolution.

Contribution

First detailed spatially resolved analysis of [CII] emission in high-redshift quasar host galaxies, linking gas dynamics to galaxy evolution models.

Findings

Decreasing [CII]-FIR ratios with increasing FIR surface brightness.

Detection of velocity gradients and turbulent gas clumps.

Host galaxy dynamical mass estimated at 2.6x10^10 Msun.

Abstract

We present ALMA observations of the [CII] 158 $\mu$m fine structure line and dust continuum emission from two quasars, SDSS J104433.04-012502.2 and SDSS J012958.51-003539.7, at z=5.78. The ALMA observations at 0.2'' resolution map the dust and gas on kpc scales. The spatially resolved emission shows a similar trend of decreasing [CII]-FIR ratios with increasing FIR surface brightnesses as was found in the infrared luminous galaxies with intense star formation. We confirm the velocity gradients of [CII] emission found previously in SDSS J0129-0035. No clear evidence of order motion is detected in SDSS J1044-0125. The velocity maps and PV-diagrams also suggest turbulent gas clumps in both objects. We tentatively detect a [CII] peak offset 4.9 kpc to the East of SDSS J1044-0125. This may be associated with an infalling companion, or node of gas outflow. All these results suggest significant dynamical evolution of the ISM in the nuclear region of these young quasar-starburst systems. We fit the velocity map of the [CII] emission from SDSS J0129-0035 with a rotating disk model. The result suggests a face-on system with an inclination angle of 16 +/- 20 degree and constrains the lower limit of the host galaxy dynamical mass to be 2.6x10^10 Msun within the [CII] emitting region. It is likely that SDSS J0129-0035, as well as other young quasars with super massive black hole masses on order of 10^7 to 10^8 Msun, falls close to the black hole and host galaxy mass relation defined by local galaxies.

Tables

Table 1. Table 1: Line and continuum parameters

	SDSS J0129$-$0035	SDSS J1044$-$0125
[C II] line flux ($\mathrm{Jy}\mathrm{km}{s}^{-1}$)	2.11$\pm$0.09	1.52$\pm$0.15
Peak ($\mathrm{Jy}{\mathrm{beam}}^{-1}\mathrm{km}{s}^{-1}$)	0.80$\pm$0.03	0.64$\pm$0.08
${\mathrm{FWHM}}_{\left[\mathrm{CII}\right]}$ ($\mathrm{km}{s}^{-1}$)	200$\pm$8	440$\pm$40
$[$C II$]$ central position (J2000)	${01}^h{29}^m{58.513}^s-{00}^{\circ }{35}^{\prime }{39.83}^{\prime \prime }$	${10}^h{44}^m{33.042}^s-{01}^{\circ }{25}^{\prime }{02.08}^{\prime \prime }$
Deconvolved FWHM [C II] size (imfit)	$\left(0^{\prime \prime }.32\pm 0^{\prime \prime }.03\right)\times \left(0^{\prime \prime }.27\pm 0^{\prime \prime }.03\right)$	$\left(0^{\prime \prime }.41\pm 0^{\prime \prime }.07\right)\times \left(0^{\prime \prime }.30\pm 0^{\prime \prime }.06\right)$
FWHM [C II] size (uvmodelfit)	$\left(0^{\prime \prime }.32\pm 0^{\prime \prime }.01\right)\times \left(0^{\prime \prime }.26\pm 0^{\prime \prime }.01\right)$	$\left(0^{\prime \prime }.32\pm 0^{\prime \prime }.02\right)\times \left(0^{\prime \prime }.26\pm 0^{\prime \prime }.03\right)$
$L_{\left[\mathrm{CII}\right]}\left(L_{\odot }\right)$	(1.96$\pm$0.08)$\times {10}^9$	(1.41$\pm$0.14)$\times {10}^9$
Continuum (mJy)	2.82$\pm$0.04	3.02$\pm$0.11
Continuum central position (J2000)	${01}^h{29}^m{58.513}^s-{00}^d{35}^{\prime }{39.84}^{\prime \prime }$	${10}^h{44}^m{33.040}^s-{01}^d{25}^{\prime }{02.10}^{\prime \prime }$
Deconvolved FWHM continuum size (imfit)	$\left(0^{\prime \prime }.18\pm 0^{\prime \prime }.01\right)\times \left(0^{\prime \prime }.16\pm 0^{\prime \prime }.01\right)$	$\left(0^{\prime \prime }.16\pm 0^{\prime \prime }.02\right)\times \left(0^{\prime \prime }.15\pm 0^{\prime \prime }.02\right)$
FWHM continuum size (uvmodelfit)	$\left(0^{\prime \prime }.172\pm 0^{\prime \prime }.003\right)\times \left(0^{\prime \prime }.158\pm 0^{\prime \prime }.005\right)$	$\left(0^{\prime \prime }.159\pm 0^{\prime \prime }.005\right)\times \left(0^{\prime \prime }.154\pm 0^{\prime \prime }.008\right)$
Optical quasar position (J2000)	${01}^h{29}^m{58.524}^s-{00}^d{35}^{\prime }{39.77}^{\prime \prime }$	${10}^h{44}^m{33.041}^s-{01}^d{25}^{\prime }{02.09}^{\prime \prime }$
Cycle 0 [C II] line fluxa ($\mathrm{Jy}\mathrm{km}{s}^{-1}$)	1.99$\pm$0.12	1.70$\pm$0.30
Cycle 0 Continuum (mJy)	2.57$\pm$0.06	3.12$\pm$0.09

Table 2. Table 2: Spectral fitting parameters for SDSS J1044-0125

Component	Center	FWHM	Flux
	$\mathrm{km}{s}^{-1}$	$\mathrm{km}{s}^{-1}$	$\mathrm{Jy}\mathrm{km}{s}^{-1}$
A	-230$\pm$18	103$\pm$40	0.21$\pm$0.09
B	-39$\pm$12	180$\pm$43	0.76$\pm$0.15
C	194$\pm$19	160$\pm$44	0.42$\pm$0.11

Full text

galaxies: active — galaxies: evolution — galaxies: high-redshift

1 Introduction

Quasars at high redshifts provide the best opportunity to study the formation and early evolution of the first supermassive black holes (SMBHs) and their host galaxies. Due to the great success of wide-field optical and near-infrared surveys in the last two decades, more than two hundred quasars have been discovered at z$\geq$5.7, an epoch close to the end of cosmic reionization (e.g. Fan et al., 2006; Mortlock et al., 2011; Venemans et al., 2013; Bañados et al., 2016, 2018; Jiang et al., 2009, 2016; Matsuoka et al., 2016, 2018a, 2018b; Chehade et al., 2018). This quasar sample covers a wide range of SMBH masses and quasar luminosities, including the most massive sources with SMBH mass above $\rm 10^{10}\,M_{\odot}$ (e.g., Fan et al., 2000; Wu et al., 2015) as well as objects with SMBH masses of a few $\rm 10^{7}$ to $\rm 10^{8}\,M_{\odot}$ (Willott et al., 2010; Jiang et al., 2009; Matsuoka et al., 2016, 2018a, 2018b) which are more similar to the common quasar population observed at lower redshifts.

The dust and gas in the host galaxies of these earliest quasars have been studied at submillimeter and millimeter wavelengths using thermal continuum emission, molecular (mostly CO), and fine structure lines, such as the [C II] 158$\mu$m emission line (Walter et al., 2004, 2009; Bertoldi et al., 2003a, b; Petric et al., 2003; Priddey et al., 2003; Robson et al., 2004; Riechers et al., 2009; Omont et al., 2013; Wang et al., 2008, 2010, 2013; Willott et al., 2015, 2017; Bañados et al., 2015; Venemans et al., 2016, 2017a, 2017b, 2017c; Decarli et al., 2017, 2018; Izumi et al., 2018). Dust continuum surveys at $\sim$0.5-1 mJy rms sensitivity have detected the most FIR-luminous objects, with far-infrared (FIR) luminosities of a few $\rm\times 10^{12}$ to $\rm 10^{13}\,L_{\odot}$ and dust masses of $\rm>10^{8}\,M_{\odot}$ (Bertoldi et al., 2003a; Priddey et al., 2003; Robson et al., 2004; Wang et al., 2007, 2008). These FIR luminous quasars are bright in CO line emission and [C II] 158$\mu$m fine structure line emission, indicating molecular gas with masses of order of $\rm 10^{10}\,M_{\odot}$ and bursts of star formation in the nuclear region, co-eval with the SMBH accretion and quasar activity (Bertoldi et al., 2003b; Maiolino et al., 2005; Carilli et al., 2007; Wang et al., 2013; Decarli et al., 2017, 2018; Venemans et al., 2016).

In recent years, the Atacama Large Millimeter/submillimeter Array (ALMA) has carried out comprehensive surveys of the [C II] 158$\mu$m fine structure line in high-z quasars. For example, Decarli et al. (2018) detected [C II] in 85% of 27 optically-selected quasars at z$>$5.94. This line is an important coolant that traces the ionized and neutral interstellar medium (ISM) and star forming activities (Herrea-Camus et al., 2015). The detection of [C II] line and dust continuum emission has revealed a wide range of star formation rates, from a few 10 $\rm M_{\odot}\,yr^{-1}$ to $\rm\geq 1000\,M_{\odot}\,yr^{-1}$ (Decarli et al., 2018; Izumi et al., 2018). The [C II]-FIR luminosity ratios of these quasar hosts range from $\rm 10^{-4}$ to a few $\rm\times 10^{-3}$ (Maiolino et al., 2005; Wang et al., 2013; Willott et al., 2017; Decarli et al., 2018), following the trend of decreasing [C II]-FIR ratio with increasing FIR luminosities found with the IR luminous star-forming systems (Malhotra et al., 2001; Stacey et al., 2010; Luhman et al., 2003; Díaz-Santos et al., 2013, 2017; Graciá-Carpio et al., 2011; Hailey-Dunsheath et al., 2010; Muñoz & Oh, 2016; Smith et al., 2017; Gullberg et al., 2018; Decarli et al., 2018; Rybak et al., 2019). In particular, the quasar hosts with high FIR luminosities of a few $\rm 10^{12}$ to $\rm 10^{13}\,L_{\odot}$ show low [C II]-FIR ratios similar to that found in the Ultraluminous infrared galaxies (ULIRGs) and submillimeter galaxies (SMGs) (Luhman et al., 2003; Díaz-Santos et al., 2013; Wang et al., 2013; Díaz-Santos et al., 2017; Rybak et al., 2019). The [C II]-FIR ratio is suggested to be related to the local conditions of the ISM; e.g. a decrease of [C II]-FIR ratio could be due to a high gas temperature (i.e., $\rm T>\Delta E/k\sim 91\,K$, $\rm\Delta$E is the energy separation of the two levels of the [C II] 158 $\mu$m transition) where the upper level of [C II] transition is saturated (Muñoz & Oh, 2016), or a high gas surface density (e.g., the typical density in compact ULIRGs and SMGs) where more gas is in molecular phase (Narayanan & Krumholz, 2017). It is interesting to see how the resolved distributions of the line and continuum surface brightnesses and emission ratios in these FIR luminous quasar hosts compare to that in the ULIRGs and SMGs (Smith et al., 2017; Gullberg et al., 2018; Rybak et al., 2019).

The [C II] 158$\mu$m and molecular CO lines are bright tracers of the gas content that could be detected from distant galaxies (Solomon & Vanden Bout, 2005; Carilli & Walter, 2013). Measurements of the line widths, profiles, and velocity maps of these emission lines reveal important information on the gas kinematics (Walter et al., 2004; Ivison et al., 2011; Maiolino et al., 2012; Wang et al., 2013; Venemans et al., 2016, 2019). Observations of the [C II] and CO emission lines from the young quasar hosts at z$\geq$5.7 at sub-arcsecond to arcsecond resolution show a range of kinematic properties, including velocity gradients of ordered motion (Willott et al., 2013, 2017; Venemans et al., 2016; Shao et al., 2017; Feruglio et al., 2018), large velocity dispersion/turbulence (e.g., SDSS J231038.88+185519.7, Feruglio et al., 2018), gas outflows (e.g., SDSS J114816.64+515250.3, Maiolino et al., 2012; Cicone et al., 2015), and very compact sources with no evidence of rotation (e.g., ULAS J112001.48+064124.3, Venemans et al., 2017c). All these submillimeter/millimeter (submm/mm) observations argue for an early phase of SMBH-galaxy co-evolution in the z$\geq$5.7 quasar sample. Deep imaging of the gas component at kpc scale is required to fully resolve the gas distribution and kinematic properties in these systems.

In this paper, we present the observations, results and conclusion for two out of three quasars from our ALMA Cycle 1 program. Our team has carried out ALMA observations at $\rm 0^{\prime\prime}.2-0^{\prime\prime}.3$ angular resolution to image the [C II] line emission from the host galaxies of three quasars, ULAS J131911.29+095051.4 (hereafter ULAS J1319+0950), SDSS J104433.04-012502.2 (hereafter SDSS J1044-0125), and SDSS J012958.51-003539.7 (hereafter SDSS J0129-0035) at z$\geq$5.7. These three objects are among the brightest dust continuum, and [C II] detections at z$\sim$6 (Wang et al., 2013), with FIR luminosities of $\rm 5\sim 10\times 10^{12}\,L_{\odot}$ 111We adopt the rest-frame wavelength range of 42.5$\mu$m to 122.5$\mu$m for FIR luminosity calculation and [C II] luminosities of $\rm 1.5\sim 4.4\times 10^{9}\,L_{\odot}$ (Wang et al., 2013), suggesting massive star formation in the quasar hosts. The three objects are different in quasar UV luminosities and SMBH masses. SDSS J1044-0125 and ULAS J1319+0950 are among the most optically luminous quasars at z$\sim$6 with rest-frame 1450$\AA$ apparent magnitudes of $\rm m_{1450}$=19.21 and 19.65, respectively (Fan et al., 2000; Mortlock et al., 2009). The SMBH mass is $\rm(2.7\pm 0.6)\times 10^{9}\,M_{\odot}$ for ULAS J1319+0950 (Shao et al., 2017) and $\rm(5.6\pm 0.6)\times 10^{9}\,M_{\odot}$ for SDSS J1044-0125 (Shen et al., 2019). SDSS J0129-0035 is much fainter in quasar emission with $\rm m_{1450}$=22.28 (Jiang et al., 2009). If Eddington accretion is assumed (Willott et al., 2010), the SMBH mass is $\rm 1.7\times 10^{8}\,M_{\odot}$. These objects provide a ideal sample for a pilot study of the ISM distribution and dynamics in the earliest SMBH-starburst systems over a range of quasar luminosities.

We report our ALMA [C II] line imaging of ULAS J1319+0950 at $\rm 0.3^{\prime\prime}$ resolution in Shao et al. (2017). The [C II] velocity map of this object shows clear velocity gradients and the spectrum appears to be broad ($\rm FWHM\sim 540\,km\,s^{-1}$) and flat at the top, suggesting that the [C II]-emitting gas is mainly distributed in an inclined rotating disk extending to a radius of $\sim$3.2 kpc. We performed dynamical modeling of the gas velocity field and measured the disk inclination angle and rotation velocity of the atomic/ionized gas in the nuclear region. The derived host galaxy dynamical mass is $\rm 1.3\times 10^{11}\,M_{\odot}$, suggesting a black hole-to-host galaxy mass ratio of $\sim 0.02$, i.e., about four times higher than the value expected from the black hole-bulge mass relation of local galaxies (Kormendy & Ho, 2013).

We here present ALMA [C II] imaging of the remaining two quasars, SDSS J0129-0035 and SDSS J1044-0125. We describe the observations in Section 2, present the results in Section 3, discuss the surface brightnesses of the [C II] and dust emission and gas kinematic properties in Section 4, and summarize the main conclusions in Section 5. We adopt a $\Lambda$CDM cosmology with $\rm H_{0}=71km\,s^{-1}\,Mpc^{-1}$, $\rm\Omega_{M}=0.27$ and $\rm\Omega_{\Lambda}=0.73$ throughout this paper (Spergel et al., 2007).

2 Observations

We carried out observations of the [C II] 158 $\mu$m emission line from SDSS J0129$-$0035 and SDSS J1044$-$0125 using the Band 7 receiver on ALMA in Cycle 1 (Program ID: 2012.1.00240.S). The data was taken in 2015 July and August in the C34-6 configuration with 32 12-m diameter antennas and a maximum baseline of 1.6 km. This results in a full width at half maximum (FWHM) synthesized beam size of $\sim$0.2*′′* using robust weighting, which corresponds to 1.2 kpc at z$\sim$5.8. We tuned the four 1.785 GHz spectral windows with one window centered on the redshifted [C II] line frequencies of the two targets and the other three windows observing the continuum. The correlator channel width is 15.625 MHz, corresponding to a velocity resolution of $\sim$17 $\rm km\,s^{-1}$. We checked the phase every 5$\sim$7 minutes on nearby calibrators, and the flux scale was calibrated on Ceres and J1058+015 with calibration uncertainties less than 10%. The total on-source time is about 80 min for SDSS J0129-0035 and 75 min for SDSS J1044-0125.

The data were calibrated using the pipeline of the Common Astronomy Software Application (CASA) Version 4.3.1. We then combined the data with that obtained from previous ALMA Cycle 0 observations (Program ID: 2011.0.00206.S, Wang et al., 2013). The ALMA Cycle 0 data alone has a lower angular resolution of 0.6*′′-0.7′′*, and 1$\sigma$ rms sensitivities of 0.35 mJy $\rm beam^{-1}$ for SDSS J0129-0035 and 0.65 mJy $\rm beam^{-1}$ for SDSS J1044-0125 in each 67 $\rm km\,s^{-1}$ channel (Wang et al., 2013). By including the Cycle 0 data, we could better sample the uv plane at different scale and improve the sensitivity at short uv-distance. As the Cycle 0 data were calibrated using a much earlier version of CASA, we corrected the weights before combining with the Cycle 1 data222https://casaguides.nrao.edu/index.php/DataWeightsAndCombination. We checked the amplitudes of visibility data from the two observations at overlapping uv-distance, and the flux scales are consistent with each other within the calibration uncertainties. We averaged, in the visibility domain, the emission in line-free channels to form a ’pseudo’ continuum dataset which we then subtracted from the uv-data to obtain a line-only dataset. We made the continuum image with the data from three line-free spectral windows and made the line image cubes with the continuum-subtracted visibility dataset. The images were cleaned using the CLEAN task in CASA using robust=0.5 for weighting. The typical rms noise is 0.20 mJy $\rm beam^{-1}$ per 17 $\rm km\,s^{-1}$ channel for the final image cube of SDSS J0129$-$0035, and 0.36 mJy $\rm beam^{-1}$ per 17 $\rm km\,s^{-1}$ channel for SDSS J1044$-$0125. The rms noise of the continuum is 21$\rm\mu Jy\,beam^{-1}$ for SDSS J0129$-$0035, and 40$\rm\mu Jy\,beam^{-1}$ for SDSS J1044$-$0125.

3 Results

For each of the two objects, we integrated the [C II] emission over the line-emitting channels to obtain velocity-integrated maps and averaged the line-free channels to obtain the continuum maps (Figure 1). We then integrated the emission within the 2$\sigma$ contour region to get total [C II] line fluxes and continuum flux densities (Table 1). For SDSS J0129-0035, the [C II] velocity-integrated map has a peak of $\rm 0.80\pm 0.03\,Jy\,beam^{-1}\,km\,s^{-1}$ (signal-to-noise ratio, $\rm SNR=23.5$), which is about 38% of the total line flux. The continuum intensity map has a peak of $\rm 1.58\pm 0.02\,mJy\,beam^{-1}$ ($\rm SNR=80$), 56% of the integrated continuum flux density. For SDSS J1044-025, the [C II] peak is $\rm 0.64\pm 0.08\,Jy\,beam^{-1}\,km\,s^{-1}$ (i.e., $\rm SNR=8.4$). This is about 45% of the total line flux. The peak of the continuum emission is $\rm 1.88\pm 0.04\,mJy\,beam^{-1}$ ($\rm SNR=47.5$), which is about 63% of the total continuum flux density. Thus, we consider the [C II] and continuum emission from both objects are resolved.

The visibilities versus uv-distance for the [C II] line and continuum are plotted in Figure 2. The data are radially averaged with a bin size of 50 k$\lambda$ for SDSS J0129-0035 and 75 k$\lambda$ for SDSS J1044-0125. The decrease of the visibility amplitudes toward longer uv-distance indicates that the line and continuum sources are resolved in the uv-plane. We performed elliptical Gaussian fitting to measure the sizes of the [C II] and continuum emission in both the image plane and the uv plane, using imfit and uvmodelfit tasks in CASA, respectively. The results, together with the integrated line fluxes and continuum flux densities, are summarized in Table 1. The uv-models are also plotted in Figure 2. The Gaussian FWHM source sizes (deconvolved from the synthesized beam) and peak positions in the image plane are consistent with that from the uv-models within the fitting uncertainties. We adopt the fitting positions and deconvolved source sizes from imfit for the discussions in the rest of this paper.

The fitting uncertainties of the [C II] peak positions are $\sim$7 mas in both RA and DEC for SDSS J0129-0035, and 25 mas for SDSS J1044-0125. The positional uncertainties for the continuum images are $<$3 mas for both objects. Thus, the peak positions of the [C II] line and continuum emission are consistent with each other within the fitting errors. We checked the coordinates of the optical quasars using the astrometric data from the Gaia mission (Lindegren et al., 2018). For SDSS J0129-0035, the quasar coordinates are RA=01h29m58.524s, and DEC=-00d35*′39.77′′* with an uncertainty of 149 mas in RA and 114 mas in DEC. Considering the coordinates uncertainties, the positions of [C II] and dust continuum peaks are marginally consistent with that of the optical quasar (Figure 1). The quasar coordinates for SDSS 1044-0125 are RA=10h44m33.041s, and DEC=-01d25*′02.09′′*, with an uncertainty of 27 mas in RA and 20 mas in DEC. This is well consistent with the positions of the [C II] and dust continuum peaks of this object (Figure 1).

We finally check the results by comparing the measurements to those from the low-resolution ALMA Cycle 0 data (Table 1, Wang et al., 2013). The [C II] line fluxes and continuum flux densities integrated over the source regions of both objects are consistent within the calibration uncertainties with the previous values. The ALMA data at $\rm 0.2^{\prime\prime}$ resolution reveal more detail on the spatial distribution and velocity field of the [CII]-emitting gas component in the quasar host galaxies.

SDSS J0129-0035. We present the intensity-weighted velocity map and velocity dispersion map of the [C II] line emission in Figure 3. We adopt the [C II] redshift of $\rm z_{[CII]}=5.7787$ as the zero velocity (Wang et al., 2013). The line is single-peaked and the velocity map shows a velocity gradient from the southeast to northwest direction. The velocity images also show two features at the edge of the [C II]-emitting region (marked as No. 1 and 2 in the velocity map, Figure 3). One has an intensity-weighted velocity of $\rm-190\,km\,s^{-1}$ at the southeast edge of the line emitting region and the other has an intensity-weighted velocity of about $\rm+50\,km\,s^{-1}$ on the northwest part. There is no strong ($\rm\geq 2\sigma$) dust emission detected in the regions of these two features. The intensity-weighted velocity dispersion map of the [C II] emission (upper right panel of Figure 3) also suggests a high velocity dispersion structure along the southeast-to-northwest direction, which seems to extend from the center to the two features. A structure with such a velocity dispersion cannot be explained with beam smearing of a rotating disk.

The PV diagrams along the direction of the two features and the perpendicular direction are plotted in the lower panel of Figure 3. The major gas component in the PV diagrams is consistent with rotation. The two features mentioned above are presented at the regions around (-0.4*′′, $\rm-190\,km\,s^{-1}$, No. 1 in Figure 3) and (+0.4′′, $\rm+50\,km\,s^{-1}$, No. 2) on the bottom left panel. It is unclear if these two components are connected to the outer part of the rotating gas disk or not. The PV-diagrams also shows tentative peaks ($\rm 3\sim 4\sigma$) with velocities of $\rm\sim\pm 200\,km\,s^{-1}$ within the central 0.2′′* (No. 3 and 4 in Figure 3). These components locate close to the center/peak of the [C II]-emitting region. However, the large velocities of $\rm\sim\pm 200\,km\,s^{-1}$ seem not follow the velocity field the rotating component. These tentative peaks should be checked with deeper observations at better resolution (e.g., $\lesssim 0.1^{\prime\prime}$), to confirm if they are real features and related to gas infow/outflows.

We plot the [C II] line spectrum integrated over the line-emitting region in Figure 4, compared to the CO (6-5) line spectrum published in Wang et al. (2011). Though the CO (6-5) line profile is poorly constrained due to the low signal-to-noise ratio, the line widths and redshifts measured with the [C II] and CO lines are consistent with each other within the uncertainties (Wang et al., 2011, 2013). Moreover, comparing to the Gaussian line profile that fit to the line spectrum, there is a tentative excess on the blue wing around $\rm-200\,km\,s^{-1}$ (arrow in Figure 4). This is likely to be related to the components saw in the velocity map and PV diagrams (e.g., No. 1 and 3).

SDSS J1044-0125. We calculate the intensity-weighted velocity and velocity dispersion maps with the $\rm>4\sigma$ pixels in each [C II] channel (Figure 5), which reveals the gas kinematics in the central high surface brightness part of the line-emitting region. The gas does show some velocity structure over this area. For example, the gas component to the northwest of the quasar position shows more positive velocities. However, the velocity map seems more complicated than what one would expect from a simple rotating disk. The PV-diagrams also suggest a very turbulent gas velocity field; the gas is distributed over a compact region of about 2.4 kpc (i.e., $\rm\pm 0.2^{\prime\prime}$) with velocities from -200 to $\rm\gtrsim 200\,km\,s^{-1}$.

The integrated [C II] line spectrum shows multiple peaks and is not well fit by a single Gaussian profile (Figure 6). We try to fit the line profile with three Gaussian components (A, B, and C). The fitting results of the line centers, FWHM widths, and fluxes are listed in Table 2. Here, we adopt the [C II] redshift of $\rm z_{[CII]}=5.7847$ from the single Gaussian fitting as the zero velocity (Wang et al. 2013). Wang et al. (2013) found a velocity offset between the Gaussian-fit line centers of the [C II] line and the CO (6-5) line for SDSS J1044-0125. With the new ALMA [C II] data in this work, we still see this offset. i.e., The CO (6-5) line overlaps only the blue part ($\rm-300-0\,km\,s^{-1}$) of the [C II] line emission with a line center at around $\rm-100\pm 44\,km\,s^{-1}$ (Wang et al., 2013) which is between the central velocities of components A and B. The mismatch between the [C II] and CO (6-5) spectra argues for complex gas components and kinematic properties in the nuclear region of SDSS J1044-0125 (see more discussion in Section 4.2).

On the velocity integrated map of the [C II] line (lower left of Figure 1 and left of Figure 7), there is a tentative peak $\rm 0.82^{\prime\prime}$ (i.e., 4.9 kpc) to the East of the major component. The spectrum of the emission from this secondary peak position shows positive signal from about -740 to $\rm 400\,km\,s^{-1}$. We integrate the emission in this velocity range and show the intensity map in the left panel of Figure 7. The secondary peak is detected at about 5$\sigma$ on the image, with a flux of $\rm 0.52\pm 0.10\,Jy\,km\,s^{-1}$. We also check the Cycle 0 data and the image at $\rm 0^{\prime\prime}.6$ resolution does show an extension toward the secondary peak location (see red contours in the upper left panel of Figure 7). A further check shows that this component is also present in a continuum-unsubtracted image in the same velocity range, as well as the image using only Cycle 1 data. No emission is detected at this location with the data from the other three continuum spectral windows. Fitting this component with a Gaussian profile yields a FWHM line width of $\rm 750\pm 190\,km\,s^{-1}$, centered at $\rm-190\pm 80\,km\,s^{-1}$. The [C II] luminosity derived with the peak intensity of this component is $\rm(4.8\pm 0.9)\times 10^{8}\,L_{\odot}$.

4 Analysis and Discussion

4.1 Surface brightnesses and ratios of the resolved [C II] and dust emission

The ALMA [C II] images of SDSS J0129-0035 and SDSS J1044-0125 probe the ISM in the nuclear regions on kpc scales in these quasar-starburst systems. As we described in the previous section, the line and continuum emission are resolved from the two objects at $\rm 0.2^{\prime\prime}$ (i.e., $\rm\sim 1.2$ kpc) resolution. The deconvolved FWHM source sizes of the [C II] line emission are 1.8$\pm$0.2 and 2.6$\pm$0.5 times larger than that of the continuum for SDSS J0129-0035 and SDSS J1044-0125, respectively (based on the results from imfit, Table 1). This is consistent with the results found in previous ALMA imaging of other starburst quasar host galaxies and SMGs (Wang et al., 2013; Díaz-Santos et al., 2016; Venemans et al., 2016, 2017c; Gullberg et al., 2018; Cooke et al., 2018; Rybak et al., 2019). In most of these systems, the [C II] emission is found to be 1.3 to $>$2 times more extended than that of the dust continuum emission (Gullberg et al., 2018; Cooke et al., 2018).

Based on the [C II] and dust continuum intensity maps, we compare the surface brightnesses as a function of radius in the nuclear region of the two systems (Figure 8). We also include ULAS J1319+0950 in the comparison, which is another FIR luminous quasar-starburst system and spatially resolved in [C II] and dust continuum emission at a similar resolution in our ALMA Cycle 1 program (Shao et al., 2017,see Section 1). Here we scale the continuum intensity map (Figure 1) with a factor of $\rm L_{FIR}/S_{con}$, where $\rm L_{FIR}$ is the FIR luminosity derived in Wang et al. (2013) assuming dust temperature of 47 K and emissivity index of $\rm\beta=1.6$ (Beelen et al., 2006), and $\rm S_{con}$ is the total continuum flux density close to the [C II] frequency in Table 1. This convert the monochromatic continuum map to a map of surface brightness of the FIR emission. We calculate the mean [C II] and FIR surface brightnesses within the rings at central radii of 0*′′, 0.1′′, 0.2′′, 0.3′′, 0.4′′, 0.5′′, 0.6′′* and a ring width of 0.1*′′. The scale of 0.1′′* corresponds to a physical size of 0.58 kpc for ULAS J1319+0950, and 0.60 kpc for SDSS J0129$-$0035 and SDSS J1044$-$0125. The symbols and error bars shown in Figure 8 represent the mean and standard deviation of the pixel values in each ring. The [C II]-to-FIR surface brightness ratios are also plotted in the right panel of Figure 8. We model the surface brightness with an exponential disk in the 2-dimensional image plane with a Sérsic index of n=1. We then convolve the model with the synthesized beam, calculate the mean surface brightnesses of the model in each ring described above, and fit these mean values to the data. The models are shown as dashed lines in Figure 8. They suggest that the [C II] and dust emission in the nuclear region of the quasar host galaxies follow an exponential light profile. The resolved surface brightness ratios between the [C II] and FIR dust emission range from 0.0001 to 0.001 within the central region where both [C II] and dust continuum are detected.

In Figure 9, we compare the resolved [C II] and FIR continuum emission at different radii in these three z$\sim$6 quasars with other systems on the [C II]-FIR ratio vs. FIR surface brightness plot. We collect samples of star forming galaxies and quasars at low and high redshifts that have continuum source size measurements from ALMA, including the GOALS sample of local luminous infrared galaxies (Díaz-Santos et al., 2013, 2017; Shangguan et al., 2019), the high-redshift SMGs (Riechers et al., 2013; Neri et al., 2014; Gullberg et al., 2018), Lyman Break Galaxies (LBG, Capak et al., 2015; Jones et al., 2017; Hashimoto et al., 2018), and the z$\sim$6 to 7 quasars (Wang et al., 2013; Venemans et al., 2016, 2017c; Díaz-Santos et al., 2016; Decarli et al., 2018). For the GOALS sample, we use the [C II]-FIR ratio extract from the central a few kpc region of the galaxies (i.e., the $\rm 9.4^{\prime\prime}\times 9.4^{\prime\prime}$ region of one PACS detector element where the [C II] and continuum emission peaks, see details in Díaz-Santos et al., 2013, 2017). This is comparable to the [C II] and dust emission region detected in the high-z samples. We re-calculate the FIR luminosity surface brightness based on the GOALs sample catalog333http://goals.ipac.caltech.edu/ and the decomposed infrared SED from Shangguan et al. (2019). For the high-redshift samples that have a source size measurements from ALMA at $\rm 0.6^{\prime\prime}\sim 0.7^{\prime\prime}$ resolution, we adopt the FIR surface brightness averaged within the half-light radius of the continuum emitting region. We also include the resolved [C II] and FIR emission from two submillimeter galaxies at z=2.9, which were imaged with ALMA at $0.16^{\prime\prime}$ resolution and have FIR luminosities and source sizes comparable to the two quasar host galaxies we studied in this work (Rybak et al., 2019).

The continuum emission from SDSS J0129$-$0035 and SDSS J1044$-$0125, as well as ULAS J1319+0950, was resolved by ALMA on 1$\sim$2 kpc scale and the peak FIR surface brightness is up to $\rm\Sigma_{FIR}\sim 2\times 10^{12}\,L_{\odot}\,kpc^{-2}$. This is comparable to the highest values found in the compact SMGs and ULIRGs. As was shown in Figure 9, at this FIR surface brightness, the [C II]-FIR ratio is typically an order of magnitude lower than that of normal star-forming galaxies. The [C II]-FIR ratios in the three quasar hosts get higher toward larger radii, and the FIR surface brightnesses decrease to a few $\rm 10^{10}\,L_{\odot}\,kpc^{-2}$. Thus the radial change of the [C II]-FIR ratios found with the three objects can be understand as a trend of descreasing [C II]-FIR ratio with increasing FIR surface brightness which was found in samples of IR luminous star forming systems (e.g., $\rm\Sigma_{FIR}\,\geq$ a few $\rm 10^{10}\,L_{\odot}\,kpc^{-2}$, Herrea-Camus et al., 2018). This trend have been termed ”[C II]-FIR deficit” and was discussed in a number of papers (Stacey et al., 2010; Hailey-Dunsheath et al., 2010; Ivison et al., 2010; Farrah et al., 2013; Gullberg et al., 2015, 2018; Díaz-Santos et al., 2016; Lutz et al., 2016; Smith et al., 2017; Herrea-Camus et al., 2018; Cooke et al., 2018).

In the compact dusty starburst systems, the high FIR surface brightness reveals high star formation rate surface density with high gas density and strong UV radiation field. As was mentioned in Section 1, in such environment, the [C II] emission could be suppressed due to several possible mechanisms, including e.g., saturation of the [C II] transition when the gas is heated with temperatures above 91 K, (Muñoz & Oh, 2016; Rybak et al., 2019), increase of dust absorption of UV photoms with high ionization parameter and a decrease of photoelectric heating efficiency (Herrea-Camus et al., 2018), or reduced amount of C+ medium with more carbon in form of CO at high cloud surface density (Narayanan & Krumholz, 2017). Though a model study of the origin of the [C II]-FIR deficit is beyound the goals of this work, we can conclude that the low [C II]-FIR ratio and high FIR surface brightness found in the central region of the FIR luminous quasar hosts suggest a similar high density and dusty ISM heated by a high-intensity radiation field as was found in the compact and intense starburst systems.

Of the three objects, SDSS J1044$-$0125 and ULAS J1319+0950 are about 10 times brighter than SDSS J0129-0035 in quasar UV luminosity (Jiang et al., 2009; Shen et al., 2019). However, their [C II]/FIR surface brightnesses and luminosity ratios are comparable. Though the central energetic AGN with a hard UV radiation field and strong X-ray emission may contribute to the ISM ionization and heating in the nuclear region (Stacey et al., 2010; Smith et al., 2017; Shangguan et al., 2018; Herrea-Camus et al., 2018), it is unlikely to be the dominant cause of the [C II]-FIR deficit in the quasar host galaxies (see also Decarli et al., 2018). As was pointed out by the spatially resolved studies of SMGs (Gullberg et al., 2018; Rybak et al., 2019) and local AGNs (Smith et al., 2017), the change of [C II]/FIR ratio with FIR surface brightness at different radii in the quasar host galaxies suggest that the [C II]/FIR deficit depends on the local conditions (radiation field, gas density,etc.), rather than the global properties of the systems.

We did not detect dust continuum emission at the [C II] peak to the east of the quasar SDSS J1044-0125. The continuum image sets a 3$\sigma$ upper limit of 0.12 mJy for the flux density at the [C II] frequency. Assuming dust temperatures of 25 to 45 K, this constrains the FIR luminosity of this component to be $\rm<2\times 10^{11}\,L_{\odot}$, and a [C II]-to-FIR luminosity ratio of $\rm>2.4\times 10^{-3}$ which is much higher than the central region of the quasar host. The high [C II]-to-FIR ratio of this component reveals ISM conditions different from that in the intense starburst nuclear region of the quasar host galaxy. Deep imaging of molecular CO and other fine structure line emission will help to measure the physical properties and excitation of the medium at this offset [C II] peak and address the nature of this component (i.e., core of mergering galaxy or node of outflow, etc. see more discussions below).

4.2 Gas kinematics in the nuclear region

The [C II]-emitting gas in SDSS J0129-0035 shows clear velocity gradients, suggesting that most of the gas is likely from a rotating disk. The ordered motion constrains the dynamical mass of the quasar host galaxy within the [C II]-emitting region (see details Section 4.3). The $\rm 0.2^{\prime\prime}$ resolution [C II] image of SDSS J0129-0035 also reveals additional gas clumps, in particular, the components in the $\leq 0.2^{\prime\prime}$ region of the PV-diagram with velocities of $\rm\sim\pm 200\,km\,s^{-1}$ (Marked as No. 3 and 4 in the lower panel of Figure 3) and the two extended features as were described in Section 3. Though the individual tentative peaks with limited ($\rm 3\sim 4\sigma$, Figure 3) SNR should be checked with deeper observations, we argue that the rotating velocity field with turbulent clumps reveals complex kinematics of the gas component in the nuclear region. Note that similar turbulent rotating gas disks were also reported with the CO(6-5) and [C II] images of SDSS J231038.88+185519.7 at z=6.0 (Feruglio et al., 2018) and the [C II] image of VIKING J030516.92$-$315056.0 at z=6.6 (Venemans et al., 2019) and was discussed in recent numerical simulations (e.g., Lupi et al., 2019); the gas around these powerful quasars show significant dynamical evolution with sub-structures, including spirals, outflow and inflow clumps, and/or truncated disk.

The [C II] map of SDSS J1044-0125 does not show a clear sign of rotation. The PV-diagrams (Figure 5) suggest a very turbulent gas velocity field in the central area; the gas within the $\rm 0.2^{\prime\prime}$ region shows velocities ranging from -200 to $\rm\gtrsim 200\,km\,s^{-1}$. In addition, the [C II] line spectrum suggest multiple peaks (Table 2). Comparison between the CO (6-5) and [C II]158$\mu$m line spectra of this object (Figure 6) shows that the CO (6-5) line does not cover the same velocity range as that of the [C II] line. The high-J CO (e.g. CO (6-5), (7-6)) and [C II] lines in other millimeter-bright quasars at z$\sim$6 usually show similar velocities and line widths (Riechers et al., 2009; Wang et al., 2013), suggesting that these lines trace similar kinematics of the ISM from the nuclear star forming region. The case of SDSS J1044-0125 suggest a multiple-phase ISM with different kinematic properties. The [C II] emission in galaxies is reported to trace the disk gas and star forming regions, as well as the warm diffuese neutral/ionized medium (Stacey et al., 2010; Pineda et al., 2014; Gullberg et al., 2015; Croxall et al., 2017; Decarli et al., 2018; Lagache et al., 2018). [C II] emission in the quasar host galaxies could also powered by the central powerful AGN (e.g., from the narrow line region, the X-ray dominated regions, etc (Stacey et al., 2010; Uzgil et al., 2016)). Asymmetric [C II] line emission that is broader than the molecular CO emission was detected the radio galaxy 3C 326N, suggest [C II] emission from a warm diffuse and turbulent molecular gas component powered by AGN-jet activity (Guillard et al, 2015). On the other hand, the CO(6-5) line emission found in the high-z quasar-starburst systems usually traces the dense molecular gas where stars are actively forming (Greve et al., 2014; Liu et al., 2015; Riechers et al., 2009; Venemans et al., 2017a). Thus, it is possible that, in SDSS J1044-0125, most of the [C II] emission covering a similar velocity range with the CO (6-5) line is from the intense star forming regions, while the additional components in the velocity range without CO detection may be from more diffuse/warm medium. This mismatch between the CO and [C II] line spectra is a curious feature and should be checked with deeper and high-resolution imaging of the molecular line emission as well as other fine structure lines to address the complex physical condition, power resource, and kinematics of the ISM in this region.

The peak to the east of SDSS J1044-0125 reveals a [C II]-emitting source with a very broad line width of $\rm FWHM=750\pm 190\,km\,s^{-1}$ at a projected distance of 4.9 kpc from the quasar (Component II in Figure 7). The nature of this offset component is unclear yet. Broad (line width up to $\rm\geq 1000\,km/s$), blue/red-shifted, and low surface brightness [C II] emission was detected as gas outflows in the young quasar hosts at high-redshift, which are evidences of strong AGN feedback (Cicone et al., 2014, 2015; Fiore et al., 2017; Bischetti et al., 2018). The offset peak we found in SDSS J1044-0125 could be a bright node of such outflowing gas component. Companion [C II]/submm continuum sources are detected close to other luminous quasars at high redshift (Decarli et al., 2017; Trakhtenbrot et al., 2017), revealing that these massive objects are in the centers of very active environments and that galaxy mergers may be a major triggering mechanism for SMBH and galaxy growth in these systems (Li et al., 2007; Narayanan et al., 2008, 2015; Volonteri et al., 2015; Watson et al., 2015; Valiante et al., 2017; Marrone et al., 2018). SDSS J1044-0125, as one of the most luminous quasars known at this redshift, may lie in a massive dark matter halo, and has already experienced several galaxy merger events (Li et al., 2007; Valiante et al., 2017; Marrone et al., 2018). Thus, it is also possible that this offset [C II] component is the dense core of a satellite galaxy that will merge with the quasar host.

In summary, the ALMA observations at 0.2*′′* resolution reveal complex gas kinematics in the nuclear region of the quasar starburst systems at z$\sim$6, e.g., rotation and/or turbulent gas clumps. The turbulent gas kninematis, as well as the multiphase ISM suggested by the [C II] and CO emission in SDSS J1044-0125, indicate rich dynamical activities of ISM powered by the central AGN and/or the intense nuclear star formation. It will be interesting to expand the high resolution imaging studies to more examples to investigate the kinematic activities of the ISM in these extreme quasar-starburst systems. The nature of the peak to the east of SDSS J1044-0125 should also be checked with deeper observations of other emission lines to constrain the ISM properties and its dynamical connection to the quasar host.

4.3 Constraints on the host galaxy dynamics of SDSS J0129-0035

We here model the gas velocity map of SDSS J0129-0035 with a simple circular rotating gas disk, following the method used in Shao et al. (2017) (see also Jones et al., 2017), to constrain the gas dynamics in this system. We use the ROTCUR task in the Groningen Image Processing System (GIPSY3; van der Hulst et al., 1992) to fit the observed velocity field. As was described in Shao et al. (2017), the parameters for the model are: (i) The sky coordinates of the rotational center of the disk. (ii) The system velocity $\rm V_{sys}$. (iii) The circular velocity at the fitting radius $\rm V_{c}(R)$. (iv) The position angle of the major axis on the receding half of the disk, anti-clockwise from the north direction, $\rm\phi$. (v) The inclination angle between the normal direction of the ring and line of sight, incl, and (vi) the azimuthal angle related to $\rm incl$, $\rm\phi$, and the disk center. We fix the disk center as the Gaussian peak center of the [C II] emission, and choose $\rm V_{sys}=0\,km\,s^{-1}$ to refer to the [C II] redshift of 5.7787 (Wang et al., 2013). We then use the data within a radius of 0.35*′′* to determine $\rm\phi$ and incl. This is approximately the $\rm 2\sigma$ contour region of the [C II] line emission (Figure 1). This fit gives $\rm\phi$ and incl as $\rm\phi=56\pm 1^{\circ}$ and $\rm incl=16\pm 20^{\circ}$, respectively, where the error bars denote the fitting errors. This suggests a gas disk that is very close to face-on with incl$\leq 36^{\circ}$.

In Wang et al. (2013), we adopted the same inclined circular disk assumption and estimated the inclination angle from the deconvolved axis ratio of the Cycle 0 [C II] map, i.e., $\rm incl=arccos{(a_{minor}/a_{major})}=56^{\circ}$. This is much larger than the value we obtained in this paper. It is likely that the two extended features (No. 1 and 2 in Figure 3) are unresolved in the low-resolution ALMA Cycle 0 image, which results in a larger major axis value, i.e., $\rm a_{major}=0.41^{\prime\prime}\pm 0.06^{\prime\prime}$ from the Cycle 0 data compared to $\rm a_{major}=0.32^{\prime\prime}\pm 0.03^{\prime\prime}$ in this work (Table 1). We cannot determine if these two features are part of the rotating disk or not, thus, we did not include those pixels in the fitting in this work.

Adopting the $\rm\phi$ and incl values above, we constrain the circular velocities at radii of 0.1*′′, 0.2′′, and 0.3′′* (Figure 10), corresponding to physical scales of 0.6 kpc, 1.2 kpc, and 1.8 kpc. If the inclination angle is fixed to 16*∘*, the circular velocity at the largest radius (i.e., 1.8 kpc) is $\rm V_{c}$=250 $\rm km\,s^{-1}$. Considering the fitting uncertainties of incl, we estimate a lower limit to the circular velocity of $\rm V_{c}$=110 $\rm km\,s^{-1}$ with $\rm incl=36^{\circ}$. An upper limit of the circular velocity cannot be constrained as incl could be very close to zero. Using $\rm V_{c}=250\,km\,s^{-1}$, we derive the dynamical mass within a radius of 1.8 kpc to be $\rm M_{dyn}=2.6\times 10^{10}\,M_{\odot}$. If the lower limit of the circular velocity of 110 $\rm km\,s^{-1}$ is adopted, the dynamical mass will descrease to $\rm 5.2\times 10^{9}\,M_{\odot}$.

Based on the quasar 1450$\AA$ magnitude from the SDSS photometry (Jiang et al., 2009) and a 1450$\AA$ bolometric correction of $\rm L_{bol}=4.2\times\lambda L_{\lambda 1450}$ (Runnoe et al., 2012), we detrive a bolometric lumiosity of $\rm 5.7\times 10^{12}\,L_{\odot}$ for SDSS J0129-0035. Assuming that the SMBH is accreting at the Eddington limit, we estimate the SMBH mass of this object to be $\rm M_{BH}=1.7^{+1.7}_{-0.9}\times 10^{8}\,M_{\odot}$. We consider an uncertainty of 0.3 dex in the SMBH mass, which represents the scatter of the Eddington ratio of the z$\sim$6 quasar sample (De Rosa et al., 2011). This result in a SMBH to host galaxy dynamical mass ratio of $\rm M_{BH}/M_{dyn}=0.007$ with $\rm V_{c}=250\,km\,s^{-1}$. If the lower limit of $\rm V_{c}=110\,km\,s^{-1}$ is adopted, the mass ratio will become 0.03. Considering the facts that i) the model has large uncertainties for face-on system and cannot accurately determine the inclination angle and circular velocity, ii) we do not reach the flat part of the rotation curve as shown in Figure 10, and iii) the outer part of the disk may not be traced by the [C II] velocity map due to low signal-to-noise ratio or is disrupted (Venemans et al., 2019; Lupi et al., 2019), the $\rm M_{dyn}$ value derived above should be considered as a lower limit of the host galaxy dynamical mass for SDSS J0129-0035. The corresponding $\rm M_{BH}/M_{dyn}$ ratio hence represents an upper limit.

We compare the SMBH mass and the lower limit of the host galaxy dynamical mass of SDSS J0129-0035 to those of other quasars at z$\sim$6 to 7 quasars and the relationship of local galaxies (Kormendy & Ho, 2013). We collect the dynamical mass measurements of the z$\sim$6 to 7 quasars from the literatures. These objects are all detected in [C II] and/or CO line emission (Wang et al., 2013, 2016; Cicone et al., 2015; Willott et al., 2015, 2017; Venemans et al., 2016; Feruglio et al., 2018; Izumi et al., 2018; Decarli et al., 2018; Wang et al., 2019) and have source size measurements from [C II] or CO data and inclination angles estimated from the minor and major axis ratios of the line emission. The two objects, SDSS J0129-0035 (this work) and ULAS J1319+0950 (Shao et al., 2017), were observed with ALMA at a much better resolution of $\rm 0.2^{\prime\prime}\sim 0.3^{\prime\prime}$. The [C II] emission is resolved and the velocity maps show clear velocity gradients. These, combined with the disk dynamical model described above, provide more reliable constraints on the disk inclination angle, circular velocity, and dynamical mass of the quasar host galaxies. We denote these two objects as blue and red stars in the plot. According to Figure 11, the black hole-to-host mass relationship of local galaxies given by Kormendy & Ho (2013) predicts a mass ratio of 0.0043 for a SMBH mass of $\rm 1.7\times 10^{8}\,M_{\odot}$. Thus, it is unlikely that the black hole-to-host mass ratio of SDSS J0129-0035 is significantly above the local relationship. This is unlike the most luminous/massive quasars at this redshift, which usually show ratios a few to $\gtrsim 10$ times higher (Venemans et al., 2016; Decarli et al., 2018; Wang et al., 2019). This agrees with the results of other [C II]-detected z$\gtrsim$6 quasars with similar or lower SMBH masses, in which the ratios are consistent with the trend defined by the local systems (Willott et al., 2010, 2017; Izumi et al., 2018). The dynamical mass constraints based on the [C II] or CO imaging of these z$\sim$6 to 7 quasars suggest that, in the early universe, the most massive SMBHs with masses of $\rm 10^{9}$ to $\rm 10^{10}\,M_{\odot}$ may grows faster than that of their host galaxies (Walter et al., 2004; Venemans et al., 2016; Wang et al., 2016; Decarli et al., 2018; Wang et al., 2019) while the less massive systems ($\rm 10^{7}\sim 10^{8}\,M_{\odot}$) are evolving more closer to the trend of local galaxies (Willott et al., 2017; Izumi et al., 2018).

5 Summary

Based on new ALMA observations we present images of the [C II] and dust continuum emission from the host galaxies of two z=5.78 quasars, SDSS J0129$-$0035 and SDSS J1044$-$0125 on $\rm 0.2^{\prime\prime}$ scales, resolving the ISM emission in the central $\rm\sim 4\,kpc$ of the hosts. The main conclusions are summarized as follows:

• •

The FWHM sizes of the [C II] emission is 1.8-2.4 kpc and that of the dust continuum $\sim$1kpc for the two bojects, constraining the size of the intense star-forming region around the accreting SMBHs. As was found in other starburst quasar host galaxies and SMGs at high redshift, The [C II] emission regions appear to be more extended compared to that of dust continuum (Gullberg et al., 2018; Cooke et al., 2018).

• •

We derive the [C II] and FIR surface brightness and emission ratios at different radii based on the resolved [C II] and continuum images. The resolved [C II]-FIR ratios are decreasing with increasing $\rm\Sigma_{FIR}$ toward the center, following the trend of the IR luminous star forming galaxies with $\rm\Sigma_{FIR}$ in the same range of a few $\rm 10^{10}\,L_{\odot}\,kpc^{-2}$ to $\rm 2\times 10^{12}\,L_{\odot}\,kpc^{-2}$ (Stacey et al., 2010; Gullberg et al., 2015, 2018; Díaz-Santos et al., 2013; Lutz et al., 2016; Herrea-Camus et al., 2018; Cooke et al., 2018; Rybak et al., 2019). The low [C II]-FIR ratios of $0.0001\sim 0.0002$ found at the peak $\rm\Sigma_{FIR}$ in the central region of these two quasar hosts suggest dense ISM with high-intensity radiation field similar to that in the compact dusty starburst systems (Díaz-Santos et al., 2017; Herrea-Camus et al., 2018; Gullberg et al., 2018; Rybak et al., 2019).

• •

The velocity field and PV diagrams of SDSS J0129-0035 reveal a combination of rotation with turbulent gas clumps; We detect clear velocity gradients in the velocity map. We also see features with velocities offsets of -190 and $\rm+50\,km\,s^{-1}$ at the outer part of the [C II] emitting region and gas clumps with velocity offsets of $\rm\pm 200\,km\,s^{-1}$ within the central 0.2*′′* region, suggesting other kinematic activities, such as inflow or outflows.

• •

The [C II]-emitting gas in the central region of SDSS J1044-0125 does not show a clear sign of rotation. The PV-diagram of the [C II] emission suggests very turbulent gas kinematics in the central area. The [C II] line spectrum shows multiple peaks, which is broader and slightly offset from the CO (6-5) line emission. In addition, we tentatively detect an offset peak, 4.9 kpc away to the east of the quasar with a very broad line width ($\rm 750\pm 190\,km\,s^{-1}$) which may relate to an infalling companion source or node of outflows.

• •

Overall, the gas kinematics on kpc scale in the nuclear region of the two starburst quasar host galaxies reveal significant dynamical evolution of the ISM, which could be powered by both AGN and the intense star formation.

• •

Based on the dynamical modeling of the [C II] velocity gradients of SDSS J0129$-$0035, we derive a disk inclination angle very close to face-on (i.e., $\rm incl=16\pm 20^{\circ}$). We can thus only put a lower limit on the host dynamical mass within the [C II]-emitting region. The result argues against an SMBH-bulge mass ratio that is significantly above (e.g., larger by a factor of 10) the value derived from the local relationship. As was discussed in previous studies of the less luminous/massive quasar population at z$\sim$6, the mass ratio of SDSS J0129-0035 is likely to be closer to the local relationship than to that found for the most luminous and massive quasars at the same redshift (Venemans et al., 2016; Willott et al., 2010; Izumi et al., 2018; Wang et al., 2019).

The ALMA images of the [C II] emission in these young quasar host galaxies at kpc resolution reveal rich features of gas kinematics in the nuclear region around the powerful AGNs. Between the two objects studied in this work, we do see more turbulent/clumpy gas in the nuclear region and an curious discrepancy between the [C II] and CO line profiles in SDSS J1044$-$0125, which hosts a more massive SMBH with higher quasar luminosity. However, it is still unclear if this is due to the feedback of the powerful AGN or nuclear star formation. If the east peak is associated with infalling companion, the tidal perturbations could also affect the gas dynamic in the nuclear region. Deeper images of the [C II] and molecular CO emission from this object are still required to confirm these features. The high-resolution imaging studies should also be expanded to a larger sample to check for significant evolution of gas kinematic properties with AGN activity, e.g., more turbulent or disrupting gas disk in more luminous AGNs. Deep imaging of the molecular gas (e.g., with ALMA) and the stellar component (e.g., with the JWST in the future) on similar scales are also urgently required to fully probe host galaxy evolution of these young quasars at the earliest epochs.

This work was supported by National Key Program for Science and Technology Research and Development (grant 2016YFA0400703). We are thankful for the supports from the National Science Foundation of China (NSFC) grants No.11721303, 11373008, 11533001, the Strategic Priority Research Program ′The Emergence of Cosmological Structures′ of the Chinese Academy of Sciences, grant No. XDB09000000, the National Key Basic Research Program of China 2014CB845700. RW acknowledge supports from the Thousand Youth Talents Program of China, the NSFC grants No. 11443002 and 11473004 D.R. acknowledges support from the National Science Foundation under grant number AST-1614213. DN was funded in part by NSF AST-1715206 and HST AR-15043.0001. The National Radio Astronomy Observatory (NRAO) is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2011.0.00206.S and #2012.1.00240.S . ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada) and NSC and ASIAA (Taiwan), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. RW thanks Dr. Matus Rybak for providing the resolved [C II] and FIR continuum measurements of their SMG sample. RW thanks Dr. Shangguan for providing FIR luminosities of the GOALS sample. rwangkiaa@pku.edu.cn. Facilities: ALMA.