Source author record

Mitsuru Kokubo

Mitsuru Kokubo appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2026arXiv

Diverse Origins of Broad H$α$ Lines in Heavily Obscured AGNs Revealed by Multi-epoch Spectroscopy

According to the classical AGN model, broad emission lines originate from the broad-line region (BLR) and are observable only when the attenuation by the dusty torus is small. However, we recently found several heavily-obscured ($A_V > 50$ mag) AGNs with broad H$α$ detections: MCG -3-34-64, UGC 5101, and Mrk 268. To investigate the origin of the observed broad line in these AGNs, we performed multi-epoch optical spectroscopic observations to search for flux variability of the broad H$α$ line. For MCG -3-34-64 and UGC 5101, no significant variability was detected, suggesting that the broad line of these AGNs may arise from sources other than the BLR. Spectral fitting analysis suggests possible large contribution of ionized outflows to the observed broad component of MCG -3-34-64, while both the outflow and scattering by polar material can explain that of UGC 5101. For Mrk 268, we detected a significant ($4.3σ$) flux variation of the broad H$α$ line by using the flux ratio of the H$α$ complex and the [SII]$λ\lambda6716$, 6731 doublet, indicating that the broad line originates directly from the BLR. The lack of significant flux variation in the optical continuum implies that the line of sight to the nucleus of Mrk 268 is mildly obscured. Our results demonstrate that the observed broad H$α$ lines in obscured AGNs likely have multiple origins. Such complexity may introduce additional uncertainties in black hole mass measurements of distant AGNs revealed by e.g., JWST.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2021arXiv

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

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.

preprint2021arXiv

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

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/.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2019arXiv

An optical search for transients lasting a few seconds

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.

preprint2019arXiv

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

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.

preprint2019arXiv

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

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.

preprint2016arXiv

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

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.

preprint2016arXiv

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

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.

preprint2016arXiv

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

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.

preprint2015arXiv

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

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.

preprint2014arXiv

Kiso Supernova Survey (KISS): Survey Strategy

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.

preprint2014arXiv

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

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.