Source author record

Ken'ichi Nomoto

Ken'ichi Nomoto 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

97works
10topics
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

97 published item(s)

preprint2022arXiv

Hydrodynamic Simulations of Electron-capture Supernovae: Progenitor and Dimension Dependence

We present neutrino-transport hydrodynamic simulations of electron-capture supernovae (ECSNe) in \texttt{FLASH} with new two-dimensional (2D) collapsing progenitor models. These progenitor models feature the 2D modelling of oxygen-flame propagation until the onset of core collapse. We perform axisymmetric simulations with 6 progenitor models that, at the time of collapse, span a range of propagating flame front radii. For comparison, we also perform a simulation with the same setup using the canonical, spherically-symmetrical progenitor model n8.8. We found that the variations in the progenitor models inherited from simulations of stellar evolution and flame propagation do not significantly alter the global properties of the neutrino-driven ECSN explosion, such as the explosion energy ($\sim1.36$-$1.48\times10^{50}$ erg) and the mass ($\sim0.017$-$0.018M_\odot$) and composition of the ejecta. Due to aspherical perturbations induced by the 2D flame, the ejecta contains a small amount ($\lesssim1.8\times10^{-3}~M_\odot$) of low-$Y_e$ ($0.35<Y_e<0.4$) component. The baryonic mass of the protoneutron star is $\sim1.34~M_\odot$ ($\sim1.357~M_\odot$) with the new (n8.8) progenitor models when simulations end at $\sim400$ ms and the discrepancy is due to updated weak-interaction rates in the progenitor evolutionary simulations. Our results reflect the nature of ECSN progenitors containing a strongly degenerate ONeMg core and suggest a standardized ECSN explosion initialized by ONeMg core collapse. Moreover, we carry out a rudimentary three-dimensional simulation and find that the explosion properties are fairly compatible with the 2D counterpart. Our paper facilitates a more thorough understanding of ECSN explosions following the ONeMg core collapse, though more three-dimensional simulations are still needed.

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

The Late-Time Light Curves of Type Ia Supernovae: Confronting Models with Observations

Type Ia supernovae (SNe Ia) play a crucial role as standardizable candles in measurements of the Hubble constant and dark energy. Increasing evidence points towards multiple possible explosion channels as the origin of normal SNe Ia, with possible systematic effects on the determination of cosmological parameters. We present, for the first time, a comprehensive comparison of publicly-available SN Ia model nucleosynthetic data with observations of late-time light curve observations of SN Ia events. These models span a wide range of white dwarf (WD) progenitor masses, metallicities, explosion channels, and numerical methodologies. We focus on the influence of $^{57}$Ni and its isobaric decay product $^{57}$Co in powering the late-time ($t > 1000$ d) light curves of SNe Ia. $^{57}$Ni and $^{57}$Co are neutron-rich relative to the more abundant radioisotope $^{56}$Ni, and are consequently a sensitive probe of neutronization at the higher densities of near-Chandrashekhar (near-$M_{\rm Ch}$) progenitor WDs. We demonstrate that observations of one SN Ia event, SN 2015F is only consistent with a sub-$M_{\rm Ch}$ WD progenitor. Observations of four other events (SN 2011fe, SN 2012cg, SN 2014J, SN2013aa) are consistent with both near-$M_{\rm Ch}$ and sub-$M_{\rm Ch}$ progenitors. Continued observations of late-time light curves of nearby SNe Ia will provide crucial information on the nature of the SN Ia progenitors.

preprint2021arXiv

EMPRESS. IV. Extremely Metal-Poor Galaxies (EMPGs) Including Very Low-Mass Primordial Systems with M*=10^4--10^5 M_sun and 2--3% (O/H)_sun: High (Fe/O) Suggestive of Metal Enrichment by Hypernovae/Pair-Instability Supernovae

We present Keck/LRIS follow-up spectroscopy for 13 photometric candidates of extremely metal poor galaxies (EMPGs) selected by a machine-learning technique applied to the deep (~26 AB mag) optical and wide-area (~500 deg^2) Subaru imaging data in the EMPRESS survey. Nine out of the 13 candidates are EMPGs with an oxygen abundance (O/H) less than ~10% solar value (O/H)_sun, and four sources are contaminants of moderately metal-rich galaxies or no emission-line objects. Notably, two out of the nine EMPGs have extremely-low stellar masses and oxygen abundances of 5*10^4--7*10^5 M_sun and 2--3% (O/H)_sun, respectively. With a sample of five EMPGs with (Fe/O) measurements, two (three) of which are taken from this study (the literature), we confirm that two EMPGs with the lowest (O/H) ratios of ~2% (O/H)_sun show high (Fe/O) ratios of ~0.1, close to the solar abundance ratio. Comparing galaxy chemical enrichment models, we find that the two EMPGs cannot be explained by a scenario of metal-poor gas accretion/episodic star-formation history due to their low (N/O) ratios. We conclude that the two EMPGs can be reproduced by an inclusion of bright hypernovae and/or hypothetical pair-instability supernovae (SNe) preferentially produced in a metal-poor environment. This conclusion implies that primordial galaxies at z~10 could have a high abundance of Fe that is not originated from Type Ia SNe with delays, and that Fe may not serve as a cosmic clock for primordial galaxies.

preprint2021arXiv

Spatially Resolved X-ray Study of Supernova Remnant G306.3$-$0.9 with Unusually High Calcium Abundance

G306.3$-$0.9 is an asymmetric Galactic supernova remnant (SNR), whose progenitor has been thought to be a Type Ia supernova (SN), but its high Ca abundance appears inconsistent with the Type Ia origin. Hoping to uncover the reason for its asymmetry and the origin of this SNR, we performed a spatially resolved X-ray spectroscopic analysis of XMM-Newton and Chandra observation data. We divided the SNR into 13 regions and analyzed the spectra using two-temperature models (0.2 keV + 1 keV). Compared to the southwestern regions, the northeastern regions have higher metal abundances and a lower gas density. This suggests that the asymmetric morphology results from the non-uniform ambient environment. We found that neither Type Ia nor core-collapse SN models can account for the abnormally high abundance ratios of Ar/Si, Ca/Si, or the shape of the abundance curve. A comparison with the Ca-rich transient models shows that G306.3$-$0.9 is likely to be the first identified Galactic "Ca-rich transient" remnant, although the theoretical production of element S is lower. We also note that the conclusion for the SNR's origin relies on the measured abundance ratios and existing nucleosynthesis models. Between two groups of Ca-rich transient explosion models, we prefer the He shell detonation for an accreting WD, rather than the merger of a white dwarf and a neutron star.

preprint2020arXiv

A Minimum Dilution Scenario for Supernovae and Consequences for Extremely Metal-Poor Stars

To date no metal-free stars have been identified by direct observations. The most common method of constraining their properties is searching the spectra of the most metal-poor stars for the chemical elements created in the first stars and their supernova. In this approach, modelled supernova yields are compared to the observed abundance patterns in extremely metal-poor stars. The method typically only uses the abundance ratios, i.e., the yields are diluted to the observed level. Following the usual assumption of spherical symmetry we compute a simple lower limit of the mass a supernova can mix with and find that it is consistent with all published simulations of early chemical enrichment in the interstellar medium. For three different cases, we demonstrate that this dilution limit can change the conclusions from the abundance fitting. There is a large discrepancy between the dilution found in simulations of SN explosions in minihaloes and the dilution assumed in many abundance fits. Limiting the dilution can significantly alter the likelihood of which supernovae are possible progenitors of observed CEMP-no stars. In particular, some of the faint, very low-yield SNe, which have been suggested as models for the abundance pattern of SMSS0313-6708, cannot explain the measured metal abundances, as their predicted metal yields are too small by two orders of magnitude. Altogether, the new dilution model presented here emphasizes the need to better understand the mixing and dilution behaviour of aspherical SNe.

preprint2020arXiv

Chemical abundances in Sgr A East: evidence for a Type Iax supernova remnant

Recent observations have shown a remarkable diversity of observational behaviors and explosion mechanisms in thermonuclear supernovae (SNe). An emerging class of peculiar thermonuclear SNe, called Type Iax, show photometric and spectroscopic behaviors distinct from normal Type Ia. Their origin remains highly controversial, but pure turbulent deflagration of white dwarfs (WDs) has been regarded as the leading formation theory. The large population of Type Iax indicates the existence of unidentified Galactic Type Iax supernova remnants (SNRs). We report evidence that SNR Sgr A East in the Galactic center resulted from a pure turbulent deflagration of a Chandrasekhar-mass carbon-oxygen WD, an explosion mechanism used for Type Iax SNe. Our X-ray spectroscopic study of Sgr A East using 3 Ms of Chandra data shows a low ratio of intermediate-mass elements to Fe and large Mn/Fe and Ni/Fe ratios. This abundance pattern does not accord with the core-collapse or normal Type Ia models. Sgr A East is thus the first Galactic SNR for which a likely Type Iax origin has been proposed and the nearest target for studying this peculiar class. We compared Sgr A East with the Fe-rich SNRs 3C 397 and W49B, which also have high Mn and Cr abundances and were claimed to result from deflagration-to-detonation explosions of Chandrasekhar-mass WDs (although with disputes). Our study shows that they have distinct abundance patterns. The X-ray spectroscopic studies of thermonuclear SNRs provide observational evidence for the theories that there are diverse explosion channels and various metal outputs for Chandrasekhar-mass WDs.

preprint2020arXiv

Electron Capture Supernovae of Super-AGB Stars: Sensitivity on Input Physics

Stars of $\sim$ 8 - 10 $M_{\odot}$ on the main-sequence form strongly electron-degenerate O+Ne+Mg core and become super-AGB stars. If such an O+Ne+Mg core grows to 1.38 $M_\odot$, electron captures on $^{20}$Ne$(e,ν_e)^{20}$F$(e,ν_e)^{20}$O take place and ignite O-Ne deflagration around the center. In this paper, we perform two-dimensional hydrodynamics simulations of the propagation of the O-Ne flame to see whether such a flame induces a collapse of the O+Ne+Mg core due to subsequent electron capture behind the flame or triggers a thermonuclear explosion. We present a series of models to explore how the outcome depends on model parameters for the central density in the range from $10^{9.80}$ to $10^{10.20}$ g cm$^{-3}$, flame structure of both centered and off-centered ignition kernels, special and general relativistic effects, turbulent flame speed formula and the treatments of laminar burning phase. We find that the O+Ne+Mg core obtained from stellar evolutionary models has a high tendency to collapse into a neutron star. We obtain the bifurcation between the electron-capture collapse and thermonuclear explosion. We discuss the implication in nucleosynthesis and the possible observational signals of this class of supernovae.

preprint2020arXiv

Explosive Nucleosynthesis in Near-Chandrasekhar-mass White Dwarf Models for Type Iax Supernovae: Dependence on Model Parameters

The recently observed diversity of Type Ia supernovae (SNe Ia) has motivated us to conduct the theoretical modeling of SNe Ia for a wide parameter range. In particular, the origin of Type Iax supernovae (SNe Iax) has been obscure. Following our earlier work on the parameter dependence of SN Ia models, we focus on SNe Iax in the present study. For a model of SNe Iax, we adopt the currently leading model of pure turbulent deflagration (PTD) of near-Chandrasekhar mass C+O white dwarfs (WDs). We carry out 2-dimensional hydrodynamical simulations of the propagation of deflagration wave, which leaves a small WD remnant behind and eject nucleosynthesis materials. We show how the explosion properties, such as nucleosynthesis and explosion energy, depend on the model parameters such as central densities and compositions of the WDs (including the hybrid WDs), and turbulent flame prescription and initial flame geometry. We extract the associated observables in our models, and compare with the recently discovered low-mass WDs with unusual surface abundance patterns and the abundance patterns of some SN remnants. We provide the nucleosynthesis yield tables for applications to stellar archaeology and galactic chemical evolution. Our results are compared with the representative models in the literature.

preprint2020arXiv

Explosive Nucleosynthesis in Sub-Chandrasekhar Mass White Dwarf Models for Type Ia Supernovae: Dependence on Model Parameters

Recent observations of Type Ia supernovae (SNe Ia) have shown diversified properties of the explosion strength, light curves and chemical composition. To investigate possible origins of such diversities in SNe Ia, we have presented multi-dimensional hydrodynamical study of explosions and associated nucleosynthesis in the near Chandrasekhar mass carbon-oxygen (CO) white dwarfs (WDs) for a wide range of parameters (Leung and Nomoto 2018 ApJ). In the present paper, we extend our wide parameter survey of models to the explosions of sub-Chandrasekhar mass CO WDs. We take the double detonation model for the explosion mechanism. The model parameters of the survey include the metallicity of $Z = 0 - 5~Z_\odot$, the CO WD mass of $M = 0.90 - 1.20~M_\odot$, and the He envelope mass of $M_{\rm He} = 0.05 - 0.20~M_\odot$. We also study how the initial He detonation configuration, such as spherical, bubble, and ring shapes, triggers the C detonation. For these parameters, we derive the minimum He envelope mass necessary to trigger the C detonation. We then examine how the explosion dynamics and associated nucleosynthesis depend on these parameters, and compare our results with the previous representative models. We compare our nucleosynthesis yields with the unusual abundance patterns of Fe-peak elements and isotopes observed in SNe Ia 2011fe, 2012cg and 2014J, as well as SN Ia remnant 3C 397 to provide constraints on their progenitors and environments. We provide the nucleosynthesis yields table of the sub-Chandrasekhar mass explosions, to discuss their roles in the galactic chemical evolution and archaeology.

preprint2020arXiv

New Type Ia supernova Yields and the Manganese and Nickel Problems in the Milky Way and Dwarf Spheroidal Galaxies

In our quest to identify the progenitors of Type Ia supernovae (SNe Ia), we first update the nucleosynthesis yields both for near-Chandrasekhar (Ch) and sub-Ch mass white dwarfs (WDs), for a wide range of metallicity, with our two-dimensional hydrodynamical code and the latest nuclear reaction rates. We then include the yields in our galactic chemical evolution code to predict the evolution of elemental abundances in the solar neighborhood and dwarf spheroidal (dSph) galaxies: Fornax, Sculptor, Sextans, and Carina. In the observations of the solar neighborhood stars, Mn shows an opposite trend to $α$ elements, showing an increase toward higher metallicities, which is very well reproduced by deflagration-detonation transition of Ch-mass WDs, but never by double detonations of sub-Ch-mass WDs alone. The problem of Ch-mass SNe Ia was the Ni over-production at high metallicities. However, we found that Ni yields of Ch-mass SNe Ia are much lower with the solar-scaled initial composition than in previous works, which keeps the predicted Ni abundance within the observational scatter. From the evolutionary trends of elemental abundances in the solar neighborhood, we conclude that the contribution of sub-Ch-mass SNe Ia in chemical enrichment is up to 25\%. In dSph galaxies, however, larger enrichment from sub-Ch-mass SNe Ia than in the solar neighborhood may be required, which causes a decrease in [(Mg, Cr, Mn, Ni)/Fe] at lower metallicities. The observed high [Mn/Fe] ratios in Sculptor and Carina may also require additional enrichment from pure deflagrations, possibly as Type Iax supernovae. Future observations of dSph stars will provide more stringent constraints on the progenitor systems and explosion mechanism of SNe Ia.

preprint2020arXiv

Pulsational Pair-instability Supernovae. II. Neutrino Signals from Pulsations and their Detection by Terrestrial Neutrino Detectors

A Pulsational Pair-instability supernova (PPISN) evolves from a massive star with a mass $\sim 80$ -- 140 $M_{\odot}$ which develops the electron-positron pair-instability after the hydrostatic He-burning in the core has finished. In [Leung et al., ApJ 887, 72 (2019)] (Paper I) we examined the evolutionary tracks and the pulsational mass loss history of this class of stars. In this paper, we analyze the thermodynamical history to explore the neutrino observables of PPISNe. We compute the neutrino light curves and spectra during pulsation. We study the detailed neutrino emission profiles of these stars. Then, we estimate the expected neutrino detection count for different terrestrial neutrino detectors including, e.g., KamLAND and Super-Kamiokande. Finally, we compare the neutrino pattern of PPISN with other types of supernovae based on a canonical 10 kt detector. The predicted neutrino signals can provide the early warning for the telescopes to trace for the early time optical signals. Implications of neutrino physics on the expected detection are discussed.

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

Accretion-Induced Collapse of Dark Matter Admixed White Dwarfs -- I: Formation of Low-mass Neutron Stars

Recently observed pulsars with masses $\sim 1.1 ~M_{\odot}$ challenge the conventional neutron star (NS) formation path by core-collapse supernova (CCSN). Using spherically symmetric hydrodynamics simulations, we follow the collapse of a massive white dwarf (WD) core triggered by electron capture, until the formation of a proto-NS (PNS). For initial WD models with the same central density, we study the effects of a static, compact dark matter (DM) admixed core on the collapse and bounce dynamics and mass of the PNS, with DM mass $\sim 0.01 ~M_{\odot}$. We show that increasing the admixed DM mass generally leads to slower collapse and smaller PNS mass, down to about 1.0 $M_{\odot}$. Our results suggest that the accretion-induced collapse of dark matter admixed white dwarfs can produce low-mass neutron stars, such as the observed low-mass pulsar J0453+1559, which cannot be obtained by conventional NS formation path by CCSN.

preprint2019arXiv

Pulsational Pair-instability Supernovae. I. Pre-collapse Evolution and Pulsational Mass Ejection

We calculate the evolution of massive stars, which undergo pulsational pair-instability (PPI) when the O-rich core is formed. The evolution from the main-sequence through the onset of PPI is calculated for stars with the initial masses of $80 - 140$ $M_{\odot}$ and metallicities of $Z = 10^{-3} - 1.0$ $Z_\odot$. Because of mass loss, $Z \leq 0.5$ $Z_\odot$ is necessary for stars to form He cores massive enough (i.e., mass $>40 ~M_\odot$) to undergo PPI. The hydrodynamical phase of evolution from PPI through the beginning of Fe core collapse is calculated for the He cores with masses of $40 - 62 ~M_\odot$ and $Z = 0$. During PPI, electron-positron pair production causes a rapid contraction of the O-rich core which triggers explosive O-burning and a pulsation of the core. We study the mass dependence of the pulsation dynamics, thermodynamics, and nucleosynthesis. The pulsations are stronger for more massive He cores and result in such a large amount of mass ejection such as $3 - 13$ $M_\odot$ for $40 - 62 ~M_\odot$ He cores. These He cores eventually undergo Fe-core collapse. The $64 ~M_\odot$ He core undergoes complete disruption and becomes a pair-instability supernova. The H-free circumstellar matter ejected around these He cores is massive enough for to explain the observed light curve of Type I (H-free) superluminous supernovae with circumstellar interaction. We also note that the mass ejection sets the maximum mass of black holes (BHs) to be $\sim 50$ $M_{\odot}$, which is consistent with the masses of BHs recently detected by VIRGO and aLIGO.

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

Disappearance of the Progenitor of Supernova iPTF13bvn

Supernova (SN) iPTF13bvn in NGC 5806 was the first Type Ib SN to have been tentatively associated with a progenitor candidate in pre-explosion images. We performed deep ultraviolet (UV) and optical Hubble Space Telescope (HST) observations of the SN site 740 days after explosion. We detect an object in the optical bands that is fainter than the pre-explosion object. This dimming is likely not produced by dust absorption in the ejecta; thus, our finding confirms the connection of the progenitor candidate with the SN. The object in our data is likely dominated by the fading SN, which implies that the pre-SN flux is mostly due to the progenitor. We compare our revised pre-SN photometry with previously proposed progenitor models. Although binary progenitors are favored, models need to be refined. In particular, to comply with our deep UV detection limit, any companion star must be less luminous than a late-O star or substantially obscured by newly formed dust. A definitive progenitor characterization will require further observations to disentangle the contribution of a much fainter SN and its environment.

preprint2016arXiv

Do Stellar Winds Prevent the Formation of Supermassive Stars by Accretion?

Supermassive stars (SMS; ~ 10^5 M_sun) formed from metal-free gas in the early Universe attract attention as progenitors of supermassive black holes observed at high redshifts. To form SMSs by accretion, central protostars must accrete at as high rates as ~ 0.1-1 M_sun/yr. Such protostars have very extended structures with bloated envelopes, like super-giant stars, and are called super-giant protostars (SGPSs). Under the assumption of hydrostatic equilibrium, SGPSs have density inverted layers, where the luminosity becomes locally super-Eddington, near the surface. If the envelope matter is allowed to flow out, however, a stellar wind could be launched and hinder the accretion growth of SGPSs before reaching the supermassive regime. We examine whether radiation-driven winds are launched from SGPSs by constructing steady and spherically symmetric wind solutions. We find that the wind velocity does not reach the escape velocity in any case considered. This is because once the temperature falls below ~ 10^4 K, the opacity plummet drastically owing to the recombination of hydrogen and the acceleration ceases suddenly. This indicates that, in realistic non-steady cases, even if outflows are launched from the surface of SGPSs, they would fall back again. Such a "wind" does not result in net mass loss and does not prevent the growth of SGPSs. In conclusion, SGPSs will grow to SMSs and eventually collapse to massive BHs of ~ 10^5 M_sun, as long as the rapid accretion is maintained.

preprint2016arXiv

Impact of New Gamow-Teller Strengths on Explosive Type Ia Supernova Nucleosynthesis

Recent experimental results have confirmed a possible reduction in the GT$_+$ strengths of pf-shell nuclei. These proton-rich nuclei are of relevance in the deflagration and explosive burning phases of Type Ia supernovae. While prior GT strengths result in nucleosynthesis predictions with a lower-than-expected electron fraction, a reduction in the GT$_+$ strength can result in an slightly increased electron fraction compared to previous shell model predictions, though the enhancement is not as large as previous enhancements in going from rates computed by Fuller, Fowler, and Newman based on an independent particle model. A shell model parametrization has been developed which more closely matches experimental GT strengths. The resultant electron-capture rates are used in nucleosynthesis calculations for carbon deflagration and explosion phases of Type Ia supernovae, and the final mass fractions are compared to those obtained using more commonly-used rates.

preprint2016arXiv

Multicolor light curves simulations of Population III core-collapse supernovae: from shock breakout to $^{56}$Co decay

The properties of the first generation of stars and their supernova (SN) explosions remains unknown due to the lack of their actual observations. Recently many transient surveys are conducted and the feasibility of the detection of supernovae (SNe) of Pop III stars is growing. In this paper we study the multicolor light curves for a number of metal-free core-collapse SN models (25-100 M$_{\odot}$) to provide the indicators for finding and identification of first generation SNe. We use mixing-fallback supernova explosion models which explain the observed abundance patterns of metal poor stars. Numerical calculations of the multicolor light curves are performed using multigroup radiation hydrodynamic code STELLA. The calculated light curves of metal-free SNe are compared with non-zero metallicity models and several observed SNe. We have found that the shock breakout characteristics, the evolution of the photosphere's velocity, the luminosity, the duration and color evolution of the plateau - all the SN phases are helpful to estimate the parameters of SN progenitor: the mass, the radius, the explosion energy and the metallicity. We conclude that the multicolor light curves can be potentially used to identify first generation SNe in the current (Subaru/HSC) and future transient surveys (LSST, JWST). They are also suitable for identification of the low-metallicity SNe in the nearby Universe (PTF, Pan-STARRS, Gaia).

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

Shock wave structure in astrophysical flows with an account of photon transfer

For an accurate treatment of the shock wave propagation in high-energy astrophysical phenomena, such as supernova shock breakouts, gamma-ray bursts and accretion disks, knowledge of radiative transfer plays a crucial role. In this paper we consider one-dimensional (1D) special relativistic radiation hydrodynamics by solving the Boltzmann equation for radiative transfer. The structure of a radiative shock is calculated for a number of shock tube problems, including strong shock waves, and relativistic- and radiation-dominated cases. Calculations are performed using an iterative technique that consistently solves the equations of relativistic hydrodynamics and relativistic comoving radiative transfer. A comparison of radiative transfer solutions with the Eddington approximation and the M1 closure is made. A qualitative analysis of moment equations for radiation is performed and the conditions for the existence of jump discontinuity for non-relativistic cases are investigated numerically.

preprint2016arXiv

The critical mass ratio of double white dwarf binaries for violent merger-induced Type Ia supernova explosions

Mergers of carbon-oxygen (CO) white dwarfs (WDs) are considered as one of the potential progenitors of type Ia supernovae (SNe Ia). Recent hydrodynamical simulations showed that the less massive (secondary) WD violently accretes onto the more massive (primary) one, carbon detonation occurs, the detonation wave propagates through the primary, and the primary finally explodes as a sub-Chandrasekhar mass SN Ia. Such an explosion mechanism is called the violent merger scenario. Based on the smoothed particle hydrodynamics (SPH) simulations of merging CO WDs, we derived more stringent critical mass ratio (qcr) leading to the violent merger scenario than the previous results. We conclude that this difference mainly comes from the differences in the initial condition, synchronously spinning of WDs or not. Using our new results, we estimated the brightness distribution of SNe Ia in the violent merger scenario and compared it with previous studies. We found that our new qcr does not significantly affect the brightness distribution. We present the direct outcome immediately following CO WD mergers for various primary masses and mass ratios. We also discussed the final fate of the central system of the bipolar planetary nebula Henize 2-428, which was recently suggested to be a double CO WD system whose total mass exceeds the Chandrasekhar-limiting mass, merging within the Hubble time. Even considering the uncertainties in the proposed binary parameters, we concluded that the final fate of this system is almost certainly a sub-Chandrasekhar mass SN Ia in the violent merger scenario.

preprint2016arXiv

The final fate of stars that ignite neon and oxygen off-center: electron capture or iron core-collapse supernova?

In the ONeMg cores of $8.8-9.5~{\rm M}_\odot$ stars, neon and oxygen burning is ignited off-center. Whether the neon-oxygen flame propagates to the center is critical to determine whether these stars undergo Fe core collapse or electron capture induced ONeMg core collapse. We present more details of stars that ignite neon and oxygen burning off-center. The neon flame is established in a similar manner to the carbon flame of super-AGB stars, albeit with a narrower flame width. The criteria for establishing a flame are able to be met if the strict Schwarzschild criterion for convective instability is adopted. Mixing across the interface of the convective shell disrupts the conditions for the propagation of the burning front and instead the shell burns as a series of inward-moving flashes. While this may not directly affect whether the burning will reach the center (as in super-AGB stars), the core is allowed to contract between each shell flash. Reduction of the electron fraction in the shell reduces the Chandrasekhar mass and the center reaches the threshold density for the URCA process to activate and steer the remaining evolution of the core. This highlights the importance of a more accurate treatment of mixing in the stellar interior for yet another important question in stellar astrophysics - determining the properties of stellar evolution and supernova progenitors at the boundary between electron capture supernova and iron core-collapse supernova.

preprint2016arXiv

The Unusual Super-Luminous Supernovae SN 2011kl and ASASSN-15lh

Two recently discovered very luminous supernovae (SNe) present stimulating cases to explore the extents of the available theoretical models. SN 2011kl represents the first detection of a supernova explosion associated with an ultra-long duration gamma ray burst. ASASSN-15lh was even claimed as the most luminous SN ever discovered, challenging the scenarios so far proposed for stellar explosions. Here we use our radiation hydrodynamics code in order to simulate magnetar powered SNe. To avoid explicitly assuming neutron star properties we adopt the magnetar luminosity and spin-down timescale as free parameters of the model. We find that the light curve (LC) of SN 2011kl is consistent with a magnetar power source, as previously proposed, but we note that some amount of 56^Ni (> 0.08 M_sun) is necessary to explain the low contrast between the LC peak and tail. For the case of ASASSN-15lh we find physically plausible magnetar parameters that reproduce the overall shape of the LC provided the progenitor mass is relatively large (a mass of the ejecta approx 6 M_sun). The ejecta hydrodynamics of this event is dominated by the magnetar input, while the effect is more moderate for SN 2011kl. We conclude that a magnetar model may be used for the interpretation of these events and that the hydrodynamic modeling is necessary to derive the properties of powerful magnetars and their progenitors.

preprint2016arXiv

Type I Superluminous Supernovae as Explosions inside Non-Hydrogen Circumstellar Envelopes

A number of Type I (hydrogenless) superluminous supernova (SLSN) events have been discovered recently. However, their nature remains debatable. One of the most promising ideas is the shock-interaction mechanism, but only simplified semi-analytical models have been applied so far. We simulate light curves for several Type I SLSN (SLSN-I) models enshrouded by dense, non-hydrogen circumstellar envelopes, using a multi-group radiation hydrodynamics code that predicts not only bolometric, but also multicolor light curves. We demonstrate that the bulk of SLSNe-I including those with relatively narrow light curves like SN 2010gx or broad ones like PTF09cnd can be explained by the interaction of the SN ejecta with he CS envelope, though the range of parameters for these models is rather wide. Moderate explosion energy ($\sim (2 - 4)\cdot 10^{51}$ ergs) is sufficient to explain both narrow and broad SLSN-I light curves, but ejected mass and envelope mass differ for those two cases. Only 5 to 10 $M_\odot$ of non-hydrogen material is needed to reproduce the light curve of SN 2010gx, while the best model for PTF09cnd is very massive: it contains almost $ 50 M_\odot $ in the CS envelope and only $ 5 M_\odot $ in the ejecta. The CS envelope for each case extends from 10 $R_\odot$ to $\sim 10^5R_\odot$ ($7\cdot 10^{15} $ cm), which is about an order of magnitude larger than typical photospheric radii of standard SNe near the maximum light. We briefly discuss possible ways to form such unusual envelopes.

preprint2015arXiv

A systematic study of carbon-oxygen white dwarf mergers: mass combinations for Type Ia supernovae

Mergers of two carbon-oxygen (CO) white dwarfs (WDs) have been considered as progenitors of Type Ia supernovae (SNe Ia). Based on smoothed particle hydrodynamics (SPH) simulations, previous studies claimed that mergers of CO WDs lead to an SN Ia explosion either in the dynamical merger phase or stationary rotating merger remnant phase. However, the mass range of CO WDs that lead to an SN Ia has not been clearly identified yet. In the present work, we perform systematic SPH merger simulations for the WD masses ranging from $0.5~M_{\odot}$ to $1.1~M_{\odot}$ with higher resolutions than the previous systematic surveys and examine whether or not carbon burning occurs dynamically or quiescently in each phase. We further study the possibility of SN Ia explosion and estimate the mass range of CO WDs that lead to an SN Ia. We found that when the both WDs are massive, i.e., in the mass range of $0.9~M_{\odot} {\le} M_{1,2} {\le} 1.1~M_{\odot}$, they can explode as an SN Ia in the merger phase. On the other hand, when the more massive WD is in the range of $0.7~M_{\odot} {\le} M_{1} {\le} 0.9~M_{\odot}$ and the total mass exceeds $1.38~M_{\odot}$, they can finally explode in the stationary rotating merger remnant phase. We estimate the contribution of CO WD mergers to the entire SN Ia rate in our galaxy to be of ${\lt} 9\%$. So, it might be difficult to explain all galactic SNe Ia by CO WD mergers.

preprint2015arXiv

Electron capture and beta-decay rates for sd-shell nuclei in stellar environments relevant to high density O-Ne-Mg cores

Electron capture and beta-decay rates for nuclear pairs in sd-shell are evaluated at high densities and high temperatures relevant to the final evolution of electron-degenerate O-Ne-Mg cores of stars with the initial masses of 8-10 solar mass. Electron capture induces a rapid contraction of the electron-degenerate O-Ne-Mg core. The outcome of rapid contraction depends on the evolutionary changes in the central density and temperature, which are determined by the competing processes of contraction, cooling, and heating. The fate of the stars are determined by these competitions, whether they end up with electron-capture supernovae or Fe core-collapse supernovae. Since the competing processes are induced by electron capture and beta-decay, the accurate weak rates are crucially important. The rates are obtained for pairs with A=20, 23, 24, 25 and 27 by shell-model calculations in sd-shell with the USDB Hamiltonian. Effects of Coulomb corrections on the rates are evaluated. The rates for pairs with A=23 and 25 are important for nuclear URCA processes that determine the cooling rate of O-Ne-Mg core, while those for pairs with A=20 and 24 are important for the core-contraction and heat generation rates in the core. We provide these nuclear rates at stellar environments in tables with fine enough meshes at various densities and temperatures for the studies of astrophysical processes sensitive to the rates. In particular, the accurate rate tables are crucially important for the final fates of not only O-Ne-Mg cores but also a wider range of stars such as C-O cores of lower mass stars.

preprint2015arXiv

Final Evolution and Delayed Explosions of Spinning White Dwarfs in Single Degenerate Models for Type Ia Supernovae

We study the occurrence of delayed SNe~Ia in the single degenerate (SD) scenario. We assume that a massive carbon-oxygen (CO) white dwarf (WD) accretes matter coming from a companion star, making it to spin at the critical rate. We assume uniform rotation due to magnetic field coupling. The carbon ignition mass for non-rotating WDs is M_{ig}^{NR} \approx 1.38 M_{\odot}; while for the case of uniformly rotating WDs it is a few percent larger (M_{ig}^{R} \approx 1.43 M_{\odot}). When accretion rate decreases, the WD begins to lose angular momentum, shrinks, and spins up; however, it does not overflow its critical rotation rate, avoiding mass shedding. Thus, angular momentum losses can lead the CO WD interior to compression and carbon ignition, which would induce an SN~Ia. The delay, largely due to the angular momentum losses timescale, may be large enough to allow the companion star to evolve to a He WD, becoming undetectable at the moment of explosion. This scenario supports the occurrence of delayed SNe~Ia if the final CO WD mass is 1.38 M_{\odot} < M < 1.43 M_{\odot}. We also find that if the delay is longer than ~3 Gyr, the WD would become too cold to explode, rather undergoing collapse.

preprint2015arXiv

Hydrodynamical evolution of merging carbon-oxygen white dwarfs: their pre-supernova structure and observational counterparts

We perform smoothed particle hydrodynamics (SPH) simulations for merging binary carbon-oxygen (CO) white dwarfs (WDs) with masses of $1.1$ and $1.0$ $M_\odot$, until the merger remnant reaches a dynamically steady state. Using these results, we assess whether the binary could induce a thermonuclear explosion, and whether the explosion could be observed as a type Ia supernova (SN Ia). We investigate three explosion mechanisms: a helium-ignition following the dynamical merger (`helium-ignited violent merger model'), a carbon-ignition (`carbon-ignited violent merger model'), and an explosion following the formation of the Chandrasekhar mass WD (`Chandrasekhar mass model'). An explosion of the helium-ignited violent merger model is possible, while we predict that the resulting SN ejecta are highly asymmetric since its companion star is fully intact at the time of the explosion. The carbon-ignited violent merger model can also lead to an explosion. However, the envelope of the exploding WD spreads out to $\sim 0.1R_\odot$; it is much larger than that inferred for SN 2011fe ($< 0.1R_\odot $) while much smaller than that for SN 2014J ($\sim 1R_\odot$). For the particular combination of the WD masses studied in this work, the Chandrasekhar mass model is not successful to lead to an SN Ia explosion. Besides these assessments, we investigate the evolution of unbound materials ejected through the merging process (`merger ejecta'), assuming a case where the SN Ia explosion is not triggered by the helium- or carbon-ignition during the merger. The merger ejecta interact with the surrounding interstellar medium, and form a shell. The shell has a bolometric luminosity of more than $2 \times 10^{35}$ ergs$^{-1}$ lasting for $\sim 2 \times 10^4$ yr. If this is the case, Milky Way should harbor about $10$ such shells at any given time.

preprint2015arXiv

Subclasses of Type Ia Supernovae as the origin of [α/Fe] ratios in dwarf spheroidal galaxies

Recent extensive observations of Type Ia Supernovae (SNe Ia) have revealed the existence of a diversity of SNe Ia, including SNe Iax. We introduce two possible channels in the single degenerate scenario: 1) double detonations in sub-Chandrasekhar (Ch) mass CO white dwarfs (WDs), where a thin He envelope is developed with relatively low accretion rates after He novae even at low metallicities, and 2) carbon deflagrations in Ch-mass possibly hybrid C+O+Ne WDs, where WD winds occur at [Fe/H] ~ -2.5 at high accretion rates. These subclasses of SNe Ia are rarer than `normal' SNe Ia and do not affect the chemical evolution in the solar neighborhood, but can be very important in metal-poor systems with stochastic star formation. In dwarf spheroidal galaxies in the Local Group, the decrease of [α/Fe] ratios at [Fe/H] ~ -2 to -1.5 can be produced depending on the star formation history. SNe Iax give high [Mn/Fe], while sub-Ch-mass SNe Ia give low [Mn/Fe], and thus a model including a mix of the two is favoured by the available observations.

preprint2015arXiv

The Progenitor of the Type IIb SN 2008ax Revisited

Hubble Space Telescope observations of the site of the supernova (SN) 2008ax obtained in 2011 and 2013 reveal that the possible progenitor object detected in pre-explosion images was in fact multiple. Four point sources are resolved in the new, higher-resolution images. We identify one of the sources with the fading SN. The other three objects are consistent with single supergiant stars. We conclude that their light contaminated the previously identified progenitor candidate. After subtraction of these stars, the progenitor appears to be significantly fainter and bluer than previously measured. Post-explosion photometry at the SN location indicates that the progenitor object has disappeared. If single, the progenitor is compatible with a supergiant star of B to mid-A spectral type, while a Wolf-Rayet (WR) star would be too luminous in the ultraviolet to account for the observations. Moreover, our hydrodynamical modelling shows the pre-explosion mass was $4-5$ $M_\odot$ and the radius was $30-50$ $R_\odot$, which is incompatible with a WR progenitor. We present a possible interacting binary progenitor computed with our evolutionary models that reproduces all the observational evidence. A companion star as luminous as an O9-B0 main-sequence star may have remained after the explosion.

preprint2014arXiv

A Blue Point Source at the Location of Supernova 2011dh

We present Hubble Space Telescope (HST) observations of the field of the Type IIb supernova (SN) 2011dh in M51 performed at ~1161 rest-frame days after explosion using the Wide Field Camera 3 and near-UV filters F225W and F336W. A star-like object is detected in both bands and the photometry indicates it has negative (F225W - F336W) color. The observed object is compatible with the companion of the now-vanished yellow supergiant progenitor predicted in interacting binary models. We consider it unlikely that the SN is undergoing strong interaction and thus estimate that it makes a small contribution to the observed flux. The possibilities of having detected an unresolved light echo or an unrelated object are briefly discussed and judged unlikely. Adopting a possible range of extinction by dust, we constrain parameters of the proposed binary system. In particular, the efficiency of mass accretion onto the binary companion must be below 50%, if no significant extinction is produced by newly formed dust. Further multiband observations are required in order to confirm the identification of the object as the companion star. If confirmed, the companion star would already be dominant in the UV/optical regime, so it would readily provide a unique opportunity to perform a detailed study of its properties.

preprint2014arXiv

Abundance Profiling of Extremely Metal-Poor Stars and Supernova Properties in the Early Universe

After the Big Bang nucleosynthesis, the first heavy element enrichment in the Universe was made by a supernova (SN) explosion of a population (Pop) III star (Pop III SN). The abundance ratios of elements produced from Pop III SNe are recorded in abundance patterns of extremely metal-poor (EMP) stars. The observations of the increasing number of EMP stars have made it possible to statistically constrain the explosion properties of Pop III SNe. We present Pop III SN models whose nucleosynthesis yields well-reproduce individually the abundance patterns of 48 such metal-poor stars as [Fe/H] $\mathrel{\rlap{\lower 4pt \hbox{$\sim$}}\raise 1pt \hbox {$<$}}-3.5$. We then derive relations between the abundance ratios of EMP stars and certain explosion properties of Pop III SNe: the higher [(C+N)/Fe] and [(C+N)/Mg] ratios correspond to the smaller ejected Fe mass and the larger compact remnant mass, respectively. Using these relations, the distributions of the abundance ratios of EMP stars are converted to those of the explosion properties of Pop III SNe. Such distributions are compared with those of the explosion properties of present day SNe: The distribution of the ejected Fe mass of Pop III SNe has the same peak as that of the resent day SNe but shows an extended tail down to $\sim10^{-2}-10^{-5}M_\odot$, and the distribution of the mass of the compact remnant of Pop III SNe is as wide as that of the present day stellar-mass black holes. Our results demonstrate the importance of large samples of EMP stars obtained by ongoing and future EMP star surveys and subsequent high-dispersion spectroscopic observations in clarifying the nature of Pop III SNe in the early Universe.

preprint2014arXiv

Abundance stratification in Type Ia Supernovae - IV: the luminous, peculiar SN 1991T

The abundance distribution of the elements in the ejecta of the peculiar, luminous Type Ia supernova (SN Ia) 1991T is obtained modelling spectra from before maximum light until a year after the explosion, with the method of "Abundance Tomography". SN 1991T is different from other slowly declining SNe Ia (e.g. SN 1999ee) in having a weaker Si II 6355 line and strong features of iron group elements before maximum. The distance to the SN is investigated along with the abundances and the density profile. The ionization transition that happens around maximum sets a strict upper limit on the luminosity. Both W7 and the WDD3 delayed detonation model are tested. WDD3 is found to provide marginally better fits. In this model the core of the ejecta is dominated by stable Fe with a mass of about 0.15 solar masses, as in most SNe Ia. The layer above is mainly 56Ni up to v~10000 km/s (~0.78 solar masses). A significant amount of 56Ni (~3 %) is located in the outer layers. A narrow layer between 10000 km/s and ~12000 km/s is dominated by intermediate mass elements (IME), ~0.18 solar masses. This is small for a SN Ia. The high luminosity and the consequently high ionization, and the high 56Ni abundance at high velocities explain the peculiar early-time spectra of SN 1991T. The outer part is mainly of oxygen, ~0.3 solar masses. Carbon lines are never detected, yielding an upper limit of 0.01 solar masses for C. The abundances obtained with the W7 density model are qualitatively similar to those of the WDD3 model. Different elements are stratified with moderate mixing, resembling a delayed detonation.

preprint2014arXiv

Detection of the Gravitational Lens Magnifying a Type Ia Supernova

Objects of known brightness, like Type Ia supernovae (SNIa), can be used to measure distances. If a massive object warps spacetime to form multiple images of a background SNIa, a direct test of cosmic expansion is also possible. However, these lensing events must first be distinguished from other rare phenomena. Recently, a supernova was found to shine much brighter than normal for its distance, which resulted in a debate: was it a new type of superluminous supernova or a normal SNIa magnified by a hidden gravitational lens? Here we report that a spectrum obtained after the supernova faded away shows the presence of a foreground galaxy--the first found to strongly magnify a SNIa. We discuss how more lensed SNIa may be found than previously predicted.

preprint2014arXiv

Dust Production Factories in the Early Universe: Formation of Carbon Grains in Red-supergiant Winds of Very Massive Population III Stars

We investigate the formation of dust in a stellar wind during the red-supergiant (RSG) phase of a very massive Population III star with the zero-age main sequence mass of 500 M_sun. We show that, in a carbon-rich wind with a constant velocity, carbon grains can form with a lognormal-like size distribution, and that all of the carbon available for dust formation finally condense into dust for wide ranges of the mass-loss rate ((0.1-3)x10^{-3} M_sun yr^{-1}) and wind velocity (1-100 km s^{-1}). We also find that the acceleration of the wind driven by newly formed dust suppresses the grain growth but still allows more than half of gas-phase carbon to be finally locked up in dust grains. These results indicate that at most 1.7 M_sun of carbon grains can form in total during the RSG phase of 500 M_sun Population III stars. Such a high dust yield could place very massive primordial stars as important sources of dust at the very early epoch of the universe if the initial mass function of Population III stars was top-heavy. We also briefly discuss a new formation scenario of carbon-rich ultra-metal-poor stars considering the feedback from very massive Population III stars.

preprint2014arXiv

Electron-capture supernovae exploding within their progenitor wind

The most massive stars on the asymptotic giant branch (AGB), so called super-AGB stars, are thought to produce supernovae (SNe) triggered by electron captures in their degenerate O+Ne+Mg cores. Super-AGB stars are expected to have slow winds with high mass-loss rates, so their wind density is high. The explosions of super-AGB stars are therefore presumed to occur in this dense wind. We provide the first synthetic light curves (LCs) for such events by exploding realistic electron-capture supernova (ecSN) progenitors within their super-AGB winds. We find that the early LC, i.e. before the recombination wave reaches the bottom of the H-rich envelope of SN ejecta (the plateau phase), is not affected by the dense wind. However, after the plateau phase, the luminosity remains higher when the super-AGB wind is taken into account. We compare our results to the historical LC of SN 1054, the progenitor of the Crab Nebula, and show that the explosion of an ecSN within an ordinary super-AGB wind can explain the LC features. We conclude that SN 1054 could have been a Type IIn SN without any extra extreme mass loss which was previously suggested to be necessary to account for its early high luminosity. We also show that our LCs match Type IIn SNe with an early plateau phase (`Type IIn-P') and suggest that they are ecSNe within super-AGB winds. Although some ecSNe can be bright in the optical spectral range due to the large progenitor radius, their X-ray luminosity from the interaction does not necessarily get as bright as other Type IIn SNe whose optical luminosities are also powered by the interaction. Thus, we suggest that optically-bright X-ray-faint Type IIn SNe can emerge from ecSNe. Optically-faint Type IIn SNe, such as SN 2008S, can also originate from ecSNe if their H-rich envelope masses are small. Some of them can be observed as `Type IIn-b' SNe due to the small H-rich envelope mass.

preprint2014arXiv

Faint Population III supernovae as the origin of the most iron-poor stars

The most iron-deficient stars in the Milky Way provide important observational constraints on the nature of astrophysical objects that have enriched the primordial gas with heavy elements from which these stars were formed. Among them, the recently discovered iron-deficient star SMSS J031300.36-670839.3 shows a remarkable chemical composition with non-detection of iron ([Fe/H]$<-7.1$) and large enhancement of carbon and magnesium relative to calcium. We investigate the supernova yields of metal-free (Population III) stars to interpret the observed abundance pattern for this star. We report that the high [C/Ca] and [C/Mg] ratios and upper limits determined for other elemental abundances are well reproduced with the yields of core-collapse supernovae (that have normal kinetic energies of explosion $E$ of $E_{51}=E/10^{51}$erg$=1$) or hypernovae ($E_{51}\geq 10$) of the Population III 25$M_{\odot}$ or 40$M_{\odot}$ stars. The best-fit model assumes that the explosion of the Population III progenitor undergoes extensive matter mixing and fallback, leaving behind a black hole remnant. In these models, Ca is produced by static/explosive O burning and incomplete Si burning in the Population III supernova/hypernova, in contrast to the suggestion that Ca is originated from the hot CNO cycle during the presupernova evolution. Among the five most iron-poor stars with [Fe/H]$<-4.5$, four carbon-enhanced stars are consistent with the faint supernova models with the ejected mass of $^{56}$Ni less than 10$^{-3}M_{\odot}$.

preprint2014arXiv

Hybrid C-O-Ne White Dwarfs as Progenitors of Diverse SNe Ia

When carbon is ignited off-center in a CO core of a super-AGB star, its burning in a convective shell tends to propagate to the center. Whether the C flame will actually be able to reach the center depends on the efficiency of extra mixing beneath the C convective shell. Whereas thermohaline mixing is too inefficient to interfere with the C-flame propagation, convective boundary mixing can prevent the C burning from reaching the center. As a result, a C-O-Ne white dwarf (WD) is formed, after the star has lost its envelope. Such a "hybrid" WD has a small CO core surrounded by a thick ONe zone. In our 1D stellar evolution computations the hybrid WD is allowed to accrete C-rich material, as if it were in a close binary system and accreted H-rich material from its companion with a sufficiently high rate at which the accreted H would be processed into He under stationary conditions, assuming that He could then be transformed into C. When the mass of the accreting WD approaches the Chandrasekhar limit, we find a series of convective Urca shell flashes associated with high abundances of 23Na and 25Mg. They are followed by off-center C ignition leading to convection that occupies almost the entire star. To model the Urca processes, we use the most recent well-resolved data for their reaction and neutrino-energy loss rates. Because of the emphasized uncertainty of the convective Urca process in our hybrid WD models of SN Ia progenitors, we consider a number of their potentially possible alternative instances for different mixing assumptions, all of which reach a phase of explosive C ignition, either off or in the center. Our hybrid SN Ia progenitor models have much lower C to O abundance ratios at the moment of the explosive C ignition than their pure CO counterparts, which may explain the observed diversity of the SNe Ia.

preprint2014arXiv

iPTF13bvn: The First Evidence of a Binary Progenitor for a Type Ib Supernova

The recent detection in archival HST images of an object at the the location of supernova (SN) iPTF13bvn may represent the first direct evidence of the progenitor of a Type Ib SN. The object's photometry was found to be compatible with a Wolf-Rayet pre-SN star mass of ~11 Msun. However, based on hydrodynamical models we show that the progenitor had a pre-SN mass of ~3.5 Msun and that it could not be larger than ~8 Msun. We propose an interacting binary system as the SN progenitor and perform evolutionary calculations that are able to self-consistently explain the light-curve shape, the absence of hydrogen, and the pre-SN photometry. We further discuss the range of allowed binary systems and predict that the remaining companion is a luminous O-type star of significantly lower flux in the optical than the pre-SN object. A future detection of such star may be possible and would provide the first robust identification of a progenitor system for a Type Ib SN.

preprint2014arXiv

Shortest Recurrence Periods of Novae

Stimulated by the recent discovery of the 1 yr recurrence period nova M31N 2008-12a, we examined the shortest recurrence periods of hydrogen shell flashes on mass-accreting white dwarfs (WDs). We discuss the mechanism that yields a finite minimum recurrence period for a given WD mass. Calculating the unstable flashes for various WD masses and mass accretion rates, we identified a shortest recurrence period of about two months for a non-rotating 1.38 M_sun WD with a mass accretion rate of 3.6 x 10^{-7} M_sun yr^{-1}. A 1 yr recurrence period is realized for very massive (> 1.3 M_sun) WDs with very high accretion rates (> 1.5 x 10^{-7} M_sun yr^{-1}). We revised our stability limit of hydrogen shell burning, which will be useful for binary evolution calculations toward Type Ia supernovae.

preprint2014arXiv

Supernova 2010as: the Lowest-Velocity Member of a Family of Flat-Velocity Type IIb Supernovae

We present extensive optical and near-infrared photometric and spectroscopic observations of the stripped-envelope (SE) supernova SN 2010as. Spectroscopic peculiarities, such as initially weak helium features and low expansion velocities with a nearly flat evolution, place this object in the small family of events previously identified as transitional Type Ib/c supernovae (SNe). There is ubiquitous evidence of hydrogen, albeit weak, in this family of SNe, indicating that they are in fact a peculiar kind of Type IIb SNe that we name "flat-velocity Type IIb". The flat velocity evolution---which occurs at different levels between 6000 and 8000 km/s for different SNe---suggests the presence of a dense shell in the ejecta. Despite the spectroscopic similarities, these objects show surprisingly diverse luminosities. We discuss the possible physical or geometrical unification picture for such diversity. Using archival HST images we associate SN 2010as with a massive cluster and derive a progenitor age of ~6 Myr, assuming a single star-formation burst, which is compatible with a Wolf-Rayet progenitor. Our hydrodynamical modelling, on the contrary, indicates the pre-explosion mass was relatively low, of ~4 M_sol. The seeming contradiction between an young age and low pre-SN mass may be solved by a massive interacting binary progenitor.

preprint2014arXiv

The origin of low [alpha/Fe] ratios in extremely metal-poor stars

We show that the low ratios of $α$ elements (Mg, Si, and Ca) to Fe recently found for a small fraction of extremely metal-poor stars can be naturally explained with the nucleosynthesis yields of core-collapse supernovae, i.e., $13-25M_\odot$ supernovae, or hypernovae. For the case without carbon enhancement, the ejected iron mass is normal, consistent with observed light curves and spectra of nearby supernovae. On the other hand, the carbon enhancement requires much smaller iron production, and the low [$α$/Fe] of carbon enhanced metal-poor stars can also be reproduced with $13-25M_\odot$ faint supernovae or faint hypernovae. Iron-peak element abundances, in particular Zn abundances, are important to put further constraints on the enrichment sources from galactic archaeology surveys.

preprint2013arXiv

Advanced burning stages and fate of 8-10 Mo stars

The stellar mass range 8<M/Mo<12 corresponds to the most massive AGB stars and the most numerous massive stars. It is host to a variety of supernova progenitors and is therefore very important for galactic chemical evolution and stellar population studies. In this paper, we study the transition from super-AGB star to massive star and find that a propagating neon-oxygen burning shell is common to both the most massive electron capture supernova (EC-SN) progenitors and the lowest mass iron-core collapse supernova (FeCCSN) progenitors. Of the models that ignite neon burning off-center, the 9.5Mo model would evolve to an FeCCSN after the neon-burning shell propagates to the center, as in previous studies. The neon-burning shell in the 8.8Mo model, however, fails to reach the center as the URCA process and an extended (0.6 Mo) region of low Ye (0.48) in the outer part of the core begin to dominate the late evolution; the model evolves to an EC-SN. This is the first study to follow the most massive EC-SN progenitors to collapse, representing an evolutionary path to EC-SN in addition to that from SAGB stars undergoing thermal pulses. We also present models of an 8.75Mo super-AGB star through its entire thermal pulse phase until electron captures on 20Ne begin at its center and of a 12Mo star up to the iron core collapse. We discuss key uncertainties and how the different pathways to collapse affect the pre-supernova structure. Finally, we compare our results to the observed neutron star mass distribution.

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

Extraordinary Magnification of the Ordinary Type Ia Supernova PS1-10afx

Recently, Chornock and co-workers announced the Pan-STARRS discovery of a transient source reaching an apparent peak luminosity of ~4x10^44 erg s^-1. We show that the spectra of this transient source are well fit by normal Type Ia supernova (SNIa) templates. The multi-band colors and light-curve shapes are also consistent with normal SNeIa at the spectroscopically determined redshift of z=1.3883; however, the observed flux is a constant factor of ~30 times too bright in each band over time as compared to the templates. At minimum, this shows that the peak luminosities inferred from the light-curve widths of some SNeIa will deviate significantly from the established, empirical relation used by cosmologists. We argue on physical grounds that the observed fluxes do not reflect an intrinsically luminous SNIa, but rather PS1-10afx is a normal SNIa whose flux has been magnified by an external source. The only known astrophysical source capable of such magnification is a gravitational lens. Given the lack of obvious lens candidates, such as galaxy clusters, in the vicinity, we further argue that the lens is a supermassive black hole or a comparatively low-mass dark matter halo. In this case, the lens continues to magnify the underlying host galaxy light. If confirmed, this discovery could impact a broad range of topics including cosmology, gamma-ray bursts, and dark matter halos.

preprint2013arXiv

Supernova Explosions of Super-Asymptotic Giant Branch Stars: Multicolor Light Curves of Electron-Capture Supernovae

An electron-capture supernova (ECSN) is a core-collapse supernova (CCSN) explosion of a super-asymptotic giant branch (SAGB) star with a main-sequence mass $M_{\rm ms}\sim7-9.5M_\odot$. The explosion takes place in accordance with core bounce and subsequent neutrino heating and is a unique example successfully produced by first-principle simulations. This allows us to derive a first self-consistent multicolor light curves of a CCSN. Adopting the explosion properties derived by the first-principle simulation, i.e., the low explosion energy of $1.5\times10^{50}$ erg and the small $^{56}$Ni mass of $2.5\times10^{-3}M_\odot$, we perform a multigroup radiation hydrodynamics calculation of ECSNe and present multicolor light curves of ECSNe of SAGB stars with various envelope mass and hydrogen abundance. We demonstrate that a shock breakout has peak luminosity of $L\sim2\times10^{44}$ erg/s and can evaporate circumstellar dust up to $R\sim10^{17}$ cm for a case of carbon dust, that plateau luminosity and plateau duration of ECSNe are $L\sim10^{42}$ erg/s and $t\sim60-100$ days, respectively, and that a plateau is followed by a tail with a luminosity drop by $\sim4$ mag. The ECSN shows a bright and short plateau that is as bright as typical Type II plateau supernovae, and a faint tail that might be influenced by spin-down luminosity of a newborn pulsar. Furthermore, the theoretical models are compared with ECSN candidates: SN 1054 and SN 2008S. We find that SN 1054 shares the characteristics of the ECSNe. For SN 2008S, we find that its faint plateau requires an ECSN model with a significantly low explosion energy of $E\sim10^{48}$ erg.

preprint2013arXiv

Type Ia Supernova Models and Progenitor Scenarios

We review some recent developments in theoretical studies on the connection between the progenitor systems of Type Ia supernovae (SNe Ia) and the explosion mechanisms. (1) DD-subCh: In the merging of double C+O white dwarfs (DD scenario), if the carbon detonation is induced near the white dwarf (WD) surface in the early dynamical phase, it could result in the (effectively) sub-Chandrasekhar mass explosion. (2) DD-Ch: If no surface C-detonation is ignited, the WD could grow until the Chandrasekhar mass is reached, but the outcome depends on whether the quiescent carbon shell burning is ignited and burns C+O into O+Ne+Mg. (3) SD-subCh: In the single degenerate (SD) scenario, if the He shell-flashes grow strong to induce a He detonation, it leads to the sub-Chandra explosion. (4) SD-Ch: If the He-shell flashes are not strong enough, they still produce interesting amount of Si and S near the surface of C+O WD before the explosion. In the Chandra mass explosion, the central density is high enough to produce electron capture elements, e.g., stable Ni-58. Observations of the emission lines of Ni in the nebular spectra provides useful diagnostics of the sub-Chandra vs. Chandra issue. The recent observations of relatively low velocity carbon near the surface of SNe Ia provide also interesting constraint on the explosion models.

preprint2012arXiv

A Binary Progenitor for the Type IIb Supernova 2011dh in M51

We perform binary stellar evolutionary calculations following the simultaneous evolution of both stars in the system to study a potential progenitor system for the Type IIb supernova 2011dh. Pre-explosion photometry as well as light-curve modeling have provided constraints on the physical properties of the progenitor system. Here we present a close binary system that is compatible with such constraints. The system is formed by stars of solar composition with 16 Msun + 10 Msun on a circular orbit with an initial period of 125 days. The primary star ends its evolution as a yellow supergiant with a mass of ~4 Msun, a final hydrogen content of ~3-5E-03 Msun and with an effective temperature and luminosity in agreement with the HST pre-explosion observations of SN 2011dh. These results are nearly insensitive to the adopted accretion efficiency factor beta. At the time of explosion, the companion star has an effective temperature of 22 to 40 thousand Kelvin, depending on the value of beta, and lies near the zero age main sequence. Considering the uncertainties in the HST pre-SN photometry the secondary star is only marginally detectable in the bluest observed band. Close binary systems, as opposed to single stars, provide a natural frame to explain the properties of SN 2011dh.

preprint2012arXiv

A High-Resolution Spectroscopic Search for the Remaining Donor for Tycho's Supernova

In this paper, we report on our analysis using Hubble Space Telescope astrometry and Keck-I HIRES spectroscopy of the central six stars of Tycho's supernova remnant (SN 1572). With these data, we measured the proper motions, radial velocities, rotational velocities, and chemical abundances of these objects. Regarding the chemical abundances, we do not confirm the unusu- ally high [Ni/Fe] ratio previously reported for Tycho-G. Rather, we find that for all metrics in all stars, none exhibit the characteristics expected from traditional SN Ia single-degenerate-scenario calculations. The only possible exception is Tycho-B, a rare, metal-poor A-type star; however, we are unable to find a suitable scenario for it. Thus, we suggest that SN 1572 cannot be explained by the standard single-degenerate model.

preprint2012arXiv

Astrophysical Particle Simulations on Heterogeneous CPU-GPU Systems

A heterogeneous CPU-GPU node is getting popular in HPC clusters. We need to rethink algorithms and optimization techniques for such system depending on the relative performance of CPU vs. GPU. In this paper, we report a performance optimized particle simulation code "OTOO", that is based on the octree method, for heterogenous systems. Main applications of OTOO are astrophysical simulations such as N-body models and the evolution of a violent merger of stars. We propose optimal task split between CPU and GPU where GPU is only used to compute the calculation of the particle force. Also, we describe optimization techniques such as control of the force accuracy, vectorized tree walk, and work partitioning among multiple GPUs. We used OTOO for modeling a merger of two white dwarf stars and found that OTOO is powerful and practical to simulate the fate of the process.

preprint2012arXiv

Can the Growth of Dust Grains in Low-Metallicity Star-Forming Clouds Affect the Formation of Metal-Poor Low-Mass Stars?

The discovery of a low-mass star with such low metallicity as < 4.5x10^{-5} Z_sun reveals the critical role of dust in the formation of extremely metal-poor stars. In this paper we explore the effect of the growth of dust grains through accretion of gaseous refractory elements in very low-metallicity pre-stellar cores on the cloud fragmentation induced by the dust emission cooling. Employing a simple model of grain growth in a gravitationally collapsing gas, we show that Fe and Si grains can grow efficiently at hydrogen densities of ~10^{10}-10^{14} cm^{-3} in the clouds with metal abundances of -5 <~ [Fe, Si/H] <~ -3. The critical metal number abundances, above which the grain growth could induce the fragmentation of the gas clouds, are estimated to be A_{crit} ~ 10^{-9}-10^{-8}, unless the initial grain radius is too large (>~ 1 um) or the sticking probability is too small (<~ 0.01). We find that even if the initial dust-to-gas mass ratio is well below the minimum value required for the dust-induced fragmentation, the grain growth increases the dust mass high enough to cause the gas fragmentation into sub-solar mass clumps. We suggest that as long as the critical metal abundance is satisfied, the grain growth could play an important role in the formation of low-mass stars with metallicity as low as 10^{-5} Z_sun.

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

Evidence for Type Ia Supernova Diversity from Ultraviolet Observations with the Hubble Space Telescope

We present ultraviolet (UV) spectroscopy and photometry of four Type Ia supernovae (SNe 2004dt, 2004ef, 2005M, and 2005cf) obtained with the UV prism of the Advanced Camera for Surveys on the Hubble Space Telescope. This dataset provides unique spectral time series down to 2000 Angstrom. Significant diversity is seen in the near maximum-light spectra (~ 2000--3500 Angstrom) for this small sample. The corresponding photometric data, together with archival data from Swift Ultraviolet/Optical Telescope observations, provide further evidence of increased dispersion in the UV emission with respect to the optical. The peak luminosities measured in uvw1/F250W are found to correlate with the B-band light-curve shape parameter dm15(B), but with much larger scatter relative to the correlation in the broad-band B band (e.g., ~0.4 mag versus ~0.2 mag for those with 0.8 < dm15 < 1.7 mag). SN 2004dt is found as an outlier of this correlation (at > 3 sigma), being brighter than normal SNe Ia such as SN 2005cf by ~0.9 mag and ~2.0 mag in the uvw1/F250W and uvm2/F220W filters, respectively. We show that different progenitor metallicity or line-expansion velocities alone cannot explain such a large discrepancy. Viewing-angle effects, such as due to an asymmetric explosion, may have a significant influence on the flux emitted in the UV region. Detailed modeling is needed to disentangle and quantify the above effects.

preprint2012arXiv

Final Fates of Rotating White Dwarfs and Their Companions in the Single Degenerate Model of Type Ia Supernovae

Taking into account the rotation of mass-accreting white dwarfs (WDs) whose masses exceed the Chandrasekhar mass, we extend our new single degenerate model for the progenitors of Type Ia supernovae (SNe Ia), accounting for two types of binary systems, those with a main sequence companion and those with a red-giant (RG) companion. We present a mass distribution of WDs exploding as SNe Ia, where the WD mass ranges from 1.38 to 2.3 Msun. These progenitor models are assigned to various types of SNe Ia. A lower mass range of WDs (1.38 Msun < M_WD <~ 1.5 Msun), which are supported by rigid rotation, correspond to normal SNe Ia. A variety of spin-down time may lead to a variation of brightness. A higher mass range of WDs (M_WD >~ 1.5 Msun), which are supported by differential rotation, correspond to brighter SNe Ia such as SN 1991T. In this case, a variety of the WD mass may lead to a variation of brightness. We also show the evolutionary states of the companion stars at SN Ia explosions and pose constraints on the unseen companions. In the WD+RG systems, in particular, most of the RG companions have evolved to helium/carbon-oxygen WDs in the spin-down phase before the SN Ia explosions. In such a case, we do not expect any prominent signature of the companion immediately before and after the explosion. We also compare our new models with the recent stringent constraints on the unseen progenitors of SNe Ia such as SN 2011fe.

preprint2012arXiv

GRB 100418A: a Long GRB without a Bright Supernova in a High-Metallicity Host Galaxy

We present results of a search for a supernova (SN) component associated with GRB 100418A at the redshift of 0.624. The field of GRB 100418A was observed with FOCAS on Subaru 8.2m telescope under a photometric condition (seeing 0.3"-0.4") on 2010 May 14 (UT). The date corresponds to 25.6 days after the burst trigger (15.8 days in the restframe). We did imaging observations in V, Rc, and Ic bands, and two hours of spectrophotometric observations. We got the resolved host galaxy image which elongated 1.6" (= 11 kpc) from north to south. No point source was detected on the host galaxy. The time variation of Rc-band magnitude shows that the afterglow of GRB 100418A has faded to Rc \sim > 24 without SN like rebrightening, when we compare our measurement to the reports in GCN circulars. We could not identify any SN feature such as broad emission-lines or bumps in our spectrum. Assuming the SN is fainter than the 3σ noise spectrum of our observation, we estimate the upper limit on the SN absolute magnitude MIc,obs > -17.2 in observer frame Ic-band. This magnitude is comparable to the faintest type Ic SNe. We also estimate host galaxy properties from the spectrum. The host galaxy of GRB 100418A is relatively massive (log M_{star}/M_{sun} = 9.54) compared to typical long GRB host galaxies, and has 12+log(O/H) = 8.75.

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 Type IIb Supernova 2011dh from a Supergiant Progenitor

A set of hydrodynamical models based on stellar evolutionary progenitors is used to study the nature of SN 2011dh. Our modeling suggests that a large progenitor star ---with R ~200 Rsun---, is needed to reproduce the early light curve of SN 2011dh. This is consistent with the suggestion that the yellow super-giant star detected at the location of the SN in deep pre-explosion images is the progenitor star. From the main peak of the bolometric light curve and expansion velocities we constrain the mass of the ejecta to be ~2 Msun, the explosion energy to be E= 6-10 x 10^50 erg, and the 56Ni mass to be approximately 0.06 Msun. The progenitor star was composed of a helium core of 3 to 4 Msun and a thin hydrogen-rich envelope of ~0.1 M_sun with a main sequence mass estimated to be in the range of 12--15 Msun. Our models rule out progenitors with helium-core masses larger than 8 Msun, which correspond to M_ZAMS > 25 Msun. This suggests that a single star evolutionary scenario for SN 2011dh is unlikely.

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.

preprint2011arXiv

A Single Degenerate Progenitor Model for Type Ia Supernovae Highly Exceeding the Chandrasekhar Mass Limit

Recent observations of Type Ia supernovae (SNe Ia) suggest that some of the progenitor white dwarfs (WDs) had masses up to 2.4-2.8 M_sun, highly exceeding the Chandrasekhar mass limit. We present a new single degenerate (SD) model for SN Ia progenitors, in which the WD mass possibly reaches 2.3-2.7 M_sun. Three binary evolution processes are incorporated; optically thick winds from mass-accreting WDs, mass-stripping from the binary companion star by the WD winds, and WDs being supported by differential rotation. The WD mass can increase by accretion up to 2.3 (2.7) M_sun from the initial value of 1.1 (1.2) M_sun, being consistent with high luminosity SNe Ia such as SN 2003fg, SN 2006gz, SN 2007if, and SN 2009dc. There are three characteristic mass ranges of exploding WDs. In an extreme massive case, differentially rotating WDs explode as an SN Ia soon after the WD mass exceeds 2.4 M_sun because of a secular instability at T/|W|\sim 0.14. For a mid mass range of M_WD=1.5-2.4 M_sun, it takes some time (spinning-down time) until carbon is ignited to induce an SN Ia explosion after the WD mass has reached maximum, because it needs a loss or redistribution of angular momentum. For a lower mass case of rigidly rotating WDs, M_WD=1.38-1.5 M_sun, the spinning-down time depends on the timescale of angular momentum loss from the WD. The difference in the spinning-down time may produce the "prompt" and "tardy" components. We also suggest the very bright super-Chandrasekhar mass SNe Ia are born in a low metallicity environment.

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

On the effects of microphysical grain properties on the yields of carbonaceous dust from type II SNe

We study the role of the unknown microphysical properties of carbonaceous dust particles in determining the amount and size distribution of carbonaceous dust condensed in type II supernova explosions. We parametrize the microphysical properties in terms of the shape factor of the grain and the sticking coefficient of gas-phase carbon atoms onto the grain surfaces. We find that the amount of dust formed is fairly independent of these properties, within the parameter range considered, though limited by the available amount of carbon atoms not locked in CO molecules. However, we find that the condensation times and size distributions of dust grains depend sensitively on the microphysical parameters, with the mass distributions being weighted toward larger effective radii for conditions considering grains with higher sticking coefficients and/or more aspherical shapes. We discuss that this leads to important consequences on the predicted extinction law of SN dust and on the survival rate of the formed grains as they pass through the reverse shock of the SN. We conclude that a more detailed understanding of the dust formation process and of the microphysical properties of each dust species needs to be achieved before robust prediction on the SN dust yields can be performed.

preprint2011arXiv

The Role of Dust in the Early Universe I : Protogalaxy Evolution

We develop one-zone galaxy formation models in the early Universe, taking into account dust formation and evolution by supernova (SN) explosions. Especially we focus on the time evolution of dust size distribution, because ${\rm H}_{2}$ formation on the dust surface plays a critical role in the star formation process in the early Universe. In the model we assume that star formation rate (SFR) is proportional to the total amount of ${\rm H}_{2}$. We consistently treat (i) the formation and size evolution of dust, (ii) the chemical reaction networks including ${\rm H}_{2}$ formation both on the surface of dust and in gas phase, and (iii) the SFR in the model. First, we find that, because of dust destruction due to both reverse and forward shocks driven by SNe, H$_{2}$ formation is more suppressed than that without dust destruction. At the galaxy age of $sim0.8 {\rm Gyr}$, for galaxy models with virial mass $M_{\rm vir}=10^{9} M_{\odot}$ and formation redshift $z_{\rm vir}=10$, the molecular fraction is 2.5 orders of magnitude less in the model with dust destruction by both shocks than that in the model without dust destruction. Second, we show that the H$_{2}$ formation rate strongly depends on the ISM density around SN progenitors. The SFR in higher ISM density is lower, since dust destruction by reverse shocks is more effective in higher ISM density. We conclude that not only the amount but also the size distribution of dust being related with the star formation activity strongly affects the evolution of galaxies in the early Universe.

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

Formation and Evolution of Dust in Type IIb Supernova with Application to the Cassiopeia A Supernova Remnant

We investigate the formation of dust grains in the ejecta of a SN IIb and their evolution in the shocked gas in the SNR by considering the uniform and power-law density structures for the CSM. Based on these calculations, we also simulate the time evolution of thermal emission from the shock-heated dust in the SNR and compare the results with the observations of Cas A SNR. We find that the total mass of dust formed in the SN IIb is as large as 0.167 M_sun but the average radius of dust is smaller than 0.01 mum and is significantly different from those in SNe II-P with the massive H-envelope. In the explosion with the small-mass H-envelope, the expanding He core undergoes little deceleration, so that the gas density in the He core is too low for large-sized grains to form. In addition, the low-mass H-envelope of the SN IIb leads to the early arrival of the reverse shock at the dust-forming region. If the CSM is more or less spherical, therefore, the newly formed grains would be completely destroyed in the relatively dense shocked gas for the CSM density of n_H > 0.1 cm^{-3}. However, the actual CSM is likely to be non-spherical, so that a part of grains could be ejected into the ISM without being shocked. We demonstrate that the time evolution of the SED by thermal dust emission is sensitive to the ambient gas density and structure that affects the passage of the reverse shock into the ejecta. Thus, the SED evolution well reflects the evolution of dust through erosion by sputtering and stochastic heating. For Cas A, we consider the CSM produced by the steady mass loss of ~8x10^{-5} M_sun/yr during the supergiant phase. Then we find the infrared SED of Cas A is reasonably reproduced by thermal emission from the newly formed dust of 0.08 M_sun, which consists of 0.008 M_sun shocked warm dust and 0.072 M_sun unshocked cold dust.

preprint2010arXiv

Production of the p-Process Nuclei in the Carbon-Deflagration Model for Type Ia Supernovae

We calculate nucleosynthesis of proton-rich isotopes in the carbon-deflagration model for Type Ia supernovae (SNe Ia). The seed abundances are obtained by calculating the s-process nucleosynthesis that is expected to occur in the repeating helium shell flashes on the carbon-oxygen (CO) white dwarf (WD) during mass accretion from a binary companion. When the deflagration wave passes through the outer layer of the CO WD, p-nuclei are produced by photodisintegration reactions on s-nuclei in a region with the peak temperature ranging from 1.9 to 3.6 x10^9 K. We confirm the sensitivity of the p-process on the initial distribution of s-nuclei. We show that the initial C/O ratio in the WD does not much affect the yield of p-nuclei. On the other hand, the abundance of 22Ne left after the s-processing has a large influence on the p-process via 22Ne(alpha,n) reaction. We find that about 50% of p-nuclides are co-produced when normalized to their solar abundances in all adopted cases of seed distribution. Mo and Ru, which are largely underproduced in Type II supernovae (SNe II), are produced more than in SNe II although they are underproduced with respect to the yield levels of other p-nuclides. The ratios between p-nuclei and iron in the ejecta are larger than the solar ratios by a factor of 1.2. We also compare the yields of oxygen, iron, and p-nuclides in SNe Ia and SNe II and suggest that SNe Ia could make a larger contribution than SNe II to the solar system content of p-nuclei.

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.

preprint2010arXiv

Supersoft X-ray Phase of Single Degenerate Type Ia Supernova Progenitors in Early Type Galaxies

In the single degenerate (SD) scenario for Type Ia supernova (SN Ia) progenitors, an accreting white dwarf (WD) is expected to undergo a supersoft X-ray source (SSS) phase. Recently, Gilfanov & Bogdan (2010, hereafter GB10) claimed that observed X-ray fluxes of early type galaxies would be too low to be consistent with the prediction of the SD scenario based on rather simple assumptions. We present realistic evolutionary models of SD systems and calculate durations of SSS phases. In most cases, accreting WDs spend a large fraction of time in the optically thick wind phase and the recurrent nova phase rather than the SSS phase. Thus the SSS phase lasts only for a few hundred thousand years. This is by a factor of ~10 shorter than those adopted by GB10 where the SN~Ia progenitor WD was assumed to spend most of its life as a SSS. The theoretical X-ray luminosity of the SSS has a large uncertainty because of the uncertain atmospheric model of mass-accreting WDs and absorption of soft X-rays by the companion star's cool wind material. We thus adopt an average of the observed fluxes of existing symbiotic SSSs, i.e., ~0.4 x 10^{36} erg s^{-1} for 0.3--0.7 keV. Using these SSS duration and soft X-ray luminosity, we show that the observed X-ray flux obtained by GB10 is rather consistent with our estimated flux in early type galaxies based on the SD scenario. This is a strong support for the SD scenario as a main-contributor of SNe Ia in early type galaxies.

preprint2009arXiv

Empirical Abundance Scaling Laws and Implications for the Gamma-Process in Core-Collapse Supernovae

Analyzing the solar system abundances, we have found two empirical abundance scaling laws concerning the p- and s-nuclei with the same atomic number. The first scaling is s/p ratios are almost constant over a wide range of the atomic number, where the p-nculei are lighter than the s-nuclei by two or four neutrons. The second scaling is p/p ratios are almost constant, where the second $p$-nuclei are lighter than the first p-nucleus by two neutrons. These scalings are a piece of evidence that most p-nuclei are dominantly synthesized by the gamma-process in supernova explosions. The scalings lead to a novel concept of "universality of gamma-process" that the s/p and p/p ratios of nuclei produced by individual gamma-processes are almost constant, respectively. We have calculated the ratios by gamma-process based on core-collapse supernova explosion models under various astrophysical conditions and found that the scalings hold for materials produced by individual gamma-processes independent of the astrophysical conditions assumed. The universality originates from three mechanisms: the shifts of the gamma-process layers to keep their peak temperature, the weak s-process in pre-supernovae, and the independence of the s/p ratios of the nuclear reactions. The results further suggest an extended universality that the s/p ratios in the gamma-process layers are not only constant but also centered on a specific value of 3. With this specific value and the first scaling, we estimate that the ratios of $s$-process abundance contributions from the AGB stars to the massive stars are almost 6.7 for the $s$-nuclei of A > 90. We find that large enhancements of s/p ratios for Ce, Er, and W are a piece of evidence that the weak s-process actually occurred before SNe.

preprint2009arXiv

Nebular Spectra and Explosion Asymmetry of Type Ia Supernovae

The spectral signatures of asymmetry in Type Ia Supernova (SN Ia) explosions are investigated, using a sample of late-time nebular spectra. First, a kinematical model is constructed for SN Ia 2003hv, which can account for the main features in its optical, Near-Infrared (NIR), and Mid-Infrared (Mid-IR) late-time spectra. It is found that an asymmetric off-center model can explain the observed characteristics of SN 2003hv. This model includes a relatively high density, Fe-rich region which displays a large velocity off-set, and a relatively low density, extended 56Ni-rich region which is more spherically distributed. The high density region consists of the inner stable Fe-Ni region and outer 56Ni-rich region. Such a distribution may be the result of a delayed-detonation explosion, in which the first deflagration produces the global asymmetry in the innermost ejecta, while the subsequent detonation can lead to the bulk spherical symmetry. This configuration, if viewed from the direction of the off-set, can consistently explain the blueshift in some of the emission lines and virtually no observed shift in other lines in SN 2003hv. For this model, we then explore the effects of different viewing angles and the implications for SNe Ia in general. The model predicts that a variation of the central wavelength, depending on the viewing angle, should be seen in some lines (e.g., [Ni II]7378), while the strongest lines (e.g., [Fe III] blend at 4700A) will not show this effect. By examining optical nebular spectra of 12 SNe Ia, we have found that such a variation indeed exists. We suggest that the global asymmetry in the innermost ejecta, as likely imprint of the deflagration flame propagation, is a generic feature of SNe Ia (abridged).

preprint2009arXiv

SN2003bg: a broad-lined Type IIb Supernova with Hydrogen

Models for the spectra and the light curve, in the photospheric as well as in the late nebular phase, are used to infer the properties of the very radio-bright, broad-lined type IIb Supernova 2003bg. Consistent fits to the light curve and the spectral evolution are obtained with an explosion that ejected ~ 4 M_sun of material with a kinetic energy of ~ 5 10^51 erg. A thin layer of hydrogen, comprising ~ 0.05 M_sun, is inferred to be present in the ejecta at the highest velocities (v >~ 9000 km/s), while a thicker helium layer, comprising ~ 1.25 M_sun, was ejected at velocities between 6500 and 9000 km/s. At lower velocities, heavier elements are present, including ~ 0.2 M_sun of 56Ni that shape the light curve and the late-time nebular spectra. These values suggest that the progenitor star had a mass of ~ 20-25 M_sun (comparable to, but maybe somewhat smaller than that of the progenitor of the XRF/SN 2008D). The rather broad-lined early spectra are the result of the presence of a small amount of material (~ 0.03 M_sun) at velocities > 0.1 c, which carries ~ 10 % of the explosion kinetic energy. No clear signatures of a highly aspherical explosion are detected.

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

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.

preprint2007arXiv

Evolution of Dust in Primordial Supernova Remnants: Can Dust Grains Formed in the Ejecta Survive and be Injected into the Early Interstellar Medium?

We investigate the evolution of dust that formed at Population III supernova (SN) explosions and its processing through the collisions with the reverse shocks resulting from the interaction of the SN ejecta with the ambient medium. In particular, we investigate the transport of the shocked dust within the SN remnant (SNR), and its effect on the chemical composition, the size distribution, and the total mass of dust surviving in SNR. We find that the evolution of the reverse shock, and hence its effect on the processing of the dust depends on the thickness of the envelope retained by the progenitor star. Furthermore, the transport and survival of the dust grains depend on their initial radius, a_{ini}, and composition: For Type II SNRs expanding into the interstellar medium (ISM) with a density of n_{H,0}=1 cm^{-3}, small grains with a_{ini} < ~ 0.05 micron are completely destroyed by sputtering in the postshock flow, while grains with a_{ini}= 0.05--0.2 micron are trapped into the dense shell behind the forward shock. Very large grains of a_{ini} > ~ 0.2 micron are ejected into the ISM without decreasing their sizes significantly. We find that the total mass fraction of dust that is destroyed by the reverse shock ranges from 0.2 to 1.0, depending on the energy of the explosion and the density of the ambient ISM. The results of our calculations have significant impact on the abundance pattern of subsequent generation of stars that form in the dense shell of primordial SNRs.

preprint2007arXiv

Spectroscopic Studies of Extremely Metal-Poor Stars with the Subaru High Dispersion Spectrograph. IV. The alpha-Element-Enhanced Metal-Poor Star BS16934-002

A detailed elemental abundance analysis has been carried out for the very metal-poor ([Fe/H]=-2.7) star BS16934-002, which was identified in our previous work as a star exhibiting large overabundances of Mg and Sc. A comparison of the abundance pattern of this star with that of the well-studied metal-poor star HD122563 indicates excesses of O, Na, Mg, Al, and Sc in BS16934-002. Of particular interest, no excess of C or N is found in this object, in contrast to CS22949-037 and CS29498-043, two previously known carbon-rich, extremely metal-poor stars with excesses of the alpha elements. No established nucleosynthesis model exists that explains the observed abundance pattern of BS16934-002. A supernova model, including mixing and fallback, assuming severe mass loss before explosion, is discussed as a candidate progenitor of BS16934-002.

preprint2007arXiv

Supernova Nucleosynthesis in Population III 13 -- 50 $M_{\odot}$ Stars and Abundance Patterns of Extremely Metal-Poor Stars

We perform hydrodynamical and nucleosynthesis calculations of core-collapse supernovae (SNe) and hypernovae (HNe) of Population (Pop) III stars. We provide new yields for the main-sequence mass of $M_{\rm MS}=13-50$ $M_{\odot}$ and the explosion energy of $E=1-40\times10^{51}$ ergs to apply for chemical evolution studies. Our HN yields based on the mixing-fallback model of explosions reproduce the observed abundance patterns of extremely metal-poor (EMP) stars ($-4< {\rm [Fe/H]}< -3$), while those of very metal-poor (VMP) stars ($-3<{\rm [Fe/H]} < -2$) are reproduced by the normal SN yields integrated over the Salpeter initial mass function. Moreover, the observed trends of abundance ratios [X/Fe] against [Fe/H] with small dispersions for the EMP stars can be reproduced as a sequence resulting from the various combination of $M_{\rm MS}$ and $E$. This is because we adopt the empirical relation that a larger amount of Fe is ejected by more massive HNe. Our results imply that the observed trends with small dispersions do not necessarily mean the rapid homogeneous mixing in the early galactic halo at [Fe/H] $< -3$, but can be reproduced by the ``inhomogeneous'' chemical evolution model. In addition, we examine how the modifications of the distributions of the electron mole fraction $Y_{\rm e}$ and the density in the presupernova models improve the agreement with observations. In this connection, we discuss possible contributions of nucleosynthesis in the neutrino-driven wind and the accretion disk.

preprint2007arXiv

Thermal Stability of White Dwarfs Accreting Hydrogen-rich Matter and Progenitors of Type Ia Supernovae

We revisit the properties of white dwarfs accreting hydrogen-rich matter by constructing steady-state models, in which hydrogen shell burning consumes hydrogen at the same rate as the white dwarf accretes it. We obtain such steady-state models for various accretion rates and white dwarf masses. We confirm that these steady models are thermally stable only when the accretion rate is higher than \sim 10^{-7} M_sun/yr. We show that recent models of ``quiescent burning'' in the ``surface hydrogen burning'' at a much wider range of accretion rates results from the too large zone mass in the outer part of the models; hydrogen burning must occur in a much thinner layer. A comparison of the positions on the HR diagram suggests that most of the luminous supersoft X-ray sources are white dwarfs accreting matter at rates high enough that the hydrogen burning shell is thermally stable. Implications on the progenitors of Type Ia supernovae are discussed.

preprint2006arXiv

Chemical Abundances in the Secondary Star of the Black Hole Binary V4641 Sagittarii (SAX J1819.3-2525)

We report on detailed spectroscopic studies performed for the secondary star in the black hole binary (micro-quasar) V4641 Sgr in order to examine its surface chemical composition and to see if its surface shows any signature of pollution by ejecta from a supernova explosion. High-resolution spectra of V4641 Sgr observed in the quiescent state in the blue-visual region are compared with those of the two bright well-studied B9 stars (14 Cyg and $ν$ Cap) observed with the same instrument. The effective temperature of V4641 Sgr (10500 $\pm$ 200 K) is estimated from the strengths of He~{\sc i} lines, while its rotational velocity, $\it v$ sin $\it i$ (95 $\pm$ 10 km s${}^{-1}$), is estimated from the profile of the Mg~{\sc ii} line at 4481 Å. We obtain abundances of 10 elements and find definite over-abundances of N (by 0.8 dex or more) and Na (by 0.8 dex) in V4641 Sgr. From line-by-line comparisons of eight other elements (C, O, Mg, Al, Si, Ti, Cr, and Fe) between V4641 Sgr and the two reference stars, we conclude that there is no apparent difference in the abundances of these elements between V4641 Sgr and the two normal late B-type stars, which have been reported to have solar abundances. An evolutionary model of a massive close binary system has been constructed to explain the abundances observed in V4641 Sgr. The model suggests that the progenitor of the black hole forming supernova was as massive as ~ 35 Msun on the main-sequence and, after becoming a ~ 10 Msun He star, underwent "dark" explosion which ejected only N and Na-rich outer layer of the He star without radioactive $^{56}$Ni.

preprint2006arXiv

Nucleosynthesis Yields of Core-Collapse Supernovae and Hypernovae, and Galactic Chemical Evolution

We present new nucleosynthesis yields as functions of the stellar mass, metallicity, and explosion energy (corresponding to normal supernovae and Hypernovae). We apply the results to the chemical evolution of the solar neighborhood. Our new yields are based on the new developments in the observational/theoretical studies of supernovae (SNe) and extremely metal-poor (EMP) stars in the halo, which have provided excellent opportunities to test the explosion models and their nucleosynthesis. We use the light curve and spectra fitting of individual SN to estimate the mass of the progenitor, explosion energy, and produced 56Ni mass. Comparison with the abundance patterns of EMP stars has made it possible to determine the model parameters of core-collapse SNe, such as mixing-fallback parameters. More specifically, we take into account the two distinct new classes of massive SNe: 1) very energetic Hypernovae, whose kinetic energy (KE) is more than 10 times the KE of normal core-collapse SNe, and 2) very faint and low energy SNe (Faint SNe). These two new classes of SNe are likely to be ``black-hole-forming'' SNe with rotating or non-rotating black holes. Nucleosynthesis in Hypernovae is characterized by larger abundance ratios (Zn,Co,V,Ti)/Fe and smaller (Mn,Cr)/Fe than normal SNe, which can explain the observed trends of these ratios in EMP stars. Nucleosynthesis in Faint SNe is characterized by a large amount of fall-back, which explains the abundance pattern of the most Fe-poor stars. These comparisons suggest that black-hole-forming SNe made important contributions to the early Galactic (and cosmic) chemical evolution.

preprint2005arXiv

GEOTAIL observation of the SGR1806-20 Giant Flare: The first 600 ms

On December 27, 2004, plasma particle detectors on the GEOTAIL spacecraft detected an extremely strong signal of hard X-ray photons from the giant flare of SGR1806-20, a magnetar candidate. While practically all gamma-ray detectors on any satellites were saturated during the first ~500 ms interval after the onset, one of the particle detectors on GEOTAIL was not saturated and provided unique measurements of the hard X-ray intensity and the profile for the first 600 ms interval with 5.48 ms time resolution. After ~50 ms from the initial rapid onset, the peak photon flux (integrated above ~50 keV) reached the order of 10^7 photons sec^{-1} cm^{-2}. Assuming a blackbody spectrum with kT=175 keV, we estimate the peak energy flux to be 21 erg sec^{-1} cm^{-2} and the fluence (for 0-600 ms) to be 2.4 erg cm^{-2}. The implied energy release comparable to the magnetic energy stored in a magnetar (~10^{47} erg) suggests an extremely efficient energy release mechanism.

preprint2003arXiv

First-generation black-hole-forming supernovae and the metal abundance pattern of a very iron-poor star

It has been proposed theoretically that the first generation of stars in the Universe (population III) would be as massive as 100 solar masses (100Mo), because of inefficient cooling of the precursor gas clouds. Recently, the most iron-deficient (but still carbon-rich) low-mass star -- HE0107-5240 -- was discovered. If this is a population III that gained its metals (elements heavier than helium) after its formation, it would challenge the theoretical picture of the formation of the first stars. Here we report that the patterns of elemental abundance in HE0107-5240 (and other extremely metal-poor stars) are in good accord with the nucleosynthesis that occurs in stars with masses of 20-130Mo when they become supernovae if, during the explosions, the ejecta undergo substantial mixing and fall-back to form massive black holes. Such supernovae have been observed. The abundance patterns are not, however, consistent with enrichment by supernovae from stars in the range 130-300 Mo. We accordingly infer that the first-generation supernovae came mostly from explosions of ~ 20-130Mo stars; some of these produced iron-poor but carbon- and oxygen-rich ejecta. Low-mass second-generation stars, like HE0107-5240, could form because the carbon and oxygen provided pathways for gas to cool.

preprint2002arXiv

Carbon Monoxide in the Type Ic SN 2000ew

We present K-band (1.9 -- 2.5 micron) spectra of the Type Ic SN 2000ew observed with IRCS on the Subaru Telescope. These data show the first detection of carbon monoxide (CO) emission in a Type Ic supernova. The detection of CO in SN 2000ew provides further evidence that molecule formation may be a common occurrence in core-collapse supernova ejecta. The spectrum also contains narrow emission lines of [Fe II] and He I probably from dense clumps of hydrogen-poor circumstellar gas surrounding SN 2000ew. Our spectrum of SN 2000ew shows no trace of an unidentified feature seen near 2.26 micron, just blueward of the CO emission, in the spectrum of SN 1987A and we discuss proposed detections of this feature in other Type II supernovae.

preprint2001arXiv

Nucleosynthesis of Zinc and Iron-Peak Elements in Pop III Type II Supernovae : Comparison with abundances of Very Metal-Poor Halo Stars

We calculate nucleosynthesis in core-collapse explosions of massive Pop III stars, and compare the results with abundances of metal-poor halo stars to constrain the parameters of Pop III supernovae. We focus on iron-peak elements and, in particular, we try to reproduce the large [Zn/Fe] observed in extremely metal-poor stars. The interesting trends of the observed ratios [Zn, Co, Mn, Cr, V/Fe] can be related to the variation of the relative mass of the complete and incomplete Si-burning regions in supernova ejecta. We find that [Zn/Fe] is larger for deeper mass-cuts, smaller neutron excess, and larger explosion energies. The large [Zn/Fe] and [O/Fe] observed in the very metal-poor halo stars suggest deep mixing of complete Si-burning material and a significant amount of fall-back in Type II supernovae. Furthermore, large explosion energies (E_51 >~ 2 for M ~ 13 Msun and E_51 >~ 20 for M >~ 20 Msun) are required to reproduce [Zn/Fe] ~ 0.5. The observed trends of the abundance ratios among the iron-peak elements are better explained with this high energy (``Hypernova'') models rather than the simple ``deep'' mass-cut effect, because the overabundance of Ni can be avoided in the hypernova models. We also present the yields of pair-instability supernova explosions of M = 130 - 300 Msun stars, and discuss that the abundance features of very metal-poor stars cannot be explained by pair-instability supernovae.

preprint1998arXiv

Dissipative Process as a Mechanism of Differentiating Internal Structures between Dwarf and Normal Elliptical Galaxies in a CDM Universe

We simulate the dynamical, chemical, and spectro-photometric evolution of dwarf and normal elliptical galaxies embedded in a dark matter halo, using a three-dimensional N-body/SPH simulation code. For the forming dwarf elliptical galaxies, supernova-driven winds propagating outwards from inside the system collide with the infalling gas and produce the super-shell in which stars are formed. The resulting stellar system forms a loosely bound virialized system due to the significant mass loss and has a large velocity dispersion and a large core. Consequently the surface brightness distribution shows an exponential profile and the color distribution shows a positive gradient such that the colors become redder away from the galaxy center in agreement with observations. On the other hand, for the normal elliptical galaxies with deep gravitational potential, the mass loss out of the system does not have a significant dynamical effect on the stellar system. The resulting surface brightness distribution has a large central concentration like de Vaucouleurs' $r^{1/4}$-profile and the color distribution shows a negative gradient as observed. Our simulation shows that different features between dwarf and normal elliptical galaxies stem from different cooling efficiencies for their respective protogalaxies in a standard CDM universe.

preprint1998arXiv

The Earliest Phases of Galaxy Evolution

In this paper we study the very early phases of the evolution of our Galaxy by means of a chemical evolution model which reproduces most of the observational constraints in the solar vicinity and in the disk. We have restricted our analysis to the solar neighborhood and present the predicted abundances of several elements (C, N, O, Mg, Si, S, Ca, Fe) over an extended range of metallicities $[Fe/H] = -4.0$ to $[Fe/H] = 0.0$ compared to previous models. We adopted the most recent yield calculations for massive stars taken from different authors (Woosley & Weaver 1995 and Thielemann et al. 1996) and compared the results with a very large sample of data, one of the largest ever used to this purpose. These data have been analysed with a new and powerful statistical method which allows us to quantify the observational spread in measured elemental abundances and obtain a more meaningful comparison with the predictions from our chemical evolution model. Our analysis shows that the ``plateau'' observed for the [$α$/Fe] ratios at low metallicities ($-3.0< [Fe/H] <-1.0$) is not perfectly constant but it shows a slope, especially for oxygen. This slope is very well reproduced by our model with both sets of yields. This is not surprising since realistic chemical evolution models, taking into account in detail stellar lifetimes, never predicted a completely flat plateau. This is due either to the fact that massive stars of different mass produce a slightly different O/Fe ratio or to the often forgotten fact that supernovae of type Ia, originating from white dwarfs, start appearing already at a galactic age of 30 million years and reach their maximum at 1 Gyr.

preprint1997arXiv

The Evolution of Dwarf Galaxies with Star Formation in Outward Propagating Super Shell

We simulate the dynamical and chemical evolution of a dwarf galaxy embedded in a dark matter halo, using a three-dimensional N-body/SPH simulation code combined with stellar population synthesis. The initial condition is adopted in accord with a $10^{10}M_\odot$ virialized sphere in a 1-sigma CDM perturbation which contains 10% baryonic mass. A supersonic spherical outflow is driven by the first star burst near the center of the galaxy and produces an expanding super shell in which stars are subsequently formed. Consecutive formation of stars in the expanding shell makes the stellar system settled with the exponential brightness profile, the positive metallicity gradient, and the inverse color gradient in agreement with observed features of dwarf galaxies. We therefore propose that the energy feedback via stellar winds and supernovae is a decisive mechanism for formation of less compact, small systems like dwarf galaxies.