Mon. Not. R. Astron. Soc. 412, 1070–1080 (2011)

doi:10.1111/j.1365-2966.2010.17960.x

Supernova dust for the extinction law in a young infrared galaxy at z ∼ 1
K. Kawara,1 H. Hirashita,2 T. Nozawa,3 T. Kozasa,4 S. Oyabu,5 Y. Matsuoka,5T. Shimizu,1 H. Sameshima1 and N. Ienaka1
1 Institute of Astronomy, the University of Tokyo, Osawa 2-21-1, Mitaka, Tokyo 181-0015, Japan
2 Institute of Astronomy and Astrophysics, Academia Sinica, PO Box 23-141, Taipei 10617, Taiwan
3 Institute for the Physics and Mathematics of the Universe, the University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa 277-8568, Japan
4 Department of Cosmosciences, Hokkaido University, Sapporo 060-0810, Japan
5 Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8602, Japan

Accepted 2010 November 2. Received 2010 November 1; in original form 2010 March 1

ABSTRACT

We apply the supernova (SN) extinction curves to reproduce the observed properties of
SST J1604+4304 which is a young infrared (IR) galaxy at z ∼ 1. The SN extinction curves
used in this work were obtained from models of unmixed ejecta of Type II supernovae for
the Salpeter initial mass function with a mass range from 8 to 30 M or 8 to 40 M . The
effect of dust distributions on the attenuation of starlight is investigated by performing the χ 2
fitting method against various dust distributions. These are the commonly used uniform dust
screen, the clumpy dust screen and the internal dust geometry. We add to these geometries
three scattering properties, namely, no scattering, isotropic scattering and forward-only scattering. Judging from the χ 2 values, we find that the uniform screen models with any scattering
property provide good approximations to the real dust geometry. Internal dust is inefficient
to attenuate starlight and thus cannot be the dominant source of the extinction. We show that
the SN extinction curves reproduce the data of SST J1604+4304 comparable to or better thanthe Calzetti extinction curve. The Milky Way extinction curve is not in satisfactory agreementwith the data unless several dusty clumps are in the line of sight. This trend may be explainedby the abundance of SN-origin dust in these galaxies; SN dust is the most abundant in theyoung IR galaxy at z ∼ 1, abundant in local starbursts and less abundant in the Galaxy. If dust
in SST J1604+4304 is dominated by SN dust, the dust production rate is ∼0.1 M per SN.
Key words: supernovae: general – dust, extinction – galaxies: evolution – galaxies: ISM –galaxies: starburst – infrared: galaxies.

1 I N T RO D U C T I O NDust is a minor constituent in the Universe , but it plays crucial roles
for evolution of various objects; for example, it cools gas clouds byradiating far-infrared (far-IR) emission, and acts as a catalyzer toform molecules. Corrections for dust extinction have significant
cosmological consequences for determining the Hubble constant,the cosmological model and the epoch of galaxy formation.
It is generally accepted that the atmospheres of asymptotic giant
branch (AGB) stars are one of the main sources of interstellar dustin the Milky Way (MW). However, AGB stars are too old to accountfor the presence of dust at z > 6 (but see Valiante et al. 2009), andsupernovae (SNe) have been recognized as a candidate of the major
source of dust in the early stage of galaxy evolution.

Stars with masses 8–40 M end up their life as Type II supernovae (SNe II)1 in ≤40 Myr since their formation (Heger et al.
2003). The formation of dust in SN ejecta was first witnessed in
SN 1987A, motivating Kozasa, Hasegawa & Nomoto (1989) andKozasa, Hasegawa & Nomoto (1991) to predict the condensed dust
composition and size. Systematic studies of dust formation in SN
ejecta have been proceeded to explore the role of dust in the formation and evolution of stars and galaxies in the early universe (Todini& Ferrara 2001; Nozawa et al. 2003; Schneider, Ferrara & Savaterra
2004): the composition of synthesized dust reflects the elemental
composition in the formation site which is affected by the type anddegree of mixing and the degree of formation efficiency of CO

1 In this paper, we simply use the term SNe II for core-collapsed SNe whose E-mail: kkawara@ioa.s.u-tokyo.ac.jp

progenitors are massive stars. We also call the dust produced by SNe II ‘SN
dust’, and the extinction curve synthesized for the SN dust ‘SN extinction
curve’.
C 2010 The Authors
C 2010 RAS
Monthly Notices of the Royal Astronomical Society Supernova dust at z ∼ 1
molecule, and the size of dust is controlled by the time evolution ofgas density and temperature depending on the type of SNe (Nozawa
et al. 2008, 2010).
Apart from the differences in the SN models and the underlying assumptions employed for dust formation calculations, thepredicted dust mass formed in the ejecta reaches 0.1 to 1 Mper SN II, which is sufficient enough to explain a vast amountof dust of >108 M discovered in host galaxies of high-redshifted
quasars [quasi-stellar objects (QSOs)] at z > 5 (e.g. Morgan &Edmunds 2003; Maiolino et al. 2006; Dwek., Galliano & Jones
2007). However, the mass and composition of dust formed in theejecta have been still controversial; the observations of nearby SNe
have claimed dust mass formed in the ejecta to be less than 10−3 M(see Kozasa et al. 2009 for a review), while the infrared and submillimetre observations of Galactic supernova remnants (SNRs) havesuggested the presence of dust reaching 0.02 to 1 M (Dunne et al.
2003, 2009; Rho et al. 2008; Gomez et al. 2009; Nozawa et al. 2010;
Temim et al. 2010). Although no information has been available on
the composition of dust formed in the ejecta except for SN 2004et(Kotak et al. 2009 ), the spectral mapping observation covering
the entire of Cas A SNR with the Spitzer Space Telescope have revealed the presence of a variety of dust species whose compositionis closely associated with the gas emission-line characteristic ofnucleosynthetic layer in the progenitor (Rho et al. 2008).
The amount, composition and size of dust supplied from SNe to
the interstellar medium (ISM) are not always the same as thosecondensed in the ejecta due to the processing in a hot plasma
produced by the reverse and forward shocks during the journeyto the ISM (Bianchi & Schneider 2007; Nozawa et al. 2007)
as well as in the hot/warm region of ISM (Jones et al. 1994;
Jones, Tielens & Hollenbach 1996; Nozawa, Kozasa & Habe 2006;
Hirashita et al. 2010). Nevertheless, Maiolino et al. (2004) foundthat the extinction curve of SDSS J1048+4647, a broad absorption line (BAL) QSO at z = 6.2 is different from those of nearby
galaxies, and is in excellent agreement with the extinction curve obtained from SN dust models by Todini & Ferrara (2001) considering
the uniformly mixed elemental composition in the ejecta (see also
Gallerani et al. 2010). Synthesizing dust extinction curves for various SN dust models by Nozawa et al. (2003), Hirashita et al. (2005)
have shown that the extinction curve observed in this BAL QSO canbe reasonably reproduced by the dust produced in SNe II without
any mixing in the ejecta. Furthermore, the extinction curve derived
from the analysis of the afterglow of γ -ray burst GRB 050904 at
z = 6.3 have suggested that SNe are the dominant source of dustin the host galaxy (Stratta et al. 2007), although Zafar et al. (2010)
have claimed no evidence of dust extinction in the afterglow. However, the observed extinction curve in the local environments may
provide less constraint on the SN dust model as well as the role ofSNe on the evolution of galaxies.
Young galaxies, in which no AGB stars evolved off the main
sequence, are obvious targets to study the properties of SN dust,
for example, extinction curves that determine the relation betweenchanges in colour and total absorption. However, it is generally
difficult to infer the amount of the dust obscuration, because different geometries of the dust spatial distribution result in different
observed colours for the same amount of dust and the same intrinsic spectrum of the underlying stellar population (e.g. Giavalisco2002). In fact, some studies of H I recombination lines in starburst
and dusty Seyfert galaxies suggest that the extinction observed in
the near-IR is generally greater than that observed by using theoptical lines (Kawara, Nishida & Phillips 1989; Puxley & Brand
1994).
C 2010 The Authors, MNRAS 412, 1070–1080
C 2010 RAS
Monthly Notices of the Royal Astronomical Society 

1071

Figure 1. UV to far-IR SED for SST J1604+4304 compared withthose of starburst-dominated ULIRGs, Arp 220(solid: local ULIRG),
MIPS J142824.0+352619(dash: z = 1.3 hyperluminous IR galaxy). Arp 220
data come from Carico et al. (1990), Rigopoulou, Lawrence & RowanRobinson (1996) and Klaas et al. (1997), and MIPS J142824.0+352619
from Desai et al. (2006) and Borys et al. (2006).

Calzetti, Kinney & Storchi-Bergmann (1994) have shown thatthe extinction laws in the MW and Large Magellanic Cloud cannot account for the ultraviolet (UV) to optical spectra of UV-brightstarburst galaxies in the local Universe in a satisfactory way. They
empirically derived the extinction law for these galaxies. Their extinction law is characterized by an overall slope which is more greyor flatter than the MW extinction law’s slope and by the absenceof the 2175 Å dust feature. In starburst galaxies, numerous SN II
explosions occur and SN-condensed dust would be injected into theISM. Thus, newly formed SN dust would produce a great impact
on the extinction law in local starbursts.
In this paper, we apply synthetic extinction curves on the unmixed
ejecta models of SNe II by Hirashita et al. (2005) to the broadband spectral energy distribution (SED) of a young ultraluminous
infrared galaxy (ULIRG) SST J1604+4304 at z ∼ 1. Comparing
with the MW and Calzetti (2001) extinction laws, we show that the
SN extinction curves fit better to the observations. We adopt H 0 =70 km s−1 Mpc−1 , m = 0.3 and  = 1 − m throughout thispaper.
2 DATA
Synthetic spectra of stellar populations reddened by dust are modelled to reproduce the data of SST J1604+4304. This is a young IR
galaxy at z = 1.135 with the characteristic of ULIRGs. The multiwavelength data from the X-ray to radio are compiled in Kawara
et al. (2010). There are no supportive evidence for active galactic
nucleus (AGN) activities in this galaxy, although the data are notdeep enough to rule out the presence of an AGN.
Fig. 1 illustrates the broad-band SED from the UV to far-IR2
along with representative SED templates of starburst-dominated
ULIRGs, namely, Arp 220 (local ULIRG) and MIPS J142824.0+
352619 (hyperluminous IR galaxy at z = 1.3). The broad-band SED
and stellar absorption lines, such as Ca II and H I lines, as well as the[O II]λ3727 line suggest that young stellar populations dominate the2 The data are taken from table 1 in Kawara et al. (2010) and the correction

for the gravitational amplification is applied.

1072

K. Kawara et al.

optical luminosity in this galaxy. The Ca II index, defined by Kawara
et al. (2010), is D(Ca II) = 0.03 ± 0.11 and the equivalent [O II] width
is EW([O II]) = 31.1 ± 4.3 Å in the rest frame which corresponds
to the observed [O II] line flux3 of 1.9 ± 0.3 × 10−17 erg cm−2 s−1 .
Kawara et al. (2010) estimated the reddening to be E(B − V) ∼ 0.8
when the Calzetti extinction law is applied. The 8–1000 μm bolometric luminosity inferred from the three Spitzer MIPS bands at the24, 70 and 160 μm is Lir = 1.78 ± 0.63 × 1012 L , corresponding
to the flux F ir = 1.00 ± 0.35 × 10−12 erg cm−2 s−1 . These values are
intrinsic to the galaxy after correcting for a gravitational amplification of 1.17 by the foreground galaxy cluster CL 1604+4304 at z =0.9. The bolometric luminosity agrees with the extinction-corrected
stellar luminosity, suggesting SST J1604+4304 can be characterized by models of pure star formation with no AGN emission. The
dust mass inferred from the far-IR data is ∼2 × 108 M , sufficient
to form a shell surrounding the galaxy and produce the observed
extinction.
The SED data at 13 photometric bands from B to 5.8 μm areused for the following analysis; the longest Spitzer IRAC 8.0-μm
band is not used, because the dust emission may be greater than thestarlight in this band corresponding to 3.7 μm in the rest frame.
3 EXTINCTION

The attenuation of starlight produced by dust in front of an extended source not only depends on the properties of dust grains,such as cross-sections, albedo and amount, but also depends on thegeometry of spatial distribution of dust. It is known that the apparent optical depth is lower than the optical depth averaged overthe source area when starlight is attenuated by a thick, clumpy dust
layer or dust is internal to the extended source. We consider the following geometries of dust distribution. For full details, the reader
may refer to Code (1973), Natta & Panagia (1984) and Calzetti et al.
(1994).
Internal dust. In the case that absorbers (dust grains) and emitters(stars) are uniformly distributed, the apparent extinction is given by
fν =

with
−τeff

τap = − ln[(1 − eDust screen. In this geometry, the absorbers are physically separated and located in a screen in front of the emitters. We adopt ananalytic expression derived by Natta & Panagia (1984) for a systemof clumps distributed according to Poisson statistics:fν =3.1 Geometry of dust distribution

0
in fν

Figure 2. The apparent optical depth τ ap is shown as a function of theeffective optical depth τ eff . The solid lines represent dust screen geometries
for N = [1, 2, 4, 8, 16, 32, 64, 128, 256, 512, ∞] from bottom to top. The
dashed line represents an internal dust geometry.

−τap

in = e,)/τeff ],(1)(2)

τeff = (1 − γ )τext

for g = 1,(3)τeff = (1 − γ )0.75 τext

for g = 0,(4)

where τ ap is the apparent optical depth of the internal dust, τ eff theeffective absorption optical depth, τ ext the sum of the optical depth
for absorption and scattering, γ the albedo, f 0ν is the unobscured,
intrinsic flux density and f ν is the reddened, observed flux density.
g ≡ cos (θ ) is the scattering asymmetry factor, where g = 1 forthe forward-only scattering and g = 0 for the isotropic scattering.
Sometimes τ eff = (1 − γ )0.5 τ ext is used for the isotropic scattering(i.e. Natta & Panagia 1984; Calzetti et al. 1994). However, τ eff =(1 − γ )0.75 τ ext approximates better to the solution by the two-stream
approximation for the geometry with a uniform distribution of theabsorbers and the emitters. For the full details, refer to Fig. A1 in
Appendix A.

3 The [O II] line flux published in Kawara et al. (2010) is incorrect.

0
sc fν

−τap
,
sc = ewithτap = N (1 − e−τeff /N ),

(5)
(6)

where N is the average number of clumps along the line of sight, andall clumps are assumed to have one and the same effective opticaldepth τ eff /N. For N → ∞, τ ap → τ eff , corresponding to the uniform
dust screen. For the dust screen which is physically distant from thesystem of stars, the effect of the dust is to remove photons from theline of sight through absorption and scattering. On the other hand,for the dust screen which is located close to the system of stars, thescattered photons can get in the line of sight. Thus, τ eff is given asτeff = τext

for no scattering,(7)

τeff = (1 − γ )τext

for scattering with g = 1,(8)τeff = (1 − γ )0.5 τext

for scattering with g = 0.

(9)Again, g = 1 is for the forward-only scattering, and g = 0 for theisotropic scattering.
Composite geometry. It is quite likely that dust grains distributewithin the system of stars as internal dust and in the screen in frontof the system. We approximate this case asfν =0
sc in fν .

(10)

Fig. 2 shows the apparent optical depth τ ap as a function of τ eff . It
is apparent that the growth of τ ap turns over at a smaller value of τ eff
in the dust screen with a smaller value of N. As shown in this figure,the extinction by dust in the screen with a small N or by internal
dust strongly depends on the geometry and tends to produce a grey
extinction, and dust in a uniform screen produces the most effective
extinction.
3.2 Extinction curves
We examine two extinction curves synthesized in ejecta of SNe II:
the empirical extinction curve derived by Calzetti (2001), and the
MW extinction curve given by Draine (2003). Fig. 3 plots theseextinction curves as a function of wavelengths, showing that theC 2010 The Authors, MNRAS 412, 1070–1080
C 2010 RAS
Monthly Notices of the Royal Astronomical Society Supernova dust at z ∼ 1

Figure 3. Synthetic SN extinction curve compared with those for the MW(Draine 2003) and local starbursts (Calzetti 2001). The curves are normalized
to the extinction (absorption + scattering) in magnitude at 0.8 μm. Note thata spike on the SN curve is produced by crystal SiO2 at 8.7 μm.

SN extinction curves are more grey or flatter than the Calzetti and
MW extinction curves from 0.1 to 1 μm.
The MW extinction curve exhibits a pronounced feature at
2175 Å. This feature is absent in the typical extinction curves of the
Small Magellanic Clouds (SMC). The dust in local starbursts andLyman break galaxies exhibits the properties of the SMC-like dust
properties, while Noll et al. (2009) claimed that about 30 per cent of
the spectra of actively star-forming, intermediate-mass galaxies at
1 < z < 2.5 exhibit broad absorption features at ∼2175 Å. To checkthe effect of the feature on the extinction law in SST J1604+4303, in
addition to the MW extinction curve by Draine (2003), we will examine the pseudo-MW extinction curve which is made by removing
the 2175-Å feature from the MW extinction curve.
The Calzetti (2001) extinction curve was derived by analysing
local starbursts and blue compact dwarf galaxies (Calzetti et al.
1994), and is widely used to, for example, derive photometric redshifts of high-redshift galaxies. Calzetti (2001) gives the extinction
curve in a narrow range from 0.12 to 2.2 μm. To analyse multiwavelength data of high-redshift objects from the UV to IR suchas given in Fig. 1, the extinction curve should be extended to longward of 2.2 μm. The functional form given by Calzetti (2001) cannot be used to extend the extinction curve, because the extinction
value sharply drops beyond 2.2 μm. The Calzetti extinction curveis almost identical to the Draine (2003)’s MW curve from 0.6 to
1.6 μm. We thus extend the Calzetti curve to longer wavelengths
by adopting the MW extinction curve for wavelengths ≥1.6 μm. In
addition, because Calzetti (2001) does not provide the albedo, we
adopt the albedo given in Draine (2003) with removing the dip in
the albedo which is attributed to the familiar 2175-Å feature. This
curve is called the Calzetti extinction curve in this paper. It is quite
likely that the Calzetti extinction curve reflects the dust composition
from the two origins, namely dust condensed in expanding atmospheres of AGB stars and ejecta of SNe II, because in local starburst
galaxies low-mass stars have already evolved off the main sequence
and reached to the dust-forming stages. Because the Calzetti curve
was empirically obtained by assuming a geometry of uniform foreground screens (i.e. f ν = e−τ f 0ν ), we apply this curve in this form
only.
The two SN extinction curves are compared with the MW andCalzetti extinction curves. The SN extinction curves were synthesized in unmixed ejecta of SNe II by Hirashita et al. (2005), based on
the dust production models by Nozawa et al. (2003). Nozawa et al.

C 2010 The Authors, MNRAS 412, 1070–1080
C 2010 RAS
Monthly Notices of the Royal Astronomical Society 

1073(2003) have investigated the dust condensation in SNe II with various progenitor masses, taking into account the time evolution of gas
temperatures in the ejecta by solving the radiative transfer equation
with energy deposition of radioactive elements. The two extreme
cases are considered for the elemental composition in the He core:
the unmixed case and the uniformly mixed case. The formation of
all possible condensations is calculated simultaneously by applying a theory of non-steady state nucleation and grain growth under
the assumption of complete formation of CO and SiO molecules.
Various dust species condense according to the elemental composition at the formation site in the unmixed ejecta; carbon in the inner
C-rich He layer, Al2 O3 , Mg silicates (Mg2 SiO4 and MgSiO3 ) in
the O-rich layer, MgO in the outer O-rich layer, SiO2 in the innerO-rich layer, Si and FeS in the Si–S layer and FeS and Fe in the
inner most Fe–S layer. In the mixed ejecta, only oxide grains such as
Al2 O3 , Mg silicates, SiO2 and Fe3 O4 condense because the elemental composition assuming complete formation of CO is oxygen richin the entire He core. In the unmixed cases, C and Si dust species
dominate the dust mass in the low-mass progenitor (Nozawa et al.
2003), and Mg2 SiO4 is getting dominant with increasing progenitor
mass, while in the mixed ejecta SiO2 and Mg2 SiO4 dust dominate
the dust mass.
Combining the optical constants in the literature with the grain
species and size distribution of SN-condensed dust derived by
Nozawa et al. (2003), and using the Mie theory, Hirashita et al.
(2005) have calculated absorption and scattering cross-sections ofhomogeneous spherical grains with various sizes condensed in
mixed and unmixed ejecta. It should be pointed out here that in
the unmixed ejecta C and Si dust grains with radius a ≥ 0.1 μm
are significantly abundant in the low-mass progenitors than in thehigh-mass progenitors. This results in the different behaviour ofthe extinction curves: the extinction curve in UV region tends tobe flatter in the unmixed cases than in the mixed cases (Hirashita
et al. 2005). The details for the adopted quantities are described in
Hirashita et al. (2005). The extinction curves in unmixed ejecta of
SNe II reasonably reproduce the extinction law of the z = 6.2 BAL
QSO found by Maiolino et al. (2004). In this paper, we use the two
extinction curves of unmixed SNe II for the Salpeter (1955) initial
mass function (IMF); one has a mass range of progenitors from 8
to 30 M , and the other from 8 to 40 M . These were obtained byinterpolating or extrapolating extinction curves that Hirashita et al.
(2005) calculated for unmixed SNe II with progenitor masses of 13,
20, 25 and 30 M .4
4 SYNTHETIC SPECTRA OF STELLARP O P U L AT I O N S
The evolutionary synthesis codes by Bruzual & Charlot (2003) areused to generate synthetic spectra of evolutionary stellar population
models. We adopt the Salpeter (1955) IMF with the lower and uppermass cut-offs of 0.1 and 100 M , respectively. It is noted that theweight of massive stars is greater in the Chabrier (2003) IMF than in
the Salpeter IMF, resulting in a mass to luminosity ratio roughly 1.6
times heavier in the Salpeter IMF. The star formation histories are
exponentially declining models with the star formation rate (SFR) ∝t−1
sfr exp (− t/tsfr ), where tsfr is the e-folding time-scale and t is the

4 In Hirashita et al. (2005), the slope of extinction curve for 25 and 30 Mprogenitors was overestimated because of the error in the calculation processes. We have replaced the results with correct ones. Nevertheless, theconclusions in Hirashita et al. (2005) are not altered by this change.

1074

K. Kawara et al.

age, i.e. the time after the onset of the first star formation. Although
Kawara et al. (2010) used the instantaneous burst models, suchmodels are not appropriate to reproduce simultaneously the strong[O II]λ3727 emission line from H II regions ionized by young, hotstars and the observed SED. The synthetic spectra are generated ina parameter range of t = 3 × 105 to 8 × 109 yr in a step of 0.2 dex
and tsfr = 106 to 1.6 × 1010 yr in a step of 0.3 dex. For the shortest
tsfr , the star formation history is near an instantaneous burst, while
for the longest tsfr , the star formation rate is almost constant overthe cosmic time. The metallicity of the synthetic spectra has a rangeof Z = 0.005–2.5 Z (i.e. six points at 0.005, 0.02, 0.2, 0.4, 1.0 and
2.5 Z ).
The extinction models discussed in Section 3 are then applied
to the synthetic spectra generated in the above procedure, yielding
template spectra which are used for the subsequent χ 2 fitting. In
the initial run, the extinction measured at 0.3 μm has a range ofτ ext (0.3 μm) = 0–10 in a step of 0.1. In case that the initial range is
too narrow to have the minimum χ 2 for a certain model, the range
and the step are widen by a factor of 3 to search for the minimum
χ 2 , and so on.
5 χ 2 T E M P L AT E - F I T T I N G M E T H O D
We search for the best set of parameters for SST J1604+4304 byusing a χ 2 template-fitting method as defined below:
2
2
2
+ χ[O
+ χCa
+ χfir2 ,χ 2 = χSED
II]
II

(11)

where the degree of freedom is ν = 12. In what follows, we define
the terms in equation (11):22=(Fobs − bFtemplate )2 /σSED
,
(12)
χSED
filters

where b is the parameter determined by fitting the template flux
density F template to the observed F obs . F obs are the photometric datameasured in 13 broad-band filters from B to 5.8 μm. The Spitzer
IRAC 8.0 μm is not used, because dust emission may be significant in this band corresponding to 3.7 μm in the rest frame. It
is noted that 5 per cent of the flux was added in the quadratic22= σ ph+ (0.05F obs )2 .
form to the photometric error σ ph , i.e. σ SED
This is necessary to avoid too small and non-realistic photometricerror.
The strong [O II]λ3727 line requires the presence of OB stars.
We thus include χ 2[O II] in the χ 2 template-fitting method. χ 2[O II] isdefined as2II][O II]
2
2
χ[O= EW[O/σ[O,(13)obs − EWtemplate
II]
II]
II]
[O II]
where EW [O
obs = 31.1 ± 4.3 Å. To derive EW template , we first obtainthe Hα luminosity by counting ionizing photons in the synthesized
spectra. We then adopt the [O II] and Hα calibration performedby Kewley, Geller & Jansen (2004). For the metallicity scaling of[O II]/Hα, we use their theoretical curves for q = 3 × 107 cm s−1 withan additional condition that [O II]/Hα should not be less than 0.1.
This additional condition is applied to avoid a negative or unrealistic
value of [O II]/Hα outside of the calibrated range of metallicity. The
calibrated [O II]/Hα has a peak value of 1.6 at Z = 0.4 Z and drops
to a value of 0.1 at Z = 0.14 Z or Z = 2.44 Z . Because Kewley
et al. (2004) estimated the uncertainty to be ∼ 30 per cent, we add
22= σ [O
this uncertainty to the observational error as σ SEDII]SED +[O II]
2(0.3 EWtemplate ) .
As discussed in Kawara et al. (2010), the D(Ca II) index, whichis the ratio of the Ca II to H I absorption lines in equivalent width, isuseful to measure the age of the relatively old stellar populations.
χ 2Ca II is defined as
22= (D(Ca II)obs − D(Ca II)template )2 /σD(Ca
,
χCa
II
II)

(14)

where σ D(Ca II) is the observational error.
The luminosity absorbed by dust should mostly be re-emitted in
the far-IR. If the energy source of the galaxy is dominated by starformation, the absorbed luminosity should match with the far-IR
luminosity. We include this effect as follows:χfir2 = (Fir − Fabsorbed )2 /σir2 ,(15)

where F ir , F absorbed and σ ir are the observed flux in 8–1000 μm, theabsorbed flux in the UV to near-IR and the observed error of F ir ,respectively. It is noted that F absorbed is calculated for a given b which
is determined by ∂χ 2SED /∂b = 0.
6 R E S U LT S
First of all, the Calzetti extinction models, which have a uniform
foreground screen, yield the best-fitting model having the reduced
χ square χ 2ν = 5.1 with χ 2[O II] = 0.0, χ 2Ca II = 0.3 and χ 2fir = 0.9.
This model has a parameter set of the optical depth τ ext (0.3 μm) =5.3, the age t = 2.0 × 108 yr with the star formation time-scale oftsfr = 6.3 × 107 yr and the metallicity in stellar atmosphere Z =1 Z , EW([O II]) = 29 Å, D(Ca II) = 0.09 and F absorbed = 6.6 ×
10−13 erg cm−2 s−1 . This is our reference to which the goodness of
our other models will be judged.
The left-hand panels of Fig. 4 plot χ 2ν as a function of τ in /(τ in +
τ sc ) for the uniform foreground screen mixed with internal dust. τ inand τ sc are τ ext for internal dust and screen dust, respectively. For
the pure internal dust geometry, τ in /(τ in + τ sc ) = 1, while the pureforeground screen geometry, τ in /(τ in + τ sc ) = 0. The scattering
properties are assumed to be the no scattering, the isotropic scattering (g = 0) or the forward-only scattering (g = 1). In case of no
scattering, forward scattering is assumed for internal dust.
The pseudo-MW extinction curve, which is made by removing
the 2175-Å feature from the MW extinction curve, has a smallerχ 2ν than the MW extinction curve, suggesting the 2175-Å feature isvery weak or unrecognizable in SST J1604+4304. Because of thelarge χ 2ν values5 relative to the Calzetti and SN extinction curves,
the MW extinction curve is unlikely regardless of the presence ofthe 2175-Å feature. The pure internal dust geometry can be ruledout regardless of the choice of extinction curves. Another difficultyof the pure internal geometry is that too much dust is required toreproduce the observed SED (i.e. τ ext = 70–200). The estimateddust mass is far greater than that observed by a factor of 10 or evenmore. The general trend is that χ 2ν becomes smaller as τ in /(τ in + τ sc )approaching to zero, indicating that the extinction due to internal
dust is a small fraction of the total extinction as already found in
highly reddened H II regions (i.e. Natta & Panagia 1984). The two
SN extinction curves have χ 2ν less than that of the Calzetti extinction
curve (i.e. χ 2ν = 5.1) for a typical range of τ in /(τ in + τ sc ) = 0.0–0.7.
The middle and right-hand panels of Fig. 4 plot χ 2ν as a functionof the number of clumps along the line of sight for the clumpy
foreground screen; pure clumpy foreground screens (i.e. τ in /(τ in +
τ sc ) = 0) are shown in the middle panels and foreground screenswith internal dust having τ in /(τ in + τ sc ) = 0.4 in the right-hand
5 The approximate values of the increment of χ 2 are 1.14, 1.75 and 2.36 forν68 per cent(1σ ), 95 per cent (2σ ) and 99.6 per cent (3σ ) confidence where
the degree of freedom is ν = 12 (Avni 1976).
C 2010 The Authors, MNRAS 412, 1070–1080
C 2010 RAS
Monthly Notices of the Royal Astronomical Society Supernova dust at z ∼ 1

1075

Figure 4. The left-hand panels show χ 2ν (ν = 12) as a function of τ in /(τ in + τ sc ) for no scattering (top), forward scattering (middle) and isotropic scattering(bottom), where τ in is the optical depth for internal dust and τ sc for screen dust. The uniform screen is assumed. In case of no scattering, forward scattering isassumed for internal dust. Solid lines represent SN dust with 8–30 M , dash lines SN dust with 8–40 M , dash–dot lines MW dust and dash–dot–dot lines
pseudo-MW dust without the 2175-Å feature. The middle panels are same as the left-hand panels, but as a function of the number of clumps in the line of sightcontained in the screen for no internal dust (i.e. τ in /(τ in + τ sc ) = 0.0). The right-hand panels are same as the left-hand panels, but as a function of the numberof clumps in the line of sight contained in the screen for τ in /(τ in + τ sc ) = 0.4. χ 2ν = 5.1 results if the Calzetti extinction curve is applied.

panels. The MW extinction curves with and without the 2175-Åfeature have larger χ 2ν than that of the Calzetti extinction curve
except for clumpy screens with the average number of clumps of
N = 4 (see the upper right-hand panel). The overall trends are that
small χ 2ν values (i.e. χ 2ν < 5.1) are found in larger N. Uniform dust
screens have been explicitly or implicitly assumed in many cases in
the literature. These figures indicate that uniform dust screens aregood approximations to the real distribution of dust. The two SN
extinction curves have χ 2ν smaller than that of the Calzetti extinction
curve for large N.
Fig. 5 gives typical features of parameters as a function ofτ in /(τ in + τ sc ), taking isotropic scattering models as an example.
C 2010 The Authors, MNRAS 412, 1070–1080
C 2010 RAS
Monthly Notices of the Royal Astronomical Society Those parameters are the total optical depth τ total (0.3 μm) = τ in +
τ sc at 0.3 μm, the metallicity Z, the time after the onset of thefirst star formation t and the e-folding time-scale tsfr . Because theeffect of scattering reduces the effective optical depth τ eff as seenin equations (7)–(9), scattering dust screens require a larger opticaldepth than no scattering screens; taking models without internaldust, τ ext (0.3 μm) ∼ 5 for the no scattering screen, τ ext (0.3 μm) ∼ 8
for the isotropic scattering and τ ext (0.3 μm) ∼ 15 for the forwardonly scattering. The optical depth of the SN 8–40 M extinction
curve requires larger optical depths than those of the other curves.
This is because albedo is greater in the SN 8–40 M curve than theothers (see Table 1). Z is smaller on the SN curves than the MW

1076

K. Kawara et al.

Figure 5. Sample illustration of the total optical depth τ in + τ sc at 0.3 μm(upper left), the metallicity Z in units of the solar value (bottom left), thetime after the onset of the first star formation t (yr) (upper right) andthe e-folding time-scale of the star formation tsfr as a function of the ratio of the optical depth of internal dust to internal+screen dust τ in /(τ in +
τ sc ). The isotropic scattering (g = 0) is assumed. Solid lines represent SN
dust with 8–30 M , dash lines SN dust with 8–40 M , dash–dot lines MW
dust and dash–dot–dot lines pseudo-MW dust without the 2175-Å feature.

curve. On the other hand, t is older on the SN curves than the MW
curves. The older t on the SN curves can be attributed to the factthat the slope of the SN curves are flatter than that of the MW curve.
To reproduce the redness of the SED, the steep MW curve requires
bluer and younger stellar populations.
Fig. 6 shows an example of model SEDs along with the observed
SED. It is noted again that the longest wavelength data (i.e. IRAC
8.0 μm) are not used for χ 2ν fitting. The sample model is obtained byassuming the SN 8–40 M extinction curve with the isotropic scattering without internal dust. The fitted parameters are τ sc (0.3 μm) =10.5, t = 3.2 × 108 yr, tsfr = 1.3 × 108 yr and Z = 1 Z . This modelhas χ 2ν = 3.7 with χ 2[O II] = 0.4, χ 2Ca II = 0.6 and χ 2fir = 2.0, providing
a smaller χ 2ν than those with the Calzetti extinction curve. Comparing with χ 2ν = 5.1, χ 2[O II] = 0.0, χ 2Ca II = 0.3 and χ 2fir = 0.9 for the
best-fitting model with the Calzetti extinction curve, the reductionin χ 2ν in this SN model is attributed to the reduction in χ 2SED , indicating that the reddening by the SN extinction curve reproduces theSED better than the Calzetti extinction curve.
7 S N II E X T I N C T I O N C U RV E
As shown so far, the SN extinction curves synthesized by Hirashita
et al. (2005) reasonably reproduce the observed characteristics of
ULIRG SST J1604+4303. The goodness of fit is comparable toor better than those with the Calzetti extinction curve. The SN
extinction curves must be useful to study young dusty galaxies not
only at high redshift but also at modest redshift (z ∼ 1). We thus
present the details of the SN extinction curves from the UV to
submillimetre.
Fig. 7 shows κ abs (λ) of the SN dust along with that of the MW
dust, where κ abs (λ) is the mass absorption cross-section in unitsof cm2 g−1 . Table 1 gives κ ext (λ) the mass extinction cross-section,
κ abs (λ) the mass absorption cross-section and albedo for the SN
dust; SN 8–30 M dust in two–four columns and SN 8–40 M
dust in five–seven columns. It is noted that the optical constant ofcrystalline SiO2 is used in the models by Hirashita et al. (2005). As
a result, the SN curves have spiky features between 9 and 21 μm.

The contribution of crystalline SiO2 to the SN extinction curve is afew per cent at most in the mid- to far-IR range except the spikes at
8.7, 12.25 and 20.5 μm. We smooth the calculated mass extinction
coefficient of Si grains around 1 μm to remove small and somehow
artificial ripple features arising from large Si grains.
As shown in Fig. 3, the SN extinction curves are flatter than theothers: RV (=AV /EB−V ) is 7.8 for the 8–30 M curve and 5.8 for the
8–40 M , while RV = 4.05 for the Calzetti curve and RV = 3.1 forthe MW curve. In other words, for a given amount of reddening, theextinction is greater in the SN curves than the others. As suggested inFig. 7, the silicate feature around 9.7 μm is much less in the SN dust
than in the MW dust. This feature could be used to distinguish the SN
extinction curve from the others by means of mid-IR observationsin the future.
Another striking feature of the SN curve is the behaviour ofκ abs (λ) from 40 to 1000 μm. Assuming a simple analytic form ofκ abs (λ) ∝ ν β , we have β 1.6 for the SN dust, while β 2 for the
MW dust. Single-temperature dust spectral distributions of the form
ν β Bν (T dust ) are widely used to fit observed SED in the far-IR and
submillimetre, where Bν (T dust ) is the Planck function. However,
fitting only a few data points results in a significant correlation
between β and T dust that can account for the data (Blain, Barnard &Chapman 2003): higher T dust with smaller β; T dust derived with the
SN dust is 13 per cent higher than that with the MW dust. Assuming
the SN dust, we have T dust = 38.3 ± 1.4 K for SST J1604+4304
which is a common characteristics of 104 IRAS bright galaxies
observed at 850 μm by Dunne et al. (2000), where the best fit was
found for T dust = 35.6 ± 4.9 K and β = 1.3 ± 0.2.
8 D U S T I N S S T J 1 6 0 4 +4 3 0 4
Using T dust = 38.3 ± 1.4 K combined with the observed 160 μm
flux density and κ abs (80 μm) = 22 cm2 g−1 in Table 1 implies thatthe dust mass M dust in SST J1604+4303 is 2.2 ± 1.1 × 108 M .
Consider a spherical shell with a radius of r and a thickness of
δr, and suppose that internal dust fills inside the sphere with r andscreen dust in the shell. The mass of the internal dust can be simply
approximated to M in = 4π (τ in /κ ext )r2 /3 and that of the screen
dust to M sc = 4π (τ sc /κ ext )r2 , thus M in /M sc = (τ in /τ sc )/3. The bestmodel for the isotropic scattering uniform screen without internal
dust has τ ext (0.3 μm) = 10.5. In this case, the inferred radius of thespherical shell is about 3.3 kpc or 0.4 arcsec which is comparable
to the optical size of SST J1604+4303.
If the star formation history of SST J1604+4303 is describedby t = 3.2 × 108 yr, tsfr = 1.3 × 108 yr and Z = 1 Z , then theluminosity to mass ratio is 7.6, implying a stellar mass of 2.4 ×
1011 M . During the last 3.2 × 108 yr, ∼2.0 × 109 SNe II explodedin this galaxy. If the dust was mainly supplied by SN explosions,then we have a dust production rate of ∼0.1 M per SN by simply
dividing the dust mass by the number of SNe II.
Although χ 2ν is reduced by adopting the synthesized SN extinction curves, the goodness of fit is far from the perfect fitting (i.e.
χ 2ν ∼ 1). This may come from various reasons. For example, themodels for dust formation in SNe II are still coarse and the contribution from the dust formed in AGB stars needs to be included,or the stellar population models are too simple and the real stellar
population may be more complex and consist of stars with different
ages with different metallicities.
Nozawa et al. (2007) take into account the effect of reverse
shock destruction and Hirashita et al. (2008) calculate the extinction curve based on their result. The extinction curves after reverse
shock destruction are generally flatter than those of Hirashita et al.

C 2010 The Authors, MNRAS 412, 1070–1080
C 2010 RAS
Monthly Notices of the Royal Astronomical Society Supernova dust at z ∼ 1

1077

Table 1. Extinction cross-sections of the SN II dust models. A machine-readable version of this table is available as Supporting Information.

λ (μm)
6.67e−02
9.12e−02
1.00e−01
1.20e−01
1.40e−01
1.60e−01
1.80e−01
1.90e−012.00e−01
2.10e−01
2.20e−012.30e−012.50e−01
3.00e−01
3.55e−01
3.64e−01
4.41e−01
4.69e−01
5.47e−01
6.17e−01
6.42e−01
6.49e−01
7.00e−01
7.48e−01
8.02e−01
8.66e−01
8.93e−01
9.00e−01
1.00e+00
1.22e+00
1.63e+00
2.19e+00
2.75e+00
3.45e+00
3.80e+00
4.75e+00
6.00e+00
7.00e+00
8.00e+00
8.28e+00
8.71e+00
9.17e+00
9.70e+00
1.13e+01
1.20e+01
1.38e+01
1.53e+01
1.61e+01
1.70e+01
1.79e+01
1.98e+01
2.08e+01
2.19e+01
2.50e+01
4.00e+01
6.00e+01
1.00e+02
1.40e+02
2.00e+02
3.50e+02
8.50e+02
1.00e+03

κ ext(cm2 g−1 )4.69e+04
4.73e+04
4.71e+04
5.07e+04
5.09e+04
4.68e+04
4.04e+04
3.87e+04
3.79e+04
3.77e+04
3.77e+04
3.77e+04
3.75e+04
3.40e+04
3.17e+04
3.14e+04
2.80e+04
2.69e+04
2.48e+04
2.21e+04
2.10e+04
2.08e+04
2.03e+04
1.95e+04
1.86e+04
1.75e+04
1.70e+04
1.68e+04
1.51e+04
1.26e+04
9.48e+03
5.83e+03
4.14e+03
2.73e+03
2.30e+03
1.66e+03
9.46e+02
8.27e+02
6.33e+02
6.21e+02
1.56e+03
6.34e+02
6.63e+02
4.69e+02
3.93e+02
2.56e+02
3.22e+02
6.18e+02
5.95e+02
2.74e+02
2.50e+02
2.90e+02
2.03e+02
1.47e+02
5.69e+01
3.30e+01
1.59e+01
1.09e+01
5.08e+00
2.03e+00
4.38e−01
3.17e−01

SN 8–30 Mκ abs(cm2 g−1 )2.34e+04
2.34e+04
2.30e+04
2.57e+04
2.43e+04
2.03e+04
1.67e+04
1.57e+04
1.50e+04
1.46e+04
1.45e+04
1.46e+04
1.44e+04
1.39e+04
1.32e+04
1.31e+04
1.20e+04
1.16e+04
1.08e+04
9.88e+03
9.66e+03
9.60e+03
9.21e+03
8.90e+03
8.55e+03
8.16e+03
7.99e+03
7.95e+03
7.36e+03
6.08e+03
4.48e+03
3.07e+03
2.42e+03
1.87e+03
1.59e+03
1.13e+03
6.35e+02
5.46e+02
4.72e+02
4.73e+02
1.25e+03
4.99e+02
5.69e+02
3.61e+02
3.25e+02
2.29e+02
2.97e+02
6.01e+02
5.84e+02
2.67e+02
2.45e+02
2.77e+02
1.99e+02
1.45e+02
5.66e+01
3.30e+01
1.59e+01
1.09e+01
5.08e+00
2.03e+00
4.38e−01
3.17e−01

C 2010 The Authors, MNRAS 412, 1070–1080
C 2010 RAS
Monthly Notices of the Royal Astronomical Society Albedo5.01e−01
5.05e−015.12e−01
4.93e−01
5.23e−01
5.67e−01
5.87e−01
5.94e−01
6.04e−01
6.12e−01
6.15e−01
6.14e−01
6.16e−01
5.92e−01
5.82e−01
5.81e−01
5.71e−01
5.69e−01
5.65e−01
5.53e−01
5.40e−01
5.39e−01
5.46e−01
5.45e−01
5.40e−01
5.34e−01
5.29e−01
5.27e−015.12e−01
5.19e−01
5.27e−01
4.74e−01
4.14e−01
3.15e−01
3.06e−01
3.21e−01
3.29e−01
3.40e−012.53e−01
2.39e−01
1.97e−012.13e−01
1.42e−012.29e−01
1.74e−01
1.02e−01
7.65e−02
2.80e−02
1.74e−02
2.68e−022.00e−02
4.49e−02
1.97e−02
1.41e−02
4.89e−03
1.61e−03
4.08e−04
1.96e−04
8.02e−05
2.25e−05
3.09e−06
2.00e−06

κ ext(cm2 g−1 )4.86e+04
4.95e+04
4.95e+04
5.23e+04
5.25e+04
4.97e+04
4.54e+04
4.29e+04
4.14e+04
4.04e+04
3.96e+04
3.88e+04
3.72e+04
3.10e+04
2.73e+04
2.68e+04
2.26e+04
2.14e+04
1.93e+04
1.70e+04
1.61e+04
1.60e+04
1.55e+04
1.49e+04
1.41e+04
1.32e+04
1.27e+04
1.26e+04
1.13e+04
9.64e+03
7.52e+03
4.82e+03
3.49e+03
2.40e+03
2.02e+03
1.53e+03
7.70e+02
7.37e+02
4.80e+02
4.70e+02
1.61e+03
5.56e+02
6.42e+02
3.86e+02
3.20e+02
1.94e+02
3.49e+02
8.32e+02
8.06e+02
3.01e+02
2.70e+02
3.08e+02
2.12e+02
1.44e+02
4.55e+01
2.68e+01
1.28e+01
8.45e+00
4.08e+00
1.61e+00
3.57e−012.65e−01

SN 8–40 Mκ abs (cm2 g−1 )Albedo2.37e+04
2.44e+04
2.44e+04
2.64e+04
2.47e+04
2.12e+04
1.61e+04
1.40e+04
1.23e+04
1.13e+04
1.08e+04
1.07e+04
1.04e+04
9.80e+03
9.00e+03
8.89e+03
8.15e+03
7.84e+03
7.30e+03
6.64e+03
6.51e+03
6.47e+03
6.19e+03
5.98e+03
5.73e+03
5.46e+03
5.36e+03
5.33e+03
4.95e+03
4.12e+03
3.03e+03
2.05e+03
1.62e+03
1.17e+03
9.94e+02
6.87e+02
3.72e+02
3.25e+02
3.10e+02
3.33e+02
1.38e+03
4.62e+02
5.71e+02
3.22e+02
2.79e+02
1.79e+02
3.38e+02
8.23e+02
7.99e+02
2.96e+02
2.67e+02
3.04e+02
2.10e+02
1.43e+02
4.53e+01
2.68e+01
1.28e+01
8.45e+00
4.08e+00
1.61e+00
3.57e−012.65e−01

5.11e−01
5.07e−01
5.07e−01
4.96e−01
5.29e−01
5.74e−01
6.46e−01
6.74e−01
7.02e−01
7.19e−01
7.26e−01
7.24e−01
7.21e−01
6.84e−01
6.70e−01
6.68e−01
6.40e−01
6.34e−01
6.22e−01
6.10e−01
5.96e−01
5.95e−01
6.00e−01
5.99e−01
5.94e−01
5.86e−01
5.79e−01
5.77e−01
5.62e−01
5.72e−01
5.96e−01
5.75e−01
5.35e−01
5.11e−01
5.08e−01
5.51e−01
5.17e−01
5.58e−01
3.53e−01
2.94e−01
1.46e−01
1.69e−01
1.10e−01
1.64e−01
1.27e−01
8.14e−02
3.37e−02
1.14e−02
8.84e−03
1.61e−02
1.17e−02
1.54e−02
1.07e−02
8.60e−03
3.96e−03
1.34e−03
3.41e−04
1.47e−04
6.63e−05
1.87e−05
2.42e−06
1.61e−06

1078

K. Kawara et al.
ample, shattering in warm ISM is occurring (Hirashita et al. 2010),
in which case the extinction curve after shattering is similar to thatused in this paper. Moreover, the extinction curve by Hirashita et al.
(2005) is supported by observations in the sense that it is consistent
with the extinction curve derived by Maiolino et al. (2004) for ahigh-z BAL QSO at z = 6.2. Note that the Maiolino et al. (2004)’s
extinction curve is also consistent with a high-z (z = 6.3) gamma-ray
burst (Stratta et al. 2007).
We have shown that a flatter extinction curve (or a larger RV )than the MW curve fits the SED of SST J1604+4304 better. A flat
extinction curve may also be realized if dust growth occurs in dense
clouds (Draine 2003). Thus, it is worth investigating whether or
not a flat extinction curve predicted from the grain growth fits theobserved SED of SST J1604+4304. This is left for future work.
9 CONCLUSIONS

Figure 6. Sample illustration of model SEDs. The illustrated model isobtained by assuming the SN 8–40 M extinction curve with the isotropic
scattering without internal dust. This model has χ 2ν = 3.7, τ sc = 10.2 at
0.3 μm, t = 5.1 × 108 yr, tsfr = 2.5 × 108 yr and Z = 1 Z . The top panel
shows the model spectrum along with the observed SED. The upper curve isthe unobscured model spectrum. The bottom panel shows χ 2ν as a functionof t and τ sc at 0.3 μm.

Figure 7. Plot of κ abs (λ) the absorption cross-section per dust mass in unitsof cm2 g−1 . Note that spiky features between 9 and 21 μm on the SN curve
are mostly attributed to crystal SiO2 and MgO; spiky features of crystal SiO2
occur at 8.7, 12.25 and 20.5 μm, while a spike MgO feature at 16.5 μm.

(2005), since small grains are selectively destroyed. In this paper, we
simply adopt the extinction curves by Hirashita et al. (2005). In
other words, the effect of reverse shock destruction is not taken
into account. Indeed, the flat extinction curves after reverse shock
destruction are not consistent with the large reddening required for
SST J1604+4304. This implies that the reverse shock destruction
may not be efficient or that the small grain production by, for ex-

We have applied the SN extinction curves to reproduce the observed
properties of SST J1604+4304 which is a young IR galaxy at z ∼1. The SN extinction curves used in this work was obtained frommodels of unmixed ejecta of SNe II for the Salpeter IMF with a mass
range from 8 to 30 M or 8 to 40 M . These curves are comparedwith the MW and Calzetti extinction curves. The data used in theχ 2 fitting methods are the fluxes in 13 photometric bands from B to
IRAC 5.8 μm, the equivalent width of the [O II]λ3727, the D(Ca II)
index and the IR flux in 8–1000 μm.
The effect of dust distributions on the attenuation of starlight wasinvestigated by performing the χ 2 fitting method with various dust
distributions. These are the commonly used uniform dust screen,
the clumpy dust screen and the internal dust geometry. In addition
to these geometries, we have investigated three scattering properties, namely, no scattering, isotropic scattering and forward-only
scattering. Judging from the χ 2 values, the uniform screen models
with any scattering property provide good approximations to thereal dust geometry. Internal dust is inefficient to attenuate starlightand thus cannot be the dominant source of the extinction.
We have found that the SN extinction curves reproduce the data of
SST J1604+4304 comparable or better than the Calzetti extinction
curve. The MW extinction curve is too steep and is not satisfactory
agreement with the data unless several dusty clumps are in the lineof sight. This trend may be explained by the abundance of SN-origin
dust in these galaxies; SN dust is the most abundant in the young
IR galaxy at z ∼ 1, abundant in local starbursts and less abundantin the Galaxy. If dust in SST J1604+4304 is dominated by SN dust,
the dust production rate is ∼0.1 M per SN.
We provide the extinction cross-sections of the SN dust alongwith albedo from UV to submillimetre wavelengths in Table 1(see Supporting Information for a machine-readable version of thistable).
AC K N OW L E D G M E N T S
We thank the anonymous referee for very helpful and constructive
comments. This work has been supported in part by Grants-in-Aid
for Specially Promoted Research on Innovative Areas (22111503),
Specially Promoted Research (20001003) and the Promotion of Science (18104003, 20340038, 22684004) from JSPS. HH is supported
by NSC grant 99-2112-M-001-006-MY3.
REFERENCES
Avni Y., 1976, ApJ, 210, 642
Bianchi S., Schneider R., 2007, MNRAS, 378, 973
C 2010 The Authors, MNRAS 412, 1070–1080
C 2010 RAS
Monthly Notices of the Royal Astronomical Society Supernova dust at z ∼ 1Blain A. W., Barnard V. E., Chapman S. C., 2003, MNRAS, 338, 733
Borys C. et al., 2006, ApJ, 636, 134
Bruzual G., Charlot S., 2003, MNRAS, 344, 1000
Calzetti D., 2001, PASP, 113, 1449
Calzetti D., Kinney A. L., Storchi-Bergmann T., 1994, ApJ, 429, 582
Carico D. P., Sanders D. B., Soifer B. T., Matthews K., Neugebauer G.,
1990, AJ, 100, 70
Chabrier G., 2003, PASP, 115, 763
Code A. D., 1973, in Greenberg J. M., Van de Hulst A. C., eds, Proc.
IAU Symp. 52, Interstellar Dust and Related Topics. Reidel, Dordrecht,
p. 505
Desai V. et al., 2006, ApJ, 641, 133
Draine B. T., 2003, ARA&A, 41, 241
Dunne L., Eales S., Edmunds M., Ivison R., Alexander P., Clements D. L.,
2000, MNRAS, 315, 115
Dunne L., Eales S., Ivison R., Morgan H., Edmunds M., 2003, Nat, 424,
285
Dunne L. et al., 2009, MNRAS, 394, 1307
Dwek E., Galliano F., Jones A. P., 2007, ApJ, 662, 927
Gallerani S. et al., 2010, A&A, 523, 85
Giavalisco M., 2002, ARA&A, 40, 579
Gomez H. L. et al., 2009, MNRAS, 397, 1621
Heger A., Fryer C. L., Woosley S. E., Langer N., Hartmann D. H., 2003,
ApJ, 591, 288
Hirashita H., Nozawa T., Kozasa T., Ishii T. T., Takeuchi T. T., 2005,
MNRAS, 357, 1077
Hirashita H., Nozawa T., Takeuchi T. T., Kozasa T., 2008, MNRAS, 384,
1725
Hirashita H., Nozawa T., Yan H., Kozasa T., 2010, MNRAS, 404,
1437
Jones A. P., Tielens A. G. G. M., Hollenbach D. J., McKee C. F., 1994, ApJ,
433, 797
Jones A. P., Tielens A. G. G. M., Hollenbach D. J., 1996, ApJ, 469, 740
Kawara K., Nishida M., Phillips M. M., 1989, ApJ, 337, 230
Kawara K. et al., 2010, MNRAS, 402, 335
Kewley L. J., Geller M. J., Jansen R. A., 2004, AJ, 127, 2002
Klaas U., Haas M., Heinrichsen I., Schulz B., 1997, A&A, 325, L21
Kotak R. et al., 2009, ApJ, 704, 306
Kozasa T., Hasegawa H., Nomoto K., 1989, ApJ, 344, 325
Kozasa T., Hasegawa H., Nomoto K., 1991, A&A, 249, 474
Kozasa T., Nozawa T., Tominaga N., Umeda H., Maeda K., Nomoto K.,
2009, in Henning Th., Grün E., Steinacker J., eds, ASP Conf. Ser.
Vol. 414, Cosmic Dust – Near and Far. Astron. Soc. Pac., San Francisco,p. 43
Maiolino R., Schneider R., Oliva E., Bianchi S., Ferrara A., Mannucci F.,
Pedani M., Roca Sogorb M., 2004, Nat, 431, 533
Maiolino R. et al., 2006, Mem. Soc. Astron. Ital., 77, 643
Morgan H. L., Edmunds M. G., 2003, MNRAS, 343, 427
Natta A., Panagia N., 1984, ApJ, 287, 228
Noll S. et al., 2009, A&A, 499, 69
Nozawa T., Kozasa T., Umeda H., Maeda K., Nomoto K., 2003, ApJ, 598,
785
Nozawa T., Kozasa T., Habe A., 2006, ApJ, 648, 435
Nozawa T., Kozasa T., Habe A., Dwek E., Umeda H., Tominaga N., Maeda
K., Nomoto K., 2007, ApJ, 666, 955
Nozawa T. et al., 2008, ApJ, 684, 1343
Nozawa T., Kozasa T., Tominaga T., Maeda K., Umeda H., Nomoto K.,
Krause O., 2010, ApJ, 713, 356
Puxley P. J., Brand P. W. J. L., 1994, MNRAS, 266, 431
Rho J. et al., 2008, ApJ, 673, 271
Rigopoulou D., Lawrence A., Rowan-Robinson M., 1996, MNRAS, 278,
1049
Salpeter E. E., 1955, ApJ, 121, 161
Schneider R., Ferrara A., Salvaterra R., 2004, MNRAS, 351, 1379
Stratta G., Maiolino R., Fiore F., D’Elia V., 2007, ApJ, 661, L9
Temim T., Slane P., Reynolds S. P., Raymond J. C., Borkowski K. J., 2010,
ApJ, 710, 309
Todini P., Ferrara A., 2001, MNRAS, 325, 726
C 2010 The Authors, MNRAS 412, 1070–1080
C 2010 RAS
Monthly Notices of the Royal Astronomical Society 

1079

Valiante R., Schneider R., Bianchi S., Andersen A. C., 2009, MNRAS, 397,
1661
Zafar T., Watson D. J., Malesani D., Vreeswijk P. M., Fynbo J. P. U., HjorthJ., Levan A. J., Michałowski M. J., 2010, A&A, 515, 94

APPENDIX A: TWO STREAM
A P P ROX I M AT I O N F O R R A D I AT I O N
TRANSFER IN A DUSTY SLAB
The two stream approximation in radiation transfer was applied toan spherically symmetric dust envelope surrounding a star by Code(1973). Here we consider a plane-parallel dust slab.
We measure distance z outward for a plane-parallel slab, thenσν
μ dIν (z, μ)= jν − κν Iν (z, μ) +p(μ, μ )Iν (z, μ) dω ,
ρ
dz
4π 4π
(A1)
where I ν (z, μ) is the specific intensity, μ = cos θ the directional
cosine, ν the frequency, κ ν the mass extinction coefficient, σ ν themass scattering coefficient, ρ the mass density of dust and jν theratio of the stellar emission per unit dust mass. The phase function
p(μ, μ ) is defined by
1 1p(μ, μ ) dω = 4π orp(μ, μ ) dμ = 1.
(A2)
2 −1
4π
By introducing the optical depth τ ν which is defined by dτ ν =−ρκ ν dz, the radiation transfer equation is reduced toγν 1dIν (τ, ν)= −jν /κν + Iν −p(μ, μ )Iν (τ, μ ) dμ , (A3)
μ
dτν
2 −1
where γ ν = σ ν /κ ν is the albedo. In what follows, we omit thesubscript ν for simplicity.
A1 Slab illuminated at the backside
Consider a slab which contains no emitters and illuminated at thebackside, i.e. jν = 0. In the two stream approximation, the specificintensity I(τ , μ) is represented by the specific intensities at two
specified directions at ±μ; I + (τ ) = I(τ , μ) and I − (τ ) = I(τ , −μ)for 0 ≤ μ ≤ 1. I + (τ ) and I − (τ ) are assumed to be independent of μ
over the upward and the downward hemispheres, respectively.
According to Code (1973), we introduce the forward and backward scattering coefficients σ + and σ − which are defined by


1 1
1 1
σ+
σ−


=
=
p(μ, μ ) dμ andp(μ, −μ ) dμ ,
σ
2 0
σ
2 0
(A4)
respectively. Note that σ + + σ − = σ . A modified asymmetry factor
g is defined byσ+ − σ−,(A5)
σnote that g = 1 for forward scattering, g = −1 for backward scattering and g = 0 for isotropic scattering. Then using the relations
σ + /σ = (1 + g)/2 and σ − /σ = (1 − g)/2, the equations (A3) is
reduced to(1 + g) +(1 − g) −dI + (τ )(A6)= I + (τ ) − γ
I (τ ) − γ
I (τ ),
μ
dτ
2
2
g=

(1 + g) −(1 − g) +dI − (τ )= I − (τ ) − γ
I (τ ) − γ
I (τ ). (A7)dτ22Here we replace the optical depth along the z-axis with the optical depth along the line of sight dτ ext = dτ /μ and introduce new
−μ

1080

K. Kawara et al.

variables x = I + − I − and y = I + + I − , the above equations become
simplydx= (1 − γ ) y
dτext

and

dy
= (1 − γ g) x.
dτext

(A8)

Then the solutions of I + (τ ext ) and I − (τ ext ) are given by
I + (τext ) =1[(β + 1)A exp (ατext ) − (β − 1)B exp (−ατext )],2(A9)

I − (τext ) =1[(β − 1)A exp (ατext ) − (β + 1)B exp (−ατext )],2(A10)

where α 2 = (1 − γ )(1 − γ g) and β = [(1 − γ g)/(1 − γ )]1/2 .
Using the boundary conditions of I − (0) = 0 and I + (τ 1ext ) = I 0 for
determining the constants A and B, where τ 1ext is the optical depthof the slab along the line of sight, then the emergent intensity I + (0)is given by
I + (0) =2βI0.
(β + 1) exp (ατ1ext ) + (β − 1) exp (−ατ1ext )(A11)

A2 Slab containing emitters
We now consider a slab that contains emitters, i.e. stars. The slab isilluminated by these internal stars only. Then the radiation transfer
equations for I ± are simply given byμdI + (τ )(1 + g) += I + (τ ) − γ
I (τ )dτ2(1 − g) −−γ
I (τ ) − ,
2

(A12)

(1 + g) −dI − (τ )= I − (τ ) − γ
I (τ )dτ2(1 − g) +(A13)I (τ ) − ,−γ2where  = jν /κν . Then x = I + − I − and y = I + + I − can be
determined by solving the following simultaneous equations:dy= (1 − γ g) x.
dτext

(A14)

The intensities I ± (τ ext ) are given by
I + (τext ) =1(β + 1) [A − E0 (τext )] exp (ατext )2
1− (β − 1) [B − E1 (τext )] exp (−ατext ) ,2

1
I (τext ) = (β − 1) [A − E0 (τext )] exp (ατext )2
1− (β + 1) [B − E1 (τext )] exp (−ατext ) ,2where τext E0 (τext ) =dτext
(τext) exp −ατext
and0 τext   dτext .
(τext) exp ατext
E1 (τext ) =I + (0) =(β + 1)E0 (τ1ext ) exp (ατ1ext ) − (β − 1)E1 (τ1ext ) exp (−ατ1ext ).
(β + 1)2 exp (ατ1ext ) − (β − 1)2 exp (−ατ1ext )(A18)

If  = const against τ ext , I + (0) → τ 1eff for τ 1ext → 0. Thus,
the apparent extinction in as defined in equation (1) isin =2β (β + 1)[exp (ατ1ext ) − 1] − (β − 1)[1 − exp (−ατ1ext )].
ατ1ext(β + 1)2 exp (ατ1ext ) − (β − 1)2 exp (−ατ1ext )(A19)
Fig. A1 compares equation (A19) with the approximations by−τeff)/τeff for τ eff = (1 − γ )τ ext , τ eff = (1 − γ )0.5 τ ext andin = (1 − eτ eff = (1 − γ )0.75 τ ext . The geometry is internal dust with γ = 0.65
and g = 0. As shown in Fig. A1, τ eff = (1 − γ )0.75 τ ext approximatesto equation (A19) better than τ eff = (1 − γ )0.5 τ ext for the isotropic

scattering.
(A15)

S U P P O RT I N G I N F O R M AT I O N
Additional Supporting Information may be found in the online version of this article:−

0

Given the boundary conditions of I − (0) = I + (τ 1ext ) = 0, theemergent intensity I + (0) is given by

2β

−μ

dx
= (1 − γ ) y − 2 anddτext

Figure A1. Plot of the apparent extinction in as a function of the opticaldepth along the line of sight τ ext . The geometry is the internal dust with γ =0.65 and g = 0 (spherical scattering). The solid line represents the solutionof equation (A19) for  = const against τ ext . The other lines represent theabsorption approximated by in = (1 − e−τeff )/τeff ; the dotted line for τ eff =(1 − γ )τ ext , the dashed line for τ eff = (1 − γ )0.5 τ ext and the dash–dot linefor τ eff = (1 − γ )0.75 τ ext .

Table 1. Extinction cross-sections of the SN II dust models.
(A16)

(A17)

Please note: Wiley-Blackwell are not responsible for the content orfunctionality of any supporting materials supplied by the authors.
Any queries (other than missing material) should be directed to thecorresponding author for the article.
This paper has been typeset from a TEX/LATEX file prepared by the author.

C 2010 The Authors, MNRAS 412, 1070–1080
C 2010 RAS
Monthly Notices of the Royal Astronomical Society 