Source author record

Masaomi Tanaka

Masaomi Tanaka 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

74works
11topics
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

74 published item(s)

preprint2022arXiv

A low-mass binary neutron star: long-term ejecta evolution and kilonovae with weak blue emission

We study the long-term evolution of ejecta formed in a binary neutron star (BNS) merger that results in a long-lived remnant NS by performing a hydrodynamics simulation with the outflow data of a numerical relativity simulation as the initial condition. At the homologously expanding phase, the total ejecta mass reaches $\approx0.1\,M_\odot$ with an average velocity of $\approx0.1\,c$ and lanthanide fraction of $\approx 0.005$. We further perform the radiative transfer simulation employing the obtained ejecta profile. We find that, contrary to a naive expectation from the large ejecta mass and low lanthanide fraction, the optical emission is not as bright as that in GW170817/AT2017gfo, while the infrared emission can be brighter. This light curve property is attributed to preferential diffusion of photons toward the equatorial direction due to the prolate ejecta morphology, large opacity contribution of Zr, Y, and lanthanides, and low specific heating rate of the ejecta. Our results suggest that these light curve features could be used as an indicator for the presence of a long-lived remnant NS. We also found that the bright optical emission broadly consistent with GW170817/AT2017gfo is realized for the case that the high-velocity ejecta components in the polar region are suppressed. These results suggest that the remnant in GW170817/AT2017gfo is unlikely to be a long-lived NS, but might have collapsed to a black hole within ${\cal O}(0.1)$ s.

preprint2022arXiv

Abundance stratification in type Ia supernovae -- VI: the peculiar slow decliner SN\,1999aa

The abundance distribution in the ejecta of the peculiar slowly declining Type Ia supernova (SN\,Ia) SN\,1999aa is obtained by modelling a time series of optical spectra. Similar to SN\,1991T, SN\,1999aa was characterised by early-time spectra dominated by \FeIII\ features and a weak \SiII\,6355\,Å line, but it exhibited a high-velocity \CaII\,H\&K line and morphed into a spectroscopically normal SN\,Ia earlier. Three explosion models are investigated, yielding comparable fits. The innermost layers are dominated by $\sim 0.3$\,\Msun\ of neutron-rich stable Fe-group elements, mostly stable iron. Above that central region lies a \Nifs-dominated shell, extending to $v \approx 11,000$ -- $12,000$\,\kms, with mass $\sim 0.65$\,\Msun. These inner layers are therefore similar to those of normal SNe\,Ia. However, the outer layers exhibit composition peculiarities similar to those of SN\,1991T: the intermediate-mass elements shell is very thin, containing only $\sim 0.2$\,\Msun, and is sharply separated from an outer oxygen-dominated shell, which includes $\sim 0.22$\,\Msun. These results imply that burning suddenly stopped in SN\,1999aa. This is a feature SN\,1999aa shares with SN\,1991T, and explain the peculiarities of both SNe, which are quite similar in nature apart from the different luminosities. The spectroscopic path from normal to SN\,1991T-like SNe\,Ia cannot be explained solely by a temperature sequence. It also involves composition layering differences, suggesting variations in the progenitor density structure or in the explosion parameters.

preprint2022arXiv

Constraints on Explosion Timescale of Core-Collapse Supernovae Based on Systematic Analysis of Light Curves

Explosion mechanism of core-collapse supernovae is not fully understood yet. In this work, we give constraints on the explosion timescale based on $^{56}$Ni synthesized by supernova explosions. First, we systematically analyze multi-band light curves of 82 stripped-envelope supernovae (SESNe) to obtain bolometric light curves, which is among the largest samples of the bolometric light curves of SESNe derived from the multi-band spectral energy distribution. We measure the decline timescale and the peak luminosity of the light curves and estimate the ejecta mass ($M_{\rm ej}$) and $^{56}$Ni mass ($M_{\rm Ni}$) to connect the observed properties with the explosion physics. We then carry out one-dimensional hydrodynamics and nucleosynthesis calculations, varying the progenitor mass and the explosion timescale. From the calculations, we show that the maximum $^{56}$Ni mass that $^{56}$Ni-powered SNe can reach is expressed as $M_{\rm Ni} \lesssim 0.2 \ M_{\rm ej}$. Comparing the results from the observations and the calculations, we show that the explosion timescale shorter than 0.3 sec explains the synthesized $^{56}$Ni mass of the majority of the SESNe.

preprint2022arXiv

Deep simultaneous limits on optical emission from FRB 20190520B by 24.4 fps observations with Tomo-e Gozen

We conduct 24.4~fps optical observations of repeating Fast Radio Burst (FRB) 20190520B using Tomo-e Gozen, a high-speed CMOS camera mounted on the Kiso 105-cm Schmidt telescope, simultaneously with radio observations carried out using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). We succeeded in the simultaneous optical observations of 11 radio bursts that FAST detected. However, no corresponding optical emission was found. The optical fluence limits as deep as 0.068 Jy ms are obtained for the individual bursts (0.029 Jy ms on the stacked data) corrected for the dust extinction in the Milky Way. The fluence limit is deeper than those obtained in the previous simultaneous observations for an optical emission with a duration $\gtrsim 0.1$ ms. Although the current limits on radio--optical spectral energy distribution (SED) of FRBs are not constraining, we show that SED models based on observed SEDs of radio variable objects such as optically detected pulsars, and a part of parameter spaces of theoretical models in which FRB optical emission is produced by inverse-Compton scattering in a pulsar magnetosphere or a strike of a magnetar blastwave into a hot wind bubble, can be ruled out once a similar fluence limit as in our observation is obtained for a bright FRB with a radio fluence $\gtrsim 5$ Jy ms.

preprint2022arXiv

Deep-learning Real/Bogus classification for the Tomo-e Gozen transient survey

We present a deep neural network Real/Bogus classifier that improves classification performance in the Tomo-e Gozen transient survey by handling label errors in the training data. In the wide-field, high-frequency transient survey with Tomo-e Gozen, the performance of conventional convolutional neural network classifier is not sufficient as about $10^6$ bogus detections appear every night. In need of a better classifier, we have developed a new two-stage training method. In this training method, label errors in the training data are first detected by normal supervised learning classification, and then they are unlabeled and used for training of semi-supervised learning. For actual observed data, the classifier with this method achieves an area under the curve (AUC) of 0.9998 and a false positive rate (FPR) of 0.0002 at true positive rate (TPR) of 0.9. This training method saves relabeling effort by humans and works better on training data with a high fraction of label errors. By implementing the developed classifier in the Tomo-e Gozen pipeline, the number of transient candidates was reduced to $\sim$40 objects per night, which is $\sim$1/130 of the previous version, while maintaining the recovery rate of real transients. This enables more efficient selection of targets for follow-up observations.

preprint2022arXiv

Lanthanide Features in Near-infrared Spectra of Kilonovae

The observations of GW170817/AT2017gfo have provided us with evidence that binary neutron star mergers are sites of $r$-process nucleosynthesis. However, the observed signatures in the spectra of GW170817/AT2017gfo have not been fully decoded especially in the near-infrared (NIR) wavelengths. In this paper, we investigate the kilonova spectra over the entire wavelength range with the aim of elemental identification. We systematically calculate the strength of bound-bound transitions by constructing a hybrid line list that is accurate for important strong transitions and complete for weak transitions. We find that the elements on the left side of the periodic table, such as Ca, Sr, Y, Zr, Ba, La, and Ce, tend to produce prominent absorption lines in the spectra. This is because such elements have a small number of valence electrons and low-lying energy levels, resulting in strong transitions. By performing self-consistent radiative transfer simulations for the entire ejecta, we find that La III and Ce III appear in the NIR spectra, which can explain the absorption features at $λ\sim 12000$-14000 A in the spectra of GW170817/AT2017gfo. The mass fractions of La and Ce are estimated to be $>2\times 10^{-6}$ and $\sim$ (1-100)$\times 10^{-5}$, respectively. An actinide element Th can also be a source of absorption as the atomic structure is analogous to that of Ce. However, we show that Th III features are less prominent in the spectra because of the denser energy levels of actinides compared to those of lanthanides.

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.

preprint2022arXiv

Opacity of the highly ionized lanthanides and the effect on the early kilonova

We investigate the effect of the presence of lanthanides (Z = 57- 71) on the kilonova at t~hours after the neutron star merger for the first time. For this purpose, we calculate the atomic structures and the opacities for selected lanthanides: Nd (Z = 60), Sm (Z = 62), and Eu (Z = 63). We consider the ionization degree up to tenth (XI), applicable for the ejecta at t ~ a few hours after the merger, when the temperature is T ~ 10^5 K. We find that the opacities for the highly ionized lanthanides are exceptionally high, reaching k_exp~1000 cm^2/g for Eu, due to the highly dense energy levels. Using the new opacity, we perform radiative transfer simulations to show that the early light curves become fainter by a (maximum) factor of four, in comparison to lanthanide-free ejecta at t~0.1 day. However, the period at which the light curves are affected is relatively brief due to the rapid time evolution of the opacity in the outermost layer of the ejecta. We predict that for a source at a distance of ~100 Mpc, UV brightness for lanthanide-rich ejecta shows a drop to ~21-22 mag at t~0.1 day and the UV peaks around t~0.2 day with a magnitude of ~19 mag. Future detection of such a kilonova by the existing UV satellite like Swift or the upcoming UV satellite ULTRASAT will provide useful constraints on the abundance in the outer ejecta and the corresponding nucleosynthesis conditions in the neutron star mergers.

preprint2022arXiv

Statistical properties of the nebular spectra of 103 stripped envelope core collapse supernovae

We present an analysis of the nebular spectra of 103 stripped envelope (SE) supernovae (SNe) collected from the literature and observed with the Subaru Telescope from 2002 to 2012, focusing on [O I] 6300, 6363. The line profile and width of [O I] are employed to infer the ejecta geometry and the expansion velocity of the inner core. These two measurements are then compared with the SN sub types, and further with the [O I]/[Ca II] ratio, which is used as an indicator of the progenitor CO core mass. Based on the best fit results of the [O I] profile, the objects are classified into different morphological groups, and we conclude that the deviation from spherical symmetry is a common feature for all types of SESNe. There is a hint (at about 1 sigma level) that the distributions of the line profile fractions are different between canonical SESNe and broad-line SNe Ic. A correlation between [O I] width and [O I]/[Ca II] ratio is discerned, indicating that the oxygen-rich material tends to expand faster for objects with a more massive CO core. Such a correlation can be utilized to constrain the relation between the progenitor mass and the kinetic energy of the explosion. Further, when [O I]/[Ca II] ratio increases, the fraction of objects with Gaussian [O I] profile increases, while those with double-peaked profile decreases. This phenomenon connects ejecta geometry and the progenitor CO core mass.

preprint2022arXiv

Tungsten vs Selenium as a potential source of kilonova nebular emission observed by Spitzer

Infrared emission lines arising from transitions between fine structure levels of heavy elements are expected to produce kilonova nebular emission. For the kilonova in GW170817, strong emission at 4.5 ${\rm μm}$ at late times was detected by the Spitzer Space Telescope but no source was detected at 3.6 ${\rm μm}$. This peculiar spectrum indicates that there exist strong line emitters around 4.5 ${\rm μm}$ and the absence of strong lines around 3.6 ${\rm μm}$. To model the spectrum we prepare a line list based on the selection rules in LS coupling from the experimentally calibrated energy levels in the NIST database. This method enables to generate the synthetic spectra with accurate line wavelengths. We find that the spectrum is sensitive to the abundance pattern whether or not the first r-process peak elements are included. In both cases, the synthetic spectra can match the observed data, leading to two possible interpretations. If the first peak elements are abundant a Se III line dominates the flux. If otherwise, W III with Os III, Rh III, and Ce IV can be the main sources. Observing nebular spectra for the future kilonovae in a wider wavelength range can provide more conclusive elemental identification.

preprint2022arXiv

Video Observations of Tiny Near-Earth Objects with Tomo-e Gozen

We report the results of video observations of tiny (diameter less than 100 m) near-Earth objects (NEOs) with Tomo-e Gozen on the Kiso 105 cm Schmidt telescope. A rotational period of a tiny asteroid reflects its dynamical history and physical properties since smaller objects are sensitive to the YORP effect. We carried out video observations of 60 tiny NEOs at 2 fps from 2018 to 2021 and successfully derived the rotational periods and axial ratios of 32 NEOs including 13 fast rotators with rotational periods less than 60 s. The fastest rotator found during our survey is 2020 HS7 with a rotational period of 2.99 s. We statistically confirmed that there is a certain number of tiny fast rotators in the NEO population, which have been missed with any previous surveys. We have discovered that the distribution of the tiny NEOs in a diameter and rotational period (D-P) diagram is truncated around a period of 10 s. The truncation with a flat-top shape is not explained well either by a realistic tensile strength of NEOs or suppression of YORP by meteoroid impacts. We propose that the dependence of the tangential YORP effect on the rotational period potentially explains the observed pattern in the D-P diagram.

preprint2021arXiv

Optical follow-up observation for GW event S190510g using Subaru/Hyper Suprime-Cam

A gravitational wave event, S190510g, which was classified as a binary-neutron-star coalescence at the time of preliminary alert, was detected by LIGO/Virgo collaboration on May 10, 2019. At 1.7 hours after the issue of its preliminary alert, we started a target-of-opportunity imaging observation in Y-band to search for its optical counterpart using the Hyper Suprime-Cam (HSC) on the Subaru Telescope. The observation covers a 118.8 deg$^2$ sky area corresponding to 11.6% confidence in the localization skymap released in the preliminary alert and 1.2% in the updated skymap. We divided the observed area into two fields based on the availability of HSC reference images. For the fields with the HSC reference images, we applied an image subtraction technique; for the fields without the HSC reference images, we sought individual HSC images by matching a catalog of observed objects with the PS1 catalog. The search depth is 22.28 mag in the former method and the limit of search depth is 21.3 mag in the latter method. Subsequently, we performed visual inspection and obtained 83 candidates using the former method and 50 candidates using the latter method. Since we have only the 1-day photometric data, we evaluated probability to be located inside the 3D skymap by estimating their distances with photometry of associated extended objects. We found three candidates are likely located inside the 3D skymap and concluded they could be an counterpart of S190510g, while most of 133 candidates were likely to be supernovae because the number density of candidates was consistent with the expected number of supernova detections. By comparing our observational depth with a light curve model of such a kilonova reproducing AT2017gfo, we show that early-deep observations with the Subaru/HSC can capture the rising phase of blue component of kilonova at the estimated distance of S190510g (~230 Mpc).

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.

preprint2020arXiv

$g$-matrix folding-model approach to reaction cross sections for scattering of Ca isotopes on a C target

We first predict the ground-state properties of Ca isotopes, using the Gogny-D1S Hartree-Fock-Bogoliubov (GHFB) with and without the angular momentum projection (AMP). We find that $^{64}$Ca is an even-dripline nucleus and $^{59}$Ca is an odd-dripline nucleus, using $A$ dependence of the one-neutron separation energy $S_{1}$ and the two-neutron separation energy, $S_{2}$. As for $S_{1}$, $S_{2}$ and the binding energies $E_{\rm B}$, our results agree with the experimental data in $^{40-58}$Ca. As other ground-state properties of $^{40-60,62,64}$Ca, we predict charge, proton, neutron, matter radii, neutron skin and deformation. As for charge radii, our results are consistent with the experimental data in $^{40-52}$Ca. For $^{48}$Ca, our results on proton, neutron, matter radii agree with the experimental data. Very lately, Tanaka et. al. measured interaction cross sections for $^{42-51}$Ca scattering on a $^{12}$C target at an incident energy per nucleon of $E_{\rm lab}=280$MeV. Secondly, we predict reaction cross sections $σ_{\rm R}$ for $^{40-60,62,64}$Ca, using a chiral $g$-matrix double-folding model (DFM). To show the reliability of the present DFM for $σ_{\rm R}$, we apply the DFM for the data on $^{12}$C scattering on $^{9}$Be, $^{12}$C, $^{27}$Al targets in $30 < E_{\rm lab} < 400 $MeV, and show that the present DFM is good in $30 < E_{\rm lab} < 100 $MeV and $250 < E_{\rm lab} < 400 $MeV. For $110 < E_{\rm lab} < 240 $MeV, our results have small errors. To improve the present DFM for $σ_{\rm R}$, we propose two prescriptions.

preprint2020arXiv

Constraint on the ejecta mass for a black hole-neutron star merger event candidate S190814bv

We derive the upper limit to the ejecta mass of S190814bv, a black hole-neutron star merger candidate, through the radiative transfer simulations for kilonovae with the realistic ejecta density profile as well as the detailed opacity and heating rate models. The limits to the ejecta mass strongly depend on the viewing angle. For the face-on observations ($\le45^\circ$), the total ejecta mass should be smaller than $0.1\,M_\odot$ for the average distance of S190814bv ($D=267$ Mpc), while larger mass is allowed for the edge-on observations. We also derive the conservative upper limits of the dynamical ejecta mass to be $0.02\,M_\odot$, $0.03\,M_\odot$, and $0.05\,M_\odot$ for the viewing angle $\le 20^\circ$, $\le 50^\circ$, and for $\le 90^\circ$, respectively. We show that the {\it iz}-band observation deeper than $22$ mag within 2 d after the GW trigger is crucial to detect the kilonova with the total ejecta mass of $0.06\,M_\odot$ at the distance of $D=300$ Mpc. We also show that a strong constraint on the NS mass-radius relation can be obtained if the future observations put the upper limit of $0.03\,M_\odot$ to the dynamical ejecta mass for a BH-NS event with the chirp mass smaller than $\lesssim 3\,M_\odot$ and effective spin larger than $\gtrsim 0.5$.

preprint2020arXiv

Energy level structure and transition data of Er$^{2+}$

The main aim of this paper is to present accurate energy levels of the ground [Xe]$4f^{12}$ and first excited [Xe]$4f^{11}5d$ configurations of Er$^{2+}$. The energy level structure of the Er$^{2+}$ ion was computed using the multiconfiguration Dirac-Hartree-Fock and relativistic configuration interaction (RCI) methods, as implemented in the GRASP2018 program package. The Breit interaction, self-energy and vacuum polarization corrections were included in the RCI computations. The zero-first-order approach was used in the computations. Energy levels with the identification in $LS$ coupling for all (399) states belonging to the [Xe]$4f^{12}$ and [Xe]$4f^{11}5d$ configurations are presented. Electric dipole (E1) transition data between the levels of these two configurations are computed. The accuracy of the these data are evaluated by studying the behaviour of the transition rates as functions of the gauge parameter as well as by evaluating the cancellation factors. The core electron correlations were studied using different strategies. Root-mean-square deviations obtained in this study for states of the ground and excited configurations from the available experimental or semi-empirical data are 649 cm$^{-1}$, and 747 cm$^{-1}$, respectively.

preprint2020arXiv

Extended calculations of energy levels and transition rates for singly ionized lanthanide elements I: Pr - Gd

Lanthanide elements play important roles as an opacity source in the ejected material from neutron star mergers. Accurate and complete atomic data are necessary to evaluate the opacities and to analyze the observed data. In this paper, we perform extended, ab-initio atomic calculations from Pr II (Z=59) to Gd II (Z=64). By using multiconfiguration Dirac-Hartree-Fock and relativistic configuration-interaction methods, implemented in the general-purpose relativistic atomic structure package GRASP2K, we calculate the energy levels and transition data of electric dipole transitions. These computations are based on strategies (with small variations) of Nd II published by Gaigalas et al. (2019). Accuracy of data is evaluated by comparing computed energy levels with the NIST database or other works. For the energy levels, we obtain the average relative accuracy of 8%, 12%, 6%, 8%, and 7% for Pr II, Pm II, Sm II, Eu II, and Gd II ions, respectively as compared with the NIST data. Accuracy of energy transfer to the wavelength as 3%, 14% and 11% for Pr II, Eu II and Gd II. Our computed E1 type transition probabilities are in good agreement with experimental values presented by other authors especially for strong transitions.

preprint2020arXiv

Follow-up Observations for IceCube-170922A: Detection of Rapid Near-Infrared Variability and Intensive Monitoring of TXS 0506+056

We present our follow-up observations to search for an electromagnetic counterpart of the IceCube high-energy neutrino, IceCube-170922A. Monitoring observations of a likely counterpart, TXS 0506+056, are also described. First, we quickly took optical and near-infrared images of 7 flat-spectrum radio sources within the IceCube error region right after the neutrino detection and found a rapid flux decline of TXS 0506+056 in Kanata/HONIR J-band data. Motivated by this discovery, intensive follow-up observations of TXS 0506+056 are continuously done, including our monitoring imaging observations, spectroscopic observations, and polarimetric observations in optical and near-infrared wavelengths. TXS 0506+056 shows a large amplitude (~1.0 mag) variability in a time scale of several days or longer, although no significant variability is detected in a time scale of a day or shorter. TXS 0506+056 also shows a bluer-when-brighter trend in optical and near-infrared wavelengths. Structure functions of variabilities are examined and indicate that TXS 0506+056 is not a special blazar in terms of optical variability. Polarization measurement results of TXS 0506+056 are also discussed.

preprint2020arXiv

Late-Phase Spectropolarimetric Observations of Superluminous Supernova SN 2017egm to Probe the Geometry of the Inner Ejecta

We present our spectropolarimetric observations of SN 2017egm, a Type I superluminous supernova (SLSN-I) in a nearby galaxy NGC 3191, with Subaru telescope at 185.0 days after the g-band maximum light. This is the first spectropolarimetric observation for SLSNe at late phases. We find that the degree of the polarization in the late phase significantly changes from that measured at the earlier phase. The spectrum at the late phase shows a strong Ca emission line and therefore we reliably estimate the interstellar polarization component assuming that the emission line is intrinsically unpolarized. By subtracting the estimated interstellar polarization, we find that the intrinsic polarization at the early phase is only ~ 0.2 %, which indicates an almost spherical photosphere, with an axial ratio ~ 1.05. The intrinsic polarization at the late phase increases to ~ 0.8 %, which corresponds to the photosphere with an axial ratio 1.2. A nearly constant position angle of the polarization suggests the inner ejecta are almost axisymmetric. By these observations, we conclude that the inner ejecta are more aspherical than the outer ejecta. It may suggest the presence of a central energy source producing aspherical inner ejecta.

preprint2020arXiv

On the possibility of GW190425 being a black hole--neutron star binary merger

We argue that the kilonova/macronova associated with the gravitational-wave event GW190425 could have been bright enough to be detected if it was caused by the merger of a low-mass black hole and a neutron star. Although tidal disruption occurs for such a low-mass black hole as is generally expected, the masses of the dynamical ejecta are limited to <~10^{-3}M_sun, which is consistent with previous work in the literature. The remnant disk could be as massive as 0.05--0.1M_sun, and the disk outflow of ~0.01--0.03M_sun is likely to be driven by viscous or magnetohydrodynamic effects. The disk outflow may not be neutron-rich enough to synthesize an abundance of lanthanide elements, even in the absence of strong neutrino emitter, if the ejection is driven on the viscous time scale of >~0.3s. If this is the case, the opacity of the disk outflow is kept moderate, and a kilonova/macronova at the distance of GW190425 reaches a detectable brightness of 20--21mag at 1day after merger for most viewing angles. If some disk activity ejects the mass within ~0.1s, instead, lanthanide-rich outflows would be launched and the detection of emission becomes challenging. Future possible detections of kilonovae/macronovae from GW190425-like systems will disfavor the prompt collapse of binary neutron stars and a non-disruptive low-mass black hole--neutron star binary associated with a small neutron-star radius, whose mass ejection is negligible. The host-galaxy distance will constrain the viewing angle and deliver further information about the mass ejection.

preprint2020arXiv

Photometric classification of HSC transients using machine learning

The advancement of technology has resulted in a rapid increase in supernova (SN) discoveries. The Subaru/Hyper Suprime-Cam (HSC) transient survey, conducted from fall 2016 through spring 2017, yielded 1824 SN candidates. This gave rise to the need for fast type classification for spectroscopic follow-up and prompted us to develop a machine learning algorithm using a deep neural network (DNN) with highway layers. This machine is trained by actual observed cadence and filter combinations such that we can directly input the observed data array into the machine without any interpretation. We tested our model with a dataset from the LSST classification challenge (Deep Drilling Field). Our classifier scores an area under the curve (AUC) of 0.996 for binary classification (SN Ia or non-SN Ia) and 95.3% accuracy for three-class classification (SN Ia, SN Ibc, or SN II). Application of our binary classification to HSC transient data yields an AUC score of 0.925. With two weeks of HSC data since the first detection, this classifier achieves 78.1% accuracy for binary classification, and the accuracy increases to 84.2% with the full dataset. This paper discusses the potential use of machine learning for SN type classification purposes.

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

Systematic Opacity Calculations for Kilonovae

Coalescence of neutron stars gives rise to kilonova, thermal emission powered by radioactive decays of freshly synthesized r-process nuclei. Although observational properties are largely affected by bound-bound opacities of r-process elements, available atomic data have been limited. In this paper, we study element-to-element variation of the opacities in the ejecta of neutron star mergers by performing systematic atomic structure calculations of r-process elements for the first time. We show that the distributions of energy levels tend to be higher as electron occupation increases for each electron shell due to the larger energy spacing caused by larger effects of spin-orbit and electron-electron interactions. As a result, elements with a fewer number of electrons in the outermost shells tend to give larger contributions to the bound-bound opacities. This implies that Fe is not representative for the opacities of light r-process elements. The average opacities for the mixture of r-process elements are found to be kappa ~ 20-30 cm^2 g^{-1} for the electron fraction of Ye < 0.20, kappa ~ 3-5 cm^2 g^{-1} for Ye = 0.25-0.35, and kappa ~ 1 cm^2 g^{-1} for Ye = 0.40 at T = 5,000-10,000 K, and they steeply decrease at lower temperature. We show that, even with the same abundance or Ye, the opacity in the ejecta changes with time by one order of magnitude from 1 to 10 days after the merger. Our radiative transfer simulations with the new opacity data confirm that ejecta with a high electron fraction (Ye >~ 0.25, with no lanthanide) are needed to explain the early, blue emission in GW170817/AT2017gfo while lanthanide-rich ejecta (with a mass fraction of lanthanides ~ 5 x 10^{-3}) reproduce the long-lasting near-infrared emission.

preprint2020arXiv

The HSC-SSP Transient Survey: Implications from Early Photometry and Rise Time of Normal Type Ia Supernovae

With a booming number of Type Ia supernovae (SNe Ia) discovered within a few days of their explosions, a fraction of SNe Ia that show luminosity excess in the early phase (early-excess SNe Ia) have been confirmed. In this article, we report early-phase observations of seven photometrically normal SNe Ia (six early detections and one deep non-detection limit) at the COSMOS field through a half-year transient survey as a part of the Hyper Suprime-Cam Subaru Strategic Program (HSC SSP). In particular, a blue light-curve excess was discovered for HSC17bmhk, a normal SN Ia with rise time longer than 18.8 days, during the first four days after the discovery. The blue early excess in optical wavelength can be explained not only by interactions with a non-degenerate companion or surrounding dense circumstellar matter but also radiation powered by radioactive decays of $^{56}$Ni at the surface of the SN ejecta. Given the growing evidence of the early-excess discoveries in normal SNe Ia that have longer rise times than the average and a similarity in the nature of the blue excess to a luminous SN Ia subclass, we infer that early excess discovered in HSC17bmhk and other normal SNe Ia are most likely attributed to radioactive $^{56}$Ni decay at the surface of the SN ejecta. In order to successfully identify normal SNe Ia with early excess similar to that of HSC17bmhk, early UV photometries or high-cadence blue-band surveys are necessary.

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

Diversity of kilonova light curves

We perform radiative transfer simulations for kilonova in various situations, including the cases of prompt collapse to a black hole from neutron-star mergers, high-velocity ejecta possibly accelerated by magnetars, and a black hole-neutron star merger. Our calculations are done employing ejecta profiles predicted by numerical-relativity simulations and a new line list for all the r-process elements. We found that (i) the optical emission for binary neutron stars promptly collapsing to a black hole would be fainter by $\gtrsim1$-$2\,{\rm mag}$ than that found in GW170817, while the infrared emission could be as bright as that in GW170817 if the post-merger ejecta is as massive as $\approx0.01\,M_\odot$; (ii) the kilonova would be brighter than that observed in GW170817 for the case that the ejecta is highly accelerated by the electromagnetic energy injection from the remnant, but it would decline rapidly and the magnitude would become fainter than in GW170817 within a few days; (iii) the optical emission from a black hole-neutron star merger ejecta could be as bright as that observed in GW170817 for the case that sufficiently large amount of matter is ejected ($\gtrsim0.02\,M_\odot$), while the infrared brightness would be brighter by $1$-$2\,{\rm mag}$ at the same time. We show that the difference in the ejecta properties would be imprinted in the differences in the peak brightness and time of peak. This indicates that we may be able to infer the type of the central engine for kilonovae by observation of the peak in the multiple band.

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

Constraining the Age and Distance of the Galactic Supernova Remnant G156.2+5.7 by H-alpha Expansion Measurements

We present deep H-alpha images of portions of the X-ray bright but optically faint Galactic supernova remnant G156.2+5.7, revealing numerous and delicately thin nonradiative filaments which mark the location of the remnant's forward shock. These new images show that these filaments have a complex structure not visible on previous lower resolution optical images. By comparing H-alpha images taken in 2004 at the McDonald Observatory and in 2015-2016 at the Kiso Observatory, we set a stringent 1-sigma upper limit of expansion to be 0.06 arcsec/yr. This proper motion, combined with a shock speed of 500 km/s inferred from X-ray spectral analyses, gives a distance of > 1.7 kpc. In addition, a simple comparison of expansion indices of several SNRs allows us to infer the age of the remnant to be a few 10,000 yr old. These estimates are more straightforward and reliable than any other previous studies, and clearly rule out a possibility that G156.2+5.7 is physically associated with part of the Taurus-Auriga cloud and dust complex at a distance of 200-300 pc.

preprint2016arXiv

J-GEM Follow-Up Observations of The Gravitational Wave Source GW151226

We report the results of optical--infrared follow-up observations of the gravitational wave (GW) event GW151226 detected by the Advanced LIGO in the framework of J-GEM (Japanese collaboration for Gravitational wave ElectroMagnetic follow-up). We performed wide-field optical imaging surveys with Kiso Wide Field Camera (KWFC), Hyper Suprime-Cam (HSC), and MOA-cam3. The KWFC survey started at 2.26 days after the GW event and covered 778 deg$^2$ centered at the high Galactic region of the skymap of GW151226. We started the HSC follow-up observations from 12 days after the event and covered an area of 63.5 deg$^2$ of the highest probability region of the northern sky with the limiting magnitudes of 24.6 and 23.8 for i band and z band, respectively. MOA-cam3 covered 145 deg$^2$ of the skymap with MOA-red filter 2.5 months after the GW alert. Total area covered by the wide-field surveys was 986.5 deg$^2$. The integrated detection probability of all the observed area was $\sim$29%. We also performed galaxy-targeted observations with six optical and near-infrared telescopes from 1.61 days after the event. Total of 238 nearby (<100 Mpc) galaxies were observed with the typical I band limiting magnitude of $\sim$19.5. We detected 13 supernova candidates with the KWFC survey, and 60 extragalactic transients with the HSC survey. Two third of the HSC transients were likely supernovae and the remaining one third were possible active galactic nuclei. With our observational campaign, we found no transients that are likely to be associated with GW151226.

preprint2016arXiv

J-GEM Follow-Up Observations to Search for an Optical Counterpart of The First Gravitational Wave Source GW150914

We present our optical follow-up observations to search for an electromagnetic counterpart of the first gravitational wave source GW150914 in the framework of the Japanese collaboration for Gravitational wave ElectroMagnetic follow-up (J-GEM), which is an observing group utilizing optical and radio telescopes in Japan, as well as those in New Zealand, China, South Africa, Chile, and Hawaii. We carried out a wide-field imaging survey with Kiso Wide Field Camera (KWFC) on the 1.05-m Kiso Schmidt telescope in Japan and a galaxy-targeted survey with Tripole5 on the B&C 61-cm telescope in New Zealand. Approximately 24 deg2 regions in total were surveyed in i-band with KWFC and 18 nearby galaxies were observed with Tripole5 in g-, r-, and i-bands 4-12 days after the gravitational wave detection. Median 5-sigma depths are i~18.9 mag for the KWFC data and g~18.9 mag, r~18.7 mag, and i~18.3 mag for the Tripole5 data. Probability for a counterpart to be in the observed area is 1.2% in the initial skymap and 0.1% in the final skymap. We do not find any transient source associated to an external galaxy with spatial offset from its center, which is consistent with the local supernova rate. We summarize future prospects and ongoing efforts to pin down electromagnetic counterparts of binary black hole mergers as well as neutron star mergers.

preprint2016arXiv

Kilonova/Macronova Emission from Compact Binary Mergers

We review current understanding of kilonova/macronova emission from compact binary mergers (mergers of two neutron stars or a neutron star and a black hole). Kilonova/macronova is optical and near-infrared emission powered by radioactive decays of r-process nuclei. Emission from the dynamical ejecta with ~0.01 Msun is likely to have a luminosity of ~10^{40}-10^{41} erg s^{-1} with a characteristic timescale of about 1 week. The spectral peak is located in red optical or near-infrared wavelengths. A subsequent accretion disk wind may provide an additional luminosity, or an earlier/bluer emission if it is not absorbed by the precedent dynamical ejecta. The detection of near-infrared excess in the afterglow of short GRB 130603B and possible optical excess in GRB 060614 supports the concept of the kilonova/macronova scenario. At 200 Mpc distance, a typical brightness of kilonova/macronova with 0.01 Msun ejecta is expected to be about 22 mag and the emission rapidly fades to >24 mag within ~10 days after the merger. Kilonova/macronova candidates can be distinguished from supernovae by (1) the faster time evolution, (2) fainter absolute magnitudes, and (3) redder colors. To effectively search for such objects, follow-up survey observations with multiple visits within <10 days and with multiple filters will be important. Since the high expansion velocity (v ~ 0.1-0.2c) is a robust outcome of compact binary mergers, the detection of smooth spectra will be the smoking gun to conclusively identify the GW source.

preprint2016arXiv

Machine-learning Selection of Optical Transients in Subaru/Hyper Suprime-Cam Survey

We present an application of machine-learning (ML) techniques to source selection in the optical transient survey data with Hyper Suprime-Cam (HSC) on the Subaru telescope. Our goal is to select real transient events accurately and in a timely manner out of a large number of false candidates, obtained with the standard difference-imaging method. We have developed the transient selector which is based on majority voting of three ML machines of AUC Boosting, Random Forest, and Deep Neural Network. We applied it to our observing runs of Subaru-HSC in 2015 May and August, and proved it to be efficient in selecting optical transients. The false positive rate was 1.0% at the true positive rate of 90% in the magnitude range of 22.0--25.0 mag for the former data. For the latter run, we successfully detected and reported ten candidates of supernovae within the same day as the observation. From these runs, we learned the following lessons: (1) the training using artificial objects is effective in filtering out false candidates, especially for faint objects, and (2) combination of ML by majority voting is advantageous.

preprint2016arXiv

Models of Kilonova/macronova emission from black hole-neutron star mergers

Black hole-neutron star mergers are among the promising gravitational-wave sources for ground-based detectors, and gravitational waves from black hole-neutron star mergers are expected to be detected in the next few years. Simultaneous detection of electromagnetic counterparts with gravitational-wave detection provides rich information about the merger events. Among the possible electromagnetic counterparts from the black hole-neutron star merger, the emission powered by the decay of radioactive r-process nuclei, so called kilonova/macronova, is one of the best targets for follow-up observation. We derive fitting formulas for the mass and the velocity of ejecta from a generic black hole-neutron star merger based on recently performed numerical relativity simulations. We combined these fitting formulas with a new semi-analytic model for a black hole-neutron star kilonova/macronova lightcurve which reproduces the results of radiation-transfer simulations. Specifically, the semi-analytic model reproduces the result of each band magnitude obtained by the radiation transfer simulations within ~1 mag. By using this semi-analytic model, we found that, at 400 Mpc, the kilonova/macronova is as bright as 22-24 mag for the cases with a small chirp mass and a high black hole spin, and >28 mag for a large chirp mass and a low black hole spin. We also apply our model to GRB130603B as an illustration, and show that a black hole-neutron star merger with a rapidly spinning black hole and a large neutron star radius is favored.

preprint2016arXiv

OISTER Optical and Near-Infrared Observations of the Super-Chandrasekhar Supernova Candidate SN 2012dn: Dust Emission from the Circumstellar Shell

We present extensively dense observations of the super-Chandrasekhar supernova (SC SN) candidate SN 2012dn from $-11$ to $+140$ days after the date of its $B$-band maximum in the optical and near-infrared (NIR) wavelengths conducted through the OISTER ToO program. The NIR light curves and color evolutions up to 35 days after the $B$-band maximum provided an excellent match with those of another SC SN 2009dc, providing a further support to the nature of SN 2012dn as a SC SN. We found that SN 2012dn exhibited strong excesses in the NIR wavelengths from $30$ days after the $B$-band maximum. The $H$ and $K_{s}$-band light curves exhibited much later maximum dates at $40$ and $70$ days after the $B$-band maximum, respectively, compared with those of normal SNe Ia. The $H$ and $K_{s}$-band light curves subtracted by those of SN 2009dc displayed plateaued evolutions, indicating a NIR echo from the surrounding dust. The distance to the inner boundary of the dust shell is limited to be $4.8 - 6.4\times10^{-2}$ pc. No emission lines were found in its early phase spectrum, suggesting that the ejecta-CSM interaction could not occur. On the other hand, we found no signature that strongly supports the scenario of dust formation. The mass loss rate of the pre-explosion system is estimated to be $10^{-6}-10^{-5}$ M$_{\odot}$ yr$^{-1}$, assuming that the wind velocity of the system is $10-100$ km~s$^{-1}$, which suggests that the progenitor of SN 2012dn could be a recurrent nova system. We conclude that the progenitor of this SC SN could be explained by the single-degenerate scenario.

preprint2016arXiv

Rapidly Rising Transients from Subaru Hyper Suprime-Cam Transient Survey

We present rapidly rising transients discovered by a high-cadence transient survey with Subaru telescope and Hyper Suprime-Cam. We discovered five transients at z=0.384-0.821 showing the rising rate faster than 1 mag per 1 day in the restframe near-ultraviolet wavelengths. The fast rising rate and brightness are the most similar to SN 2010aq and PS1-13arp, for which the ultraviolet emission within a few days after the shock breakout was detected. The lower limit of the event rate of rapidly rising transients is ~9 % of core-collapse supernova rates, assuming a duration of rapid rise to be 1 day. We show that the light curves of the three faint objects agree with the cooling envelope emission from the explosion of red supergiants. The other two luminous objects are, however, brighter and faster than the cooling envelope emission. We interpret these two objects to be the shock breakout from dense wind with the mass loss rate of ~10^{-3} Msun yr^{-1}, as also proposed for PS1-13arp. This mass loss rate is higher than that typically observed for red supergiants. The event rate of these luminous objects is >~1 % of core-collapse supernova rate, and thus, our study implies that more than ~1 % of massive stars can experience an intensive mass loss at a few years before the explosion.

preprint2016arXiv

The Macronova in GRB 050709 and the GRB/macronova connection

GRB 050709 was the first short Gamma-ray Burst (sGRB) with an identified optical afterglow. In this work we report a re-analysis of the publicly available data of this event and the discovery of a Li-Paczynski macronova/kilonova that dominates the optical/IR signal at t $>$ 2.5 days. Such a signal would arise from ~0.05 M$_\odot$ r-process material launched by a compact binary merger. The implied mass ejection supports the suggestion that compact binary mergers are significant and possibly main sites of heavy r-process nucleosynthesis. We have re-analyzed all afterglow data from nearby short and hybrid GRBs. A statistical study of sGRB/macronova connection reveals that macronova may have taken place in all these GRBs though the fraction as low as 0.18 cannot be ruled out. The identification of two of the three macronova candidates in the I-band implies a more promising detection prospect for ground-based surveys.

preprint2015arXiv

Kepler's Supernova: An Overluminous Type Ia Event Interacting with a Massive Circumstellar Medium at a Very Late Phase

We have analyzed XMM-Newton, Chandra, and Suzaku observations of Kepler's supernova remnant (SNR) to investigate the properties of both the SN ejecta and the circumstellar medium (CSM). For comparison, we have also analyzed two similarly-aged, ejecta-dominated SNRs: Tycho's SNR, thought to be the remnant of a typical Type Ia SN, and SNR 0509-67.5 in the Large Magellanic Cloud, thought to be the remnant of an overluminous Type Ia SN. By simply comparing the X-ray spectra, we find that line intensity ratios of iron-group elements (IGE) to intermediate-mass elements (IME) for Kepler's SNR and SNR 0509-67.5 are much higher than those for Tycho's SNR. We therefore argue that Kepler is the product of an overluminous Type Ia SN. This inference is supported by our spectral modeling, which reveals the IGE and IME masses respectively to be ~0.95 M_sun and ~0.12 M_sun (Kepler's SNR), ~0.75 M_sun and ~0.34 M_sun (SNR 0509-67.5), and ~0.35 M_sun and ~0.70 M_sun (Tycho's SNR). We find that the CSM component in Kepler's SNR consists of tenuous diffuse gas (~0.3 M_sun) present throughout the entire remnant, plus dense knots (~0.035 M_sun). Both of these components have an elevated N abundance (N/H ~ 4 times the solar value), suggesting that they originate from CNO-processed material from the progenitor system. The mass of the diffuse CSM allows us to infer the pre-SN mass-loss rate to be ~1.5e-5 (V_w/10 km/s) M_sun/yr, in general agreement with results from recent hydrodynamical simulations. Since the dense knots have slow proper motions and relatively small ionization timescales, they were likely located a few pc away from the progenitor system. Therefore, we argue that Kepler's SN was an overluminous event that started to interact with massive CSM a few hundred years after the explosion. This supports the possible link between overluminous SNe and the so-called "Ia-CSM" SNe.

preprint2015arXiv

OISTER Optical and Near-Infrared Observations of Type Iax Supernova 2012Z

We report observations of the Type Iax supernova (SN Iax) 2012Z at optical and near-infrared wavelengths from immediately after the explosion until $\sim$ $260$ days after the maximum luminosity using the Optical and Infrared Synergetic Telescopes for Education and Research (OISTER) Target-of-Opportunity (ToO) program and the Subaru telescope. We found that the near-infrared (NIR) light curve evolutions and color evolutions are similar to those of SNe Iax 2005hk and 2008ha. The NIR absolute magnitudes ($M_{J}\sim-18.1$ mag and $M_{H}\sim-18.3$ mag) and the rate of decline of the light curve ($Δ$ $m_{15}$($B$)$=1.6 \pm 0.1$ mag) are very similar to those of SN 2005hk ($M_{J}\sim-17.7$ mag, $M_{H}\sim$$-18.0$ mag, and $Δ$ $m_{15}$($B$)$\sim1.6$ mag), yet differ significantly from SNe 2008ha and 2010ae ($M_{J}\sim-14 - -15$ mag and $Δ$ $m_{15}$($B$)$\sim2.4-2.7$ mag). The estimated rise time is $12.0 \pm 3.0$ days, which is significantly shorter than that of SN 2005hk or any other Ia SNe. The rapid rise indicates that the $^{56}$Ni distribution may extend into the outer layer or that the effective opacity may be lower than that in normal SNe Ia. The late-phase spectrum exhibits broader emission lines than those of SN 2005hk by a factor of 6--8. Such high velocities of the emission lines indicate that the density profile of the inner ejecta extends more than that of SN 2005hk. We argue that the most favored explosion scenario is a `failed deflagration' model, although the pulsational delayed detonations is not excluded.

preprint2015arXiv

Radioactive decay products in neutron star merger ejecta: heating efficiency and $γ$-ray emission

The radioactive decay of the freshly synthesized $r$-process nuclei ejected in compact binary mergers power optical/infrared macronovae (kilonovae) that follow these events. The light curves depend critically on the energy partition among the different products of the radioactive decay and this plays an important role in estimates of the amount of ejected $r$-process elements from a given observed signal. We study the energy partition and $γ$-ray emission of the radioactive decay. We show that $20$-$50\%$ of the total radioactive energy is released in $γ$-rays on timescales from hours to a month. The number of emitted $γ$-rays per unit energy interval has roughly a flat spectrum between a few dozen keV and $1$ MeV so that most of this energy is carried by $\sim 1$ MeV $γ$-rays. However at the peak of macronova emission the optical depth of the $γ$-rays is $\sim 0.02$ and most of the $γ$-rays escape. The loss of these $γ$-rays reduces the heat deposition into the ejecta and hence reduces the expected macronova signals if those are lanthanides dominated. This implies that the ejected mass is larger by a factor of $2$-$3$ than what was previously estimated. Spontaneous fission heats up the ejecta and the heating rate can increase if a sufficient amount of transuranic nuclei are synthesized. Direct measurements of these escaping $γ$-rays may provide the ultimate proof for the macronova mechanisms and an identification of the $r$-process nucleosynthesis sites. However, the chances to detect these signals are slim with current X-ray and $γ$-ray missions. New detectors, more sensitive by at least a factor of ten, are needed for a realistic detection rate.

preprint2014arXiv

Discovery of Dramatic Optical Variability in SDSS J1100+4421: A Peculiar Radio-Loud Narrow-Line Seyfert 1 Galaxy?

We present our discovery of dramatic variability in SDSS J1100+4421 by the high-cadence transient survey Kiso Supernova Survey (KISS). The source brightened in the optical by at least a factor of three within about half a day. Spectroscopic observations suggest that this object is likely a narrow-line Seyfert 1 galaxy (NLS1) at z=0.840, however with unusually strong narrow emission lines. The estimated black hole mass of ~ 10^7 Msun implies bolometric nuclear luminosity close to the Eddington limit. SDSS J1100+4421 is also extremely radio-loud, with a radio loudness parameter of R ~ 4 x 10^2 - 3 x 10^3, which implies the presence of relativistic jets. Rapid and large-amplitude optical variability of the target, reminiscent of that found in a few radio- and gamma-ray loud NLS1s, is therefore produced most likely in a blazar-like core. The 1.4 GHz radio image of the source shows an extended structure with a linear size of about 100 kpc. If SDSS J1100+4421 is a genuine NLS1, as suggested here, this radio structure would then be the largest ever discovered in this type of active galaxies.

preprint2014arXiv

Early-phase photometry and spectroscopy of transitional Type Ia SN 2012ht: Direct constraint on the rise time

We report photometric and spectroscopic observations of the nearby Type Ia Supernova (SN Ia) 2012ht from $-15.8$ days to $+49.1$ days after $B$-band maximum. The decline rate of the light curve is $Δm_{15}$($B$)$=1.39~\pm~0.05$ mag, which is intermediate between normal and subluminous SNe Ia, and similar to that of the `transitional' Type Ia SN 2004eo. The spectral line profiles also closely resemble those of SN 2004eo. We were able to observe SN 2012ht at very early phase, when it was still rising and was about three magnitudes fainter than at the peak. The rise time to the $B$-band maximum is estimated to be $17.6 \pm 0.5$ days and the time of the explosion is MJD $56277.98 \pm 0.13$. SN 2012ht is the first transitional SN Ia whose rise time is directly measured without using light curve templates, and the fifth SN Ia overall. This rise time is consistent with those of the other four SNe within the measurement error, even including the extremely early detection of SN 2013dy. The rising part of the light curve can be fitted by a quadratic function, and shows no sign of a shock-heating component due to the interaction of the ejecta with a companion star. The rise time is significantly longer than that inferred for subluminous SNe such as SN 1991bg, which suggests that a progenitor and/or explosion mechanism of transitional SNe Ia are more similar to normal SNe Ia rather than subluminous SNe Ia.

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

Radiation hydrodynamics simulations of wide-angle outflows from super-critical accretion disks around black holes

By performing two-dimensional radiation hydrodynamics simulations with large computational domain of 5000 Schwarzschild radius, we revealed that wide-angle outflow is launched via the radiation force from the super-critical accretion flows around black holes. The angular size of the outflow, of which the radial velocity (v_r) is over the escape velocity (v_esc), increases with an increase of the distance from the black hole. As a result, the mass is blown away with speed of v_r > v_esc in all direction except for the very vicinity of the equatorial plane, theta=0-85^circ, where theta is the polar angle. The mass ejected from the outer boundary per unit time by the outflow is larger than the mass accretion rate onto the black hole, ~150L_Edd/c^2, where L_Edd and c are the Eddington luminosity and the speed of light. Kinetic power of such wide-angle high-velocity outflow is comparable to the photon luminosity and is a few times larger than the Eddington luminosity. This corresponds to ~10^39-10^40 erg/s for the stellar mass black holes. Our model consistent with the observations of shock excited bubbles observed in some ultra-luminous X-ray sources (ULXs), supporting a hypothesis that ULXs are powered by the super-critical accretion onto stellar mass black holes.

preprint2014arXiv

Random walks and effective optical depth in relativistic flow

We investigate the random walk process in relativistic flow. In the relativistic flow, photon propagation is concentrated in the directions of the flow velocity due to relativistic beaming effect. We show that, in the pure scattering case, the number of scatterings is proportional to the size parameter $ξ\equiv L/l_0$ if the flow velocity $β\equiv v/c$ satisfies $β/Γ\gg ξ^{-1}$, while it is proportional to $ξ^2$ if $β/Γ\ll ξ^{-1}$ where $L$ and $l_0$ are the size of the system in the observer frame and the mean free path in the comoving frame, respectively. We also examine the photon propagation in the scattering and absorptive medium. We find that, if the optical depth for absorption $τ_{\rm a}$ is considerably smaller than the optical depth for scattering $τ_{\rm s}$ ($τ_{\rm a}/τ_{\rm s} \ll 1$) and the flow velocity satisfies $β\gg \sqrt{2τ_{\rm a}/τ_{\rm s}}$, the effective optical depth is approximated by $τ_*\simeqτ_{\rm a}(1+β)/β$. Furthermore, we perform Monte Carlo simulations of radiative transfer and compare the results with the analytic expression for the number of scattering. The analytic expression is consistent with the results of the numerical simulations. The expression derived in this Letter can be used to estimate the photon production site in relativistic phenomena, e.g., gamma-ray burst and active galactic nuclei.

preprint2014arXiv

The red supergiant and supernova rate problems: implications for core-collapse supernova physics

Mapping supernovae to their progenitors is fundamental to understanding the collapse of massive stars. We investigate the red supergiant problem, which concerns why red supergiants with masses $\sim16$-$30 M_\odot$ have not been identified as progenitors of Type IIP supernovae, and the supernova rate problem, which concerns why the observed cosmic supernova rate is smaller than the observed cosmic star formation rate. We find key physics to solving these in the compactness parameter, which characterizes the density structure of the progenitor. If massive stars with compactness above $ξ_{2.5} \sim 0.2$ fail to produce canonical supernovae, (i) stars in the mass range $16$-$30 M_\odot$ populate an island of stars that have high $ξ_{2.5}$ and do not produce canonical supernovae, and (ii) the fraction of such stars is consistent with the missing fraction of supernovae relative to star formation. We support this scenario with a series of two- and three-dimensional radiation hydrodynamics core-collapse simulations. Using more than 300 progenitors covering initial masses $10.8$-$75 M_\odot$ and three initial metallicities, we show that high compactness is conducive to failed explosions. We then argue that a critical compactness of $\sim 0.2$ as the divide between successful and failed explosions is consistent with state-of-the-art three-dimensional core-collapse simulations. Our study implies that numerical simulations of core collapse need not produce robust explosions in a significant fraction of compact massive star initial conditions.

preprint2013arXiv

Detectability of High-Redshift Superluminous Supernovae with Upcoming Optical and Near-Infrared Surveys - II. Beyond z=6

Observational identification of the first stars is one of the great challenges in the modern astronomy. Although a single first star is too faint to be detected, supernova explosions of the first stars can be bright enough. An important question is whether such supernovae can be detected in the limited observational area with realistic observational resources. We perform detailed simulations to study the detectability of superluminous supernovae (SLSNe) at high redshifts, using the observationally-calibrated star formation rate density and supernova occurrence rate. We show that a 100 deg^2 survey with the limiting magnitude of 26 mag in near-infrared wavelengths will be able to discover about 10 SLSNe at z > 10. If the survey is extended to 200 deg^2 with 27 mag depth, about 10 SLSNe can be discovered at z > 15. We emphasize that the observations at > 3 um are crucial to detect and select SLSNe at z > 10. Our simulations are also applied to the planned survey with Euclid, WFIRST, and WISH. These surveys will be able to detect about 1000, 400, and 3000 SLSNe up to z ~ 5, 7, and 12, respectively. We conclude that detection of SLSNe at z > 10 is in fact achievable in the near future.

preprint2013arXiv

Early UV/Optical Emission of The Type Ib SN 2008D

We propose an alternative explanation for the post-breakout emission of SN 2008D associated with the X-ray transient 080109. Observations of this object show a very small contrast of 0.35 dex between the light-curve minimum occurring soon after the breakout, and the main luminosity peak that is due to radioactive heating of the ejecta. Hydrodynamical models show that the cooling of a shocked Wolf-Rayet star leads to a much greater difference (> 0.9 dex). Our proposed scenario is that of a jet produced during the explosion which deposits 56Ni-rich material in the outer layers of the ejecta. The presence of high-velocity radioactive material allows us to reproduce the complete luminosity evolution of the object. Without outer 56Ni it could be possible to reproduce the early emission purely from cooling of the shocked envelope by assuming a larger progenitor than a Wolf-Rayet star, but that would require an initial density structure significantly different from what is predicted by stellar evolution models. Analytic models of the cooling phase have been proposed reproduce the early emission of SN 2008D with an extended progenitor. However, we found that the models are valid only until 1.5 days after the explosion where only two data of SN 2008D are available. We also discuss the possibility of the interaction of the ejecta with a binary companion, based on published analytic expressions. However, the binary separation required to fit the early emission should be < 3 Rsun which is too small for a system containing two massive stars.

preprint2013arXiv

Progenitor models of the electromagnetic transient associated with the short GRB 130603B

An electromagnetic transient powered by the radioactive decay of r-process elements, a so-called kilonova/macronova, is one of the possible observable consequences of compact binary mergers including at east one neutron star. Recent observations strongly suggest the first discovery of the electromagnetic transient, which is associated with the short GRB 130603B. We explore a possible progenitor of this event combining the numerical-relativity simulations and radiative transfer simulations of the dynamical ejecta of binary neutron star and black hole - neutron star mergers. We show that the ejecta models within a realistic parameter range consistently reproduce the observed near-infrared excess. We also show that the soft equation of state models for binary neutron star mergers and the stiff equation of state models for black hole - neutron star mergers are favored to reproduce the observed luminosity.

preprint2013arXiv

Radiative Transfer Simulations for Neutron Star Merger Ejecta

The merger of binary neutron stars (NSs) is among the most promising gravitational wave (GW) sources. Next-generation GW detectors are expected to detect signals from the NS merger within 200 Mpc. Detection of electromagnetic wave (EM) counterpart is crucial to understand the nature of GW sources. Among possible EM emission from the NS merger, emission powered by radioactive r-process nuclei is one of the best targets for follow-up observations. However, prediction so far does not take into account detailed r-process element abundances in the ejecta. We perform radiative transfer simulations for the NS merger ejecta including all the r-process elements from Ga to U for the first time. We show that the opacity in the NS merger ejecta is about kappa = 10 cm^2 g^{-1}, which is higher than that of Fe-rich Type Ia supernova ejecta by a factor of ~ 100. As a result, the emission is fainter and longer than previously expected. The spectra are almost featureless due to the high expansion velocity and bound-bound transitions of many different r-process elements. We demonstrate that the emission is brighter for a higher mass ratio of two NSs and a softer equation of states adopted in the merger simulations. Because of the red color of the emission, follow-up observations in red optical and near-infrared (NIR) wavelengths will be the most efficient. At 200 Mpc, expected brightness of the emission is i = 22 - 25 AB mag, z = 21 - 23 AB mag, and 21 - 24 AB mag in NIR JHK bands. Thus, observations with wide-field 4m- and 8m-class optical telescopes and wide-field NIR space telescopes are necessary. We also argue that the emission powered by radioactive energy can be detected in the afterglow of nearby short gamma-ray bursts.

preprint2013arXiv

Radioactively Powered Emission from Black Hole-Neutron Star Mergers

Detection of electromagnetic counterparts of gravitational wave (GW) sources is important to unveil the nature of compact binary coalescences. We perform three-dimensional, time-dependent, multi-frequency radiative transfer simulations for radioactively powered emission from the ejecta of black hole (BH) - neutron star (NS) mergers. Depending on the BH to NS mass ratio, spin of the BH, and equations of state of dense matter, BH-NS mergers can eject more material than NS-NS mergers. In such cases, radioactively powered emission from the BH-NS merger ejecta can be more luminous than that from NS-NS mergers. We show that, in spite of the expected larger distances to BH-NS merger events, observed brightness of BH-NS mergers can be comparable to or even higher than that of NS-NS mergers. We find that, when the tidally disrupted BH-NS merger ejecta are confined to a small solid angle, the emission from BH-NS merger ejecta tends to be bluer than that from NS-NS merger ejecta for a given total luminosity. Thanks to this property, we might be able to distinguish BH-NS merger events from NS-NS merger events by multi-band observations of the radioactively powered emission. In addition to the GW observations, such electromagnetic observations can potentially provide independent information on the nature of compact binary coalescences.

preprint2013arXiv

The very energetic, broad-lined type Ic Supernova 2010ah (PTF10bzf) in the context of GRB/SNe

SN2010ah, a very broad-lined type Ic SN discovered by the Palomar Transient Factory, was interesting because of its relatively high luminosity and the high velocity of the absorption lines, which was comparable to that of GRB/SNe, suggesting a high explosion kinetic energy. However, no GRB was detected in association with the SN. Here, the properties of SN2010ah are determined with higher accuracy than previous studies through modelling. New Subaru telescope photometry is presented. A bolometric light curve is constructed taking advantage of the spectral similarity with SN1998bw. Radiation transport tools are used to reproduce the spectra and the light curve. The results thus obtained regarding ejecta mass, composition and kinetic energy are then used to compute a synthetic light curve. This is in reasonable agreement with the early bolometric light curve of SN2010ah, but a high abundance of 56Ni at high velocity is required to reproduce the early rise, while a dense inner core must be used to reproduce the slow decline at late phases. The high-velocity 56Ni cannot have been located on our line of sight, which may be indirect evidence for an off-axis, aspherical explosion. The main properties of SN2010ah are: ejected mass ~ 3 Mo; kinetic energy ~10^52 erg, M(56Ni) ~ 0.25 Mo. The mass located at v >~ 0.1c is ~0.2 Mo. Although these values, in particular the kinetic energy, are quite large for a SN Ic, they are all smaller (especially the ejecta mass) than those typical of GRB/SNe. This confirms the tendency for these quantities to correlate, and suggests that there are minimum requirements for a GRB/SN, which SN2010ah may not meet although it comes quite close. Depending on whether a neutron star or a black hole was formed following core collapse, SN2010ah was the explosion of a CO core of ~ 5-6 Mo, pointing to a progenitor mass of ~24 - 28 Mo.

preprint2012arXiv

A Search for Infrared Emission from Core-Collapse Supernovae at the Transitional Phase

Most of the observational studies of supernova (SN) explosions are limited to early phases (< a few yr after the explosion) of extragalactic SNe and observations of SN remnants (> 100 yr) in our Galaxy or very nearby galaxies. SNe at the epoch between these two, which we call "transitional" phase, have not been explored in detail except for several extragalactic SNe including SN 1987A in the Large Magellanic Cloud. We present theoretical predictions for the infrared (IR) dust emissions by several mechanisms; emission from dust formed in the SN ejecta, light echo by circumstellar and interstellar dust, and emission from shocked circumstellar dust. We search for IR emission from 6 core-collapse SNe at the transitional phase in the nearby galaxies NGC 1313, NGC 6946, and M101 by using the data taken with the AKARI satellite and Spitzer. Among 6 targets, we detect the emission from SN 1978K in NGC 1313. SN 1978K is associated with 1.3 x 10^{-3} Msun of silicate dust. We show that, among several mechanisms, the shocked circumstellar dust is the most probable emission source to explain the IR emission observed for CSM-rich SN 1978K. IR emission from the other 5 objects is not detected. Our current observations are sensitive to IR luminosity of > 10^{38} erg s^{-1}, and the non-detection of SN 1962M excludes the existence of the shocked circumstellar dust for a high gas mass-loss rate of sim 10^{-4} Msun yr^{-1}. Observations of SNe at the transitional phase with future IR satellites will fill the gap of IR observations of SNe with the age of 10-100 years, and give a new opportunity to study the circumstellar and interstellar environments of the progenitor, and possibly dust formation in SNe.

preprint2012arXiv

A theoretical color-velocity correlation for supernovae associated with gamma-ray bursts

We carry out the first multi-dimensional radiative transfer calculations to simultaneously compute synthetic spectra and light curves for models of supernovae driven by fast bipolar outflows. These allow us to make self-consistent predictions for the orientation dependence of both colour evolution and spectral features. We compare models with different degrees of asphericity and metallicity and find significant observable consequences of both. In aspherical models, we find spectral and light curve features that vary systematically with observer orientation. In particular, we find that the early phase light curves are brighter and bluer when viewed close to the polar axis but that the peak flux is highest for equatorial (off-axis) inclinations. Spectral line features also depend systematically on observer orientation, including the velocity of the SiII 6355A line. Consequently, our models predict a correlation between line velocity and color that could assist the identification of supernovae associated with off-axis jet-driven explosions. The amplitude and range of this correlation depends on the degree of asphericity, the metallicity and the epoch of observation but we find that it is always present and acts in the same direction.

preprint2012arXiv

Detectability of High-Redshift Superluminous Supernovae with Upcoming Optical and Near-Infrared Surveys

Observations of high-redshift supernovae (SNe) open a novel opportunity to study the massive star population in the early Universe. We study the detectability of superluminous SNe with upcoming optical and near-infrared (NIR) surveys. Our calculations are based on the cosmic star formation history, the SN occurence rate, the characteristic colour and the light curve of the SNe that are all calibrated by available observations. We show that 15-150 SNe up to z ~ 4 will be discovered by the proposed Subaru/Hyper Suprime-Cam deep survey: 30 deg^2 survey with 24.5 AB mag depth in z-band for 3 months. With its ultra-deep layer (3.5 deg^2 with 25.6 AB mag depth in z-band for 4 months), the highest redshift can be extended to z ~ 5. We further explore the detectability by upcoming NIR survey utilizing future satellites such as Euclid, WFIRST, and WISH. The wide-field NIR surveys are very efficient to detect high-redshift SNe. With a hypothetical deep NIR survey for 100 deg^2 with 26 AB mag depth at 1-4 um, at least ~ 50 SNe will be discovered at z>3 in half a year. The number of the detected SNe can place a strong constraint on the stellar initial mass function or its slope especially at the high-mass end. Superluminous SNe at high redshifts can be distinguished from other types of SNe by the long time-scale of their light curves in the observer's frame, the optical colours redder than other core-collapse SNe and the NIR colours redder than any other types of SNe.

preprint2012arXiv

Elemental Abundances in the Possible Type Ia Supernova Remnant G344.7-0.1

Recent studies on the Galactic supernova remnant (SNR) G344.7-0.1 have commonly claimed its origin to be a core-collapse supernova (SN) explosion, based on its highly asymmetric morphology and/or proximity to a star forming region. In this paper, however, we present an X-ray spectroscopic study of this SNR using Suzaku, which is supportive of a Type Ia origin. Strong K-shell emission from lowly ionized Fe has clearly been detected, and its origin is determined, for the first time, to be the Fe-rich SN ejecta. The abundance pattern is highly consistent with that expected for a somewhat-evolved Type Ia SNR. It is suggested, therefore, that the X-ray point-like source CXOU J170357.8-414302 located at the SNR's geometrical center is not associated with the SNR but is likely to be a foreground object. Our result further indicates that G344.7-0.1 is the first possible Type Ia SNR categorized as a member of the so-called "mixed-morphology" class. In addition, we have detected emission from He-like Al at ~1.6 keV, the first clear detection of this element in the spectrum of an extended X-ray source. The possible enhancement of the Al/Mg abundance ratio from the solar value suggests that the ambient interstellar medium has a relatively high metallicity (not less than 10% of the solar value), if this SNR has indeed a Type Ia origin. We also report marginal detection of Cr and Mn, although the measured fluxes have large statistical and systematic uncertainties.

preprint2012arXiv

Light Curve Modeling of Superluminous Supernova 2006gy: Collision between Supernova Ejecta and Dense Circumstellar Medium

We show model light curves of superluminous supernova 2006gy on the assumption that the supernova is powered by the collision of supernova ejecta and its dense circumstellar medium. The initial conditions are constructed based on the shock breakout condition, assuming that the circumstellar medium is dense enough to cause the shock breakout within it. We perform a set of numerical light curve calculations by using a one-dimensional multigroup radiation hydrodynamics code STELLA. We succeeded in reproducing the overall features of the early light curve of SN 2006gy with the circumstellar medium whose mass is about 15 Msun (the average mass-loss rate ~ 0.1 Msun/yr). Thus, the progenitor of SN 2006gy is likely a very massive star. The density profile of the circumstellar medium is not well constrained by the light curve modeling only, but our modeling disfavors the circumstellar medium formed by steady mass loss. The ejecta mass is estimated to be comparable to or less than 15 Msun and the explosion energy is expected to be more than 4e51 erg. No 56Ni is required to explain the early light curve. We find that the multidimensional effect, e.g., the Rayleigh-Taylor instability, which is expected to take place in the cool dense shell between the supernova ejecta and the dense circumstellar medium, is important in understanding supernovae powered by the shock interaction. We also show the evolution of the optical and near-infrared model light curves of high-redshift superluminous supernovae. They can be potentially used to identify SN 2006gy-like superluminous supernovae in the future optical and near-infrared transient surveys.

preprint2012arXiv

Super-Chandrasekhar-Mass Light Curve Models for the Highly Luminous Type Ia Supernova 2009dc

Several highly luminous Type Ia supernovae (SNe Ia) have been discovered. Their high luminosities are difficult to explain with the thermonuclear explosions of the Chandrasekhar-mass white dwarfs (WDs). In the present study, we estimate the progenitor mass of SN 2009dc, one of the extremely luminous SNe Ia, using the hydrodynamical models as follows. Explosion models of super-Chandrasekhar-mass (super-Ch-mass) WDs are constructed, and multi-color light curves (LCs) are calculated. The comparison between our calculations and the observations of SN 2009dc suggests that the exploding WD has a super-Ch mass of 2.2-2.4 solar masses, producing 1.2-1.4 solar masses of Ni-56, if the extinction by its host galaxy is negligible. If the extinction is significant, the exploding WD is as massive as \sim2.8 solar masses, and \sim1.8 solar masses of Ni-56 is necessary to account for the observations. Whether the host-galaxy extinction is significant or not, the progenitor WD must have a thick carbon-oxygen layer in the outermost zone (20-30% of the WD mass), which explains the observed low expansion velocity of the ejecta and the presence of carbon. Our estimate on the mass of the progenitor WD, especially for the extinction-corrected case, is challenging to the current scenarios of SNe Ia. Implications on the progenitor scenarios are also discussed.

preprint2012arXiv

The Highly Energetic Expansion of SN2010bh Associated with GRB 100316D

We present the spectroscopic and photometric evolution of the nearby (z = 0.059) spectroscopically confirmed type Ic supernova, SN 2010bh, associated with the soft, long-duration gamma-ray burst (X-ray flash) GRB 100316D. Intensive follow-up observations of SN 2010bh were performed at the ESO Very Large Telescope (VLT) using the X-shooter and FORS2 instruments. Owing to the detailed temporal coverage and the extended wavelength range (3000--24800 A), we obtained an unprecedentedly rich spectral sequence among the hypernovae, making SN 2010bh one of the best studied representatives of this SN class. We find that SN 2010bh has a more rapid rise to maximum brightness (8.0 +/- 1.0 rest-frame days) and a fainter absolute peak luminosity (L_bol~3e42 erg/s) than previously observed SN events associated with GRBs. Our estimate of the ejected (56)Ni mass is 0.12 +/- 0.02 Msun. From the broad spectral features we measure expansion velocities up to 47,000 km/s, higher than those of SNe 1998bw (GRB 980425) and 2006aj (GRB 060218). Helium absorption lines He I lambda5876 and He I 1.083 microm, blueshifted by ~20,000--30,000 km/s and ~28,000--38,000 km/s, respectively, may be present in the optical spectra. However, the lack of coverage of the He I 2.058 microm line prevents us from confirming such identifications. The nebular spectrum, taken at ~186 days after the explosion, shows a broad but faint [O I] emission at 6340 A. The light-curve shape and photospheric expansion velocities of SN 2010bh suggest that we witnessed a highly energetic explosion with a small ejected mass (E_k ~ 1e52 erg and M_ej ~ 3 Msun). The observed properties of SN 2010bh further extend the heterogeneity of the class of GRB supernovae.

preprint2012arXiv

Three-Dimensional Explosion Geometry of Stripped-Envelope Core-Collapse Supernovae. I. Spectropolarimetric Observations

We study the multi-dimensional geometry of supernova (SN) explosions by means of spectropolarimetric observations of stripped-envelope SNe, i.e., SNe without a H-rich layer. We perform spectropolarimetric observations of 2 stripped-envelope SNe, the Type Ib SN 2009jf and the Type Ic SN 2009mi. Both objects show non-zero polarization at the wavelength of the strong lines. They also show a loop in the Stokes Q-U diagram, which indicates a non-axisymmetric, three-dimensional ion distribution in the ejecta. We show that five out of six stripped-envelope SNe which have been observed spectropolarimetrically so far show such a loop. This implies that a three-dimensional geometry is common in stripped-envelope SNe. We find that stronger lines tend to show higher polarization. This effect is not related to the geometry, and must be corrected to compare the polarization of different lines or different objects. Even after the correction, however, there remains a dispersion of polarization degree among different objects. Such a dispersion might be caused by three-dimensional clumpy ion distributions viewed from different directions.

preprint2012arXiv

Unusual Long and Luminous Optical Transient in the Subaru Deep Field

We present observations of SDF-05M05, an unusual optical transient discovered in the Subaru Deep Field (SDF). The duration of the transient is > ~800 d in the observer frame, and the maximum brightness during observation reached approximately 23 mag in the i' and z' bands. The faint host galaxy is clearly identified in all 5 optical bands of the deep SDF images. The photometric redshift of the host yields z~0.6 and the corresponding absolute magnitude at maximum is ~-20. This implies that this event shone with an absolute magnitude brighter than -19 mag for approximately 300 d in the rest frame, which is significantly longer than a typical supernova and ultra-luminous supernova. The total radiated energy during our observation was 1x10^51 erg. The light curves and color evolution are marginally consistent with some of luminous IIn supernova. We suggest that the transient may be a unique and peculiar supernova at intermediate redshift.

preprint2011arXiv

Formation of Dust in the Ejecta of Type Ia Supernovae

We investigate the formation of dust grains in the ejecta of Type Ia supernovae (SNe Ia), adopting the carbon-deflagration W7 model. In the calculations of dust formation, we apply the nucleation and grain growth theory and consider the two cases with and without formation of CO and SiO molecules. The results of the calculations show that for the sticking probability of alpha_j=1, C, silicate, Si, and FeS grains can condense at early times of ~100--300 days after the explosion, whereas Fe and SiC grains cannot form substantially. Due to the low gas density in SNe Ia with no H-envelope, the average radii of the newly formed grains are generally below 0.01 micron, being much smaller than those in Type II-P SNe. This supports our previous conclusion that the radius of dust formed in the ejecta is smaller in SNe with less massive envelopes. The total dust mass ranges from 3x10^{-4} M_sun to 0.2 M_sun for alpha_j=0.1--1, depending on whether or not CO and SiO molecules are formed. We also estimate the optical depths and thermal emission by the newly formed dust and compare to the relevant observations of SNe Ia. We find that the formation of C grains in SNe Ia is suppressed to be consistent with the observational constraints. This implies that energetic photons and electrons heavily depress the formation efficiency of C grains or that the outermost C-O layer of SNe Ia is almost fully burned. Finally, we perform the calculations of dust destruction in the SN remnants and find that dust grains formed in the ejecta of SNe Ia are almost completely destroyed in the shocked gas before being injected into the interstellar medium. This indicates that SNe Ia are unlikely to be the major sources of interstellar dust.

preprint2011arXiv

Unburned Material in the Ejecta of Type Ia Supernovae

(Abridged) The presence of unburned material in the ejecta of normal Type Ia supernovae (SNeIa) is investigated using early-time spectroscopy obtained by the CSP. The tell-tale signature of pristine material from a C+O white dwarf progenitor star is the presence of carbon. The most prominent carbon lines in optical spectra of SNeIa are due to C II. We find that at least 30% of the objects in the sample show an absorption at ~6300 AA which is attributed to C II 6580. An association to H_alpha is considered to be unlikely. These findings imply a larger incidence of carbon in SNeIa ejecta than previously noted. We show how observational biases and physical conditions may hide these C II lines, and account for the scarcity of previous carbon detections. This relatively large frequency of carbon detections has crucial implications on our understanding of the explosive process. The observations indicate that carbon is present merely ~1000 km s^-1 above the bulk of Si II. Spectral modeling shows the detections are consistent with a mass of carbon of 10^-3--10^-2 M_sun. This would imply substantial mixing, which may be related to asymmetries of the flame propagation or clumps of unburned material along the line of sight. However, the uniformity of the relation between C II and Si II velocities is not consistent with the latter scenario. The observational properties of SNeIa with and without carbon signatures are compared. A trend toward bluer color and lower luminosity is found for objects with carbon.

preprint2010arXiv

A Core-Collapse Supernova Model for the Extremely Luminous Type Ic Supernova 2007bi: An Alternative to the Pair-Instability Supernova Model

We present a core-collapse supernova model for the extremely luminous Type Ic supernova 2007bi. By performing numerical calculations of hydrodynamics, nucleosynthesis, and radiation transport, we find that SN 2007bi is consistent with the core-collapse supernova explosion of a 43 Msun carbon and oxygen core obtained from the evolution of a progenitor star with a main sequence mass of 100 Msun and metallicity of Z = Zsun/200, from which its hydrogen and helium envelopes are artificially stripped. The ejecta mass and the ejecta kinetic energy of the models are 40 Msun and 3.6*10^{52} erg. The ejected 56Ni mass is as large as 6.1 Msun, which results from the explosive nucleosynthesis with large explosion energy. We also confirm that SN 2007bi is consistent with a pair-instability supernova model as has recently been claimed. We show that the earlier light curve data can discriminate between the models for such luminous supernovae.

preprint2010arXiv

Abundance stratification in Type Ia supernovae - III. The normal SN 2003du

The element abundance distributions in the ejecta of Type Ia supernova (SN) is studied by modelling a time series of optical spectra of SN 2003du until ~1 year after the explosion. Since SN 2003du is a very normal Type Ia SN both photometrically and spectroscopically, the abundance distribution derived for it can be considered as representative of normal Type Ia SNe. We find that the innermost layers are dominated by stable Fe-group elements, with a total mass of ~ 0.2 Msun, which are synthesized through electron capture. Above the core of stable elements there are thick 56Ni-rich layers. The total mass of 56Ni is 0.65 Msun. The Si- and S-rich layers are located above the 56Ni-rich layers. The dominant element in the outermost layers (M_r > 1.1 Msun, v > 13000 km/s) is O, with a small amount of Si. Little unburned C remains, with an upper limit of 0.016 Msun. The element distributions in the ejecta are moderately mixed, but not fully mixed as seen in three-dimensional deflagration models.

preprint2010arXiv

Fallback Supernovae: A Possible Origin of Peculiar Supernovae with Extremely Low Explosion Energies

We perform hydrodynamical calculations of core-collapse supernovae (SNe) with low explosion energies. These SNe do not have enough energy to eject the whole progenitor and most of the progenitor falls back to the central remnant. We show that such fallback SNe can have a variety of light curves (LCs) but their photospheric velocities can only have some limited values with lower limits. We also perform calculations of nucleosynthesis and LCs of several fallback SN model, and find that a fallback SN from the progenitor with a main-sequence mass of 13 Msun can account for the properties of the peculiar Type Ia supernova SN 2008ha. The kinetic energy and ejecta mass of the model are 1.2*10^{48} erg and 0.074 Msun, respectively, and the ejected 56Ni mass is 0.003 Msun. Thus, SN 2008ha can be a core-collapse SN with a large amount of fallback. We also suggest that SN 2008ha could have been accompanied with long gamma-ray bursts and long gamma-ray bursts without associated SNe may be accompanied with very faint SNe with significant amount of fallback which are similar to SN 2008ha.

preprint2010arXiv

Spectropolarimetry of Extremely Luminous Type Ia Supernova 2009dc: Nearly Spherical Explosion of Super-Chandrasekhar Mass White Dwarf

We present the first spectropolarimetric observations of a candidate of super-Chandrasekhar mass Type Ia supernova (SN): SN 2009dc. The observations were performed at 5.6 and 89.5 days after the B-band maximum. The data taken at the later epoch are used to determine the interstellar polarization. Continuum polarization is found to be small (<0.3 %), indicating that the explosion is nearly spherically symmetric. This fact suggests that a very aspherical explosion is not a likely scenario for SN 2009dc. Polarization at the Si II and Ca II lines clearly shows a loop in the Q-U plane, indicating a non-axisymmetric, clumpy distribution of intermediate-mass elements. The degree of line polarization at the Si and Ca line is moderate (0.5% +- 0.1% and 0.7% +- 0.1%, respectively), but it is higher than expected from the trend of other Type Ia SNe. This may suggest that there are thick enough, clumpy Si-rich layers above the thick 56Ni-rich layers (>~ 1.2 Msun). The observed spectropolarimetric properties, combined with the photometric and spectroscopic properties, suggest that the progenitor of SN 2009dc has a super-Chandrasekhar mass, and that the explosion geometry is globally spherically symmetric, with clumpy distribution of intermediate-mass elements.

preprint2009arXiv

Early Spectral Evolution of the Rapidly Expanding Type Ia SN 2006X

We present optical spectroscopic and photometric observations of Type Ia supernova (SN) 2006X from --10 to +91 days after the $B$-band maximum. This SN exhibits one of the highest expansion velocity ever published for SNe Ia. At premaximum phases, the spectra show strong and broad features of intermediate-mass elements such as Si, S, Ca, and Mg, while the O{\sc i}$λ$7773 line is weak. The extremely high velocities of Si{\sc ii} and S{\sc ii} lines and the weak O{\sc i} line suggest that an intense nucleosynthesis might take place in the outer layers, favoring a delayed detonation model. Interestingly, Si{\sc ii}$λ$5972 feature is quite shallow, resulting in an unusually low depth ratio of Si{\sc ii}$λ$5972 to $λ$6355, $\cal R$(Si{\sc ii}). The low $\cal R$(Si{\sc ii}) is usually interpreted as a high photospheric temperature. However, the weak Si{\sc iii}$λ$4560 line suggests a low temperature, in contradiction to the low $\cal R$(Si{\sc ii}). This could imply that the Si{\sc ii}$λ$5972 line might be contaminated by underlying emission. We propose that $\cal R$(Si{\sc ii}) may not be a good temperature indicator for rapidly expanding SNe Ia at premaximum phases.

preprint2009arXiv

Spectropolarimetry of the Unique Type Ib Supernova 2005bf: Larger Asymmetry Revealed by Later-Phase Data

We present an optical spectropolarimetric observation of the unique Type Ib supernova (SN) 2005bf at 8 days after the second maximum. The data, combined with the polarization spectrum taken at 6 days before the second maximum (Maund et al. 2007a), enable us to closely examine the intrinsic properties of the SN. The polarization percentage is smaller at the later epoch over a wide wavelength range, while the position angle is similar at the two epochs. We find that an assumption of complete depolarization of strong lines at the emission peak is not necessarily correct. The intrinsic polarization of the SN is larger, and thus, the ejecta of SN 2005bf would be more asymmetric than previously expected. The axis ratio of the photosphere projected on the sky deviates from unity by at least 20 %. If the position angle of interstellar polarization is aligned with the spiral arm of the host galaxy, the deviation is larger than 25 %. The line polarization at the He I, Ca II and Fe II lines is also smaller at the later epoch. The changes in the position angle across these lines, which were observed at the earlier epoch, are still marginally present at the later epoch. The similar polarimetric behavior suggests that the distributions of these ions are correlated. Properties of polarization, as well as the light curve and the spectra both at photospheric and nebular phases, can be explained by an aspherical, possibly unipolar explosion of a WN star in which the blob of 56Ni penetrates C+O core and stops within the He layer.

preprint2009arXiv

The broad line type Ic supernova SN 2007ru: Adding to the diversity of type Ic supernovae

Photometric and spectral evolution of the Type Ic supernova SN 2007ru until around 210 days after maximum are presented. The spectra show broad spectral features due to very high expansion velocity, normally seen in hypernovae. The photospheric velocity is higher than other normal Type Ic supernovae. It is lower than SN 1998bw at $\sim$ 8 days after the explosion, but is comparable at later epochs. The light curve evolution of SN 2007ru indicates a fast rise time of 8$\pm$3 days to $B$ band maximum and post-maximum decline more rapid than other broad-line Type Ic supernovae. With an absolute $V$ magnitude of -19.06, SN 2007ru is comparable in brightness with SN 1998bw and lies at the brighter end of the observed Type Ic supernovae. The ejected mass of \Nifs is estimated to be $\sim0.4\Msun$. The fast rise and decline of the light curve and the high expansion velocity suggest that SN 2007ru is an explosion with a high kinetic energy/ejecta mass ratio ($E_{\rm K}/M_{\rm {ej}}$). This adds to the diversity of Type Ic supernovae. Although the early phase spectra are most similar to those of broad-line SN 2003jd, the [OI] line profile in the nebular spectrum of SN 2007ru shows the singly-peaked profile, in contrast to the doubly-peaked profile in SN 2003jd. The singly-peaked profile, together with the high luminosity and the high expansion velocity, may suggest that SN 2007ru could be an aspherical explosion viewed from the polar direction. Estimated oxygen abundance 12 + log(O/H) of $\sim$8.8 indicates that SN 2007ru occurred in a region with nearly solar metallicity.

preprint2009arXiv

The Extremely Luminous Supernova 2006gy at Late Phase: Detection of Optical Emission from Supernova

We performed optical spectroscopy and photometry of SN 2006gy at late time, ~400 days after the explosion, with the Subaru/FOCAS in a good seeing condition. We found that the SN faded by ~3 mag from ~200 to ~400 days after the explosion (i.e., by ~5 mag from peak to ~400 days) in R band. The overall light curve is marginally consistent with the 56Ni heating model, although the flattening around 200 days suggests the optical flux declined more steeply between ~200 and ~400 days. The late time spectrum was quite peculiar among all types of SNe. It showed many intermediate width (~2000 km/s FWHM) emission lines, e.g., [Fe II], [Ca II], and Ca II. The absence of the broad [O I] 6300, 6364 line and weakness of [Fe II] and [Ca II] lines compared with Ca II IR triplet would be explained by a moderately high electron density in the line emitting region. This high density assumption seems to be consistent with the large amount of ejecta and low expansion velocity of SN 2006gy. The H-alpha line luminosity was as small as ~1x10^39 erg/s, being comparable with those of normal Type II SNe at similar epochs. Our observation indicates that the strong CSM interaction had almost finished by ~400 days. If the late time optical flux is purely powered by radioactive decay, at least M_Ni ~ 3 M_sun should be produced at the SN explosion. In the late phase spectrum, there were several unusual emission lines at 7400--8800 AA and some of them might be due to Ti or Ni synthesized at the explosion. (abridged)

preprint2009arXiv

Type Ib Supernova 2008D associated with the Luminous X-ray Transient 080109: An Energetic Explosion of a Massive Helium Star

We present a theoretical model for supernova (SN) 2008D associated with the luminous X-ray transient 080109. The bolometric light curve and optical spectra of the SN are modelled based on the progenitor models and the explosion models obtained from hydrodynamic/nucleosynthetic calculations. We find that SN 2008D is a more energetic explosion than normal core-collapse supernovae, with an ejecta mass of Mej = 5.3 +- 1.0 Msun and a kinetic energy of E = 6.0 +- 2.5 x 10^{51} erg. The progenitor star of the SN has a 6-8 Msun He core with essentially no H envelope (< 5 x 10^{-4} Msun) prior to the explosion. The main-sequence mass of the progenitor is estimated to be Mms =20-25 Msun, with additional systematic uncertainties due to convection, mass loss, rotation, and binary effects. These properties are intermediate between those of normal SNe and hypernovae associated with gamma-ray bursts. The mass of the central remnant is estimated as 1.6 - 1.8 Msun, which is near the boundary between neutron star and black hole formation.

preprint2008arXiv

Multi-Dimensional Simulations of Radiative Transfer in Aspherical Core-Collapse Supernovae

We study optical radiation of aspherical supernovae (SNe) and present an approach to verify the asphericity of SNe with optical observations of extragalactic SNe. For this purpose, we have developed a multi-dimensional Monte-Carlo radiative transfer code, SAMURAI (SupernovA MUlti-dimensional RAdIative transfer code). The code can compute the optical light curve and spectra both at early phases (<~ 40 days after the explosion) and late phases (~ 1 year after the explosion), based on hydrodynamic and nucleosynthetic models. We show that all the optical observations of SN 1998bw (associated with GRB 980425) are consistent with polar-viewed radiation of the aspherical explosion model with kinetic energy 20 x 10^{51} ergs. Properties of off-axis hypernovae are also discussed briefly.

preprint2008arXiv

Spectral luminosity indicators in SNe Ia - Understanding the R(SiII) line strength ratio and beyond

SNe Ia are good distance indicators because the shape of their light curves, which can be measured independently of distance, varies smoothly with luminosity. This suggests that SNe Ia are a single family of events. Similar correlations are observed between luminosity and spectral properties. In particular, the ratio of the strengths of the SiII λ5972 and λ6355 lines, known as R(SiII), was suggested as a potential luminosity indicator. Here, the physical reasons for the observed correlation are investigated. A Monte-Carlo code is used to construct a sequence of synthetic spectra resembling those of SNe with different luminosities near B maximum. The influence of abundances and of ionisation and excitation conditions on the synthetic spectral features is investigated. The ratio R(SiII) depends ssentially on the strength of SiII λ5972, because SiII λ6355 is saturated. In less luminous objects, SiII λ5972 is stronger because of a rapidly increasing SiII/SiIII ratio. Thus, the correlation between R(SiII) and luminosity is the effect of ionisation balance. The SiII λ5972 line itself may be the best spectroscopic luminosity indicator for SNe Ia, but all indicators discussed show scatter which may be related to abundance distributions.

preprint2008arXiv

Subaru and Keck Observations of the Peculiar Type Ia Supernova 2006gz at Late Phases

Recently, a few peculiar Type Ia supernovae (SNe) that show exceptionally large peak luminosity have been discovered. Their luminosity requires more than 1 Msun of 56Ni ejected during the explosion, suggesting that they might have originated from super-Chandrasekhar mass white dwarfs. However, the nature of these objects is not yet well understood. In particular, no data have been taken at late phases, about one year after the explosion. We report on Subaru and Keck optical spectroscopic and photometric observations of the SN Ia 2006gz, which had been classified as being one of these "overluminous" SNe Ia. The late-time behavior is distinctly different from that of normal SNe Ia, reinforcing the argument that SN 2006gz belongs to a different subclass than normal SNe Ia. However, the peculiar features found at late times are not readily connected to a large amount of 56Ni; the SN is faint, and it lacks [Fe II] and [Fe III] emission. If the bulk of the radioactive energy escapes the SN ejecta as visual light, as is the case in normal SNe Ia, the mass of 56Ni does not exceed ~ 0.3 Msun. We discuss several possibilities to remedy the problem. With the limited observations, however, we are unable to conclusively identify which process is responsible. An interesting possibility is that the bulk of the emission might be shifted to longer wavelengths, unlike the case in other SNe Ia, which might be related to dense C-rich regions as indicated by the early-phase data. Alternatively, it might be the case that SN 2006gz, though peculiar, was actually not substantially overluminous at early times.