SBS~0335-052E is a young star-forming dwarf galaxy with a total stellar mass of $M_{*} \lesssim 10^{8}~M_{\odot}$ and an extremely low metallicity ($Z \sim 1/40~Z_{\odot}$), which has long been considered to be devoid of an active galactic nucleus (AGN). Here we report the detection of temporal flux variability of SBS~0335-052E in near-infrared (NIR) 3-4\ ${\rm \mu}$m bands on timescales of several years, showing dimming and brightening of up to 50\% over 14~years, based on archival data from the Wide-field Infrared Survey Explorer. Our spectral energy distribution (SED) fitting of archival ultraviolet (UV)-NIR photometry, including AGN SED models, indicates that the variable NIR emission arises from an edge-on AGN dust torus. The UV-optical emission from the accretion disk is obscured and does not reach us, leading to the dominance of the host galaxy's young stellar population in the UV-optical wavelengths. This analysis favors the presence of a Compton-thick, heavily obscured AGN in SBS~0335-052E, consistent with its observed X-ray weakness. From the SED fitting, we estimate an AGN bolometric luminosity of $L_{\rm bol} = 1.2\times10^{43}\ {\rm erg\ s^{-1}}$, which implies a black hole mass of $M_{\rm BH} \simeq 10^{5}\ M_\odot$ if the AGN is accreting at the Eddington limit. If confirmed, SBS~0335-052E would be the least massive galaxy known to host an AGN, likely harboring an intermediate-mass black hole.

Understanding the rapid formation of supermassive black holes (SMBHs) in the early universe requires insight into stellar mass growth in host galaxies. Here, we present NIRSpec rest-frame optical spectra and NIRCam imaging from JWST of two galaxies at z>6, both hosting moderate-luminosity quasars. These galaxies exhibit Balmer absorption lines, similar to low-redshift post-starburst galaxies. Our analyses of the medium-resolution spectra and multiband photometry show bulk of the stellar mass (log (M_* / M_sun) > 10.6) formed in starburst episodes at redshift 9 and 7. One of the galaxies shows a clear Balmer break and lacks spatially resolved H alpha emission. It falls well below the star formation main sequence at z = 6, indicating quiescence. The other is transitioning to quiescence; together, these massive galaxies are among the most distant post-starburst systems known. The blueshifted wings of the quasar [O III] emission lines suggest quasar-driven outflow possibly influencing star formation. Direct stellar velocity dispersion measurements reveal one galaxy follows the local black hole mass-sigma_* relation while the other is overmassive. The existence of massive post-starburst galaxies hosting billion-solar-mass BHs in short-lived quasar phases suggests SMBHs and host galaxies played a major role in each other's rapid early formation.

We show the Eddington ratio distributions of supermassive black holes at a wide redshift range (0 < z < 8) obtained with a semi-analytic model of galaxy formation. The distribution is broadly consistent with observational estimates at low redshift. We find that the growth rate of black holes at higher redshift is more likely to exceed the Eddington limit because the typical gas fraction of the host galaxies is higher at higher redshift. We also find that the super- Eddington growth is more common for less massive supermassive black holes, supporting an idea that supermassive black holes have been formed via super-Eddington accretion. These results indicate the "slowing down" of cosmic growth of supermassive black holes: the growth of supermassive black holes with a higher Eddington ratio peaks at higher redshift. We also show the effect of the sample selection on the shape of the Eddington ratio distribution functions and find that shallower observations will miss active galactic nuclei with not only the smaller but also higher Eddington ratios.

We apply a new X-ray clumpy torus model called XCLUMPY (Tanimoto et al. 2019), where the clump distribution is assumed to be the same as in the infrared clumpy torus model (CLUMPY) by Nenkova et al. 2008, to the broadband X-ray spectra of type-1 active galactic nuclei (AGNs) for the first time. We analyze the archival data of IC 4329A and NGC 7469 observed with NuSTAR/Suzaku and NuSTAR/XMM-Newton, respectively, whose infrared spectra were studied with CLUMPY by Ichikawa et al. 2015 and optical extinctions ($A_{\rm V}$) of the tori were estimated. We consider two models, invoking (Model 1) a relativistic reflection component from the accretion disk and (Model 2) a partial absorber. Assuming that the narrow Fe K$\alpha$ emission line at 6.4 keV originates from the torus, we separate the contribution of the torus reflection components in the total spectra. Our models yield equatorial hydrogen column densities of the tori to be $N^{\rm Equ}_{\rm H} =$ (0.53--1.43) $\times 10^{23}~\rm cm^{-2}$ and $N^{\rm Equ}_{\rm H} =$ (0.84--1.43) $\times 10^{24}~\rm cm^{-2}$, for IC 4329A and NGC 7469, respectively. We find that the $N_{\rm H}/A_{\rm V}$ ratios in the tori are by factors of 25--68 (IC 4329A) and 2.4--3.9 (NGC 7469) smaller than that in the Galactic interstellar medium (ISM). These results suggest that a non-negligible fraction of AGNs are 'dust-rich' compared with the Galactic ISM, as opposite to the general trend previously reported in many obscured AGNs.

We present deep Subaru/FOCAS spectra for two extreme emission line galaxies (EELGs) at $z\sim 1$ with strong {\sc[Oiii]}$\lambda$5007 emission lines, exhibiting equivalent widths (EWs) of $2905^{+946}_{-578}$ \AA\ and $2000^{+188}_{-159}$ \AA, comparable to those of EELGs at high redshifts that are now routinely identified with JWST spectroscopy. Adding a similarly large {\sc [Oiii]} EW ($2508^{+1487}_{-689}$ \AA) EELG found at $z\sim 2$ in the JWST CEERS survey to our sample, we explore for the physical origins of the large {\sc [Oiii]} EWs of these three galaxies with the Subaru spectra and various public data including JWST/NIRSpec, NIRCam, and MIRI data. While there are no clear signatures of AGN identified by the optical line diagnostics, we find that two out of two galaxies covered by the MIRI data show strong near-infrared excess in the spectral energy distributions (SEDs) indicating obscured AGN. Because none of the three galaxies show clear broad H$\beta$ lines, the upper limits on the flux ratios of broad-H$\beta$ to {\sc [Oiii]} lines are small, $\lesssim 0.15$ that are comparable with Seyfert $1.8-2.0$ galaxies. We conduct \texttt{Cloudy} modeling with the stellar and AGN incident spectra, allowing a wide range of parameters including metallicities and ionization parameters. We find that the large {\sc [Oiii]} EWs are not self-consistently reproduced by the spectra of stars or unobscured AGN, but obscured AGN that efficiently produces O$^{++}$ ionizing photons with weak nuclear and stellar continua that are consistent with the SED shapes.

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 &lt; M_{1450} &lt; -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 present the stellar population and ionized-gas outflow properties of ultra-luminous IR galaxies (ULIRGs) at $z=$ 0.1-1.0, which are selected from AKARI FIR all-sky survey. We construct a catalog of 1077 ULIRGs to examine feedback effect after major mergers. 202 out of the 1077 ULIRGs are spectroscopically identified by SDSS and Subaru/FOCAS observations. Thanks to deeper depth and higher resolution of AKARI compared to the previous IRAS survey, and reliable identification from WISE MIR pointing, the sample is unique in identifying optically-faint (i$\sim$20) IR-bright galaxies, which could be missed in previous surveys. A self-consistent spectrum-SED decomposition method, which constrains stellar population properties in SED modeling based on spectral fitting results, has been employed for 149 ULIRGs whose optical continua are dominated by host galaxies. They are massive galaxies ($M_{\rm star}\sim10^{11}$-$10^{12}$ M$_{\odot}$), associated with intense star formation activities (SFR $\sim$ 200-2000 M$_{\odot}$ yr$^{-1}$). The sample covers a range of AGN bolometric luminosity of $10^{10}$-$10^{13}$ L$_{\odot}$, and the outflow velocity measured from [OIII] 5007A line shows a correlation with AGN luminosity. Eight galaxies show extremely fast outflows with velocity up to 1500-2000 km s$^{-1}$. However, the co-existence of vigorous starbursts and strong outflows suggests the star formation has not been quenched during the ULIRG phase. By deriving stellar mass and mass fraction of young stellar population, we find no significant discrepancies between stellar properties of ULIRGs with weak and powerful AGNs. The results are not consistent with the merger-induced evolutionary scenario, which predicts that SF-dominated ULIRGs show smaller stellar mass and younger stellar populations compared to AGN-dominated ULIRGs.

Which galaxies in the general population turn into active galactic nuclei (AGNs) is a keystone of galaxy formation and evolution. Thanks to SRG/eROSITA's contiguous 140 square degree pilot survey field, we constructed a large, complete, and unbiased soft X-ray flux-limited (F_X&gt;6.5\times 10^{-15} erg s$^{-1}$ cm$^{-2}$) AGN sample at low redshift, $0.05

We show the significance of the super-Eddington accretion for the cosmic growth of supermassive black holes (SMBHs) with a semi-analytical model for galaxy and black hole evolution. The model explains various observed properties of galaxies and active galactic nuclei at a wide redshift range. By tracing the growth history of individual SMBHs, we find that the fraction of the SMBH mass acquired during the super-Eddington accretion phases to the total SMBH mass becomes larger for less massive black holes and at higher redshift. Even at z = 0, SMBHs with > 1e+9 Msun have acquired more than 50% of their mass by super-Eddington accretions, which is apparently inconsistent with classical Soltan's argument. However, the mass-weighted radiation efficiency of SMBHs with > 1e+8 Msun obtained with our model, is about 0.08 at z = 0, which is consistent with Soltan's argument within the observational uncertainties. We, therefore, conclude that Soltan's argument cannot reject the possibility that SMBHs are grown mainly by super-Eddington accretions.

We present a study of the molecular gas in five closely-spaced (R_{\perp}&lt;20 kpc) dual quasars ($L_{\rm bol}\gtrsim10^{44}~\mathrm{erg~s}^{-1}$) at redshifts $0.4

We present the latest results of a semi-analytic model of galaxy formation, "New Numerical Galaxy Catalogue", which is combined with large cosmological N-body simulations. This model can reproduce statistical properties of galaxies at z < 6.0. We focus on the properties of active galactic nuclei (AGNs) and supermassive black holes, especially on the accretion timescale onto black holes. We find that the number density of AGNs at z < 1.5 and at hard X-ray luminosity 10^{ 44 }< erg/s is underestimated compared with recent observational estimates when we assume the exponentially decreasing accretion rate and the accretion timescale which is proportional to the dynamical time of the host halo or the bulge, as is often assumed in semi-analytic models. We show that to solve this discrepancy, the accretion timescale of such less luminous AGNs instead should be a function of the black hole mass and the accreted gas mass. This timescale can be obtained from a phenomenological modelling of the gas angular momentum loss in the circumnuclear torus and/or the accretion disc. Such models predict a longer accretion timescale for less luminous AGNs at z < 1.0 than bright QSOs whose accretion timescale would be 10^{ 7-8 } yr. With this newly introduced accretion timescale, our model can explain the observed luminosity functions of AGNs at z < 6.0.

We compiled the X-ray and soft gamma-ray observations of the Galactic black hole binary XTE J1859+226 in the 1999--2000 outburst from RXTE, ASCA, BeppoSAX and CGRO. Throughout systematic spectral analysis using a two-component model consisting of a multi-temperature accretion disk plus a fraction of its flux convolved with an empirical Comptonized powerlaw component, we found that the innermost radius ($r_{\rm in}$) and temperature (Tin) of the disk are very variable with time in the rising phase of soft X-ray flux where Type-A/B/C low-frequency quasi-periodic oscillations (QPOs) were found. After this phase, $r_{\rm in}$ remains constant at around 60 km assuming a distance of 8 kpc and an inclination angle of 67$^{\circ}$, and Tin smoothly decays with time. The constant $r_{\rm in}$ suggests a presence of the innermost stable circular orbit (ISCO), with $r_{\rm in}$ repeatedly moving closer and farther away from the ISCO in the rising phase. Both disk parameters are remarkably correlated with independently analyzed timing properties such as QPO frequency and rms variability. Type-A/B QPOs are seen only when $r_{\rm in}$ is close to the ISCO, while Type-C are seen when $r_{\rm in}$ is truncated and the frequency changes with a relation of $r^{-1.0}_{\rm in}$, supporting that Type-C QPOs occur at the inner edge of the truncated disk. Accurate determinations of the frequency--$r_{\rm in}$ relation for various objects should be a powerful tool to discriminate plausible Type-C QPO models. Furthermore, we suggest that jet ejection events may occur when $r_{\rm in}$ rapidly approaches to the ISCO, along with rapid changes of the disk flux, the rms variability and the hardness ratio. A rapid shrinkage of $r_{\rm in}$ down to the ISCO can be a useful index as a precursor of radio flares for triggering Target-of-Opportunity observations and would provide constraints on jet launching mechanisms.

The relationship between quasars and their galaxy environment is important for understanding the evolution of galaxies and supermassive black holes, but it is not fully understood. We perform a wide and deep exploration of the environment of quasars at 0.4 &lt; z &lt; 1.0 using the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) survey. We investigate the environment of the 1,912 spectroscopically selected quasars from the Sloan Digital Sky Survey (SDSS), using photometrically selected galaxies from the HSC-SSP data, over an area of 505 deg${^{2}}$. The quasar environment is compared to the environment of matched galaxies with similar stellar mass and redshift. We employ the $k$-nearest neighbor method to define the local galaxy number density for both the quasars and the matched galaxies at a scale of a few hundred kpc. As a result, we find that the number density of galaxies around SDSS quasars is lower than that of the matched galaxies by $\sim$11--$20\%$. We also investigate possible correlations between the local galaxy number densities and the quasar properties such as black hole mass and Eddington ratio. As a result, no correlation is found between the local galaxy number densities and these properties of quasars. These results suggest that the quasar activity is not triggered by the high number density of surrounding galaxies at the scale of a few hundred kpc.

High-$z$ radio galaxies (HzRGs) are considered important objects for understanding the formation and evolution of massive galaxies in the early universe. However, till date, detailed studies of the stellar population of HzRGs such as the star-formation history have been scarce. Therefore, this study conducted a new survey to establish a less-biased sample of HzRGs and consequently investigate their properties. We utilized a sample of $g$-dropout Lyman break galaxies (LBGs) obtained from an optical wide and deep imaging survey made by Subaru Hyper Suprime-Cam (HSC). Based on the cross-matching of this LBG sample with the VLA FIRST radio survey data, we constructed a photometric sample of high-redshift radio galaxies (HzRGs) at $z \sim 4$ for $\sim$560 deg$^2$ survey field. Consequently, we identified 146 HzRG candidates. To analyze the characteristics of these candidates, we focus on objects exhibiting the near-infrared photometry of VIKING or UKIDSS and the mid-infrared photometry of unWISE (28 objects). The results indicate that 7 objects exhibit SEDs consistent with galaxies at $z \sim 4$. The HzRG candidates have very large stellar masses with $\sim 4.2 \times 10^{11} M_{\odot}$ on average. This stellar mass is similar to that of previously discovered USS HzRGs at $z \sim 4$, though our sample is affected by a sample selection bias that selects only HzRGs with $M_{\star} > 10^{11} M_{\odot}$. Further, the SEDs of those HzRG candidates suggest a past fast quenching with a rough timescale of $\sim$0.1 Gyr, as evidenced from the rest-frame UVJ diagram.

HSC 120505.09-000027.9 (J1205$-$0000) is one of the highest redshift ($z=6.72$) dust-reddened quasars (red quasars) known to date. We present an improved analysis of Atacama Large Millimeter/submillimeter Array data of the [CII] $158\ \rm{\mu m}$ line and the underlying rest-frame far-infrared (FIR) continuum emission, previously reported in Izumi et al. (2021a), toward J1205$-$0000. Red quasars are thought to be a transitional phase from an obscured starburst to a luminous blue quasar, in some cases associated with massive outflows driven by the active galactic nucleus (AGN). J1205$-$0000 has a high FIR luminosity, $L_{\mathrm{FIR}}=2.5\times 10^{12}\ L_{\odot}$ and a total IR luminosity of $L_{\mathrm{TIR}}=3.5\times 10^{12}\ L_{\odot}$, corresponding to a star formation rate (SFR) of $\sim 528\ M_{\odot}\ \mathrm{yr}^{-1}$. With the [CII]-based dynamical mass of$\sim 1 \times 10^{11}~M_\odot$, we conclude that J1205$-$0000 is hosted by a starburst galaxy. In contradiction to Izumi et al. (2021a), our improved analysis shows no hint of a broad component in the [CII] line spectrum. Thus there is no evidence for a host galaxy-scale fast [CII] outflow, despite the fact that J1205$-$0000 has fast nuclear ionized outflows seen in the rest-frame UV. We explore several scenarios for this discrepancy (e.g., early phase of AGN feedback, reliability of the [CII] line as a tracer of outflows), and we claim that it is still too early to conclude that there is no significant negative AGN feedback on star formation in this red quasar.

Dust-obscured galaxies (DOGs) are considered to be in a co-evolution phase, with the associated active galactic nuclei (AGN) obscured by dust and gas. Although the DOGs are thought to harbor rapidly growing SMBHs, their X-ray statistical properties, crucial for understanding the properties of obscuring gas as well as the accretion disk state and the hot electron corona around the SMBHs, remain unexplored due to the combination of the low number density of DOGs and the lack of X-ray surveys achieving both of the wide-area and uniformly high-sensitivity observations. We construct a sample of X-ray-detected DOGs in the eROSITA Final Equatorial Depth Survey (eFEDS) field and examine their X-ray statistical properties. By using Subaru/HSC SSP, VIKING, and WISE all-sky surveys, our results reveal the discovery of 5738 IR-bright DOGs in the footprint covered by both of the eFEDS and VIKING surveys (60 deg^2), with 65 objects identified as X-ray-detected DOGs (eFEDS-DOGs). Among them, 41 eFEDS-DOGs show a power-law slope in the near to mid-IR bands (power-law DOGs), indicating dust-obscured AGN. The hydrogen column density (N_H) suggests that eFEDS-DOGs cover even unobscured AGN, spanning 10^20 < N_H <= 10^23. On the other hand, the majority of IR-bright DOGs are not detected by eROSITA, suggesting that most IR-bright DOGs are heavily obscured by dust and gas with N_H > 10^23. Therefore, eFEDS-DOGs, discovered thanks to the wide-area survey by eROSITA, are newly found populations showing less obscured phases among the lifetime of DOGs. Additionally, some eFEDS-DOGs exhibit deviations, down to nearly 1.0 dex below the monochromatic luminosity at 6 micron versus absorption-corrected intrinsic X-ray luminosity between 0.5-2 keV relation, suggesting that it may signal high Eddington ratios reaching the Eddington limit.

We present ionizing spectra estimated at 13.6--100 eV for ten dwarf galaxies with strong high ionization lines of He {\sc {ii}}$\lambda$4686 and [Ne {\sc{v}}]$\lambda$3426 ([Ne {\sc{iv}}]$\lambda$2424) at $z=0$ ($z=8$) that are identified in our Keck/LRIS spectroscopy and the literature (the JWST ERO program). With the flux ratios of these high ionization lines and $>10$ low-ionization lines of hydrogen, helium, oxygen, neon, and sulfur, we determine ionizing spectra consisting of stellar and non-thermal power-law radiation by photoionization modeling with free parameters of nebular properties including metallicity and ionization parameter, cancelling out abundance ratio differences. We find that all of the observed flux ratios are well reproduced by the photoinization models with the power law index $\alpha_{\rm EUV}$ of $\alpha_{\rm EUV}\sim (-1)-0$ and the luminosity $L_{\rm EUV}$ of $L_{\rm EUV}\sim 10^{40}-10^{42}$ erg s$^{-1}$ at $\sim 55-100$ eV for six galaxies, while four galaxies include large systematics in $\alpha_{\rm EUV}$ caused by stellar radiation contamination. We then compare $\alpha_{\rm EUV}$ and $L_{\rm EUV}$ of these six galaxies with those predicted by the black hole (BH) accretion disk models, and find that $\alpha_{\rm EUV}$ and $L_{\rm EUV}$ are similar to those of the intermediate mass black holes (IMBHs) in BH accretion disk models {albeit with possibilities of the other scenarios.} Confirming these results with a known IMBH having a mass $M_{\rm BH}$ of $M_{\rm BH}=10^{5.75} \ M_\odot$, we find that four local galaxies and one $z=7.665$ galaxy have ionizing spectra consistent with those of IMBHs with $M_{\rm BH} \sim 10^3-10^5 \ M_\odot$.

We examine the clustering of quasars over a wide luminosity range, by utilizing 901 quasars at $\overline{z}_{\rm phot}\sim3.8$ with $-24.73

We present ALMA [C II] 158 $\mu$m line and far-infrared (FIR) continuum emission observations toward HSC J120505.09$-$000027.9 (J1205$-$0000) at $z = 6.72$with the beam size of$\sim 0''.8 \times 0''.5$(or 4.1 kpc$\times$ 2.6 kpc), the most distant red quasar known to date. Red quasars are modestly reddened by dust, and are thought to be in rapid transition from an obscured starburst to an unobscured normal quasar, driven by powerful active galactic nucleus (AGN) feedback which blows out a cocoon of interstellar medium (ISM). The FIR continuum of J1205$-$0000 is bright, with an estimated luminosity of $L_{\rm FIR} \sim 3 \times 10^{12}~L_\odot$. The [C II] line emission is extended on scales of $r \sim 5$ kpc, greater than the FIR continuum. The line profiles at the extended regions are complex and broad (FWHM $\sim 630-780$ km s$^{-1}$). Although it is not practical to identify the nature of this extended structure, possible explanations include (i) companion/merging galaxies and (ii) massive AGN-driven outflows. For the case of (i), the companions are modestly star-forming ($\sim 10~M_\odot$ yr$^{-1}$), but are not detected by our Subaru optical observations ($y_{\rm AB,5\sigma} = 24.4$ mag). For the case of (ii), our lower-limit to the cold neutral outflow rate is $\sim 100~M_\odot$ yr$^{-1}$. The outflow kinetic energy and momentum are both much smaller than what predicted in energy-conserving wind models, suggesting that the AGN feedback in this quasar is not capable of completely suppressing its star formation.

We present ALMA [CII] 158 $\mu$m line and underlying far-infrared (FIR) continuum emission observations ($0''.70 \times 0''.56$ resolution) toward HSC J124353.93$+$010038.5 (J1243$+$0100) at $z = 7.07$, the only low-luminosity ($M_{\rm 1450} > -25$ mag) quasar currently known at $z > 7$. The FIR continuum is bright (1.52 mJy) and resolved with a total luminosity of $L_{\rm FIR} = 3.5 \times 10^{12}~L_\odot$. The spatially extended component is responsible for $\sim 40\%$ of the emission. The area-integrated [CII] spectrum shows a broad wing (${\rm FWHM} = 997$ km s$^{-1}$, $L_{\rm [CII]} = 1.2 \times 10^9~L_\odot$) as well as a bright core (${\rm FWHM} = 235$km s$^{-1}$, $L_{\rm [CII]} = 1.9 \times 10^9~L_\odot$). This wing is the first detection of a galactic-scale quasar-driven outflow (atomic outflow rate $> 447~M_\odot$ yr$^{-1}$) at $z > 7$. The estimated large mass loading factor of the total outflow (e.g., $\gtrsim 9$ relative to the [CII]-based SFR) suggests that this outflow will soon quench the star-formation of the host. The core gas dynamics are governed by rotation, with a rotation curve suggestive of a compact bulge ($\sim 3.3 \times 10^{10}~M_\odot$), although it is not yet spatially resolved. Finally, we found that J1243$+$0100 has a black hole mass-to-dynamical mass ratio (and -to-bulge mass ratio) of $\sim 0.4\%$ ($\sim 1\%$), consistent with the local value within uncertainties. Our results therefore suggest that the black hole-host co-evolution relation is already in place at $z \sim 7$ for this object.