Researcher profile

Masaomi Tanaka

Masaomi Tanaka contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
26works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

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