Formation and evolution of dustin hydrogen-poor supernovae
Takaya Nozawa(IPMU, Univ. of Tokyo)
Collaborators;
T. Kozasa (Hokkaido Univ.), N. Tominaga (Konan Univ.),
K. Maeda (IPMU), H. Umeda (Univ. of Tokyo),
K. Nomoto (IPMU/Univ. of Tokyo)

1-1. Introduction (1)SNe are the major sources of interstellar dust?
・ mass estimate of dust condensing in the ejecta➔ Mdust < 10-3 Msun from observations of SNe(e.g., Miekle et al. 2007; Kotak et al. 2009)➔ Mdust = 0.1-1 Msun from theoretical studies(Todini & Ferrara 2001; Nozawa et al. 2003)

It is essential to reveal the composition, size, and massof dust by comparing theoretical model with observation－ a limited number of observations of dust forming-SNe－ a lack of sophisticated radiative transport model1-2. Introduction (2)・ Comparison of models with IR observations of SNRs
－ erosion by sputtering and stochastic heating－ envelope-poor SNe are more favorable(For type II-P SNe, the reverse shock collides with the He
core after a few thousand years)

・ Dust formation and evolution in SN with no envelope－ ~30 % of CCSNe explode as SNe Ib/c and SN IIb(Prieto et al. 2009; Smartt et al. 2008; Boissier & Prantzos 2009)・ investigate the composition, size, and mass of dust
formed in hydrogen-deficient SNe・ apply the results to the IR observations of young SNR
such as Cas A

2-1. Calculation of dust formation○ nucleation and grain growth theory (Kozasa & Hasegawa 1987)
steady-state nucleation rate

grain growth rate

・ key species:a gas species with the least collision frequency among
reactants
・ sticking probability; α=1
・ Tdust = Tgas (dust temperature is the same as that of gas)

2-2. Model of Dust formation calculation○ 1-D models of SNe➔ ☓ clumpy structure and asymmetry explosion○ formation of molecules (CO and SiO)

－ complete－ taking account of formation and destruction(formation efficiency of molcules is 0.01-0.1)○ time evolution of gas temperature－ following a light curve calculation－ adopting adiabatic approximation○ mixing of elements within the He-core－ unmixed case (original onion-like composition)－ uniformly mixed case

2-3. Model of Type IIb SN○ SN IIb model (SN1993J-like model)－ Meje = 2.94 Msun
MH-env = 0.08 Msun
Mstar = 18 Msun
－ E51 = 1－ M(56Ni) = 0.07 Msun

2-4. Composition and mass of dust formed
Mass of dust formed

Total mass of dust formed :0.167 Msun in SN IIb
0.1-1 Msun in SN II-P

condensation time・ different species of dust cancondense in different layers

・ condensation time:300-700 days

2-5. Radius of dust formed in the ejecta
average radius

SN II-PSN IIb

SN Ib
Grain radius➔ > 0.01 μm for SN IIP➔ < 0.01 μm for SN IIb
Dust grains formed in Hdeficient SNe are small

Nozawa et al. (2008)

3-1. Dust destruction calculation○ time evolution of shock wave
－ numerically (Nozawa et al. 2007)－ semi-analytically (Bianchi & Schneider 2007)○ deceleration of dust due to the gas drag－ inversely proportional to grain size and bulk density○ erosion of dust by sputtering－ da / dt ～ 10-6 nH μm yr-1 cm3 for T > 106 K○ ambient gas density－ nH = 0.1, 1.0 and 10.0 cm-3・ contribution from ISM dust?
・ unshocked ejecta-dust?

3-2. Calculation of dust evolution in SNRs
○ Model of calculations(Nozawa et al. 2006, 2007)

・ejecta model－ hydrodynamic model fordust formation calculation
・CSM gas density－ nH = 1.0 and 10.0 cm-3
－ nH(r) ∝ M / (4πr2 vw) g/cm-3(M = 2x10-5 Msun/yr)
・treating dust as a test particle－ erosion by sputtering－ deceleration by gas drag－ collsional heating3-3. Evolution of dust in Type IIb SNR
n=1.0 /cc

n=10.0 /cc

3-4. Time evolution of dust mass
n=1.0 /cc

n=10.0 /cc

Mdust ~ 10-4 Msun at 105 yr

Mdust = 0 Msun at 105 yr

Newly formed dust grains in the ejecta are completely
destroyed in the shocked gas within the SNR
Core-collapse SNe with low-mass outer envelope cannot
be main sources of dust

3-5. Time evolution of IR thermal emission (1)n=1.0 /cc

Data: Hines et al. (2004)
red: with SH
green: without SH

3-6. Time evolution of IR thermal emission (2)n=10.0 /cc

Data: Hines et al. (2004)
red: with SH
green: without SH

3-7. Contribution from unshocked dustdM/dt = 2x10-5 Msun/yr

Md,warm ~ 0.006 Msun,Md,cool ~ 0.08 Msun
dM/dt = 8x10-5 Msun/yr
ρx4

4-1. Dust formation calculation for SN Ia○ Type Ia SN modelW7 model (C-deflagration)(Thielemann et al. 1986)－ Mpr = 1.38 Msun

－ E51 = 1
－ M(56Ni) = 0.6 Msun

○ Dust formation theory－ non-steady nucleation andgrain growth theory(Nozawa et al. 2003)－ onion-like composition－ sticking probability; αs = 1

4-2. Results of dust formation calculation
Condensation time

Average radius of grain

Condensation time of dust : 100-300 days
Average radius of dust : < 0.01 μm

4-3. Mass of dust formed in SN Ia
Mass of dust formed
Optical depth at 300 days

Total mass of dust
Mdust =0.12 Msun

τ(0.55) ~ 100τ(0.55) ~ 60 by C grains
τ(0.55) ~ 35 by Si and FeS➔ too highSummary1) The radius of dust formed in the ejecta of Type IIb SN isquite small (＜0.01 μm) because of low ejecta density
2) Small dust grains formed in Type IIb SN cannot survive
destruction in the shocked gas within the SNR3) Model of dust destruction and heating in Type IIb SNR
to reproduce the observed SED of Cas A isMd,warm = 0.005-0.007 Msun, Md,cool = 0.07-0.09 Msun
dM/dt = 6-8x10-5 Msun/yr
4) IR SED reflects the destruction and stochastic heating➔ properties (size and composition) of dust➔ density structure of circumstellar medium