Showing 1–39 of 39 results for author: Kokubo, M

Challenging the AGN scenario for JWST/NIRSpec broad H$α$ emitters/Little Red Dots in light of non-detection of NIRCam photometric variability and X-ray

Authors: Mitsuru Kokubo, Yuichi Harikane

Abstract: JWST has uncovered a substantial population of high-$z$ ($z \gtrsim 4$) galaxies exhibiting broad H$α$ emission line with a Full Width at Half Maximum exceeding $1,000$~km~s$^{-1}$. This population includes a subset known as 'Little Red Dots', characterized by their compact morphology and extremely red rest-frame optical colors. If all of these broad H$α$ emitters were attributed to type $1-1.9$ A… ▽ More JWST has uncovered a substantial population of high-$z$ ($z \gtrsim 4$) galaxies exhibiting broad H$α$ emission line with a Full Width at Half Maximum exceeding $1,000$~km~s$^{-1}$. This population includes a subset known as 'Little Red Dots', characterized by their compact morphology and extremely red rest-frame optical colors. If all of these broad H$α$ emitters were attributed to type $1-1.9$ Active Galactic Nuclei (AGNs), it would imply a significantly higher number density of low-luminosity AGNs than extrapolated from that of more luminous AGNs. Here, we have examined the rest-frame ultraviolet (UV)-optical flux variability of five JWST broad H$α$ emitters using multi-epoch, multi-band JWST/NIRCam imaging data. The rest-frame temporal sampling interval of the NIRCam data ($\sim 400-500$~days$/(1+z)$) is comparable to typical variability timescales of AGNs with black hole (BH) masses of $M_{\text{BH}} \sim 10^{7}~M_{\odot}$; thus, the flux variations should be detectable if AGNs were present. However, no measurable flux variation over the rest-frame wavelength range of $λ_{\text{rest}} \sim 1,500-9,000$Å has been detected, placing stringent upper limits on the variability amplitudes. This result, combined with the X-ray faintness confirmed by the ultra-deep {\it Chandra} data, indicates that, under the AGN scenario, we need to postulate peculiar Compton-thick broad-line AGNs with either ($a$) an intrinsically non-variable AGN disk continuum, ($b$) a host galaxy-dominated continuum, or ($c$) scattering-dominated AGN emission. Alternatively, ($d$) they could be non-AGNs where the broad-line emission originates from unusually fast and dense/low-metallicity star-formation-driven outflows or inelastic Raman scattering of stellar UV continua by neutral hydrogen atoms. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Comments: 33 pages, 12 figures, 6 tables, submitted to ApJ

Updated Picture of the Active Galactic Nuclei with Dusty/Dust-free Gas Structures and Effects of the Radiation Pressure

Authors: Shoichiro Mizukoshi, Takeo Minezaki, Hiroaki Sameshima, Mitsuru Kokubo, Hirofumi Noda, Taiki Kawamuro, Satoshi Yamada, Takashi Horiuchi

Abstract: This study investigates the properties of two gas structures of X-ray selected active galactic nuclei (AGNs), that is, dusty and dust-free gas components, by separating them with the line-of-sight dust extinction ($A_V$) and the neutral gas column density ($N_{\mathrm{H}}$). The typical column density of the dusty and dust-free gas differs depending on the Seyfert type, indicating that both struct… ▽ More This study investigates the properties of two gas structures of X-ray selected active galactic nuclei (AGNs), that is, dusty and dust-free gas components, by separating them with the line-of-sight dust extinction ($A_V$) and the neutral gas column density ($N_{\mathrm{H}}$). The typical column density of the dusty and dust-free gas differs depending on the Seyfert type, indicating that both structures have anisotropic column density distributions. The number of targets with the dusty gas column density ($N_{\mathrm{H,d}}$) of $\log N_{\mathrm{H,d}}\ [\mathrm{cm^{-2}}]>23$ is much smaller than that with the same column density of the dust-free gas. This result indicates that the optically-thick part of the dusty gas structure is very thin. There are very few targets with a larger Eddington ratio ($f_{\mathrm{Edd}}$) than the effective Eddington limit of the dusty gas and the covering factor of the dusty gas with $22\leq \log N_{\mathrm{H,d}}\ [\mathrm{cm^{-2}}]<24$ exhibits a clear drop at the effective Eddington limit. These results support the scenario wherein the covering factor of the dusty torus decreases in a high Eddington ratio owing to the radiation-driven dusty gas outflow. The covering factor of the dust-free gas with the column density ($N_{\mathrm{H,df}}$) of $22\leq \log N_{\mathrm{H,df}}\ [\mathrm{cm^{-2}}]<24$ similarly exhibits the decrease in high Eddington ratio, although it may be owing to the dust-free gas outflow driven by certain other mechanisms than the radiation pressure. Finally, we propose an updated picture of the AGN gas structure based on our results and the literature. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 16 pages, 9 figures, accepted for publication in MNRAS

Optical Variability of Blazars in the Tomo-e Gozen Northern Sky Transient Survey

Authors: TianFang Zhang, Mamoru Doi, Mitsuru Kokubo, Shigeyuki Sako, Ryou Ohsawa, Nozomu Tominaga, Masaomi Tanaka, Yasushi Fukazawa, Hidenori Takahashi, Noriaki Arima, Naoto Kobayashi, Ko Arimatsu, Shin-ichiro Okumura, Sohei Kondo, Toshihiro Kasuga, Yuki Mori, Yuu Niino

Abstract: We studied the optical variability of 241 BL Lacs and 83 flat-spectrum radio quasars (FSRQ) from the 4LAC catalog using data from the Tomo-e Gozen Northern Sky Transient Survey, with $\sim$ 50 epochs per blazar on average. We excluded blazars whose optical variability may be underestimated due to the influence of their host galaxy, based on their optical luminosity ($L_O$). FSRQs with $γ$-ray phot… ▽ More We studied the optical variability of 241 BL Lacs and 83 flat-spectrum radio quasars (FSRQ) from the 4LAC catalog using data from the Tomo-e Gozen Northern Sky Transient Survey, with $\sim$ 50 epochs per blazar on average. We excluded blazars whose optical variability may be underestimated due to the influence of their host galaxy, based on their optical luminosity ($L_O$). FSRQs with $γ$-ray photon index greater than 2.6 exhibit very low optical variability, and their distribution of standard deviation of repeated photometry is significantly different from that of the other FSRQs (KS test P value equal to $5 \times 10^{-6}$ ). Among a sample of blazars at any particular cosmological epoch, those with lower $γ$-ray luminosity ($L_γ$) tend to have lower optical variability, and those FSRQs with $γ$-ray photon index greater than 2.6 tend to have low $L_γ$. We also measured the structure function of optical variability and found that the amplitude of the structure function for FSRQs is higher than previously measured and higher than that of BL Lacs at multiple time lags. Additionally, the amplitude of the structure function of FSRQs with high $γ$-ray photon index is significantly lower than that of FSRQs with low $γ$-ray photon index. The structure function of FSRQs of high $γ$-ray photon index shows a characteristic timescale of more than 10 days, which may be the variability timescale of the accretion disk. In summary, we infer that the optical component of FSRQs with high $γ$-ray photon index may be dominated by the accretion disk. △ Less

Submitted 20 May, 2024; originally announced May 2024.

Comments: 22 pages, 13 figures

X-ray Winds In Nearby-to-distant Galaxies (X-WING) - I: Legacy Surveys of Galaxies with Ultrafast Outflows and Warm Absorbers in $z \sim 0$-$4$

Authors: Satoshi Yamada, Taiki Kawamuro, Misaki Mizumoto, Claudio Ricci, Shoji Ogawa, Hirofumi Noda, Yoshihiro Ueda, Teruaki Enoto, Mitsuru Kokubo, Takeo Minezaki, Hiroaki Sameshima, Takashi Horiuchi, Shoichiro Mizukoshi

Abstract: As an inaugural investigation under the X-ray Winds In Nearby-to-distant Galaxies (X-WING) program, we assembled a dataset comprising 132 active galactic nuclei (AGNs) spanning redshifts $z \sim 0$-$4$ characterized by blueshifted absorption lines indicative of X-ray winds. Through an exhaustive review of previous research, we compiled the outflow parameters for 573 X-ray winds, encompassing key a… ▽ More As an inaugural investigation under the X-ray Winds In Nearby-to-distant Galaxies (X-WING) program, we assembled a dataset comprising 132 active galactic nuclei (AGNs) spanning redshifts $z \sim 0$-$4$ characterized by blueshifted absorption lines indicative of X-ray winds. Through an exhaustive review of previous research, we compiled the outflow parameters for 573 X-ray winds, encompassing key attributes such as outflow velocities ($V_{\rm out}$), ionization parameters ($ξ$), and hydrogen column densities. By leveraging the parameters $V_{\rm out}$ and $ξ$, we systematically categorized the winds into three distinct groups: ultrafast outflows (UFOs), low-ionization parameter (low-IP) UFOs, and warm absorbers. Strikingly, a discernible absence of linear correlations in the outflow parameters, coupled with distributions approaching instrumental detection limits, was observed. Another notable finding was the identification of a velocity gap around $V_{\rm out} \sim 10,000~{\rm km~s^{-1}}$. This gap was particularly evident in the winds detected via absorption lines within the $\lesssim$2 keV band, indicating disparate origins for low-IP UFOs and warm absorbers. In cases involving Fe XXV/Fe XXVI lines, where the gap might be attributed to potential confusion between emission/absorption lines and the Fe K-edge, the possibility of UFOs and galactic-scale warm absorbers being disconnected is considered. An examination of the outflow and dust sublimation radii revealed a distinction: UFOs appear to consist of dust-free material, whereas warm absorbers likely comprise dusty gas. From 2024, the X-Ray Imaging and Spectroscopy Mission (XRISM) is poised to alleviate observational biases, providing insights into the authenticity of the identified gap, a pivotal question in comprehending AGN feedback from UFOs. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: 34 pages, 4 tables, 7 figures, accepted in ApJS. The CSV files of Tables B2-B4 are available from the URL: http://www.kusastro.kyoto-u.ac.jp/~styamada/yamada2024bApJS/XWING-PaperI_tables.zip

The relationship of SMBHs and host galaxies at $z<4$ in the deep optical variability-selected AGN sample in the COSMOS field

Authors: Atsushi Hoshi, Toru Yamada, Mitsuru Kokubo, Yoshiki Matsuoka, Tohru Nagao

Abstract: We present the study on the relationship between SMBHs and their host galaxies using our variability-selected AGN sample ($i_\mathrm{AB} \leq 25.9,\ z \leq 4.5$) constructed from the HSC-SSP Ultra-deep survey in the COSMOS field. We estimated the BH mass ($M_\mathrm{BH}=10^{5.5-10}\ M_{\odot}$) based on the single-epoch virial method and the total stellar mass (… ▽ More We present the study on the relationship between SMBHs and their host galaxies using our variability-selected AGN sample ($i_\mathrm{AB} \leq 25.9,\ z \leq 4.5$) constructed from the HSC-SSP Ultra-deep survey in the COSMOS field. We estimated the BH mass ($M_\mathrm{BH}=10^{5.5-10}\ M_{\odot}$) based on the single-epoch virial method and the total stellar mass ($M_\mathrm{star}=10^{10-12}\ M_{\odot}$) by separating the AGN component with SED fitting. We found that the redshift evolution of the BH-stellar mass ratio ($M_\mathrm{BH}/M_\mathrm{star}$) depends on the $M_\mathrm{BH}$ which is caused by the no significant correlation between $M_\mathrm{BH}$ and $M_\mathrm{star}$. Variable AGNs with massive SMBHs ($M_\mathrm{BH}>10^{9}\ M_{\odot}$) at $1.5<z<3$ show considerably higher BH-stellar mass ratios ($>\sim1\%$) than the BH-bulge ratios ($M_\mathrm{BH}/M_\mathrm{bulge}$) observed in the local universe for the same BH range. This implies that there is a typical growth path of massive SMBHs which is faster than the formation of the bulge component as final products seen in the present day. For the low-mass SMBHs ($M_\mathrm{BH}<10^{8}\ M_{\odot}$) at $0.5<z<3$, on the other hand, variable AGNs show the similar BH-stellar mass ratios with the local objects ($\sim 0.1\%$) but smaller than those observed at $z > 4$. We interpret that host galaxies harboring less massive SMBHs at intermediate redshift have already acquired sufficient stellar mass, although high-z galaxies are still in the early stage of galaxy formation relative to those at the intermediate/local universe. △ Less

Submitted 21 April, 2024; originally announced April 2024.

Comments: Accepted for publication in ApJ, 14 pages, 5 figures

Journal ref: 2024, ApJ, 969, 11

Rayleigh and Raman scattering cross-sections and phase matrices of the ground-state hydrogen atom, and their astrophysical implications

Authors: Mitsuru Kokubo

Abstract: We present explicit expressions for Rayleigh and Raman scattering cross-sections and phase matrices of the ground $1s$ state hydrogen atom based on the Kramers-Heisenberg dispersion formula. The Rayleigh scattering leaves the hydrogen atom in the ground-state while the Raman scattering leaves the hydrogen atom in either $ns$ ($n\geq2$; $s$-branch) or $nd$ ($n\geq3$; $d$-branch) excited state, and… ▽ More We present explicit expressions for Rayleigh and Raman scattering cross-sections and phase matrices of the ground $1s$ state hydrogen atom based on the Kramers-Heisenberg dispersion formula. The Rayleigh scattering leaves the hydrogen atom in the ground-state while the Raman scattering leaves the hydrogen atom in either $ns$ ($n\geq2$; $s$-branch) or $nd$ ($n\geq3$; $d$-branch) excited state, and the Raman scattering converts incident ultraviolet (UV) photons around the Lyman resonance lines into optical-infrared (IR) photons. We show that this Raman wavelength conversion of incident flat UV continuum in dense hydrogen gas with a column density of $N_{\text{H}} > 10^{21}~\text{cm}^{-2}$ can produce broad emission features centred at Balmer, Paschen, and higher-level lines, which would mimic Doppler-broadened hydrogen lines with the velocity width of $\gtrsim 1,000~\text{km}~\text{s}^{-1}$ that could be misinterpreted as signatures of Active Galactic Nuclei, supernovae, or fast stellar winds. We show that the phase matrix of the Rayleigh and Raman $s$-branch scatterings is identical to that of the Thomson scattering while the Raman $d$-branch scattering is more isotropic, thus the Paschen and higher-level Raman features are depolarized compared to the Balmer features due to the flux contribution from the Raman $d$-branch. We argue that observations of the line widths, line flux ratios, and linear polarization of multiple optical/IR hydrogen lines are crucial to discriminate between the Raman-scattered broad emission features and Doppler-broadened emission lines. △ Less

Submitted 12 July, 2024; v1 submitted 9 August, 2023; originally announced August 2023.

Comments: 32 pages, 8 figures, 2 tables, accepted for publication in MNRAS (the journal version is more abbreviated than the arXiv version)

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 529, Issue 3, pp.2131-2149, April 2024

Active Massive Black Hole Found in the Young Star-Forming Dwarf Galaxy SBS 0335-052E

Authors: Shun Hatano, Masami Ouchi, Kimihiko Nakajima, Toshihiro Kawaguchi, Mitsuru Kokubo, Satoshi Kikuta, Nozomu Tominaga, Yi Xu, Kuria Watanabe, Yuichi Harikane, Yuki Isobe, Akinori Matsumoto, Moka Nishigaki, Yoshiaki Ono, Masato Onodera, Yuma Sugahara, Hiroya Umeda, Yechi Zhang

Abstract: Since the late 20th century, young star-forming dwarf galaxies with extremely low heavy-element abundances have been identified in the local universe\cite{1972ApJ...173...25S, 1990Natur.343..238I, 2022ApJS..262....3N}. Such a population of galaxies, represented by SBS 0335-052E, is intensively studied as a laboratory of galaxy formation, and interpreted as galaxies with hard emission and energetic… ▽ More Since the late 20th century, young star-forming dwarf galaxies with extremely low heavy-element abundances have been identified in the local universe\cite{1972ApJ...173...25S, 1990Natur.343..238I, 2022ApJS..262....3N}. Such a population of galaxies, represented by SBS 0335-052E, is intensively studied as a laboratory of galaxy formation, and interpreted as galaxies with hard emission and energetic outflows driven by massive stars\cite{1990Natur.343..238I, 2004ApJ...606..213T, 2009AA...503...61I}. Here we report the temporal flux variability of SBS 0335-052E in $3-4 \ μ$m bands on timescale of months to years with dimming and brightening up to 50\% over 12 years. This is a clear signature of dust torus emission of an active massive black hole in SBS 0335-052E. The deep optical spectrum reveals a very broad component with $1.24\ (\pm 0.01) \times 10^{4}$ km s$^{-1}$ in H$α$ emission, suggesting emission originated from broad line regions around the massive black hole. The black hole mass is estimated to be $\sim 100$ thousand solar masses from the spectral energy distribution. This is the first example of the massive black hole existence in a long-considered young star-forming dwarf galaxy. △ Less

Submitted 7 April, 2023; originally announced April 2023.

Comments: 32 pages, 6 figures, 3 tables

eROSITA Final Equatorial-Depth Survey (eFEDS): eFEDS X-ray view of WERGS Radio Galaxies selected by the Subaru/HSC and VLA/FIRST survey

Authors: Kohei Ichikawa, Takuji Yamashita, Andrea Merloni, Junyao Li, Teng Liu, Mara Salvato, Masayuki Akiyama, Riccardo Arcodia, Tom Dwelly, Xiaoyang Chen, Masatoshi Imanishi, Kohei Inayoshi, Toshihiro Kawaguchi, Taiki Kawamuro, Mitsuru Kokubo, Yoshiki Matsuoka, Tohru Nagao, Malte Schramm, Hyewon Suh, Masayuki Tanaka, Yoshiki Toba, Yoshihiro Ueda

Abstract: We construct the eROSITA X-ray catalog of radio galaxies discovered by the WERGS survey that is made by the cross-matching of the wide-area Subaru/HSC optical survey and VLA/FIRST 1.4 GHz radio survey. We find 393 eROSITA detected radio galaxies in the 0.5--2 keV band in the eFEDS field covering 140~deg$^2$. Thanks to the wide and medium depth eFEDS X-ray survey, the sample contains the rare and m… ▽ More We construct the eROSITA X-ray catalog of radio galaxies discovered by the WERGS survey that is made by the cross-matching of the wide-area Subaru/HSC optical survey and VLA/FIRST 1.4 GHz radio survey. We find 393 eROSITA detected radio galaxies in the 0.5--2 keV band in the eFEDS field covering 140~deg$^2$. Thanks to the wide and medium depth eFEDS X-ray survey, the sample contains the rare and most X-ray luminous radio galaxies above the knee of the X-ray luminosity function, spanning 44<log L(0.5-2keV,abs)<46.5 at $1<z<4$. Based on the X-ray properties obtained by the spectral fitting, 37 sources show obscured AGN signature with $\log (N_\mathrm{H}/\mathrm{cm}^{-2})>22$. Those obscured and radio AGN reside in $0.4<z<3.2$, indicating that they are obscured counterparts of the radio-loud quasar, which are missed in the previous optical quasar surveys. By combining radio and X-ray luminosities, the jet production efficiency $η_\mathrm{jet}$ is investigated, and we find 14 sources with extremely high jet production efficiency at $η_\mathrm{jet}\approx1$. This high $η_\mathrm{jet}$ value might be a result of 1) the decreased radiation efficiency of $η_\mathrm{rad}<0.1$ due to the low accretion rate for those sources and/or 2) the boosting due to the decline of $L_\mathrm{bol}$ by a factor of 10--100 by keeping $P_\mathrm{jet}$ constant in the previous Myr, indicating the experience of the AGN feedback. Finally, inferring the BH masses from the stellar-mass, we find that X-ray luminous sources show the excess of the radio emission with respect to the value estimated from the fundamental plane. Such radio emission excess cannot be explained by the Doppler booming alone, and therefore disk-jet connection of X-ray luminous eFEDS-WERGS is fundamentally different from the conventional fundamental plane which mainly covers low accretion regime. △ Less

Submitted 15 February, 2023; originally announced February 2023.

Comments: 20 pages, 11 figures, accepted for publication in A&A

Journal ref: A&A 672, A171 (2023)

Narrow Fe-K$α$ Reverberation Mapping Unveils the Deactivated Broad-Line Region in a Changing-Look Active Galactic Nucleus

Authors: Hirofumi Noda, Taisei Mineta, Takeo Minezaki, Hiroaki Sameshima, Mitsuru Kokubo, Taiki Kawamuro, Satoshi Yamada, Takashi Horiuchi, Hironori Matsumoto, Makoto Watanabe, Kumiko Morihana, Yoichi Itoh, Koji S. Kawabata, Yasushi Fukazawa

Abstract: "Changing-look active galactic nuclei" (CLAGNs) are known to change their apparent types between types 1 and 2, usually accompanied by a drastic change in their luminosity on timescales of years. However, it is still unclear whether materials in broad-line regions (BLRs) in CLAGNs appear and disappear during the type-transition or remain at the same location while the line production is simply act… ▽ More "Changing-look active galactic nuclei" (CLAGNs) are known to change their apparent types between types 1 and 2, usually accompanied by a drastic change in their luminosity on timescales of years. However, it is still unclear whether materials in broad-line regions (BLRs) in CLAGNs appear and disappear during the type-transition or remain at the same location while the line production is simply activated or deactivated. Here we present our X-ray-optical monitoring results of a CLAGN, NGC 3516, by Suzaku, Swift, and ground telescopes, with our primary focus on the narrow Fe-K$α$ emission line, which is an effective probe of the BLR materials. We detected significant variations of the narrow Fe-K$α$ line on a timescale of tens of days during the type-2 (faint) phase in 2013-2014, and conducted "narrow Fe-K$α$ reverberation mapping," comparing its flux variation with those of the X-ray continuum from a corona and $B$-band continuum from an accretion disk. We derived, as a result, a time lag of $10.1^{+5.8}_{-5.6}$ days ($1σ$ errors) for the Fe-K$α$ line behind the continuum, which is consistent with the location of the BLR determined in optical spectroscopic reverberation mapping during the type-1 (bright) phase. This finding shows that the BLR materials remained at the same location without emitting optical broad-lines during the type-2 phase. Considering the drastic decrease of the radiation during the type-transition, our result is possibly inconsistent with the hotly-discussed formation models of the BLR which propose that the radiative pressure from an accretion disk should be the main driving force. △ Less

Submitted 5 December, 2022; originally announced December 2022.

Comments: 15 pages, 7 figures, 3 tables, accepted for the publication in ApJ

H1821+643: The most X-ray and infrared luminous AGN in the Swift/BAT survey in the process of rapid stellar and supermassive black hole mass assembly

Authors: Hikaru Fukuchi, Kohei Ichikawa, Masayuki Akiyama, Claudio Ricci, Sunmyon Chon, Mitsuru Kokubo, Ang Liu, Takuya Hashimoto, Takuma Izumi

Abstract: H1821+643 is the most X-ray luminous non-beamed AGN of $L_\mathrm{14-150 keV}= 5.2\times 10^{45}$ erg s$^{-1}$ in the Swift/BAT ultra-hard X-ray survey and it is also a hyper-luminous infrared (IR) galaxy $L_\mathrm{IR} = 10^{13.2} L_\odot$ residing in the center of a massive galaxy cluster, which is a unique environment achieving the rapid mass assembly of black holes (BH) and host galaxies in th… ▽ More H1821+643 is the most X-ray luminous non-beamed AGN of $L_\mathrm{14-150 keV}= 5.2\times 10^{45}$ erg s$^{-1}$ in the Swift/BAT ultra-hard X-ray survey and it is also a hyper-luminous infrared (IR) galaxy $L_\mathrm{IR} = 10^{13.2} L_\odot$ residing in the center of a massive galaxy cluster, which is a unique environment achieving the rapid mass assembly of black holes (BH) and host galaxies in the local universe. We decompose the X-ray to IR spectral energy distribution (SED) into the AGN and starburst component using the SED fitting tool CIGALE-2022.0 and show that H1821+643 consumes a large amount of cold gas ($\dot{M}_\mathrm{con}$) with star-formation rate of $\log ( \mathrm{SFR}/M_{\odot}~\mathrm{yr}^{-1}) = 3.01 \pm 0.04$ and BH accretion rate of $\log (\dot{M}_\mathrm{BH}/M_{\odot}~\mathrm{yr}^{-1}) = 1.20 \pm 0.05$. This high $\dot{M}_\mathrm{con}$ is larger than the cooling rate ($\dot{M}_\mathrm{cool}$) of the intra-cluster medium (ICM), $\dot{M}_\mathrm{con}/\dot{M}_\mathrm{cool} \gtrsim 1$, which is one to two order magnitude higher than the typical value of other systems, indicating that H1821 provides the unique and extreme environment of rapid gas consumption. We also show that H1821+643 has an efficient cooling path achieving from $10^7$ K to $10^2$ K thanks to [OIII] 63 $μ\mathrm{m}$, which is a main coolant in low temperature range ($10^4$ K to $10^2$ K) with a cooling rate of $\dot{M}_{\mathrm{cool}}=3.2\times 10^5\ M_{\odot}\mathrm{~yr^{-1}}$, and the star-forming region extends over 40 kpc scale. △ Less

Submitted 19 September, 2022; originally announced September 2022.

Comments: 20 pages, 8 figures, accepted for publication in ApJ

Measurement of AGN dust extinction based on the near-infrared flux variability of WISE data

Authors: Shoichiro Mizukoshi, Takeo Minezaki, Shoichi Tsunetsugu, Atsuhiro Yoshida, Hiroaki Sameshima, Mitsuru Kokubo, Hirofumi Noda

Abstract: We present the measurement of the line-of-sight extinction of the dusty torus for a large number of obscured active galactic nuclei (AGNs) based on the reddening of the colour of the variable flux component in near-infrared (NIR) wavelengths. We collected long-term monitoring data by $\textit{Wide-field Infrared Survey Explorer (WISE)}$ for 513 local AGNs catalogued by the $\mathit{Swift/}$BAT AGN… ▽ More We present the measurement of the line-of-sight extinction of the dusty torus for a large number of obscured active galactic nuclei (AGNs) based on the reddening of the colour of the variable flux component in near-infrared (NIR) wavelengths. We collected long-term monitoring data by $\textit{Wide-field Infrared Survey Explorer (WISE)}$ for 513 local AGNs catalogued by the $\mathit{Swift/}$BAT AGN Spectroscopic Survey (BASS) and found that the multi-epoch NIR flux data in two different bands (WISE $W1$ and $W2$) are tightly correlated for more than 90% of the targets. The flux variation gradient (FVG) in the $W1$ and $W2$ bands was derived by applying linear regression analysis, and we reported that those for unobscured AGNs fall in a relatively narrow range, whereas those for obscured AGNs are distributed in a redder and broader range. The AGN's line-of-sight dust extinction ($A_V$) is calculated using the amount of the reddening in the FVG and is compared with the neutral hydrogen column density ($N_{\rm{}H}$) of the BASS catalogue. We found that the $N_{\rm{}H}/A_V$ ratios of obscured AGNs are greater than those of the Galactic diffuse interstellar medium (ISM) and are distributed with a large scatter by at most two orders of magnitude. Furthermore, we found that the lower envelope of the $N_{\rm{}H}/A_V$ of obscured AGNs is comparable to the Galactic diffuse ISM. These properties of the $N_{\rm{}H}/A_V$ can be explained by increase in the $N_{\rm{}H}$ attributed to the dust-free gas clouds covering the line of sight in the broad-line region. △ Less

Submitted 16 September, 2022; originally announced September 2022.

Comments: 11 pages, 7 figures, published in MNRAS

Journal ref: MNRAS, 516, 2876 (2022)

Finding of a population of active galactic nuclei showing a significant luminosity decline in the past $\sim 10^{3-4}$ yrs

Authors: J. Pflugradt, K. Ichikawa, M. Akiyama, M. Kokubo, B. Vijarnwannaluk, H. Noda, X. Chen

Abstract: Recent observations have revealed an interesting active galactic nuclei (AGN) subclass that shows strong activity at large scales ($\sim1$ kpc) but weaker at small scales ($<10$ pc), suggesting a strong change in the mass accretion rate of the central engine in the past $10^{3-4}$ yr. We systematically search for such declining or fading AGN by cross-matching the SDSS type-1 AGN catalog at… ▽ More Recent observations have revealed an interesting active galactic nuclei (AGN) subclass that shows strong activity at large scales ($\sim1$ kpc) but weaker at small scales ($<10$ pc), suggesting a strong change in the mass accretion rate of the central engine in the past $10^{3-4}$ yr. We systematically search for such declining or fading AGN by cross-matching the SDSS type-1 AGN catalog at $z<0.4$, covering the [OIII]$\lambda5007$ emission line which is a tracer for the narrow-line region (NLR) emission, with the WISE mid-infrared (MIR) catalog covering the emissions from the dusty tori. Out of the 7,653 sources, we found 57 AGN whose bolometric luminosities estimated from the MIR band are at least one order of magnitude fainter than those estimated from the [OIII]$\lambda5007$ emission line. This luminosity declining AGN candidate population shows four important properties: 1) the past AGN activity estimated from the [OIII]$\lambda5007$ line reaches around the Eddington-limit, 2) more than 30% of the luminosity declining AGN candidates show a large absolute variability of $ΔW1 > 0.45$ mag in the previous $\sim10$ yr at the WISE 3.4 $μ$m band, 3) the median ratio of $\log$([NII]$\lambda6584/ \mathrm{H}α\lambda6563)=-0.52$, suggesting a lower gas metallicity and/or higher ionization parameter compared to other AGN populations. 4) the second epoch spectra of the population indicate a spectral type change for 15% of the sources. This population provides insights on the possible connection between the luminosity decline which started $\sim10^{3-4}$ yr ago and the decline in the recent $10$ yr. △ Less

Submitted 25 August, 2022; originally announced August 2022.

Comments: 21 pages, 10 figures, 2 tables (1 MRT), accepted for publication in ApJ

MUSSES2020J: The Earliest Discovery of a Fast Blue Ultraluminous Transient at Redshift 1.063

Authors: Ji-an Jiang, Naoki Yasuda, Keiichi Maeda, Nozomu Tominaga, Mamoru Doi, Željko Ivezić, Peter Yoachim, Kohki Uno, Takashi J. Moriya, Brajesh Kumar, Yen-Chen Pan, Masayuki Tanaka, Masaomi Tanaka, Ken'ichi Nomoto, Saurabh W. Jha, Pilar Ruiz-Lapuente, David Jones, Toshikazu Shigeyama, Nao Suzuki, Mitsuru Kokubo, Hisanori Furusawa, Satoshi Miyazaki, Andrew J. Connolly, D. K. Sahu, G. C. Anupama

Abstract: In this Letter, we report the discovery of an ultraluminous fast-evolving transient in rest-frame UV wavelengths, MUSSES2020J, soon after its occurrence by using the Hyper Suprime-Cam (HSC) mounted on the 8.2 m Subaru telescope. The rise time of about 5 days with an extremely high UV peak luminosity shares similarities to a handful of fast blue optical transients whose peak luminosities are compar… ▽ More In this Letter, we report the discovery of an ultraluminous fast-evolving transient in rest-frame UV wavelengths, MUSSES2020J, soon after its occurrence by using the Hyper Suprime-Cam (HSC) mounted on the 8.2 m Subaru telescope. The rise time of about 5 days with an extremely high UV peak luminosity shares similarities to a handful of fast blue optical transients whose peak luminosities are comparable with the most luminous supernovae while their timescales are significantly shorter (hereafter "fast blue ultraluminous transient," FBUT). In addition, MUSSES2020J is located near the center of a normal low-mass galaxy at a redshift of 1.063, suggesting a possible connection between the energy source of MUSSES2020J and the central part of the host galaxy. Possible physical mechanisms powering this extreme transient such as a wind-driven tidal disruption event and an interaction between supernova and circumstellar material are qualitatively discussed based on the first multiband early-phase light curve of FBUTs, although whether the scenarios can quantitatively explain the early photometric behavior of MUSSES2020J requires systematical theoretical investigations. Thanks to the ultrahigh luminosity in UV and blue optical wavelengths of these extreme transients, a promising number of FBUTs from the local to the high-z universe can be discovered through deep wide-field optical surveys in the near future. △ Less

Submitted 10 June, 2022; v1 submitted 30 May, 2022; originally announced May 2022.

Comments: 14 pages, 5 figures, accepted for publication in The Astrophysical Journal Letters

Journal ref: ApJL 933, L36 (2022)

Third Data Release of the Hyper Suprime-Cam Subaru Strategic Program

Authors: Hiroaki Aihara, Yusra AlSayyad, Makoto Ando, Robert Armstrong, James Bosch, Eiichi Egami, Hisanori Furusawa, Junko Furusawa, Sumiko Harasawa, Yuichi Harikane, Bau-Ching Hsieh, Hiroyuki Ikeda, Kei Ito, Ikuru Iwata, Tadayuki Kodama, Michitaro Koike, Mitsuru Kokubo, Yutaka Komiyama, Xiangchong Li, Yongming Liang, Yen-Ting Lin, Robert H. Lupton, Nate B Lust, Lauren A. MacArthur, Ken Mawatari , et al. (42 additional authors not shown)

Abstract: The paper presents the third data release of Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP), a wide-field multi-band imaging survey with the Subaru 8.2m telescope. HSC-SSP has three survey layers (Wide, Deep, and UltraDeep) with different area coverages and depths, designed to address a wide array of astrophysical questions. This third release from HSC-SSP includes data from 278 nights of ob… ▽ More The paper presents the third data release of Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP), a wide-field multi-band imaging survey with the Subaru 8.2m telescope. HSC-SSP has three survey layers (Wide, Deep, and UltraDeep) with different area coverages and depths, designed to address a wide array of astrophysical questions. This third release from HSC-SSP includes data from 278 nights of observing time and covers about 670 square degrees in all five broad-band filters at the full depth ($\sim26$~mag at $5σ$) in the Wide layer. If we include partially observed area, the release covers 1,470 square degrees. The Deep and UltraDeep layers have $\sim80\%$ of the originally planned integration times, and are considered done, as we have slightly changed the observing strategy in order to compensate for various time losses. There are a number of updates in the image processing pipeline. Of particular importance is the change in the sky subtraction algorithm; we subtract the sky on small scales before the detection and measurement stages, which has significantly reduced false detections. Thanks to this and other updates, the overall quality of the processed data has improved since the previous release. However, there are limitations in the data (for example, the pipeline is not optimized for crowded fields), and we encourage the user to check the quality assurance plots as well as a list of known issues before exploiting the data. The data release website is https://hsc-release.mtk.nao.ac.jp/. △ Less

Submitted 30 August, 2021; originally announced August 2021.

Comments: 25 pages, 19 figures, submitted to PASJ. Data available at https://hsc-release.mtk.nao.ac.jp/

Reverberation Mapping of Two Luminous Quasars: the Broad-line Region Structure and Black Hole Mass

Authors: Sha-Sha Li, Sen Yang, Zi-Xu Yang, Yong-Jie Chen, Yu-Yang Songsheng, He-Zhen Liu, Pu Du, Bin Luo, Zhe Yu, Chen Hu, Bo-Wei Jiang, Dong-Wei Bao, Wei-Jian Guo, Zhi-Xiang Zhang, Yan-Rong Li, Ming Xiao, Kai-Xing Lu, Luis C. Ho, Jing-Min Bai, Wei-Hao Bian, Jesús Aceituno, Takeo Minezaki, Mitsuru Kokubo, Jian-Min Wang

Abstract: We report the results of a multi-year spectroscopic and photometric monitoring campaign of two luminous quasars, PG~0923+201 and PG~1001+291, both located at the high-luminosity end of the broad-line region (BLR) size-luminosity relation with optical luminosities above $10^{45}~{\rm erg~s^{-1}}$. PG~0923+201 is for the first time monitored, and PG~1001+291 was previously monitored but our campaign… ▽ More We report the results of a multi-year spectroscopic and photometric monitoring campaign of two luminous quasars, PG~0923+201 and PG~1001+291, both located at the high-luminosity end of the broad-line region (BLR) size-luminosity relation with optical luminosities above $10^{45}~{\rm erg~s^{-1}}$. PG~0923+201 is for the first time monitored, and PG~1001+291 was previously monitored but our campaign has a much longer temporal baseline. We detect time lags of variations of the broad H$β$, H$γ$, Fe {\sc ii} lines with respect to those of the 5100~Å continuum. The velocity-resolved delay map of H$β$ in PG~0923+201 indicates a complicated structure with a mix of Keplerian disk-like motion and outflow, and the map of H$β$ in PG~1001+291 shows a signature of Keplerian disk-like motion. Assuming a virial factor of $f_{\rm BLR}=1$ and FWHM line widths, we measure the black hole mass to be $118_{-16}^{+11}\times 10^7 M_{\odot}$ for PG~0923+201 and $3.33_{-0.54}^{+0.62}\times 10^7 M_{\odot}$ for PG~1001+291. Their respective accretion rates are estimated to be $0.21_{-0.07}^{+0.06} \times L_{\rm Edd}\,c^{-2}$ and $679_{-227}^{+259}\times L_{\rm Edd}\,c^{-2}$, indicating that PG~0923+201 is a sub-Eddington accretor and PG~1001+291 is a super-Eddington accretor. While the H$β$ time lag of PG~0923+201 agrees with the size-luminosity relation, the time lag of PG~1001+291 shows a significant deviation, confirming that in high-luminosity AGN the BLR size depends on both luminosity and Eddington ratio. Black hole mass estimates from single AGN spectra will be over-estimated at high luminosities and redshifts if this effect is not taken into account. △ Less

Submitted 10 June, 2021; originally announced June 2021.

Comments: 21 pages, 14 figures, accepted

Journal ref: 2021, The Astrophysical Journal

Multiple giant eruptions and X-ray emission in the recoiling AGN/LBV candidate SDSS1133

Authors: Mitsuru Kokubo

Abstract: We present a comprehensive analysis of 20 years worth of multi-color photometric light curves, multi-epoch optical spectra, and X-ray data of an off-nuclear variable object SDSS1133 in Mrk 177 at $z=0.0079$. The UV-optical light curves reveal that SDSS1133 experienced four outbursts in 2001, 2014, 2019, and 2021. The persistent UV-optical luminosity in the non-outbursting state is $\sim 10^{41}$ e… ▽ More We present a comprehensive analysis of 20 years worth of multi-color photometric light curves, multi-epoch optical spectra, and X-ray data of an off-nuclear variable object SDSS1133 in Mrk 177 at $z=0.0079$. The UV-optical light curves reveal that SDSS1133 experienced four outbursts in 2001, 2014, 2019, and 2021. The persistent UV-optical luminosity in the non-outbursting state is $\sim 10^{41}$ erg/s with small-scale flux variations, and peak luminosities during the outbursts reach $\sim 10^{42}$ erg/s. The optical spectra exhibit enduring broad hydrogen Balmer P-Cygni profiles with the absorption minimum at $\sim -2,000$ km/s, indicating the presence of fast moving ejecta. {\it Chandra} detected weak X-ray emission at a $0.3-10$ keV luminosity of $L_{X} = 4 \times 10^{38}$ erg/s after the 2019 outburst. These lines of evidence suggests that SDSS1133 is an extreme luminous blue variable (LBV) star experiencing multiple giant eruptions with interactions of the ejected shell with different shells and/or circumstellar medium (CSM), and disfavors the recoiling Active Galactic Nuclei (AGN) scenario suggested in the literature. We suggest that pulsational pair-instability may provide a viable explanation for the multiple energetic eruptions in SDSS1133. If the current activity of SDSS1133 is a precursor of a supernova explosion, we may be able to observe a few additional giant eruptions and then the terminal supernova explosion or collapse to a massive black hole in future observations. △ Less

Submitted 16 June, 2022; v1 submitted 19 January, 2021; originally announced January 2021.

Comments: Accepted for publication in MNRAS. 30 pages, 17 figures, 5 tables. Table 1 and optical spectra are available online as ancillary files for this preprint

Constraints on the rate of supernovae lasting for more than a year from Subaru/Hyper Suprime-Cam

Authors: Takashi J. Moriya, Ji-an Jiang, Naoki Yasuda, Mitsuru Kokubo, Kojiro Kawana, Keiichi Maeda, Yen-Chen Pan, Robert M. Quimby, Nao Suzuki, Ichiro Takahashi, Masaomi Tanaka, Nozomu Tominaga, Ken'ichi Nomoto, Jeff Cooke, Lluis Galbany, Santiago Gonzalez-Gaitan, Chien-Hsiu Lee, Giuliano Pignata

Abstract: Some supernovae such as pair-instability supernovae are predicted to have the duration of more than a year in the observer frame. To constrain the rates of supernovae lasting for more than a year, we conducted a long-term deep transient survey using Hyper Suprime-Cam (HSC) on the 8.2m Subaru telescope. HSC is a wide-field (a 1.75 deg2 field-of-view) camera and it can efficiently conduct transient… ▽ More Some supernovae such as pair-instability supernovae are predicted to have the duration of more than a year in the observer frame. To constrain the rates of supernovae lasting for more than a year, we conducted a long-term deep transient survey using Hyper Suprime-Cam (HSC) on the 8.2m Subaru telescope. HSC is a wide-field (a 1.75 deg2 field-of-view) camera and it can efficiently conduct transient surveys. We observed the same 1.75 deg2 field repeatedly using the g, r, i, and z band filters with the typical depth of 26 mag for 4 seasons (from late 2016 to early 2020). Using these data, we searched for transients lasting for more than a year. Two supernovae were detected in 2 continuous seasons, one supernova was detected in 3 continuous seasons, but no transients lasted for all 4 seasons searched. The discovery rate of supernovae lasting for more than a year with the typical limiting magnitudes of 26 mag is constrained to be 1.4^{+1.3}_{-0.7}(stat.)^{+0.2}_{-0.3}(sys.) events deg-2 yr-1. All the long-lasting supernovae we found are likely Type IIn supernovae and our results indicate that about 40% of Type IIn supernovae have long-lasting light curves. No plausible pair-instability supernova candidates lasting for more than a year are discovered. By comparing the survey results and survey simulations, we constrain the luminous pair-instability supernova rate up to z ~ 3 should be of the order of 100 Gpc-3 yr-1 at most, which is 0.01 - 0.1 per cent of the core-collapse supernova rate. △ Less

Submitted 30 November, 2020; originally announced December 2020.

Comments: 18 pages, 9 figures, 6 tables, accepted by The Astrophysical Journal

Journal ref: The Astrophysical Journal, Volume 908, Issue 2, id.249, 13 pp. (2021)

Relationship between Radar Cross Section and Optical Magnitude based on Radar and Optical Simultaneous Observations of Faint Meteors

Authors: Ryou Ohsawa, Akira Hirota, Kohei Morita, Shinsuke Abe, Daniel Kastinen, Johan Kero, Csilla Szasz, Yasunori Fujiwara, Takuji Nakamura, Koji Nishimura, Shigeyuki Sako, Jun-ichi Watanabe, Tsutomu Aoki, Noriaki Arima, Ko Arimatsu, Mamoru Doi, Makoto Ichiki, Shiro Ikeda, Yoshifusa Ita, Toshihiro Kasuga, Naoto Kobayashi, Mitsuru Kokubo, Masahiro Konishi, Hiroyuki Maehara, Takashi Miyata , et al. (18 additional authors not shown)

Abstract: Radar and optical simultaneous observations of meteors are important to understand the size distribution of the interplanetary dust. However, faint meteors detected by high power large aperture radar observations, which are typically as faint as 10 mag. in optical, have not been detected until recently in optical observations, mainly due to insufficient sensitivity of the optical observations. In… ▽ More Radar and optical simultaneous observations of meteors are important to understand the size distribution of the interplanetary dust. However, faint meteors detected by high power large aperture radar observations, which are typically as faint as 10 mag. in optical, have not been detected until recently in optical observations, mainly due to insufficient sensitivity of the optical observations. In this paper, two radar and optical simultaneous observations were organized. The first observation was carried out in 2009 to 2010 using Middle and Upper Atmosphere Radar (MU radar) and an image-intensified CCD camera. The second observation was carried out in 2018 using the MU radar and a mosaic CMOS camera, Tomo-e Gozen, mounted on the 1.05-m Kiso Schmidt Telescope. In total, 331 simultaneous meteors were detected. The relationship between radar cross sections and optical V-band magnitudes was well approximated by a linear function. A transformation function from the radar cross section to the V-band magnitude was derived for sporadic meteors. The transformation function was applied to about 150,000 meteors detected by the MU radar in 2009--2015, large part of which are sporadic, and a luminosity function was derived in the magnitude range of $-1.5$ to $9.5$ mag. The luminosity function was well approximated by a single power-law function with the population index of $r = 3.52{\pm}0.12$. The present observation indicates that the MU radar has capability to detect interplanetary dust of $10^{-5}$ to $10^{0}$ g in mass as meteors. △ Less

Submitted 5 February, 2021; v1 submitted 20 August, 2020; originally announced August 2020.

Comments: Accepted for publication in Planetary & Space Science. 20 pages, 8 figures

Journal ref: Planetary and Space Science 194 (2020) 10511

Dust Reverberation Mapping of Type 2 AGN NGC 2110 Realized with X-ray and 3-5 $μ$m IR monitoring

Authors: Hirofumi Noda, Taiki Kawamuro, Mitsuru Kokubo, Takeo Minezaki

Abstract: The dust reverberation mapping is one of powerful methods to investigate the structure of the dusty tori in AGNs, and it has been performed on more than a hundred type 1 AGNs. However, no clear results have been reported on type 2 AGNs because their strong optical-UV extinction completely hides their accretion disc emission. Here we focus on an X-ray-bright type 2 AGN, NGC 2110, and utilize 2-20 k… ▽ More The dust reverberation mapping is one of powerful methods to investigate the structure of the dusty tori in AGNs, and it has been performed on more than a hundred type 1 AGNs. However, no clear results have been reported on type 2 AGNs because their strong optical-UV extinction completely hides their accretion disc emission. Here we focus on an X-ray-bright type 2 AGN, NGC 2110, and utilize 2-20 keV X-ray variation monitored by MAXI to trace disc emission, instead of optical-UV variation. Comparing it with light curves in the WISE infrared (IR) W1 band ($λ= 3.4$ $μ$m) and W2 band ($λ= 4.6$ $μ$m) with cross-correlation analyses, we found candidates of the dust reverberation time lag at $\sim60$ days, $\sim130$ days, and $\sim1250$ days between the X-ray flux variation and those of the IR bands. By examining the best-fitting X-ray and IR light curves with the derived time lags, we found that the time lag of $\sim130$ days is most favoured. With this time lag, the relation between the time lag and luminosity of NGC 2110 is consistent with those in type 1 AGNs, suggesting that the dust reverberation in NGC 2110 mainly originates in hot dust in the torus innermost region, the same as in type 1 AGNs. As demonstrated by the present study, X-ray and IR simultaneous monitoring can be a promising tool to perform the dust reverberation mapping on type 2 AGNs. △ Less

Submitted 14 May, 2020; originally announced May 2020.

Comments: 10 pages, 7 figures, 1 table, accepted for the publication in MNRAS

Properties of AGN Multiband Optical Variability in the HSC SSP Transient Survey

Authors: Yuki Kimura, Toru Yamada, Mitsuru Kokubo, Naoki Yasuda, Tomoki Morokuma, Tohru Nagao, Yoshiki Matsuoka

Abstract: We study variability of active galactic nuclei (AGNs) by using the deep optical multiband photometry data obtained from the Hyper Suprime-Cam Subaru Strategic Program (HSC SSP) survey in the COSMOS field. The images analyzed here were taken with 8, 10, 13, and 15 epochs over three years in the $g$, $r$, $i$, and $z$ bands, respectively. We identified 491 robust variable AGN candidates, down to… ▽ More We study variability of active galactic nuclei (AGNs) by using the deep optical multiband photometry data obtained from the Hyper Suprime-Cam Subaru Strategic Program (HSC SSP) survey in the COSMOS field. The images analyzed here were taken with 8, 10, 13, and 15 epochs over three years in the $g$, $r$, $i$, and $z$ bands, respectively. We identified 491 robust variable AGN candidates, down to $i=25$ mag and with redshift up to $4.26$. Ninety percent of the variability-selected AGNs are individually identified with the X-ray sources detected in the Chandra COSMOS Legacy survey. We investigate their properties in variability by using structure function analysis and find that the structure function for low-luminosity AGNs ($L_{\mathrm{bol}}\lesssim10^{45}$ erg s$^{-1}$) shows a positive correlation with luminosity, which is the opposite trend for the luminous quasars. This trend is likely to be caused by larger contribution of the host galaxy light for lower-luminosity AGNs. Using the model templates of galaxy spectra, we evaluate the amount of host galaxy contribution to the structure function analysis and find that dominance of the young stellar population is needed to explain the observed luminosity dependence. This suggests that low-luminosity AGNs at $0.8\lesssim z\lesssim1.8$ are predominantly hosted in star-forming galaxies. The X-ray stacking analysis reveals the significant emission from the individually X-ray undetected AGNs in our variability-selected sample. The stacked samples show very large hardness ratios in their stacked X-ray spectrum, which suggests that these optically variable sources have large soft X-ray absorption by dust-free gas. △ Less

Submitted 19 April, 2020; originally announced April 2020.

Comments: 27 pages, 22 figures. Accepted for publication in ApJ

Rapidly Evolving Transients from the Hyper Suprime-Cam SSP Transient Survey

Authors: Yusuke Tampo, Masaomi Tanaka, Keiichi Maeda, Naoki Yasuda, Nozomu Tominaga, Ji-an Jiang, Takashi J. Moriya, Tomoki Morokuma, Nao Suzuki, Ichiro Takahashi, Mitsuru Kokubo, Kojiro Kawana

Abstract: Rapidly evolving transients form a new class of transients which show shorter timescales of the light curves than those of typical core-collapse and thermonuclear supernovae. We performed a systematic search for rapidly evolving transients using the deep data taken with the Hyper Suprime-Cam Subaru Strategic Program Transient Survey. By measuring the timescales of the light curves of 1824 transien… ▽ More Rapidly evolving transients form a new class of transients which show shorter timescales of the light curves than those of typical core-collapse and thermonuclear supernovae. We performed a systematic search for rapidly evolving transients using the deep data taken with the Hyper Suprime-Cam Subaru Strategic Program Transient Survey. By measuring the timescales of the light curves of 1824 transients, we identified 5 rapidly evolving transients. Our samples are found in a wide range of redshifts (0.3 $\le$ z $\le$ 1.5) and peak absolute magnitudes ($-$17 $\ge$ $M_i$ $\ge$ $-$20). The properties of the light curves are similar to those of the previously discovered rapidly evolving transients. They show a relatively blue spectral energy distribution, with the best-fit blackbody of 8,000 - 18,000 K. We show that some of the transients require power sources other than the radioactive decays of $^{56}$Ni because of their high peak luminosities and short timescales. The host galaxies of all the samples are star-forming galaxies, suggesting a massive star origin for the rapidly evolving transients. The event rate is roughly estimated to be $\sim$4,000 events yr$^{-1}$ Gpc$^{-3}$, which is about 1 $\%$ of core-collapse supernovae. △ Less

Submitted 3 March, 2020; originally announced March 2020.

Comments: 19 pages, 13 figures, accepted for publication in ApJ

Dust Destruction by Charging: A Possible Origin of Grey Extinction Curves of Active Galactic Nuclei

Authors: Ryo Tazaki, Kohei Ichikawa, Mitsuru Kokubo

Abstract: Observed extinction curves of active galactic nuclei (AGNs) are significantly different from those observed in the Milky Way. The observations require preferential removal of small grains at the AGN environment; however, the physics for this remains unclear. In this paper, we propose that dust destruction by charging, or Coulomb explosion, may be responsible for AGN extinction curves. Harsh AGN ra… ▽ More Observed extinction curves of active galactic nuclei (AGNs) are significantly different from those observed in the Milky Way. The observations require preferential removal of small grains at the AGN environment; however, the physics for this remains unclear. In this paper, we propose that dust destruction by charging, or Coulomb explosion, may be responsible for AGN extinction curves. Harsh AGN radiation makes a dust grain highly charged through photoelectric emission, and grain fission via the Coulomb explosion occurs when the electrostatic tensile stress of a charge grain exceeds its tensile strength. We show that the Coulomb explosion can preferentially remove both small silicate and graphite grains and successfully reproduce both flat extinction curves and the absence of 2175Å~bump. △ Less

Submitted 19 February, 2020; originally announced February 2020.

Comments: 8 pages, 6 Figures; Accepted for publication in ApJ

Tracing the co-evolution path of super massive black holes and spheroids with AKARI-selected ultra-luminous IR galaxies at intermediate redshifts

Authors: Xiao-Yang Chen, Masayuki Akiyama, Kohei Ichikawa, Hirofumi Noda, Yoshiki Toba, Issei Yamamura, Toshihiro Kawaguchi, Abdurro'uf, Mitsuru Kokubo

Abstract: 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 resolu… ▽ More 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. △ Less

Submitted 13 July, 2020; v1 submitted 11 November, 2019; originally announced November 2019.

Comments: 28 pages, 25 figures, 7 tables; accepted by ApJ

An optical search for transients lasting a few seconds

Authors: Michael W. Richmond, Masaomi Tanaka, Tomoki Morokuma, Shigeyuki Sako, Ryou Ohsawa, Noriaki Arima, Nozomu Tominaga, Mamoru Doi, Tsutomu Aoki, Ko Arimatsu, Makoto Ichiki, Shiro Ikeda, Yoshifusa Ita, Toshihiro Kasuga, Koji S. Kawabata, Hideyo Kawakita, Naoto Kobayashi, Mitsuru Kokubo, Masahiro Konishi, Hiroyuki Maehara, Hiroyuki Mito, Takashi Miyata, Yuki Mori, Mikio Morii, Kentaro Motohara , et al. (15 additional authors not shown)

Abstract: Using a prototype of the Tomo-e Gozen wide-field CMOS mosaic camera, we acquire wide-field optical images at a cadence of 2 Hz and search them for transient sources of duration 1.5 to 11.5 seconds. Over the course of eight nights, our survey encompasses the equivalent of roughly two days on one square degree, to a fluence equivalent to a limiting magnitude about $V = 15.6$ in a 1-second exposure.… ▽ More Using a prototype of the Tomo-e Gozen wide-field CMOS mosaic camera, we acquire wide-field optical images at a cadence of 2 Hz and search them for transient sources of duration 1.5 to 11.5 seconds. Over the course of eight nights, our survey encompasses the equivalent of roughly two days on one square degree, to a fluence equivalent to a limiting magnitude about $V = 15.6$ in a 1-second exposure. After examining by eye the candidates identified by a software pipeline, we find no sources which meet all our criteria. We compute upper limits to the rate of optical transients consistent with our survey, and compare those to the rates expected and observed for representative sources of ephemeral optical light. △ Less

Submitted 24 October, 2019; originally announced October 2019.

Comments: Accepted for publication in PASJ, 26 pages, 10 figures. The posted version here is our original manuscript, before any changes in response to the referee's comments; to see a copy of the version which was accepted for publication, go to http://spiff.rit.edu/richmond/pub/tomoe_trans.pdf

New constraint on the atmosphere of (50000) Quaoar from a stellar occultation

Authors: Ko Arimatsu, Ryou Ohsawa, George L. Hashimoto, Seitaro Urakawa, Jun Takahashi, Miyako Tozuka, Yoichi Itoh, Misato Yamashita, Fumihiko Usui, Tsutomu Aoki, Noriaki Arima, Mamoru Doi, Makoto Ichiki, Shiro Ikeda, Yoshifusa Ita, Toshihiro Kasuga, Naoto Kobayashi, Mitsuru Kokubo, Masahiro Konishi, Hiroyuki Maehara, Noriyuki Matsunaga, Takashi Miyata, Mikio Morii, Tomoki Morokuma, Kentaro Motohara , et al. (13 additional authors not shown)

Abstract: We report observations of a stellar occultation by the classical Kuiper belt object (50000) Quaoar occurred on 28 June 2019. A single-chord high-cadence (2 Hz) photometry dataset was obtained with the Tomo-e Gozen CMOS camera mounted on the 1.05 m Schmidt telescope at Kiso Observatory. The obtained ingress and egress data do not show any indication of atmospheric refraction and allow to set new… ▽ More We report observations of a stellar occultation by the classical Kuiper belt object (50000) Quaoar occurred on 28 June 2019. A single-chord high-cadence (2 Hz) photometry dataset was obtained with the Tomo-e Gozen CMOS camera mounted on the 1.05 m Schmidt telescope at Kiso Observatory. The obtained ingress and egress data do not show any indication of atmospheric refraction and allow to set new $1σ$ and $3σ$ upper limits of 6 and 16 nbar, respectively, for the surface pressure of a pure methane atmosphere. These upper limits are lower than the saturation vapor pressure of methane at Quaoar's expected mean surface temperature ($T \sim 44$ K) and imply the absence of a $\sim$10 nbar-level global atmosphere formed by methane ice on Quaoar's surface. △ Less

Submitted 22 October, 2019; originally announced October 2019.

Comments: 12 pages, 3 figures, accepted for publication in the Astronomical Journal

Rapid Luminosity Decline and Subsequent Reformation of the Innermost Dust Distribution in the Changing-look AGN Mrk 590

Authors: Mitsuru Kokubo, Takeo Minezaki

Abstract: We examine the long-term optical/near-infrared (NIR) flux variability of a "changing-look" active galactic nucleus (AGN) Mrk 590 between 1998 and 2007. Multi-band multi-epoch optical/NIR photometry data from the SDSS Stripe 82 database and the Multicolor Active Galactic Nuclei Monitoring (MAGNUM) project reveal that Mrk 590 experienced a sudden luminosity decrease during the period from 2000 to 20… ▽ More We examine the long-term optical/near-infrared (NIR) flux variability of a "changing-look" active galactic nucleus (AGN) Mrk 590 between 1998 and 2007. Multi-band multi-epoch optical/NIR photometry data from the SDSS Stripe 82 database and the Multicolor Active Galactic Nuclei Monitoring (MAGNUM) project reveal that Mrk 590 experienced a sudden luminosity decrease during the period from 2000 to 2001. Detection of dust reverberation lag signals between $V$- and $K$-band light curves obtained by the MAGNUM project during the faint state in $2003-2007$ suggests that the dust torus innermost radius $R_\text{dust}$ of Mrk 590 had become very small [$R_\text{dust} \simeq 32$~ light-days (lt-days)] by the year 2004 according to the aforementioned significant decrease in AGN luminosity. The $R_\text{dust}$ in the faint state is comparable to the H$β$ broad line region (BLR) radius of $R_{\text{H}β, \text{BLR}} \simeq 26$ lt-days measured by previous reverberation mapping observations during the bright state of Mrk 590 in $1990-1996$. These observations indicate that the innermost radius of the dust torus in Mrk 590 decreased rapidly after the AGN ultraviolet-optical luminosity drop, and that the replenishment time scale of the innermost dust distribution is less than 4 years, which is much shorter than the free fall time scale of BLR gas or dust clouds. We suggest that rapid replenishment of the innermost dust distribution can be accomplished either by new dust formation in radiatively-cooled BLR gas clouds or by new dust formation in the disk atmosphere and subsequent vertical wind from the dusty disk as a result of radiation pressure. △ Less

Submitted 2 December, 2019; v1 submitted 18 April, 2019; originally announced April 2019.

Comments: 21 pages, 8 figures, 4 tables, accepted for publication in MNRAS 2019 November 29

Discovery of a strong ionized-gas outflow in an AKARI-selected Ultra-luminous Infrared Galaxy at z = 0.5

Authors: Xiaoyang Chen, Masayuki Akiyama, Hirofumi Noda, Abdurrouf, Yoshiki Toba, Issei Yamamura, Toshihiro Kawaguchi, Mitsuru Kokubo, Kohei Ichikawa

Abstract: In order to construct a sample of ultra-luminous infrared galaxies (ULIRGs, with infrared luminosity, $L_{\rm IR} > 10^{12}$ L$_{\odot}$) at 0.5 < z < 1, we are conducting an optical follow-up program for bright 90-$μ$m FIR sources with a faint optical (i < 20) counterpart selected in the AKARI Far-Infrared Surveyor (FIS) Bright Source catalog (Ver.2). AKARI-FIS-V2 J0916248+073034, identified as a… ▽ More In order to construct a sample of ultra-luminous infrared galaxies (ULIRGs, with infrared luminosity, $L_{\rm IR} > 10^{12}$ L$_{\odot}$) at 0.5 < z < 1, we are conducting an optical follow-up program for bright 90-$μ$m FIR sources with a faint optical (i < 20) counterpart selected in the AKARI Far-Infrared Surveyor (FIS) Bright Source catalog (Ver.2). AKARI-FIS-V2 J0916248+073034, identified as a ULIRG at z = 0.49 in the spectroscopic follow-up observation, indicates signatures of an extremely strong outflow in its emission line profiles. Its [OIII] 5007Å emission line shows FWHM of 1830 km s$^{-1}$ and velocity shift of -770 km s$^{-1}$ in relative to the stellar absorption lines. Furthermore, low-ionization [OII] 3726Å 3729Å doublet also shows large FWHM of 910 km s$^{-1}$ and velocity shift of -380 km s$^{-1}$. After the removal of an unresolved nuclear component, the long-slit spectroscopy 2D image possibly shows that the outflow extends to radius of 4 kpc. The mass outflow and energy ejection rates are estimated to be 500 M$_{\odot}$ yr$^{-1}$ and $4\times10^{44}$ erg s$^{-1}$, respectively, which imply that the outflow is among the most powerful ones observed in ULIRGs and QSOs at 0.3 < z < 1.6. The co-existence of the strong outflow and intense star formation (star formation rate of 990 M$_{\odot}$ yr$^{-1}$) indicates that the feedback of the strong outflow has not severely affect the star-forming region of the galaxy. △ Less

Submitted 16 January, 2019; originally announced January 2019.

Comments: 20 pages, 18 figures, accepted to be published in PASJ

A Long-Duration Luminous Type IIn Supernova KISS15s: Strong Recombination Lines from the Inhomogeneous Ejecta-CSM Interaction Region and Hot Dust Emission from Newly Formed Dust

Authors: Mitsuru Kokubo, Kazuma Mitsuda, Tomoki Morokuma, Nozomu Tominaga, Masaomi Tanaka, Takashi J. Moriya, Peter Yoachim, Željko Ivezić, Shigeyuki Sako, Mamoru Doi

Abstract: We report the discovery of an SN1988Z-like type IIn supernova KISS15s found in a low-mass star-forming galaxy at redshift z=0.038 during the course of the Kiso Supernova Survey (KISS). KISS15s shows long-duration optical continuum and emission line light curves, indicating that KISS15s is powered by a continuous interaction between the expanding ejecta and dense circumstellar medium (CSM). The H… ▽ More We report the discovery of an SN1988Z-like type IIn supernova KISS15s found in a low-mass star-forming galaxy at redshift z=0.038 during the course of the Kiso Supernova Survey (KISS). KISS15s shows long-duration optical continuum and emission line light curves, indicating that KISS15s is powered by a continuous interaction between the expanding ejecta and dense circumstellar medium (CSM). The H$α$ emission line profile can be decomposed into four Gaussians of narrow, intermediate, blue-shifted intermediate, and broad velocity width components, with a full width at half maximum of $\lesssim 100$, $\sim 2,000$, and $\sim 14,000$ km s${}^{-1}$ for the narrow, intermediate, and broad components, respectively. The presence of the blue-shifted intermediate component, of which the line-of-sight velocity relative to the systemic velocity is about $-5,000$ km s${}^{-1}$, suggests that the ejecta-CSM interaction region has an inhomogeneous morphology and anisotropic expansion velocity. We found that KISS15s shows increasing infrared continuum emission, which can be interpreted as hot dust thermal emission of $T \sim 1,200$ K from newly formed dust in a cool, dense shell in the ejecta-CSM interaction region. The progenitor mass-loss rate, inferred from bolometric luminosity, is $\dot{M} \sim 0.4 M_{\odot} \text{yr}^{-1} (v_{w}/40 \text{km}~\text{s}^{-1})$, where $v_{w}$ is the progenitor's stellar wind velocity. This implies that the progenitor of KISS15s was a red supergiant star or a luminous blue variable that had experienced a large mass-loss in the centuries before the explosion. △ Less

Submitted 16 January, 2019; originally announced January 2019.

Comments: 43 pages, 29 figures, and 12 tables, accepted for publication in the Astrophysical Journal (ApJ) 2019 January 15

Luminosity Function of Faint Sporadic Meteors measured with a Wide-Field CMOS mosaic camera Tomo-e PM

Authors: Ryou Ohsawa, Shigeyuki Sako, Yuki Sarugaku, Fumihiko Usui, Takafumi Ootsubo, Yasunori Fujiwara, Mikiya Sato, Toshihiro Kasuga, Ko Arimatsu, Jun-ichi Watanabe, Mamoru Doi, Naoto Kobayashi, Hidenori Takahashi, Kentaro Motohara, Tomoki Morokuma, Masahiro Konishi, Tsutomu Aoki, Takao Soyano, Ken'ichi Tarusawa, Yuki Mori, Yoshikazu Nakada, Makoto Ichiki, Noriaki Arima, Yuto Kojima, Masahiro Morita , et al. (12 additional authors not shown)

Abstract: Imaging observations of faint meteors were carried out on April 11 and 14, 2016 with a wide-field CMOS mosaic camera, Tomo-e PM, mounted on the 105-cm Schmidt telescope at Kiso Observatory, the University of Tokyo. Tomo-e PM, which is a prototype model of Tomo-e Gozen, can monitor a sky of ${\sim}1.98\,\mathrm{deg^2}$ at 2\,Hz. The numbers of detected meteors are 1514 and 706 on April 11 and 14, r… ▽ More Imaging observations of faint meteors were carried out on April 11 and 14, 2016 with a wide-field CMOS mosaic camera, Tomo-e PM, mounted on the 105-cm Schmidt telescope at Kiso Observatory, the University of Tokyo. Tomo-e PM, which is a prototype model of Tomo-e Gozen, can monitor a sky of ${\sim}1.98\,\mathrm{deg^2}$ at 2\,Hz. The numbers of detected meteors are 1514 and 706 on April 11 and 14, respectively. The detected meteors are attributed to sporadic meteors. Their absolute magnitudes range from $+4$ to $+10\,\mathrm{mag}$ in the $V$-band, corresponding to about $8.3{\times}10^{-2}$ to $3.3{\times}10^{-4}\,\mathrm{g}$ in mass. The present magnitude distributions we obtained are well explained by a single power-law luminosity function with a slope parameter $r = 3.1{\pm}0.4$ and a meteor rate $\log_{10}N_0 = -5.5{\pm}0.5$. The results demonstrate a high performance of telescopic observations with a wide-field video camera to constrain the luminosity function of faint meteors. The performance of Tomo-e Gozen is about two times higher than that of Tomo-e PM. A survey with Tomo-e Gozen will provide a more robust measurement of the luminosity function. △ Less

Submitted 24 September, 2018; originally announced September 2018.

Comments: 18 pages, 13 figures, accepted for publication in Planetary and Space Science

Constraints on accretion disk size in the massive type 1 quasar PG 2308+098 from optical continuum reverberation lags

Authors: Mitsuru Kokubo

Abstract: Two years' worth of u-, g-, r-, i-, and z-band optical light curves were obtained for the massive type 1 quasar PG 2308+098 at z=0.433 using the 1.05-m Kiso Schmidt telescope/Kiso Wide Field Camera, and inter-band time lags of the light curves were measured. Wavelength-dependent continuum reverberation lag signals of several tens of days relative to the u-band were detected at g-, r-, i-, and z-ba… ▽ More Two years' worth of u-, g-, r-, i-, and z-band optical light curves were obtained for the massive type 1 quasar PG 2308+098 at z=0.433 using the 1.05-m Kiso Schmidt telescope/Kiso Wide Field Camera, and inter-band time lags of the light curves were measured. Wavelength-dependent continuum reverberation lag signals of several tens of days relative to the u-band were detected at g-, r-, i-, and z-bands, where the longer wavelength bands showed larger lags. From the wavelength-dependent lags, and assuming the standard disk temperature radial profile $T \propto R_{\rm disk}^{-3/4}$ and an X-ray/far-ultraviolet reprocessing picture, a constraint on the radius of the accretion disk responsible for the rest-frame 2500 Å disk continuum emission was derived as $R_{\rm disk} = 9.46^{+0.29}_{-3.12}$ light-days. The derived disk size is slightly (1.2-1.8 times) larger than the theoretical disk size of $R_{\rm disk} = 5.46$ light-days predicted from the black hole mass ($M_{\rm BH}$) and Eddington ratio estimates of PG 2308+098. This result is roughly in accordance with previous studies of lower mass active galactic nuclei (AGNs), where measured disk sizes have been found to be larger than the standard disk model predictions by a factor of $\sim 3$; however, the disk size discrepancy is more modest in PG 2308+098. By compiling literature values of the disk size constraints from continuum reverberation and gravitational microlensing observations for AGNs/quasars, we show that the $M_{\rm BH}$ dependence of $R_{\rm disk}$ is weaker than that expected from the standard disk model. These observations suggest that the standard Shakura-Sunyaev accretion disk theory has limitations in describing AGN/quasar accretion disks. △ Less

Submitted 30 July, 2018; originally announced July 2018.

Comments: 15 pages, 8 figures and 2 tables, accepted for publication in Publications of the Astronomical Society of Japan (PASJ) 2018 July 30

OISTER Optical and Near-Infrared Monitoring Observations of a Peculiar Radio-Loud Active Galactic Nucleus SDSS J110006.07+442144.3

Authors: Tomoki Morokuma, Masaomi Tanaka, Yasuyuki T. Tanaka, Ryosuke Itoh, Nozomu Tominaga, Poshak Gandhi, Elena Pian, Paolo Mazzali, Kouji Ohta, Emiko Matsumoto, Takumi Shibata, Hinako Akimoto, Hiroshi Akitaya, Gamal B. Ali, Tsutomu Aoki, Mamoru Doi, Nana Ebisuda, Ahmed Essam, Kenta Fujisawa, Hideo Fukushima, Shuhei Goda, Yuya Gouda, Hidekazu Hanayama, Yasuhito Hashiba, Osamu Hashimoto , et al. (58 additional authors not shown)

Abstract: We present monitoring campaign observations at optical and near-infrared (NIR) wavelengths for a radio-loud active galactic nucleus (AGN) at z=0.840, SDSS~J110006.07+442144.3 (hereafter, J1100+4421), which was identified during a flare phase in late February, 2014. The campaigns consist of three intensive observing runs from the discovery to March, 2015, mostly within the scheme of the OISTER coll… ▽ More We present monitoring campaign observations at optical and near-infrared (NIR) wavelengths for a radio-loud active galactic nucleus (AGN) at z=0.840, SDSS~J110006.07+442144.3 (hereafter, J1100+4421), which was identified during a flare phase in late February, 2014. The campaigns consist of three intensive observing runs from the discovery to March, 2015, mostly within the scheme of the OISTER collaboration. Optical-NIR light curves and simultaneous spectral energy distributions (SEDs) are obtained. Our measurements show the strongest brightening in March, 2015. We found that the optical-NIR SEDs of J1100+4421 show an almost steady shape despite the large and rapid intranight variability. This constant SED shape is confirmed to extend to $\sim5~μ$m in the observed frame using the archival WISE data. Given the lack of absorption lines and the steep power-law spectrum of $α_ν\sim-1.4$, where $f_ν\proptoν^{α_ν}$, synchrotron radiation by a relativistic jet with no or small contributions from the host galaxy and the accretion disk seems most plausible as an optical-NIR emission mechanism. The steep optical-NIR spectral shape and the large amplitude of variability are consistent with this object being a low $ν_{\rm{peak}}$ jet-dominated AGN. In addition, sub-arcsec resolution optical imaging data taken with Subaru Hyper Suprime-Cam does not show a clear extended component and the spatial scales are significantly smaller than the large extensions detected at radio wavelengths. The optical spectrum of a possible faint companion galaxy does not show any emission lines at the same redshift and hence a merging hypothesis for this AGN-related activity is not supported by our observations. △ Less

Submitted 17 July, 2017; originally announced July 2017.

Comments: 22 pages, 7 figures, 2 tables. Accepted for publication in PASJ

H$α$ intensity map of the repeating fast radio burst FRB 121102 host galaxy from Subaru/Kyoto 3DII AO-assisted optical integral-field spectroscopy

Authors: Mitsuru Kokubo, Kazuma Mitsuda, Hajime Sugai, Shinobu Ozaki, Yosuke Minowa, Takashi Hattori, Yutaka Hayano, Kazuya Matsubayashi, Atsushi Shimono, Shigeyuki Sako, Mamoru Doi

Abstract: We present the H$α$ intensity map of the host galaxy of the repeating fast radio burst FRB 121102 at a redshift of z=0.193 obtained with the AO-assisted Kyoto 3DII optical integral-field unit mounted on the 8.2-m Subaru Telescope. We detected a compact H$α$-emitting (i.e., star-forming) region in the galaxy, which has a much smaller angular size [$< 0".57$ (1.9 kpc) at full width at half maximum (… ▽ More We present the H$α$ intensity map of the host galaxy of the repeating fast radio burst FRB 121102 at a redshift of z=0.193 obtained with the AO-assisted Kyoto 3DII optical integral-field unit mounted on the 8.2-m Subaru Telescope. We detected a compact H$α$-emitting (i.e., star-forming) region in the galaxy, which has a much smaller angular size [$< 0".57$ (1.9 kpc) at full width at half maximum (FWHM)] than the extended stellar continuum emission region determined by the Gemini/GMOS z'-band image [$\simeq 1".4$ (4.6 kpc) at FWHM with ellipticity b/a=0.45]. The spatial offset between the centroid of the H$α$ emission region and the position of the radio bursts is $0".08 \pm 0".02$ ($0.26 \pm 0.07$ kpc), indicating that FRB 121102 is located within the star-forming region. This close spatial association of FRB 121102 with the star-forming region is consistent with expectations from young pulsar/magnetar models for FRB 121102, and it also suggests that the observed H$α$ emission region can make a major dispersion measure (DM) contribution to the host galaxy DM component of FRB 121102. Nevertheless, the largest possible value of the DM contribution from the H$α$ emission region inferred from our observations still requires a significant amount of ionized baryons in intergalactic medium (the so-called `missing' baryons) as the DM source of FRB 121102, and we obtain a 90\% confidence level lower limit on the cosmic baryon density in the intergalactic medium in the low-redshift universe as $Ω_{IGM} > 0.012$. △ Less

Submitted 21 June, 2017; v1 submitted 12 May, 2017; originally announced May 2017.

Comments: 11 pages, 5 figures, and 1 table, accepted for publication in the Astrophysical Journal 2017 June 20

Constraints on the optical polarization source in the luminous non-blazar quasar 3C 323.1 (PG 1545+210) from the photometric and polarimetric variability

Authors: Mitsuru Kokubo

Abstract: We examine the optical photometric and polarimetric variability of the luminous type 1 non-blazar quasar 3C 323.1 (PG 1545+210). Two optical spectro-polarimetric measurements taken during the periods 1996$-$98 and 2003 combined with a $V$-band imaging polarimetric measurement taken in 2002 reveal that (1) as noted in the literature, the polarization of 3C 323.1 is confined only to the continuum em… ▽ More We examine the optical photometric and polarimetric variability of the luminous type 1 non-blazar quasar 3C 323.1 (PG 1545+210). Two optical spectro-polarimetric measurements taken during the periods 1996$-$98 and 2003 combined with a $V$-band imaging polarimetric measurement taken in 2002 reveal that (1) as noted in the literature, the polarization of 3C 323.1 is confined only to the continuum emission, that is, the emission from the broad line region is unpolarized; (2) the polarized flux spectra show evidence of a time-variable broad absorption feature in the wavelength range of the Balmer continuum and other recombination lines; (3) weak variability in the polarization position angle ($PA$) of $\sim$ 4 deg over a time-scale of 4$-$6 years is observed; and (4) the V-band total flux and the polarized flux show highly correlated variability over a time-scale of one year. Taking the above-mentioned photometric and polarimetric variability properties and the results from previous studies into consideration, we propose a geometrical model for the polarization source in 3C 323.1, in which an equatorial absorbing region and an axi-asymmetric equatorial electron-scattering region are assumed to be located between the accretion disc and the broad line region. The scattering/absorbing regions can perhaps be attributed to the accretion disc wind or flared disc surface, but further polarimetric monitoring observations for 3C~323.1 and other quasars with continuum-confined polarization are needed to probe the true physical origins of these regions. △ Less

Submitted 13 January, 2017; originally announced January 2017.

Comments: 16 pages, 8 figures and 2 tables, accepted for publication in MNRAS 2017 January 10

X-ray and Optical Correlation of Type I Seyfert NGC 3516 Studied with Suzaku and Japanese Ground-Based Telescopes

Authors: Hirofumi Noda, Takeo Minezaki, Makoto Watanabe, Mitsuru Kokubo, Kenji Kawaguchi, Ryosuke Itoh, Kumiko Morihana, Yoshihiko Saito, Hikaru Nakao, Masataka Imai, Yuki Moritani, Katsutoshi Takaki, Miho Kawabata, Tatsuya Nakaoka, Makoto Uemura, Koji Kawabata, Michitoshi Yoshida, Akira Arai, Yuhei Takagi, Tomoki Morokuma, Mamoru Doi, Yoichi Itoh, Shin'ya Yamada, Kazuhiro Nakazawa, Yasushi Fukazawa , et al. (1 additional authors not shown)

Abstract: From 2013 April to 2014 April, we performed an X-ray and optical simultaneous monitoring of the type 1.5 Seyfert galaxy NGC 3516. It employed Suzaku, and 5 Japanese ground-based telescopes, the Pirka, Kiso Schmidt, Nayuta, MITSuME, and the Kanata telescopes. The Suzaku observations were conducted seven times with various intervals ranging from days, weeks, to months, with an exposure of $\sim50$ k… ▽ More From 2013 April to 2014 April, we performed an X-ray and optical simultaneous monitoring of the type 1.5 Seyfert galaxy NGC 3516. It employed Suzaku, and 5 Japanese ground-based telescopes, the Pirka, Kiso Schmidt, Nayuta, MITSuME, and the Kanata telescopes. The Suzaku observations were conducted seven times with various intervals ranging from days, weeks, to months, with an exposure of $\sim50$ ksec each. The optical $B$-band observations not only covered those of Suzaku almost simultaneously, but also followed the source as frequently as possible. As a result, NGC 3516 was found in its faint phase with the 2-10 keV flux of $0.21-2.70 \times 10^{-11}$ erg s$^{-1}$ cm$^{-2}$. The 2-45 keV X-ray spectra were composed of a dominant variable hard power-law continuum with a photon index of $\sim1.7$, and a non-relativistic reflection component with a prominent Fe-K$α$ emission line. Producing the $B$-band light curve by differential image photometry, we found that the $B$-band flux changed by $\sim2.7 \times 10^{-11}$ erg s$^{-1}$ cm$^{-2}$, which is comparable to the X-ray variation, and detected a significant flux correlation between the hard power-law component in X-rays and the $B$-band radiation, for the first time in NGC 3516. By examining their correlation, we found that the X-ray flux preceded that of $B$ band by $2.0^{+0.7}_{-0.6}$ days ($1σ$ error). Although this result supports the X-ray reprocessing model, the derived lag is too large to be explained by the standard view which assumes a "lamppost"-type X-ray illuminator located near a standard accretion disk. Our results are better explained by assuming a hot accretion flow and a truncated disk. △ Less

Submitted 25 May, 2016; originally announced May 2016.

Comments: 16 pages, 12 figures, accepted for publication in The Astrophysical Journal

The relationship between variable and polarized optical spectral components of luminous type 1 non-blazar quasars

Authors: Mitsuru Kokubo

Abstract: Optical spectropolarimetry carried out by Kishimoto et al. (2004) has shown that several luminous type 1 quasars show a strong decrease of the polarized continuum flux in the rest-frame near-UV wavelengths of $λ<4000$Å. In the literature, this spectral feature is interpreted as evidence of the broadened hydrogen Balmer absorption edge imprinted in the accretion disk thermal emission due to the dis… ▽ More Optical spectropolarimetry carried out by Kishimoto et al. (2004) has shown that several luminous type 1 quasars show a strong decrease of the polarized continuum flux in the rest-frame near-UV wavelengths of $λ<4000$Å. In the literature, this spectral feature is interpreted as evidence of the broadened hydrogen Balmer absorption edge imprinted in the accretion disk thermal emission due to the disk atmospheric opacity effect. On the other hand, the quasar flux variability studies have shown that the variable continuum component in UV-optical spectra of quasars, which is considered to be a good indicator of the intrinsic spectral shape of the accretion disk emission, generally have significantly flat spectral shape throughout the near-UV to optical spectral range. To examine whether the disk continuum spectral shapes revealed as the polarized flux and as the variable component spectra are consistent with each other, we carry out multi-band photometric monitoring observations for a sample of four polarization-decreasing quasars of Kishimoto et al. (2004) (4C09.72, 3C323.1, Ton 202, and B2 1208+32) to derive the variable component spectra and compare the spectral shape of them with that of the polarized flux spectra. Contrary to expectation, we confirm that the two spectral components of these quasars have totally different spectral shape in that the variable component spectra are significantly bluer compared to the polarized flux spectra. This discrepancy in the spectral shape may imply either (1) the decrease of polarization degree in the rest-frame UV wavelengths is not indicating the Balmer absorption edge feature but is induced by some unknown (de)polarization mechanisms, or (2) the UV-optical flux variability is occurring preferentially at the hot inner radii of the accretion disk and thus the variable component spectra do not reflect the whole accretion disk emission. △ Less

Submitted 15 April, 2016; originally announced April 2016.

Comments: 12 pages, 5 figures and 3 tables, accepted for publication in Publications of the Astronomical Society of Japan (PASJ) 2016 April 13

An Effective Selection Method for Low-Mass Active Black Holes and First Spectroscopic Identification

Authors: Tomoki Morokuma, Nozomu Tominaga, Masaomi Tanaka, Naoki Yasuda, Hisanori Furusawa, Yuki Taniguchi, Takahiro Kato, Ji-an Jiang, Tohru Nagao, Hanindyo Kuncarayakti, Kana Morokuma-Matsui, Hiroyuki Ikeda, Sergei Blinnikov, Ken'ichi Nomoto, Mitsuru Kokubo, Mamoru Doi

Abstract: We present a new method to effectively select objects which may be low-mass active black holes (BHs) at galaxy centers using high-cadence optical imaging data, and our first spectroscopic identification of an active 2.7x10^6 Msun BH at z=0.164. This active BH was originally selected due to its rapid optical variability, from a few hours to a day, based on Subaru Hyper Suprime-Cam~(HSC) g-band imag… ▽ More We present a new method to effectively select objects which may be low-mass active black holes (BHs) at galaxy centers using high-cadence optical imaging data, and our first spectroscopic identification of an active 2.7x10^6 Msun BH at z=0.164. This active BH was originally selected due to its rapid optical variability, from a few hours to a day, based on Subaru Hyper Suprime-Cam~(HSC) g-band imaging data taken with 1-hour cadence. Broad and narrow H-alpha and many other emission lines are detected in our optical spectra taken with Subaru FOCAS, and the BH mass is measured via the broad H-alpha emission line width (1,880 km s^{-1}) and luminosity (4.2x10^{40} erg s^{-1}) after careful correction for the atmospheric absorption around 7,580-7,720A. We measure the Eddington ratio to be as low as 0.05, considerably smaller than those in a previous SDSS sample with similar BH mass and redshift, which indicates one of the strong potentials of our Subaru survey. The g-r color and morphology of the extended component indicate that the host galaxy is a star-forming galaxy. We also show effectiveness of our variability selection for low-mass active BHs. △ Less

Submitted 7 March, 2016; originally announced March 2016.

Comments: 11 pages, 5 figures, 2 tables. accepted for publication in PASJ

Constraints on the temperature inhomogeneity in quasar accretion discs from the ultraviolet-optical spectral variability

Authors: Mitsuru Kokubo

Abstract: The physical mechanisms of the quasar ultraviolet (UV)-optical variability are not well understood despite the long history of observations. Recently, Dexter & Agol presented a model of quasar UV-optical variability, which assumes large local temperature fluctuations in the quasar accretion discs. This inhomogeneous accretion disc model is claimed to describe not only the single-band variability a… ▽ More The physical mechanisms of the quasar ultraviolet (UV)-optical variability are not well understood despite the long history of observations. Recently, Dexter & Agol presented a model of quasar UV-optical variability, which assumes large local temperature fluctuations in the quasar accretion discs. This inhomogeneous accretion disc model is claimed to describe not only the single-band variability amplitude, but also microlensing size constraints and the quasar composite spectral shape. In this work, we examine the validity of the inhomogeneous accretion disc model in the light of quasar UV-optical spectral variability by using five-band multi-epoch light curves for nearly 9 000 quasars in the Sloan Digital Sky Survey (SDSS) Stripe 82 region. By comparing the values of the intrinsic scatter $σ_{\text{int}}$ of the two-band magnitude-magnitude plots for the SDSS quasar light curves and for the simulated light curves, we show that Dexter & Agol's inhomogeneous accretion disc model cannot explain the tight inter-band correlation often observed in the SDSS quasar light curves. This result leads us to conclude that the local temperature fluctuations in the accretion discs are not the main driver of the several years' UV-optical variability of quasars, and consequently, that the assumption that the quasar accretion discs have large localized temperature fluctuations is not preferred from the viewpoint of the UV-optical spectral variability. △ Less

Submitted 27 February, 2015; originally announced March 2015.

Comments: 14 pages, 7 figures and 2 tables, accepted for publication in MNRAS 2015 February 4

Kiso Supernova Survey (KISS): Survey Strategy

Authors: Tomoki Morokuma, Nozomu Tominaga, Masaomi Tanaka, Kensho Mori, Emiko Matsumoto, Yuki Kikuchi, Takumi Shibata, Shigeyuki Sako, Tsutomu Aoki, Mamoru Doi, Naoto Kobayashi, Hiroyuki Maehara, Noriyuki Matsunaga, Hiroyuki Mito, Takashi Miyata, Yoshikazu Nakada, Takao Soyano, Ken'ichi Tarusawa, Satoshi Miyazaki, Fumiaki Nakata, Norio Okada, Yuki Sarugaku, Michael W. Richmond, Hiroshi Akitaya, Greg Aldering , et al. (36 additional authors not shown)

Abstract: The Kiso Supernova Survey (KISS) is a high-cadence optical wide-field supernova (SN) survey. The primary goal of the survey is to catch the very early light of a SN, during the shock breakout phase. Detection of SN shock breakouts combined with multi-band photometry obtained with other facilities would provide detailed physical information on the progenitor stars of SNe. The survey is performed us… ▽ More The Kiso Supernova Survey (KISS) is a high-cadence optical wide-field supernova (SN) survey. The primary goal of the survey is to catch the very early light of a SN, during the shock breakout phase. Detection of SN shock breakouts combined with multi-band photometry obtained with other facilities would provide detailed physical information on the progenitor stars of SNe. The survey is performed using a 2.2x2.2 deg field-of-view instrument on the 1.05-m Kiso Schmidt telescope, the Kiso Wide Field Camera (KWFC). We take a three-minute exposure in g-band once every hour in our survey, reaching magnitude g~20-21. About 100 nights of telescope time per year have been spent on the survey since April 2012. The number of the shock breakout detections is estimated to be of order of 1 during our 3-year project. This paper summarizes the KISS project including the KWFC observing setup, the survey strategy, the data reduction system, and CBET-reported SNe discovered so far by KISS. △ Less

Submitted 3 September, 2014; originally announced September 2014.

Comments: 13 pages, 10 figures, 3 tables. Publications of the Astronomical Society of Japan (PASJ) in press

Statistical Properties of Multi-epoch Spectral Variability of SDSS Stripe82 Quasars

Authors: Mitsuru Kokubo, Tomoki Morokuma, Takeo Minezaki, Mamoru Doi, Toshihiro Kawaguchi, Hiroaki Sameshima, Shintaro Koshida

Abstract: We investigate the UV-optical (longward of Ly$α$ 1216Å) spectral variability of nearly 9000 quasars ($0<z<4$) using multi-epoch photometric data within the SDSS Stripe 82 region. The regression slope in the flux-flux space of a quasar light curve directly measures the color of the flux difference spectrum, then the spectral shape of the flux difference spectra can be derived by taking a careful lo… ▽ More We investigate the UV-optical (longward of Ly$α$ 1216Å) spectral variability of nearly 9000 quasars ($0<z<4$) using multi-epoch photometric data within the SDSS Stripe 82 region. The regression slope in the flux-flux space of a quasar light curve directly measures the color of the flux difference spectrum, then the spectral shape of the flux difference spectra can be derived by taking a careful look at the redshift dependence of the regression slopes. First, we confirm that the observed quasar spectrum becomes bluer when the quasar becomes brighter. We infer the spectral index of the composite difference spectrum as $α_ν^{\text{dif}}\sim +1/3$ (in the form of $f_ν\propto ν^{α_ν}$), which is significantly bluer than that of the composite spectrum $α_ν^{\text{com}}\sim -0.5$. We also show that the continuum variability cannot be explained by the accretion disk models with varying mass accretion rate. Second, we examine the effects of broad emission line variability on the color-redshift space. The variability of the "Small Blue Bump" is extensively discussed. We show that the low-ionization lines of MgII and FeII are less variable compared to Balmer emission lines and high-ionization lines, and the Balmer continuum is the dominant variable source around $\sim 3000$Å. These results are compared with previous studies, and the physical mechanisms of the variability of the continuum and emission lines are discussed. △ Less

Submitted 12 February, 2014; v1 submitted 20 January, 2014; originally announced January 2014.

Comments: 20 pages, 17 figures and 2 tables, matches version published in ApJ with a typo corrected in caption of Figure 14 (α_ν---> α_λ)

Journal ref: The Astrophysical Journal 783 (2014) 46