Source author record

Toshikazu Shigeyama

Toshikazu Shigeyama 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

48works
9topics
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

48 published item(s)

preprint2026arXiv

Delayed Wind Onset in Pa 30, the Remnant of Type Iax SN 1181

Pa 30 is the recently identified remnant of the historical supernova SN 1181, likely a Type Iax event, and a nebula surrounding the central white dwarf launching a fast wind ($\sim10^9~\cm~\s^{-1}$) is observed in optical and infrared bands. X-ray observations show that this wind collides with the surrounding material and produces a termination shock, and the observed extent of the shock indicates that the wind started blowing centuries after 1181 A.D. rather than immediately after the SN explosion. We propose that the wind is triggered by delayed ignition of fallback carbon-rich material on the WD surface and investigate the conditions that reproduce such delayed ignition. We show that producing delays of several centuries requires a relatively hot post-explosion WD core with a temperature $T_c \simeq 6\times10^8~\mathrm{K}$. This supports the pure-deflagration progenitor scenario for Type Iax SN 1181, which implies the presence of a He star companion inside Pa~30; we also discuss why such a potential He star has not been detected and its prospects for discovery by future observations.

preprint2022arXiv

CHIPS: Complete History of Interaction-Powered Supernovae

We present the public release of the Complete History of Interaction-Powered Supernovae (CHIPS) code, suited to model a variety of transients that arise from interaction with a dense circumstellar medium (CSM). Contrary to existing modellings which mostly attach the CSM by hand, CHIPS self-consistently simulates both the creation of the CSM from mass eruption of massive stars prior to core-collapse, and the subsequent supernova light curve. We demonstrate the performance of CHIPS by presenting examples of the density profiles of the CSM and the light curves. We show that the gross light curve properties of putative interaction-powered transients, such as Type IIn supernovae, rapidly evolving transients and recently discovered fast blue optical transients, can be comprehensively explained with the output of CHIPS.

preprint2022arXiv

Eruption of the Envelope of Massive Stars by Energy Injection with Finite Duration

A significant fraction of supernovae show signatures of dense circumstellar material (CSM). While multiple scenarios for creating a dense CSM exist, mass eruption due to injection of energy at the base of the outer envelope is a likely possibility. We carry out radiation hydrodynamical simulations of eruptive mass loss from a typical red supergiant progenitor with initial mass of $15\ M_\odot$, for the first time focusing on the timescale of the injection as well as energy. We find that not only sufficient injection energy but also sufficient rate of energy injection per unit time, $L_{\rm{min}} \sim 8\times 10^{40}$ erg s$^{-1}$ in this particular model, is required for eruption of unbound CSM. This result suggests that the energy injection rate needs to be greater than the binding energy of the envelope divided by the dynamical timescale for the eruption. The density profile of the resulting CSM, whose shape was analytically and numerically predicted in the limit of instantaneous energy injection, similarly holds for a finite injection timescale. We discuss our findings in the framework of proposed mass outburst scenarios, specifically wave-driven outbursts and common envelope ejection.

preprint2022arXiv

Fast optical flares from M dwarfs detected by a one-second-cadence survey with Tomo-e Gozen

We report a one-second-cadence wide-field survey for M-dwarf flares using the Tomo-e Gozen camera mounted on the Kiso Schmidt telescope. We detect 22 flares from M3-M5 dwarfs with rise times and amplitudes ranging from $5\, \mathrm{sec} \lesssim t_\mathrm{rise} \lesssim 100\,\mathrm{sec}$ and $0.5 \lesssim ΔF/F_{\star} \lesssim 20$, respectively. The flare light curves mostly show steeper rises and shallower decays than those obtained from the Kepler one-minute cadence data and tend to have flat peak structures. Assuming a blackbody spectrum with temperatures of $9,000-15,000\,\mathrm{K}$, the peak luminosities and bolometric energies are estimated to be $10^{29}\,\mathrm{erg\,sec^{-1}} \lesssim L_\mathrm{peak} \lesssim 10^{31}\,\mathrm{erg\,sec^{-1}}$ and $10^{31}\,\mathrm{erg} \lesssim E_{\rm bol} \lesssim 10^{34}\,\mathrm{erg}$, which constitutes the bright end of fast optical flares for M dwarfs. We confirm that more than 90\% of the host stars of the detected flares are magnetically active based on their H$α$ emission line intensities obtained by LAMOST. The estimated occurrence rate of the detected flares is $\sim 0.7$ per day per an active star, indicating they are common in magnetically active M dwarfs. We argue that the flare light curves can be explained by the chromospheric compression model; the rise time is broadly consistent with the Alfvén transit time of a magnetic loop with a length scale of $l_\mathrm{loop} \sim 10^4\,\mathrm{km}$ and a field strength of $1,000\,\mathrm{G}$, while the decay time is likely determined by the radiative cooling of the compressed chromosphere down to near the photosphere with a temperature of $\gtrsim 10,000\,\mathrm{K}$. These flares from M dwarfs could be a major contamination source for a future search of fast optical transients of unknown types.

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

Stationary accretion flow with nuclear burning

We present a series of numerical solutions of spherically symmetric stationary flows with nuclear burning accreted by a neutron star (or black hole). We consider the accretion of matter composed of carbon and oxygen, which mimics the flow after a neutron star is engulfed by a CO star or CO core of a massive star. It is found that there are two types of transonic solutions depending on the accretion rate. The flow with a small accretion rate reaches the center (or the surface of the central object) at supersonic speeds. The other type with a large accretion rate has another sonic point inside the transonic point and the flow truncates at the sonic point. The critical accretion rate dividing these two types is derived as a function of the mass of the central object and the specific enthalpy in the ambient matter. We discuss implications from the solutions for a new mechanism of super-Chandrasekhar type Ia supernovae and type Icn supernovae.

preprint2022arXiv

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

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

preprint2021arXiv

Comparison between the first and second mass eruptions from progenitors of Type IIn supernovae

Some massive stars experience episodic and intense mass loss phases with fluctuations in the luminosity. Ejected material forms circumstellar matter around the star, and the subsequent core collapse results in a Type IIn supernova that is characterized by interaction between supernova ejecta and circumstellar matter. The energy source that triggers these mass eruptions and dynamics of the outflow have not been clearly explained. Moreover, the mass eruption itself can alter the density structure of the envelope and affect the dynamics of the subsequent mass eruption if these events are repeated. A large amount of observational evidence suggests multiple mass eruptions prior to core collapse. We investigate the density structure of the envelope altered by the first mass eruption and the nature of the subsequent second mass eruption event in comparison with the first event. We deposited extra energy at the bottom of the hydrogen envelope of 15$M_\odot$ stars twice and calculated the time evolution by radiation hydrodynamical simulation code. We did not deal with the origin of the energy source, but focused on the dynamics of repeated mass eruptions from a single massive star. There are significant differences between the first and second mass eruptions in terms of the luminosity and the color. The second eruption leads to a redder burst event in which the associated brightening phase lasts longer than the first. The amount of ejected matter is different even with the same deposited energy in the first and second event, but the difference depends on the density structure of the star. Upcoming high cadence and deep transient surveys will provide us a lot of pre-supernova activities, and some of which might show multi-peaked light curves. These should be interpreted taking the effect of density structure altered by the preceding outburst events into consideration.

preprint2021arXiv

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

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

preprint2020arXiv

Fast Luminous Blue Transients in the Reionization Era and Beyond

To determine the epoch of reionization precisely and to reveal the property of inhomogeneous reionization are some of the most important topics of modern cosmology. Existing methods to investigate reionization which use cosmic microwave background, Ly$α$ emitters, quasars, or gamma ray bursts, have difficulties in terms of accuracy or event rate. We propose that recently discovered fast luminous blue transients (FLBTs) have a potential as a novel probe of reionization. We study the detectability of FLBTs at the epoch of reionization with upcoming WFIRST Wide-Field Instruments (WFI), using a star formation rate derived from galaxy observations and an event rate of FLBTs proportional to the star formation rate. We find that if FLBTs occur at a rate of 1\% of the core-collapse supernova rate, 2 (0.3) FLBTs per year per deg$^2$ at $z>6$ ($z>8$) can be detected by a survey with a limiting magnitude of 26.5 mag in the near-infrared band and a cadence of 10 days. We conclude that the WFIRST supernova deep survey can detect $\sim20$ FLBTs at the epoch of reionization in the near future.

preprint2020arXiv

Intermediate Luminosity Red Transients by Black Holes Born from Erupting Massive Stars

We consider black hole formation in failed supernovae when a dense circumstellar medium (CSM) is present around the massive star progenitor. By utilizing radiation hydrodynamical simulations, we calculate the mass ejection of blue supergiants and Wolf-Rayet stars in the collapsing phase and the radiative shock occurring between the ejecta and the ambient CSM. We find that the resultant emission is redder and dimmer than normal supernovae (bolometric luminosity of $\sim 10^{40}-10^{41}\ {\rm erg\ s^{-1}}$, effective temperature of $\sim 5\times 10^3$ K, and timescale of 10-100 days) and shows a characteristic power-law decay, which may comprise a fraction of intermediate luminosity red transients (ILRTs) including AT 2017be. In addition to searching for the progenitor star in the archival data, we encourage X-ray follow-up observations of such ILRTs $\sim$ 1-10 yr after the collapse, targeting the fallback accretion disk.

preprint2020arXiv

Radiation hydrodynamical simulations of eruptive mass Loss from progenitors of Type Ibn/IIn supernovae

Observations suggest that some massive stars experience violent and eruptive mass loss associated with significant brightening that cannot be explained by hydrostatic stellar models. This event seemingly forms dense circumstellar matter (CSM). The mechanism of eruptive mass loss has not been fully explained. We focus on the fact that the timescale of nuclear burning gets shorter than the dynamical timescale of the envelope a few years before core collapse for some massive stars. To reveal the properties of the eruptive mass loss, we investigate its relation to the energy injection at the bottom of the envelope supplied by nuclear burning taking place inside the core. In this study, we do not specify the actual mechanism for transporting energy from the site of nuclear burning to the bottom of the envelope. Instead, we parameterize the amount of injected energy and the injection time and try to extract information on these parameters from comparisons with observations. We carried out 1-D radiation hydrodynamical simulations for progenitors of red, yellow, and blue supergiants, and Wolf-Rayet stars. We calculated the evolution of the progenitors with a public stellar evolution code. We obtain the light curve associated with the eruption, the amount of ejected mass, and the CSM distribution at the time of core-collapse. The energy injection at the bottom of the envelope of a massive star within a period shorter than the dynamical timescale of the envelope could reproduce some observed optical outbursts prior to the core-collapse and form the CSM, which can power an interaction supernova (SN) classified as type IIn.

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

A numerical light curve model for interaction-powered supernovae

We construct a numerical light curve model for interaction-powered supernovae that arise from an interaction between the ejecta and the circumstellar matter (CSM). In order to resolve the shocked region of an interaction-powered supernova, we solve the fluid equations and radiative transfer equation assuming the steady states in the rest frames of the reverse and forward shocks at each time step. Then we numerically solve the radiative transfer equation and the energy equation in the CSM with the thus obtained radiative flux from the forward shock as a radiation source. We also compare results of our models with observational data of two supernovae 2005kj and 2005ip classified as type IIn and discuss the validity of our assumptions. We conclude that our model can predict physical parameters associated with supernova ejecta and the CSM from the observed features of the light curve as long as the CSM is sufficiently dense. Furthermore, we found that the absorption of radiation in the CSM is an important factor to calculate the luminosity.

preprint2019arXiv

An optical search for transients lasting a few seconds

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

preprint2019arXiv

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

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

preprint2019arXiv

The optically thick rotating magnetic wind from a massive white dwarf merger product

WD J005311 is a newly identified white dwarf (WD) in a mid-infrared nebula. The spectroscopic observation indicates the existence of a neon-enriched carbon/oxygen wind with a terminal velocity of $v_{\infty,\rm obs}\sim 16,000\,\rm km\,s^{-1}$ and a mass loss rate of $\dot M_{\rm obs}\sim 3.5\times 10^{-6}\,M_\odot$ yr$^{-1}$. Here we consistently explain the properties of WD J005311 using a newly constructed wind solution, where the optically thick outflow is launched from the carbon burning shell on an oxygen-neon core and accelerated by the rotating magnetic field to become supersonic and unbound well below the photosphere. Our model implies that WD J005311 has a mass of $M_* \sim 1.1\mbox{-}1.3\,M_\odot$, a magnetic field of $B_* \sim (2\mbox{-}5)\times 10^7\,\rm G$, and a spin angular frequency of $Ω\sim 0.2\mbox{-}0.5 \,\rm s^{-1}$. The large magnetic field and fast spin support the carbon-oxygen WD merger origin. WD J005311 will neither explode as a type Ia supernova nor collapse into a neutron star. If the wind continues to blow another few kyr, WD J005311 will spin down significantly and join to the known sequence of slowly-rotating magnetic WDs. Otherwise it may appear as a fast-spinning magnetic WD and could be a new high energy source.

preprint2016arXiv

2D radiaition-hydrodynamic simulations of supernova shock breakout in bipolar explosions of a blue supergiant progenitor

A two-dimensional special relativistic radiation-hydrodynamics code is developed and applied to numerical simulations of supernova shock breakout in bipolar explosions of a blue supergiant. Our calculations successfully simulate the dynamical evolution of a blast wave in the star and its emergence from the surface. Results of the model with spherical energy deposition show a good agreement with previous simulations. Furthermore, we calculate several models with bipolar energy deposition and compare their results with the spherically symmetric model. The bolometric light curves of the shock breakout emission are calculated by a ray-tracing method. Our radiation-hydrodynamic models indicate that the early part of the shock breakout emission can be used to probe the geometry of the blast wave produced as a result of the gravitational collapse of the iron core.

preprint2016arXiv

Contribution of Neutron Star Mergers to the R-process Chemical Evolution in the Hierarchical Galaxy Formation

The main astronomical source of r-process elements has not yet been identified. One plausible site is neutron star mergers (NSMs), but from perspective of the Galactic chemical evolution, it has been pointed out that NSMs cannot reproduce the observed r-process abundance distribution of metal-poor stars at [Fe/H] $< -3$. Recently, Tsujimoto & Shigeyama (2014) pointed out that NSM ejecta can spread into much larger volume than ejecta from a supernova. We re-examine the enrichment of r-process elements by NSMs considering this difference in propagation using the chemical evolution model under the hierarchical galaxy formation. The observed r-process enhanced stars around [Fe/H] $\sim -3$ are reproduced if the star formation efficiency is lower for low-mass galaxies under a realistic delay time distribution for NSMs. We show that a significant fraction of NSM ejecta escape from its host proto-galaxy to pollute intergalactic matter and other proto-galaxies. The propagation of r-process elements over proto-galaxies changes the abundance distribution at [Fe/H] $< -3$ and obtains distribution compatible with observations of the Milky Way halo stars. In particular, the pre-enrichment of intergalactic medium explains the observed scarcity of EMP stars without Ba and abundance distribution of r-process elements at [Fe/H] $\lesssim -3.5$.

preprint2016arXiv

Hydrodynamical interaction of mildly relativistic ejecta with an ambient medium

Hydrodynamical interaction of spherical ejecta freely expanding at mildly relativistic speeds into an ambient cold medium is studied in semi-analytical and numerical ways to investigate how ejecta produced in energetic stellar explosions dissipate their kinetic energy through the interaction with the surrounding medium. We especially focus on the case in which the circumstellar medium is well represented by a steady wind at a constant mass-loss rate having been ejected from the stellar surface prior to the explosion. As a result of the hydrodynamical interaction, the ejecta and circumstellar medium are swept by the reverse and forward shocks, leading to the formation of a geometrically thin shell. We present a semi-analytical model describing the dynamical evolution of the shell and compare the results with numerical simulations. The shell can give rise to bright emission as it gradually becomes transparent to photons. while it is optically thick. We develop an emission model for the expected emission from the optically thick shell, in which photons in the shell gradually diffuse out to the interstellar space. Then, we investigate the possibility that radiation powered by the hydrodynamical interaction is the origin of an underluminous class of gamma-ray bursts.

preprint2015arXiv

Chemical Feature of Eu abundance in the Draco dwarf spheroidal galaxy

Chemical abundance of r-process elements in nearby dwarf spheroidal (dSph) galaxies is a powerful tool to probe the site of r-process since their small-mass scale can sort out individual events producing r-process elements. A merger of binary neutron stars is a promising candidate of this site. In faint, or less massive dSph galaxies such as the Draco, a few binary neutron star mergers are expected to have occurred at most over the whole past. We have measured chemical abundances including Eu and Ba of three red giants in the Draco dSph by Subaru/HDS observation. The Eu detection for one star with [Fe/H]=-1.45 confirms a broadly constant [Eu/H] of ~-1.3 for stars with [Fe/H]>-2. This feature is shared by other dSphs with similar masses, i.e., the Sculptor and the Carina, and suggests that neutron star merger is the origin of r-process elements in terms of its rarity. In addition, two very metal-poor stars with [Fe/H]=-2.12 and -2.51 are found to exhibit very low Eu abundances such as [Eu/H]<-2 with an implication of a sudden increase of Eu abundance by more than 0.7 dex at [Fe/H] ~ -2.2 in the Draco dSph. In addition, the detection of Ba abundances for these stars suggests that the r-process enrichment initiated no later than the time when only a few % of stars in the present-day Draco dSph was formed. Though an identification of the origin of an early Eu production inside the Draco dSph should be awaited until more abundance data of stars with [Fe/H]<-2 in the Draco as well as other faint dSphs become available, the implied early emergence of Eu production event might be reconciled with the presence of extremely metal-poor stars enriched by r-process elements in the Galactic halo.

preprint2015arXiv

On the Brightness of Surviving Companions in Type Ia Supernova Remnants

The progenitor systems for type Ia supernovae are still controversial. One of the methods to test the proposed scenario for the progenitor systems is to identify companions that are supposed to survive according to the so-called single degenerate scenario. These companions might be affected by supernova ejecta. We present several numerical simulations of surviving red-giant companions whose envelopes were stripped and heated. We find that red-giants with less-massive helium cores ($\lesssim0.30\,M_{\odot}$) can be so faint after the supernovae that we cannot detect them. In addition, we apply the results to the case of SNR 0509-67.5, and put constraints on the helium core mass, envelope stripping, and energy injection under the single degenerate scenario for type Ia supernovae.

preprint2015arXiv

Revealing progenitors of type Ia supernovae from their light curves and spectra

In the single degenerate (SD) scenario of type Ia supernovae (SNe Ia), the collision of the ejecta with its companion results in stripping hydrogen rich matter from the companion star. This hydrogen rich matter might leave its trace in the light curves and/or spectra. In this paper, we perform radiation hydrodynamical simulations of this collision for three binary systems. As a result, we find that the emission from the shock-heated region is not as strong as in the previous study. This weak emission, however, may be a result of our underestimate of the coupling between the gas and radiation in the shock interaction. Therefore, though our results suggest that the observed early light curves of SNe Ia can not rule out binary systems with a short separation as the progenitor system, more elaborate numerical studies will be needed to reach a fair conclusion. Alternatively, our results indicate that the feature observed in the early phase of a recent type Ia SN 2014J might result from interaction of the ejecta with a red giant companion star. We also discuss the dependence of spectral features of H alpha and Si II absorption lines on viewing angles and investigate whether they can constrain the event rate of the SD progenitor.

preprint2015arXiv

Spherical explosion with central energy source

We present a novel semi-analytic solution describing the propagation of a spherical blast wave driven by a central energy source. The initial density profile has a power-law function of the distance from the center and the energy is injected only into the central region at a rate given by a power-law function of time. This solution is composed of three regions separated by the contact surface and the shock front. The innermost region is assumed to be uniform and the outside of the contact surface includes the shocked matter described by self-similar solutions. We analytically derive the applicable range of parameters of this solution from requirements to satisfy boundary conditions. We perform numerical simulations for flows with various values of parameters, some of which reside out of the thus derived applicable range, and compare with the semi-analytic solutions.

preprint2014arXiv

An Investigation into Surface Temperature Distributions of High-B Pulsars

Bearing in mind the application to high-magnetic-field (high-B) radio pulsars, we investigate two-dimensional (2D) thermal evolutions of neutron stars (NSs). We pay particular attention to the influence of different equilibrium configurations on the surface temperature distributions. The equilibrium configurations are constructed in a systematic manner, in which both toroidal and poloidal magnetic fields are determined self-consistently with the inclusion of general relativistic effects. To solve the 2D heat transfer inside the NS interior out to the crust, we have developed an implicit code based on a finite-difference scheme that deals with anisotropic thermal conductivity and relevant cooling processes in the context of a standard cooling scenario. In agreement with previous studies, the surface temperatures near the pole become higher than those in the vicinity of the equator as a result of anisotropic heat transfer. Our results show that the ratio of the highest to the lowest surface temperatures changes maximally by one order of magnitude, depending on the equilibrium configurations. Despite such difference, we find that the area of such hot and cold spots is so small that the simulated X-ray spectrum could be well reproduced by a single temperature blackbody fitting.

preprint2014arXiv

Dynamical Evolution of an Ultra-relativistic Fireball Colliding with a Freely Expanding Gas

We investigate the hydrodynamical evolution of an ultra-relativistic fireball colliding with a freely expanding gas. The hydrodynamical interaction of the fireball and the gas results in the formation of a geometrically thin shell. We study the dynamical evolution of the shell by an analytical way and perform a numerical simulation equipped with an adaptive mesh refinement to investigate the internal structure of the shell. The shocked gas can give rise to bright emission in the X-ray and gamma-ray energy range. We propose that the breakout emission from the forward shock and the photospheric emission from the reverse-shocked fireball contribute to early gamma-ray emission from gamma-ray bursts.

preprint2014arXiv

Enrichment history of r-process elements shaped by a merger of neutron star pairs

The origin of r-process elements remains unidentified and still puzzles us. The recent discovery of evidence for the ejection of r-process elements from a short-duration gamma-ray burst singled out neutron star mergers (NSMs) as their origin. In contrast, core-collapse supernovae are ruled out as the main origin of heavy r-process elements (A>110) by recent numerical simulations. However, the properties characterizing NSM events - their rarity and high yield of r-process elements per event - have been claimed to be incompatible with the observed stellar records on r-process elements in the Galaxy. We add to this picture with our results, which show that the observed constant [r-process/H] ratio in faint dwarf galaxies and one star unusually rich in r-process in the Sculptor galaxy agree well with this rarity of NSM events. Furthermore, we found that a large scatter in the abundance ratios of r-process elements to iron in the Galactic halo can be reproduced by a scheme that incorporates an assembly of various protogalactic fragments, in each of which r-process elements supplied by NSMs pervade the whole fragment while supernovae distribute heavy elements only inside the regions swept up by the blast waves. Our results demonstrate that NSMs occurring at Galactic rate of 12-23 per Myr are the main site of r-process elements, and we predict the detection of gravitational waves from NSMs at a high rate with upcoming advanced detectors.

preprint2014arXiv

Origin of Ne Emission Line of Very Luminous Soft X-ray Transient MAXI J0158$-$744

We investigate the mechanism to reproduce notable spectral features at the ignition phase of nova explosion observed for a super-Eddington X-ray transient source MAXI J0158$-$744 in the Small Magellanic Cloud. These are a strong Ne IX emission line at 0.92 keV with a large equivalent width of $0.32^{+0.21}_{-0.11}$ keV and the absence of Ne X line at 1.02 keV. In this paper, we calculate the radiative transfer using a Monte Carlo code, taking into account the line blanketing effect due to transitions of N, O, Ne, Mg and Al ions in an accelerating wind emanating from a white dwarf with a structure based on a spherically symmetric stationary model. We found that the strong Ne IX line can be reproduced if the mass fraction of Ne is enhanced to $10^{-3}$ or more and that of O is reduced to $\sim5\times10^{-9}$ or less and that the absence of other lines including Ne X ions at higher energies can be also reproduced by the line blanketing effect. This enhancement of the Ne mass fraction indicates that the ejecta are enriched by Ne dredged up from the surface of the white dwarf composed of O, Ne, and Mg rather than C and O, as already pointed out in the previous work. We argue that the CNO cycle driving this nova explosion converted most of C and O into N and thus reduced the O mass fraction.

preprint2014arXiv

Signatures of A Companion Star in Type Ia Supernovae

While type Ia Supernovae (SNe Ia) have been used as precise cosmological distance indicators, their progenitor systems remain unresolved. One of the key questions is if there is a non-degenerate companion star at the time of a thermonuclear explosion of a white dwarf (WD). In this paper, we investigate if an interaction between the SN ejecta and the companion star may result in observable footprints around the maximum brightness and thereafter, by performing multi-dimensional radiation transfer simulations based on hydrodynamic simulations of the interaction. We find that such systems result in variations in various observational characteristics due to different viewing directions, while the predicted behaviors (redder and fainter for the companion direction) are opposite to what were suggested by the previous study. The variations are generally modest and within observed scatters. However, the model predicts trends between some observables different from observationally derived, thus a large sample of SNe Ia with small calibration errors may be used to constrain the existence of such a companion star. The variations in different colors in optical band passes can be mimicked by external extinctions, thus such an effect could be a source of a scatter in the peak luminosity and derived distance. After the peak, hydrogen-rich materials expelled from the companion will manifest themselves in hydrogen lines. H$α$ is however extremely difficult to identify. Alternatively, we find that P$_β$ in post-maximum near-infrared spectra can potentially provide powerful diagnostics.

preprint2014arXiv

The origins of light and heavy r-process elements identified by chemical tagging of metal-poor stars

Growing interests in neutron star (NS) mergers as the origin of r-process elements have sprouted since the discovery of evidence for the ejection of these elements from a short-duration gamma-ray burst. The hypothesis of a NS merger origin is reinforced by a theoretical update of nucleosynthesis in NS mergers successful in yielding r-process nuclides with A>130. On the other hand, whether the origin of light r-process elements are associated with nucleosynthesis in NS merger events remains unclear. We find a signature of nucleosynthesis in NS mergers from peculiar chemical abundances of stars belonging to the Galactic globular cluster M15. This finding combined with the recent nucleosynthesis results implies a potential diversity of nucleosynthesis in NS mergers. Based on these considerations, we are successful in the interpretation of an observed correlation between [light r-process/Eu] and [Eu/Fe] among Galactic halo stars and accordingly narrow down the role of supernova nucleosynthesis in the r-process production site. We conclude that the tight correlation by a large fraction of halo stars is attributable to the fact that core-collapse supernovae produce light r-process elements while heavy r-process elements such as Eu and Ba are produced by NS mergers. On the other hand, stars in the outlier, composed of r-enhanced stars ([Eu/Fe]>+1) such as CS22892-052, were exclusively enriched by matter ejected by a subclass of NS mergers that is inclined to be massive and consist of both light and heavy r-process nuclides.

preprint2013arXiv

Early Thermal X-ray Emission from Long Gamma-ray Bursts and Their Circumstellar Environments

We performed a series of hydrodynamical calculations of an ultra-relativistic jet propagating through a massive star and the circumstellar matter to investigate the interaction between the ejecta and the circumstellar matter. We succeed in distinguishing two qualitatively different cases in which the ejecta are shocked and adiabatically cool. To examine whether the cocoon expanding at subrelativistic speeds emits any observable signal, we calculate expected photospheric emission from the cocoon. It is found that the emission can explain early thermal X-ray emission recently found in some long gamma-ray bursts. The result implies that the difference of the circumstellar environment of long gamma-ray bursts can be probed by observing their early thermal X-ray emission.

preprint2013arXiv

Generation of High-Energy Photons at Ultra-Relativistic Shock Breakout in Supernovae

We present theoretical expectations for non-thermal emission due to the bulk Comptonization at the ultra-relativistic shock breakout. We calculate the transfer of photons emitted from the shocked matter with a Monte Carlo code fully taking into account special relativity. As a hydrodynamical model, we use a self-similar solution of Nakayama & Shigeyama (2005). Our calculations reveal that the spectral shape exhibits a double peak or a single peak depending on the shock temperature at the shock breakout. If it is significantly smaller than the rest energy of an electron, the spectrum has a double peak. We also display a few example of light curves, and estimate the total radiation energy. In comparison with observations of gamma-ray bursts, a part of the higher energy component in the spectra and the total energy can be reproduced by some parameter sets. Meanwhile, the lower energy counterpart in the Band function is not reproduced by our results and the duration time seems too short to represent the entire event of a gamma-ray burst. Therefore the subsequent phase will constitute the lower energy part in the spectrum.

preprint2013arXiv

On the Subclasses in Swift Long Gamma-Ray Bursts: A Clue to Different Central Engines

Analyzing light curves of a complete sample of bright Swift long gamma-ray bursts (LGRBs) of which the peak photon fluxes constructed with the bin width of 1 second in the Swift 15-350 keV energy band exceed 2.6 photons cm$^{-2}$s$^{-1}$, we confirm that there do exist the third class in GRBs in addition to short and long GRBs. Being different from previous works based on the duration, fluence, etc. our classification method is based on two properties both quantified with light curve shapes of the prompt emission: the Absolute Deviation from the Constant Luminosity of their cumulative light curve $ADCL$, and the ratio of the mean counts to the maximum counts $\bar C/C_{\rm max}$. These are independent of the distance and the jet opening angle. A cluster analysis via the Gaussian mixture model detects three subclasses: one consisting of LGRBs with small $ADCL$ and large $\bar C/C_{\rm max}$ values referred to as Type I, one with large $ADCL$ and large $\bar C/C_{\rm max}$ referred to as Type II, and one with intermediate $ADCL$ and small $\bar C/C_{\rm max}$, which is composed of contaminating short GRBs with the extended emission. This result is reinforced by different temporal and spectral indices of their X-ray afterglows. The difference is prominent in the temporal index of the steep decay phase in particular: the indices for Type I LGRBs distribute between $-6$ and $-3$ while those for Type II LGRBs between $-3$ and $-2$. From these properties, we propose a possible scenario with different central engines: an accreting black hole and a magnetar.

preprint2013arXiv

Universal scaling law in long gamma-ray bursts

Overwhelming diversity of long gamma-ray bursts (LGRBs), discovered after the launch of {\it Swift} satellite, is a major obstacle to LGRB studies. Recently, it is shown that the prompt emission of LGRBs is classified into three subclasses: Type I, Type II LGRBs populating separate fundamental planes in a 3D space defined by the peak luminosity, the duration, and the spectral peak energy, and outliers not belonging to either of the planes. Here we show that Type I LGRBs (LGRBs I) exhibit different shapes of light curves from Type II LGRBs (LGRBs II). Furthermore, we demonstrate that this classification has uncovered a new scaling law in the light curves of LGRBs II spanning 8 orders of magnitude from the prompt to late X-ray afterglow emission. The scaled light curve has four distinct phases. The first phase has a characteristic time scale while the subsequent three phases exhibit power law behaviors with different exponents. We discuss its possible interpretation in terms of the emission from an optically thick fireball propagating in the cricumstellar matter at relativistic speeds and argue that the observed four phases correspond to its hydrodynamical phases. Our classification scheme can pin down the intrinsic luminosities of LGRBs II through the scaling law from a sample of polymorphic GRBs. Further refinement of this scheme and scaling law will make a subclass of LGRBs a standard candle as reliable and accurate as Type Ia supernovae in the more distant universe than supernovae can reach.

preprint2012arXiv

Diversity of Type Ia Supernovae Imprinted in Chemical Abundances

A time delay of Type Ia supernova (SN Ia) explosions hinders the imprint of their nucleosynthesis on stellar abundances. However, some occasional cases give birth to stars that avoid enrichment of their chemical compositions by massive stars and thereby exhibit a SN Ia-like elemental feature including a very low [Mg/Fe] (~-1). We highlight the elemental feature of Fe-group elements for two low-Mg/Fe objects detected in nearby galaxies, and propose the presence of a class of SNe Ia that yield the low abundance ratios of [Cr,Mn,Ni/Fe]. Our novel models of chemical evolution reveal that our proposed class of SNe Ia (slow SNe Ia) is associated with ones exploding on a long timescale after their stellar birth, and gives a significant impact on the chemical enrichment in the Large Magellanic Cloud (LMC). In the Galaxy, on the other hand, this effect is unseen due to the overwhelming enrichment by the major class of SNe Ia that explode promptly (prompt SNe Ia) and eject a large amount of Fe-group elements. This nicely explains the different [Cr,Mn,Ni/Fe] features between the two galaxies as well as the puzzling feature seen in the LMC stars exhibiting very low Ca but normal Mg abundances. Furthermore, the corresponding channel of slow SN Ia is exemplified by performing detailed nucleosynthesis calculations in the scheme of SNe Ia resulting from a 0.8+0.6 solar mass white dwarf merger.

preprint2012arXiv

Effects of magnetic fields on the propagation of nuclear flames in magnetic white dwarfs

We investigate effects of the magnetic field on the propagation of laminar flames of nuclear reactions taking place in white dwarfs (WDs) with the mass close to the Chandrasekhar limit. We calculate the velocities of laminar flames parallel and perpendicular to uniform magnetic fields as eigenvalues of steady solutions for magnetic hydrodynamical equations. As a result, we find that even when the magnetic pressure does not dominate the entire pressure it is possible for the magnetic field to suppress the flame propagation through the thermal conduction. Above the critical magnetic field, the flame velocity decreases with increasing magnetic field strength as $v \sim B^{-1}$. In media with densities of $10^{7}, 10^{8}, \,\mathrm{and}\,10^{9} \unit{g\,cm^{-3}}$, the critical magnetic fields are orders of $\sim 10^{10}, 10^{11}, \,\mathrm{and}\,10^{12} \unit{G}$, respectively.

preprint2012arXiv

Feedback effects of aspherical supernovae explosions on galaxies

We investigate how explosions of aspherical supernovae (A-SNe) can influence star formation histories and chemical evolution of dwarf galaxies by using a new chemodynamical model. We mainly present the numerical results of two comparative models so that the A-SN feedback effects on galaxies can be more clearly seen. SNe originating from stars with masses larger than 30M_sun are A-SNe in the "ASN" model whereas all SNe are spherical ones (S-SNe) in the "SSN" model. Each S-SN and A-SN are assumed to release feedback energy of 10^{51} erg and 10^{52} erg, respectively, and chemical yields and feedback energy of A-SN ejecta depend on angles between the axis of symmetry and the ejection directions. We find that star formation can become at least by a factor of ~3 lower in the ASN model in comparison with the SSN one owing to the more energetic feedback of A-SNe. As a result of this, chemical evolution can proceed very slowly in the ASN model. A-SN feedback effects can play a significant role in the formation of giant gaseous holes and energetic gaseous outflow and unique chemical abundances (e.g., high [Mg/Ca]). Based on these results, we provide a number of implications of the A-SN feedback effects on galaxy formation and evolution.

preprint2012arXiv

Type-Ia Supernova Remnant Shell At $Z=3.5$ Seen In The Three Sightlines Toward The Gravitationally Lensed Qso B1422+231

Using the Subaru 8.2m Telescope with an IRCS Echelle spectrograph, we obtained high-resolution (R=10,000) near-infrared (1.01-1.38 μm) spectra of images A and B of the gravitationally lensed QSO B1422+231 (z=3.628) consisting of four known lensed images. We detected MgII absorption lines at z=3.54, which show a large variance of column densities (~ 0.3 dex) and velocities (~ 10 km/s) between the sightlines A and B with a projected separation of only 8.4h_{70}^{-1} pc at the redshift. This is the smallest spatial structure of the high-z gas clouds ever detected after Rauch et al. found a 20-pc scale structure for the same z=3.54 absorption system using optical spectra of images A and C. The observed systematic variances imply that the system is an expanding shell as originally suggested by Rauch et al. By combining the data for three sightlines, we managed to constrain the radius and expansion velocity of the shell (~ 50-100 pc, 130 km/s), concluding that the shell is truly a supernova remnant (SNR) rather than other types of shell objects, such as a giant HII region. We also detected strong FeII absorption lines for this system, but with much broader Doppler width than that of α-element lines. We suggest that this FeII absorption line originates in a localized FeII-rich gas cloud that is not completely mixed with plowed ambient interstellar gas clouds showing other α-element low-ion absorption lines. Along with the Fe richness, we conclude that the SNR is produced by an SNIa explosion.

preprint2010arXiv

Boron Synthesis in Type Ic Supernovae

We investigate the neutrino-process in an energetic Type Ic supernova (SN Ic) and resultant productions of the light elements including boron and its stable isotopes. SN Ic is a very unique boron source because it can produce boron through not only spallation reactions as discussed in Nakamura & Shigeyama (2004) but also the neutrino-process. The neutrino-process is considered to occur in core-collapse supernovae and previous studies were limited to Type II supernovae (SNe II). Although the progenitor star of an SN Ic does not posses a He envelope so that 7Li production via the neutrino-process is unlikely, 11B can be produced in the C-rich layers. We demonstrate a hydrodynamic simulation of SN Ic explosion and estimate the amounts of the light elements produced via the neutrino-process for the first time, and also the subsequent spallation reactions between the outermost layers of compact SN Ic progenitor and the ambient medium. We find that the neutrino-process in the current SN Ic model produces a significant amount of 11B, which is diluted by 10B from spallation reactions to get closer to B isotopic ratios observed in meteorites. We also confirm that high-temperature mu- and tau-neutrinos and their anti-neutrinos, reasonably suggested from the compact structure of SN Ic progenitors, enhance the light element production through the neutral-current reactions, which may imply an important role of SNe Ic in the Galactic chemical evolution.

preprint2010arXiv

Compton degradation of gamma-ray line emission from radioactive isotopes in the classical nova V2491 Cygni

To account for the non-thermal emission from the classical nova V2491 Cygni, we perform a series of numerical calculations of radiative transfer of γ-ray photons from the radioactive isotope $^{22}$Na in the matter ejected from a white dwarf. Using a simple wind model for the dynamical evolution of the ejecta and a monte-carlo code, we calculate radiative transfer of the γ-ray photons in the ejecta. Repeated scattering of the γ-ray photons by electrons in the ejecta, i.e., Compton degradation, results in an extremely flat spectrum in the hard X-ray range, which successfully reproduces the observed spectrum of the X-ray emission from V2491 Cygni. The amount of the isotope $^{22}$Na synthesized in the ejecta is required to be 3\times 10^{-5} M_\odot to account for the flux of the hard X-ray emission. Our model indicates that the ejecta become transparent to the γ-ray photons within several tens days. Using the results, we briefly discuss the detectability of the γ-ray line emission by the {\it INTEGRAL} gamma-ray observatory and the {\it Fermi} Gamma-ray Space Telescope.

preprint2010arXiv

Multiple Main Sequence of Globular Clusters as a Result of Inhomogeneous Big Bang Nucleosynthesis

A new mechanism for enhancing the helium abundance in the blue main sequence stars of omega Centauri and NGC 2808 is investigated. We suggest that helium enhancement was caused by the inhomogeneous big bang nucleosynthesis. Regions with extremely high baryon-to-photon ratios are assumed to be caused by the baryogenesis. Its mass scale is also assumed to be 10^6 Msun. An example of the mechanisms to realize these two things was already proposed as the Affleck-Dine baryogenesis. As the baryon-to-photon ratio becomes larger, the primordial helium abundance is enhanced. We calculated the big bang nucleosynthesis and found that there exists a parameter region yielding enough helium to account for the split of the main sequence in the aforementioned globular clusters while keeping the abundances of other elements compatible with observations. Our mechanism predicts that heavy elements with the mass number of around 100 is enhanced in the blue main sequence stars. We estimate the timescales of diffusion of the enhanced helium and mass accretion in several stages after the nucleosynthesis to investigate whether these processes diminish the enhancement of helium. We found that the diffusion does not influence the helium content. A cloud with a sufficiently large baryon-to-photon ratio to account for the multiple main sequence collapsed immediately after the recombination. Subsequently, the cloud accreted the ambient matter with the normal helium content. If the star formation occurred both in the collapsed core and the accreted envelope, then the resultant star cluster has a double main sequence.

preprint2010arXiv

Non-thermal photon production via bulk comptonization at supernova shock breakout

We investigate roles of the bulk comptonization process in the propagation of thermal photons emitted at the shock breakout of a supernova explosion. We use a hydrodynamical model based on a self-similar solution for the shock breakout. The propagation of photons is treated by a Monte-Carlo technique. Results of the simulations successfully reproduce the power-law spectrum observed by X-ray observations for XRF 080109/SN 2008D, if a mildly relativistic shock velocity is assumed. Effects of some radiative processes, photoionization, radiative recombination, and free-free absorption on the propagation of emitted photons are also investigated. It is found that none of them hardly changes the spectrum regarding the progenitor stars of type Ib or Ic supernovae. Light curves calculated under the assumption of a spherical explosion indicate that the progenitor radius is required to be $/sim 10^{13}$ cm.

preprint2010arXiv

Probing Explosion Geometry of Core-collapse Supernovae with Light Curves of the Shock Breakout

We investigate effects of aspherical energy deposition in core-collapse supernovae on the light curve of the supernova shock breakout. We performed two-dimensional hydrodynamical calculations of an aspherical supernova explosion to obtain the time when a shock wave generated in the stellar interior reaches the stellar surface in each radial direction. Using results of the calculations, light curves during the shock breakout are derived in an approximate way. We show that the light curve during the shock breakout can be a strong indicator of aspherical properties of core-collapse supernovae.

preprint2009arXiv

A conservative scheme for the relativistic Vlasov-Maxwell system

A new scheme for numerical integration of the 1D2V relativistic Vlasov-Maxwell system is proposed. Assuming that all particles in a cell of the phase space move with the same velocity as that of the particle located at the center of the cell at the beginning of each time step, we successfully integrate the system with no artificial loss of particles. Furthermore, splitting the equations into advection and interaction parts, the method conserves the sum of the kinetic energy of particles and the electromagnetic energy. Three test problems, the gyration of particles, the Weibel instability, and the wakefield acceleration, are solved by using our scheme. We confirm that our scheme can reproduce analytical results of the problems. Though we deal with the 1D2V relativistic Vlasov-Maxwell system, our method can be applied to the 2D3V and 3D3V cases.

preprint2009arXiv

Detailed Analysis of Filamentary Structure in the Weibel Instability

We present results of a 2D3V kinetic Vlasov simulation of the Weibel instability. The kinetic Vlasov simulation allows us to investigate the velocity distribution of dilute plasmas, in which the effect of collisions between particles is negligible, and has the advantage that the accuracy of the calculated velocity distribution does not depend on the density of plasmas at each point in the physical space. We succeed in reproducing some features of the Weibel instability shown by other simulations, for example, the exponentially growing phase, the saturation of the magnetic field strength, the formation of filamentary structure, and the coalescence of the filaments. Especially, we concentrate on the behavior of the filaments after the saturation of the magnetic field strength and find that there is a kind of quasi-equilibrium states before the coalescence occurs. Furthermore, it is found that an analytical solution for stationary states of the 2D3V Vlasov-Maxwell system can reproduce some dominant features of the quasi-equilibrium, e.g, the configuration of the magnetic field and the velocity distribution at each point. The analytical expression could give a plausible model for the transition layer of a collisionless shock where a strong magnetic field generated by the Weibel instability provides an effective dissipation process instead of collisions between particles.

preprint2007arXiv

Searching for a Companion Star of Tycho's Type Ia Supernova with Optical Spectroscopic Observations

We report our first results of photometric and spectroscopic observations for Tycho's supernova remnant (SNR Tycho) to search for the companion star of a type Ia supernova (SN Ia). From photometric observations using Suprime-Cam on the Subaru Telescope, we have picked up stars brighter than 22 mag (in $V$-band) for spectroscopy, which are located within a circular region with the radius of 30" around the center of SNR Tycho. If the ejecta of young supernova remnants, such as SNR Tycho, have a sufficient amount of Fe I, we should be able to detect absorption lines at 3720 Åand 3860 Åassociated with transitions from the ground state of Fe I in the spectrum of the companion star. To identify the companion star of a SN Ia using these characteristic absorption lines of Fe I, we made optical low-resolution spectroscopy of their targets using FOCAS on the Subaru Telescope. In our spectroscopic observations, we obtained spectra of 17 stars in the SNR Tycho region and compare them with template stellar spectra. We detect significant absorption lines from two stars at 3720 Å. Since widths of their absorption lines are broad, it is likely that the detected absorptions are due to Fe I in the expanding ejecta of SNR Tycho. However, none of stars exhibits a clear red wing in the observed profiles of the absorption, though a star in the background of the SNR should show it. Hence, we suggest another interpretation that the detected absorption lines might be caused by the peculiarity of stars. A star named Tycho(E) has the absorption line at 3720 Åand its projected position is close to the center of SNR Tycho. Based on our observations, Tycho(E) is a new candidate as the companion star of Tycho's supernova.

preprint2006arXiv

The Origin of Carbon-Enhancement and Initial Mass Function of Extremely Metal-Poor Stars in the Galactic Halo

It is known that the carbon-enhanced, extremely metal-poor (CEMP) stars constitute a substantial proportion in the extremely metal-poor (EMP) stars of the Galactic Halo, by far larger than CH stars in Population II stars. We investigate their origin with taking into account an additional evolutionary path to the surface carbon-enrichment, triggered by hydrogen engulfment by the helium flash convection, in EMP stars of $[Fe/H] \lesssim -2.5$. This process is distinct from the third dredge-up operating in more metal-rich stars and also in EMP stars. In binary systems of EMP stars, the secondary stars become CEMP stars through mass transfer from the primary stars of low and intermediate masses, which have developed the surface carbon-enhancement. Our binary scenario can predict the variations in the abundances not only for carbon but also for nitrogen and s-process elements and reasonably explain the observed properties such as the stellar distributions with respect to the carbon abundances, the binary periods, and the evolutionary stages. Furthermore, from the observed frequencies of CEMP stars with and without s-process element enhancement, we demonstrate that the initial mass function of EMP stars need to give the mean mass $~10\msun$ under the reasonable assumptions on the distributions of orbital separations and mass ratio of binary components. This also indicates that the currently observed EMP stars were exclusively born as the secondary members of binaries, making up $\sim 10%$ remnants of EMP binary systems of mass $~10^8\msun$ in total; in addition to CEMP stars with white dwarf companions, a significant fraction of them have experienced supernova explosions of their companions. We discuss the implications of the present results in relation to the formation of Galactic halo.