Dust in Supernovae and Supernova Remnants I: Formation Scenarios

Supernovae are considered as prime sources of dust in space. Observations of local supernovae over the past couple of decades have detected the presence of dust in supernova ejecta. The reddening of the high redshift quasars also indicate the presence of large masses of dust in early galaxies. Considering the top heavy IMF in the early galaxies, supernovae are assumed to be the major contributor to these large amounts of dust. However, the composition and morphology of dust grains formed in a supernova ejecta is yet to be understood with clarity. Moreover, the dust masses inferred from observations in mid-infrared and submillimeter wavelength regimes differ by two orders of magnitude or more. Therefore, the mechanism responsible for the synthesis of molecules and dust in such environments plays a crucial role in studying the evolution of cosmic dust in galaxies. This review summarises our current knowledge of dust formation in supernova ejecta and tries to quantify the role of supernovae as dust producers in a galaxy.


Introduction
The solid grains in space which may vary in size from the dimension of a few molecules to the order of a few microns, are broadly classified as cosmic dust. Cosmic dust can be of diverse chemical types, which includes carbon (amorphous and graphite), silicates (olivines, pyroxenes, both amorphous and crystalline), metal oxides (alumina, fayalite, magnesia, spinel), silica, pure metals (iron, magnesium, silicon), metal sulphides, carbides and some other species yet to be determined with certainty (Molster et al. 2010;Cherchneff 2013). Dust accounts for about 1-2% of the total baryonic mass in the universe. On the other hand, it is responsible for up to 50% of the total radiation from some galaxies, which rightfully justifies its importance. These dust grains are efficient in absorbing and scattering UV, optical and near-infrared (IR) radiations, and re-emitting the absorbed energy in the mid-IR and submillimeter (submm) wavelengths (Martin 1978;Li 2005). Owing to their continuum opacities, the dust grains absorb and emit over a large swathe of the electromagnetic spectrum (Draine 2003). Therefore, dust in space profoundly influences the spectra and colours of stars and galaxies. Moreover, cosmic dust acts as the basic building block for the planets in protoplanetary disks around main-sequence (MS) stellar systems (Brownlee 1985;Dwek 2006). The branches of science related to the study of cosmic dust has flourished over last four decades (Grün et al. 2001). It has brought together scientists (astronomers, meteoriticists, physicists, chemists, mineralogists and many more) from various disciplines and expertise to focus on the development of the field.
The earliest discoveries of molecule and dust clouds, leading to interstellar reddening were inferred in the first part of the 20th century (Barnard 1919;Trumpler 1930). However, only after the advent of IR astronomy in the 1960s, we have made a fast progress towards the current state of the art. The infrared region of the electromagnetic spectrum ranges from 700 nm (the red edge of the visible spectrum) to approximately 1 mm in wavelengths and is further subdivided into near-IR, mid-IR and far-IR regions. The biggest challenge to the development of IR astronomy, is that only a small part of the IR spectrum can actually be observed from ground-based facilities, owing to the absorption by earth's atmosphere, mainly by water vapour. In order to overcome the impediment of atmospheric shielding, the following schemes have been adopted: (a) to set up telescopes on the top of high altitude mountain peaks, (b) to raise the telescope above the atmosphere by helium balloons, (c) to attach the telescopes in jet transport aeroplanes, and (d) to launch IR satellite telescopes (since 1983) (Walker 2000). Later on, the ground-based telescope facilities with IR capabilities have also been developed. Some of the most important telescopes, operating in the near-IR, mid-IR, far-IR and submm, launched in the 21st century are the Hubble Space Telescope (optical and near-IR), Spitzer space telescope (3-180 µm), AKARI (2-200 µm), Herschel space observatory (55-672 µm) (until 2013) and WISE (3-25 µm). The Atacama Large Millimeter/submillimeter Array (ALMA), situated high on the Chajnantor plateau in Chile, is a high resolution millimeter-submillimeter telescope, which can study radiations from the coldest objects in the universe. The James Webb Space Telescope, under preparation for launch in 2018, is designed to boost IR astronomy with its unprecedented resolution and sensitivity in 0.6-27 µm range. The energy radiated in the IR is largely thermal in nature. Therefore, IR photometry and spectroscopy have become the crucial tools to study the dust continuum from various sources, vibrational transitions in molecules, signatures from exoplanets and the high redshift universe.
Using the newly developed techniques, confirmed evidence for the presence of cosmic dust has been established through various studies over the past couple of decades. Having said that, it is equally necessary at this point to investigate the potential sources of origin Fig. 1 Hertzsprung-Russell diagram showing luminosity-temperature relation of stars at their evolutionary phases as function of main sequence initial mass. The potential stellar sources of cosmic dust are marked with circles. The core-collapse supernova (SN), the final fate of the RSG (red super giant) phase as shown in the diagram, is a major source of dust formation for dust present in space. Cosmic dust is made up of large clusters of chemical compounds. Hence its synthesis requires chemical processes, which are sensitive to the local gas densities and temperatures. The chemical processes that are prevalent at moderately high temperatures ∼1000-2000 K and high gas densities govern the dust formation scenario. Such physical conditions are mainly encountered in evolved stellar outflows or winds and ejecta of novae and supernovae (Cherchneff 2010(Cherchneff , 2013Kozasa et al. 2009). However, the required temperature and densities in these evolving circumstellar environments persist only for a short timescale making the entire mechanism more complex.
The evolved stellar environments responsible for dust formation in galaxies are presented in Fig. 1 with their respective position on the Hertzsprung-Russell diagram. The quiescent mass outflows of the low mass (<8 M ) stars in their Asymptotic Giant Branch (AGB) phase are thought to be major dust producers in galaxies (Whittet 2003;Gobrecht et al. 2016). Core-collapse SNe, the ultimate fate of massive (8-30 M MS mass) stars, marked in the figure (not in the scale of the H-R diagram), are a prime source for dust formation (Wooden et al. 1993;Sugerman et al. 2006;Dwek 2006;Kozasa et al. 2009;Cherchneff 2013;Todini and Ferrara 2001;Szalai et al. 2011). Additionally, periodic mass ejections in massive Luminous Blue Variable (LBV) stars (Hillier et al. 2001;Ferland et al. 2005) and mass-loss events in the Wolf-Rayet stars (Harries et al. 2004) are also associated to dust formation in the galaxies. However, AGB stars and CCSNe are the known to be the main stellar sources. In case of supernova ejecta, where the stellar core has already collapsed, the physical conditions evolve very fast. Therefore, the entire process of dust synthesis occurs in less than a decade. Synthesis of dust in a suitable environment proceeds in two steps: the nucleation and the condensation (Cherchneff and Dwek 2009;Sarangi and Cherchneff 2013). Out of all the potential sources, this chapter aims to discuss and quantify the contribution of supernovae as dust producer in the galaxies.
Dust in dense clouds with a visual extinction exceeding ∼3 mag is characterized by solid cores surrounded by ice mantles. The presence of ice mantles in the dense interstellar clouds compared to the small dust present in the diffuse ISM confirms the growth of grains in the clouds. The grains grow in mass over the lifetime of the molecular clouds. For significant growth of dust grains in the ISM, the necessary conditions are: (a) the presence of preexisting grain cores on which accretion takes place, (b) the accretion timescales must be shorter than the lifetime of the clouds (Dwek and Cherchneff 2011). Interstellar dust grains are subject to various turbulent forces over the span of its lifetime, which either reduce the sizes of the grains or completely destroy them (Dwek 1998;Jones 2004). According to theoretical models, for dust grains to survive in the ISM, a continuous source of newly formed in the ejecta, with a bottom heavy size distribution. However, that being the case, the newly formed dust grains are almost likely to get completely destroyed by the reverse shock in course of their injection into the ISM .
Type Ib/c: These SNe types are characterised by the absence of H and Si II lines and the presence of He I in case of Type Ib, where as in Type Ic He I lines are also absent (Barbon et al. 1999;Matheson et al. 2000b). These objects are identified as core collapse SNe of massive stars (>25 M ) which have been stripped of their hydrogen envelope. They show evidence of shock interactions with dense circumstellar environments (Chevalier 1982) at early times after explosion. Presence of dust continuum in the spectra of Type Ib SN 1990I was first reported by Elmhamdi et al. (2004). Further, peculiar Type Ib SN2006jc has shown dust signatures as early as day 50 post explosion (Nozawa et al. 2008;Di Carlo et al. 2008), which originate from the dense shell formed in the post-shock CSM. Mid-IR continuum recorded by Spitzer from Type Ic SN2005at and SN2007gr confirm the presence of dust in the environment (Kankare et al. 2014). The source of the emission can be either an echo or the expanding ejecta, it is yet to be understood with clarity.
Type IIb: Type IIb's are core collapse supernovae which have early time signatures similar to Type II's (presence of prominent H lines) and a late time spectra similar to Type Ib's (Matheson et al. 2000a). The decline of the lightcurves in Type IIb SNe are much faster compared to other Type II, which implies a smaller mass of the ejecta. Due to the absence of the outer H-envelope after explosion, the ejecta expands faster and the typical gas densities are lower (Nozawa et al. 2010) when compared to a Type IIP. Therefore, Type IIb supernovae do not form dust as efficiently as Type IIP's Cherchneff 2014, 2016). IR and submm observations of galactic SN remnant Cassiopeia A of Type IIb have revealed the presence of dust ranging between 10 −4 M (Douvion et al. 2001;Dunne et al. 2003) to 0.08 M (Barlow et al. 2010;).

Type IIP:
The SNe Type II-P (Plateau) is the most abundant type of core-collapse supernovae in the local universe (Barbon et al. 1979). It differs from Type I's by the presence of prominent H lines in the spectra. The light curve is characterised by a flat plateau type nature during the first few months after explosion. Massive stars in the 8-25 M range can result to Type II-P supernovae at the end of their lives. The presence of a large H-envelope as the outer shell of the ejecta makes the outflowing ejecta massive and slow moving. Moreover, the densities of the expanding He core remain ∼10 2 /10 3 higher than Type IIb's. Based on the degree of stripping of the envelope, SNe can be arranged in sequence as Type IIP-IIL-IIb-Ib-Ic (Nomoto et al. 1995) in the ascending order. Also, this is the order of importance in terms of dust synthesis. Therefore, the ejecta of Type IIP SNe is the most efficient and copious source of dust among all SNe types (Szalai and Vinkó 2013;Sarangi and Cherchneff 2015;Cherchneff 2013;) and our discussion in this chapter will be largely based on Type IIP's.
Type IIn: SNe Type IIn, introduced by Schlegel (1990), are characterised by the presence of narrow (∼100 km s −1 ), intermediate (1-4 × 10 3 km s −1 ) and broad (10-15 × 10 3 km s −1 ) velocity width components. The narrow component originates from the slow moving CSM before being traversed by the supernova shock. The source of intermediate width component is the post-shock region of the CSM, where as the broad component arises from the fast expanding SN ejecta . Observations indicate early interactions between the ejecta and the circumstellar medium (Smith et al. 2012), which lead to a velocity shear followed by shocks. Many recent studies have reported the presence of warm dust in such environments (Fassia et al. 2001;Mauerhan et al. 2013;Smith et al. 2012;Gall et al. 2014).

Observations of Dust Formation in Core-Collapse Supernovae
Using modern high-resolution telescopes, a lot of important information can be derived from a single observation and a general trend for dust formation in supernovae can be studied. The formation of warm dust in the ejecta is characterised by the following: (a) an excess in the mid-IR region of the spectral energy distribution due to thermal emission (Roche et al. 1991;Bouchet and Danziger 1993), (b) a sharp decline in the optical light curve , (c) a progressive and systematic blue-shift of emission line profiles when the receding part of the ejecta is increasingly blocked by newly formed dust ), (d) diminishing line emission of elements and molecules compared to the adjacent continuum, owing to depletion in dust (Bouchet and Danziger 1993;). The first two cases might also arise either due to an IR echo from pre-existing dust in the circumstellar material (Smith et al. 2012) or from the escape of radioactive luminosity due to decrease in optical depth in the ejecta. But the occurrence of three or all of the four phenomena concurrently, clearly indicates the formation of new dust in the ejecta.
There are several local SNe that observers have closely monitored over the past few decades. We shall now discuss the findings from some important SNe with regard to their history of dust formation and its evolution in the ejecta.

Dust Formation in the Ejecta of SN 1987A
The explosion of a core-collapse SN was detected in the Large Magellanic Cloud, only 50 kpc away on the 23rd of February, 1987 (Kunkel et al. 1987). Since its detection, studies of SN 1987A have led to many discoveries of SNe and SNRs.

Historical Observations of Dust in SN 1987A
One of the major discoveries associated with the SN 1987A is the detection of newly formed dust ). From about 615 days, thermal emission of dust has been detected at near and mid-infrared (IR) with an increase in contrast and a shift towards longer wavelength in time (Fig. 3). The reported dust mass was 10 −4 M at day 775 with the temperature of 300-400 K. However, the possibility of larger dust mass was suggested because dust could be optically thick (Wooden et al. 1993). At almost the same time (∼530 days), the optical light curve declined at a much steeper rate than the predicted light curve accounting for the decay of 56 Co and 57 Co (Fig. 4: Whitelock et al. 1989). It has been interpreted that, the dust grains that are formed, are causing the extinction in optical wavelengths . Furthermore, asymmetric line-profile has also been detected simultaneously.  modelled the blue-shifted line profile absorbed by the dust that has condensed in the ejecta.
It had been suggested that dust could be present in clumps instead of being smoothly distributed across the ejecta Wooden et al. 1993). If a large mass had been evenly distributed in the ejecta, UV and optical energy would have been totally absorbed resulting in a black out (Hoyle and Wickramasinghe 1970;Dwek 1988). In order to avoid such black out, dust must be present in clumps so that UV and optical light can escape between the clumps. If dust is located in clumps, and if individual clumps are assumed to be optically thick, it is difficult to estimate the actual dust mass. It possibly suggests that the dust masses as large as 0.5 M  or 1 M (Dwek 1988). Whether it is feasible to accommodate such a large amount of dust in early days still remains a matter of  Wesson et al. (2014), based on observational data from Wooden et al. (1993), Harvey et al. (1989), Bouchet and Danziger (1993), Dwek et al. (1992), Matsuura et al. (2011) Fig. 4 The light curve of SN 1987A. The measured luminosity (squares) more sharply declined at about day 530, than the predicted light-curve from 56 Co and 57 Co decay. This further decline was attributed to the newly formed dust in the ejecta (Whitelock et al. 1989) debate (Ercolano et al. 2007;Wesson et al. 2014;Bevan and Barlow 2016;Dwek and Arendt 2015).
Since its detection, thermal emission of dust had been monitored at mid-infrared up to day 1,316 (Suntzeff and Bouchet 1990;Dwek et al. 1992). The following study was at day 6000 with nearly a 13 year gap in between the observations (Bouchet et al. 2004). By day ∼1150, the temperature had dropped to 150 K with an inferred dust mass of 5 × 10 −4 M (Dwek et al. 1992).
Ten years later, mid-infrared observations found dust thermal emission from the circumstellar ring (Bouchet et al. 2006). The ring composed of the material that were expelled from the progenitor star when it was in the red supergiant phase about 20,000 years ago (McCray 1993). The fast (∼2000 km s −1 ) blast wave from the explosion has collided into the slowly expanding (∼100 km −1 ) circumstellar ring. The collisions and subsequent shocks yield X-radiation heating dust grains with 300-400 K radiation at mid-infrared. The infrared luminosity of the ring dust continued to increase concurrently with a constant increase of the X-ray brightness, suggesting destruction of dust grains in the ring ). However, since day 8,500, infrared brightness started declining, while X-ray luminosity still increased. That is almost the time when the knots around the known ring structure started to brighten up, so that SN waves appeared to have overtaken the main part of the existing ring .

Detection of a Large Mass of Dust in the Ejecta
In 2010, the HERSCHEL SPACE OBSERVATORY discovered far-infrared emission towards SN 1987A (Fig. 5: Matsuura et al. 2011). The inferred dust mass was 0.4-0.7 M with the temperature of ∼20-30 K. The only place that can form such a large mass of dust would be the ejecta, where a large mass of refractory elements are present. The confirmation of the presence of dust in the ejecta was made by ALMA, soon after the HERSCHEL observations. ALMA high angular resolution image (angular resolution of 0.3 ) enabled to resolve SN 1987A system at submillimeter wavelengths, detecting thermal emission from cold dust towards the ejecta , with a mass of ∼0.5 M (Matsuura et al. 2015). Detailed analysis of ALMA data by  suggested that there could be an extra component beyond dust at millimetre wavelength, which can be associated with the pulsar winds.
Importantly, the ejecta's expansion velocity of SN 1987A might appear to be too fast for dust formation. Typically dust grains were predicted to be destroyed by SN winds with velocities exceeding 100 km s −1 (Jones et al. 1994). On the other hand, the velocities of SN 1987A ejecta are typically 2000 km s −1 (full width of half maximum (FWHM); McCray 1993) much larger than the dust destruction velocities. Dust grains can form and survive in such a high velocity environment, because in the rest frame of dust grains, the velocity is much smaller. The expansion velocities reach 2000 km s −1 a few months after the explosion, and following that, the velocities were measured to be constant (McCray 1993). Dust grains have formed only after these expansion velocities settled. Therefore currently in SN1987A, the dust grains are expanding at much faster speed without being impacted by shocks.

Look Back at the Historical Measurements of Dust in SN 1987A
The detection of a large mass (∼0.5 M ) of ejecta dust at ∼ day 8500 reignites the debate if such a large mass of dust is present from early days, starting from day 650. The inferred  (Matsuura et al. 2011). Spitzer 8 µm in blue, Spitzer 24 µm in red and Herschel 250 µm in green. LMC field stars are seen in blue, while SN 1987A, which contain significant mass of cold dust, is seen as a green point source in this image dust mass in early days was 10 −4 M at day 775, while measurements at late day (∼ day 10,000) suggests ∼0.5 M . There are two possibilities that can resolve the disparity. The first possibility is that inferred dust masses in early days were underestimated because the fitting assumed optically thin dust, but the emission was actually optically thick and a large mass of dust was present from the early days itself (Dwek and Arendt 2015). Alternative possibility is that the dust mass has increased over time by growing dust grains (Wesson et al. 2014;Bevan and Barlow 2016). Using currently available observational data, it is not possible to disentangle these two cases.

Molecules in SN 1987A
SN 1987A is also the first SN where molecule formation was reported. The emission of CO ν = 2 bands at 2.3 µm started appearing possibly 112-136 days after the explosion, with much more recognisable features from day 192 (Spyromilio et al. 1988;Catchpole et al. 1988). CO ν = 2 bands at 2.3 µm and ν = 1 bands at 4.6 µm continued to be detected until about day 600 and then disappeared from infrared spectra (Wooden et al. 1993).
More than 20 years after, CO was detected at millimetre wavelength after high sensitivity millimetre telescope, the ALMA, started operation in 2011 (Kamenetzky et al. 2013;. The analysis of these rotational bands suggests an excitation temperature of 20-50 K, which is much lower than the estimated temperature of 1200-4000 K in early days (up to day ∼600; Liu et al. 1992;Liu and Dalgarno 1995). It has been perceived that initially the gas temperature is high enough to excite vibrational-rotational bands. When the gas temperature drops, CO population is mainly in rotational transitions.
Not only CO but also detections of a few other molecules have been reported in SN 1987A. SiO fundamental bands at 4 and 8 µm appeared at about day 200 (Rank et al. 1988;Aitken et al. 1988;Roche et al. 1991;Meikle et al. 1989) and SiO rotational lines at millimetre wavelength were detected at day 9,173 onwards (Kamenetzky et al. 2013;. Potential detections of CO + , CS and H + 3 were also reported (Meikle et al. 1989;Miller et al. 1992) though not yet well established (Cherchneff and Dwek 2009). In later epoch, vibrationally excited H 2 at 2 µm was found since day 6,832 , and HCO + and SO at millimetre wavelength at day 10,037 .
The ALMA offers angular resolution of 50 miliarcsec images of CO and SiO in the ejecta (Abellán et al. 2017) as shown in Fig. 6. It is found that the CO and SiO clumps in the and  images overlaid on the HUBBLE SPACE TELESCOPE Hα image of the circumstellar ring of SN1987A. The figure is published by Abellán et al. (2017) and the HST image taken from  ejecta forms torus shape. These molecules fill in clumps. Three-dimensional models of SN explosions predict that Rayleigh-Taylor instabilities, triggered by the explosion, break the ejecta material into clumps (Wongwathanarat et al. 2015). After over 20 years of evolution in the post-explosion era, high resolution of the ALMA enabled us to resolve these clumps.

Cassiopeia A-Forward Shocks, Reverse Shocks and Unshocked Dust
Cassiopeia (Cas) A is the remnant of Type IIb SN, with an approximate age of 335 years Krause et al. 2008). IRAS (Infrared Astronomical Satellite), an all sky survey mission launched in 1983, detected dust in the SNR, Cas A (Dwek et al. 1987). Together with SN 1987A, detection of dust in this young SNR was one of the milestone discoveries which confirmed that dust can be formed in core collapse SNe.

SN Dust or ISM Dust?
The total dust mass in Cas A is largely debated. This is because there is heavy extinction of ISM dust towards Cas A, and evaluation of foreground and background ISM dust could affect Cas A dust mass estimates. The reported dust mass from IRAS was 7 × 10 −3 M (Dwek et al. 1987). Using SCUBA/JCMT, a dust mass of 0.2-4 M was reported in submillimeter wavelength (Dunne et al. 2003). Fitting of the 8-850 µm observations including SPITZER MIPS, indicates a presence of 3 × 10 −3 M of dust (Hines et al. 2004). AKARI and BLAST 65-500 µm observations deduced a 0.003 and 10 M of hot (99 K) and cold (17 K) dust ). The analysis of HERSCHEL observations resulted in 0.075 M of dust with a temperature of ∼35 K ( Barlow et al. 2010). Developing SN template analysis on SPITZER and HERSCHEL found 0.04 M of warm SN dust and <0.1 M of cold dust, whose origin was unclear. It can either be SN ejecta or the ISM (Arendt et al. 2014). After subtracting ISM dust templates from the SED of SPITZER and HERSCHEL, 0.4-0.7 M of cold SN dust was still present in Cas A in the unshocked region ).

Dust Compositions
Another important finding in Cas A was the spatially varying dust compositions distributed within the SNR and the type of dust compositions correlated to specific atomic lines (Lagage Fig. 7 One of the SPITZER spectra obtained in Cas A. Strong 21 µm feature was detected, which was attributed to Mg protosilicate, amorphous SiO 2 and FeO (Rho et al. 2008(Rho et al. ) et al. 1996Rho et al. 2008;Arendt et al. 2014). The spatially resolved spectra of Cas A showed that the spectra with the '21 µm feature' are accompanied with a strong [Ar II] line, while the spectra with featureless dust emission is correlated to the [Si II] line ( Fig. 7; Arendt et al. 2014;Rho et al. 2008). The '21 µm feature' could be attributed to a combination of Mg protosilicate (Arendt et al. 1999), amorphous SiO 2 and FeO (Rho et al. 2008). This is an unique case where a specific dust feature has been identified in SNRs, with recent finding in another Galactic SNR, G54.1+0.3 (Temim et al. 2017).

CO Molecules in Cas A-Survived Reverse Shocks or Reformed in Shocked Region
CO molecules have also been detected in Cas A. Filter measurements at 2 and 4 µm suggested the presence of CO emission at the reverse shocked region of Cas A (Rho et al. 2009) and was later spectroscopically confirmed by . One possibility is that the CO formed in early days might have survived reverse shocks. Alternative possibility is that CO is re-formed. Wallström et al. (2013) detected high J lines in Cas A. The wide range of temperatures (400-200 K) and high column density (5 × 10 17 cm −2 ) suggests that these CO molecules are from the dense regions of the post-shocked gas, where CO molecules were reformed. Chemical models suggest that molecules can re-form in the post-shock region, but not dust (Biscaro and Cherchneff 2014).

Dust in the Crab Nebula
The SN explosion of the Crab Nebula was recorded in 1054. Unlike Cas A, whose bright optical and X-ray nebula is collisionally heated by reverse shocks (Gotthelf et al. 2001), the major heating source of the Crab Nebula is synchrotron radiation (Hester 2008). Dust found in the Crab Nebula was also heated by the synchrotron radiation (Owen and Barlow 2015).
The Crab Nebula is composed of many filament like structures. The wind from the pulsar compresses the expanding ejecta to form these dense filaments. Not only gas but also dust is associated with these filaments (Fig. 8: Fesen and Blair 1990;Gomez et al. 2012;Grenman et al. 2017).
There has been a range of dust masses, that is inferred from observation studies of the Crab Nebula. Marsden et al. (1984) reported thermal emission of dust from IRAS flux measurements at 25, 60 and 100 µm. The dust mass was inferred to be between 5 × 10 −3 -0.03 M .  analysed the INFRARED SPACE OBSERVATORY (ISO)'s 60, 100 and 170 µm images, together with the 850 µm synchrotron image from James Clerk Maxwell Telescope (JCMT) and estimated a dust mass of 3 × 10 −3 -0.07 M with temperatures of 45-50 K. Temim et al. (2006) used SPITZER 24 and 70 µm flux, together with ISO's photometric data, resulting in the estimated dust mass of 4 × 10 −3 -1.5 × 10 −2 M with a temperature of 74 K. Adding the photometric measurements at longer wavelength of 70-500 µm from HERSCHEL SPACE OBSERVATORY, Gomez et al. (2012) found that the dust emission from the Crab Nebula is much broader than a single temperature black body/grey body, which previous authors have used. Using a two component black body, Gomez et al. (2012) found that while their hot component has similar dust temperature and mass as those derived by Temim et al. (2006). Gomez et al. (2012) found that the additional cold component, emitting at longer wavelengths, is at a mass of 0.11-0.24 M with temperatures of 28-34 K. However, the question is if mid-infrared emission from warm dust can hide emission from cold dust, which dominates the total dust mass.
In order to fit the broad emission, Temim and Dwek (2013) invoked temperature gradient of dust grains because synchrotron radiations can make small grains hotter than larger grains due to difference in heat capacities and absorption coefficients. The resultant dust mass is comparable to the estimates from two black bodies (Gomez et al. 2012), as far as the same dust emissivities are chosen. Temim and Dwek (2013)'s model assumed that all dust grains are located at a fixed distance. Owen and Barlow (2015) advanced the distance issue, by applying full radiative transfer modelling to shell structure. Their fitted results on both ionised lines and dust emissions show that the dust mass ranges between 0.17-1.5 M depending on the morphologies of the dust species and their emissivities.

Molecules in the Filaments of the Crab Nebula
Vibrationally excited molecular hydrogen lines have been detected in the Crab Nebula (Graham et al. 1990), and the H 2 emission originates in the filaments (Loh et al. 2011). Loh et al. (2012 estimated the excitation temperature of H 2 to be 2000-3000 K, which is most likely collisionally excited. The first detection of noble molecular gas in space was made in the Crab Nebula, by the detection of 618 and 1235 GHz lines of 36 ArH + found in the HERSCHEL Space Observatory's spectra . This is an unusual case because majority of molecules had been found first in interstellar medium (ISM) or circumstellar medium of evolved stars and later identified in SNRs. On the other hand, this molecule was reported in a SNR first and then identified later in the ISM (Schilke et al. 2014). Argon is a noble gas, i.e. nonreactive, but its ionised state Ar + is extremely reactive and its reaction with H 2 leads to the formation of ArH + .

Galactic and Magellanic SNRs
Infrared surveys enabled to detect and analyse Galactic and Magellanic clouds SNRs with possibilities to analyse the dust within. Launched in 1983, the IRAS is the first whole sky survey at 12, 25, 60 and 100 µm and detections of 51 and 11 SNRs were reported from the survey (Arendt 1989;Saken et al. 1992). These achieved systematic studies of infrared emission from SNRs.
Launched in 2003, the SPITZER SPACE TELESCOPE used Galactic plane surveys (|b| < 1 • and 10 • < |l| < 65 • ; e.g. Benjamin et al. 2003) with two different instruments, IRAC and MIPS. These survey data were used to detect infrared emission from SNRs. Compared to the pixel scale of 2' of the IRAS map used by Arendt (1989), ∼2.5 arcsec resolution at 3.5-8.0 µm, 6 arcsec at 24 µm and 20 arcsec at 70 µm are great improvements, which helped in resolving ISM cirrus confusion towards the SNRs.  found 17 detections at 3.5-8.0 µm and colour analysis suggested that four of the SNRs appeared to have infrared emission associated with polycyclic aromatic hydrocarbons (PAHs). The emission of the remaining thirteen SNe were found to be due to molecular lines (H 2 and CO) or atomic lines ([Fe II], [Ar II] and Brα) caused by shocks. Pinheiro Gonçalves et al. (2011) searched for counterparts of Green's SNR catalogue (Green 2004) and reported the detection of 39 SNRs at 24 and 70 µm. The mid-infrared emissions of these SNRs have various origins that include dust and lines.
SPITZER's sensitivities and angular resolutions allowed us to study SNRs beyond the Galaxy such as those in the Magellanic Clouds. That includes dedicated studies of individual SNRs, 1E 0102.2−7219 (Stanimirović et al. 2005;Sandstrom et al. 2009) and the studies (Seok et al. 2013;Lakicevic et al. 2015) using the Magellanic Clouds survey . Seok et al. (2008) also studied Magellanic Clouds SNRs using the AKARI infrared whole sky survey. Correlation between infrared emission and X-ray emission suggests that 24 and 70 µm flux of the SNRs are collisionally heated (Seok et al. 2015) with the possibility of dust destruction by sputtering (Williams et al. 2006).

Extra-Galactic SNe
The detections of ejecta dust in SNe beyond the magellanic clouds have been reported by finding asymmetric structures in the optical emission lines (e.g. Elmhamdi et al. 2003b). After the explosion, the ejecta is characterised by broad emission lines and at some stage the blue wing of the broad ejecta lines might be affected by extinction. This extinction indicates the formation of dust in the ejecta. Additionally, sharp decrease in the optical light curve can be found (Elmhamdi et al. 2004). These two methods have been effective to detect ejecta dust in extragalactic SNe.
Presence of high sensitivity instruments in Spitzer Space Telescope and the development of mid-infrared instruments on eight meter telescopes enabled to detect mid-infrared thermal emission of dust from extragalactic SNe (e.g. Sugerman et al. 2006). Detections of midinfrared dust emission help to quantify the dust masses and temperatures. Typical ejecta dust masses reported in SNe younger than day 10,000 are between 10 −6 -10 −3 M (e.g. Sugerman et al. 2006;Meikle et al. 2007;Kotak et al. 2006;Szalai et al. 2016) and it has been suggested that the dust masses increase in time (Gall et al. 2011(Gall et al. , 2014. Because the measurements of late phase dust masses are limited to a few SNe (e.g. SN 1987A and1995N;Matsuura et al. 2015; van Dyk 2013) with a large gap in epochs, it  Owen and Barlow (2015). Publication of HST image?
is not yet clear whether the dust masses in the early phase is actually such a small quantity or a large mass of dust was missed due to being optically thick (Sect. 2.1.3).
The dust thermal emission can be associated with ejecta dust, light echo, circumstellar dust or dust from circumstellar material. Often the timing of the appearance together with optical line extinction or optical light curves are used in order to identify the origin of dust emission.

Molecules as Precursor to Dust Formation?
By exploiting extra-galactic SNe, the number of SNe where molecules and dust have been detected has increased. Chemical models predict, the time sequence of formation of molecules in SNe is followed by dust formation. The binding energy of CO is one of the highest among molecules found in astronomical environment, and formation of CO starts about day ∼200 with a sharp increase in mass up to day 500, followed by a gradual increase up to day 800 (Sect. 6.2.2: Sarangi and Cherchneff 2013). After that, the CO mass does not increase so much. Whereas the dust condensation starts at about day 200, followed by silica formation at day 400. We test this scenario by looking into the past observations of molecules, mainly CO, and dust.
First, we test link between molecule and dust formation. So far, molecular detections were reported in the spectra of twelve SNe (Summarised in Table 1; extracted from Matsuura 2017). CO detections suggested in photometry measurement only are excluded from the analysis. SN 1987A is included among the twelve, but not SNRs such as Cas A and the Crab Nebula. Figure 9 shows the statistics of dust detections out of twelve SNe. We divide the sample into two parts-six SNe before the Spitzer mission and six SNe after the Spitzer mission. The latter includes SN 2002hh, whose explosion was recorded before Spitzer launch in 2003, but has been observed with Spitzer, and SN 2013eg, which was type Ic and detected after Spitzer cryogenic phase. Dust detections at late epoch were reported at day 6067 (Bouchet et al. 2004) and at day 8500 and 9090 (Matsuura et al. 2011(Matsuura et al. , 2015 c Detection of dust thermal emission was reported (Barlow et al. 2005). Most likely it is associated with progenitor dust, but potentially ejecta dust References: Out of the twelve SNe with molecular detections, dust detection was reported in six SNe (Fig. 9). Three of the no-dust detected SNe are of Type Ic and Ib/c SNe, which have fewer dust detections reported so far (Matsuura 2017). Other SNe are all Type II, mainly Type IIP, with one IIn and unclassified subclass (Table 1). From Type II, dust detection was reported in six SNe out of nine, making the succession of molecule formation followed by dust formation very likely. In particular, after the launch of Spitzer, all Type II SNe with molecular detections ended with dust detections, so this scenario seems very likely.
Next, we test the timing of molecule and dust formation/detections. We assumed the reported days of CO detections and ejecta dust detections from Matsuura (2017) and the histogram of these days are plotted in Fig. 10. CO detections from one SN at different epochs are counted individually. In these panels, CO and ejecta dust detections in late epoch (>1000 days) of SN 1987A are not plotted.
The histogram (Fig. 10) shows that CO detections have been reported about 100-300 days. These are mainly detections at 2.3 µm and 4.6 µm CO fundamental and the firstovertone band. After that period, CO emission fades, making it difficult to detect at nearinfrared wavelengths. The second highest peak at 225-275 days is due to five separate CO observations of SN 1987A, distorting the frequency of CO detection rate during this period. While the ejecta gas cools down, the gas can not excite vibrational level of CO, so that nearinfrared vibrational CO bands disappear by ∼600 days. Instead, CO rotational lines appear at millimetre wavelength ranges and is at least reported by day ∼9000 (Kamenetzky et al. 2013;.
While CO detections have a gentle peak between 100-300 days, any such peak is not found in the timings of the ejecta dust detections. In Fig. 10, records of first ejecta dust detection per SN are plotted in blue and these reported first detections spread across 100-800 days. Dust can be associated with ejecta, light echo, circumstellar and ISM dust, and it is not always easy to disentangle these cases. Nevertheless, typically the highest chance of the ejecta dust detection appears to be between 600-700 days.

Indirect Evidence of Supernova Dust
Apart from direct observations of dust signatures in the mid-IR and submillimeter wavelengths in different local SNe as discussed above, there are also some indirect yet strong evidences of dust in SNe, which helps us understand and constrain the mechanism of dust formation in such environments. Among them, the dust case of high redshift galaxies (Maiolino et al. 2004;Dwek and Cherchneff 2011) and the isotopic anomalies in presolar grains from meteorites (Clayton and Nittler 2004;Ott 2001) provide us with a lot of important information on SNe dust, as explained below.

The High Red-Shift Galaxies
The reddening of the background quasars and Lyα systems at high redshift indicates the presence of dust in those galaxies (Pei et al. 1991;Pettini et al. 1994;Hines et al. 2006). Large masses of dust of the order 1-5 × 10 8 M is estimated to be present in Sloan Digital Sky Survey (SDSS) J1148+5251, a hyperluminous quasar at z = 6.4 (Bertoldi et al. 2003;Dwek et al. 2007). Figure 11 points to the high redshift galaxies on the cosmic time scale where presence of large dust masses is reported. The presence of large masses of dust plays pivotal role in the cooling of the primeval large gas clouds, eventually leading to the formation of next generation of stars (Cherchneff 2013). The Population III stars present in these galaxies are predominantly massive and a large section is assumed to be in binary form. The sources of such large masses of dust and the formation mechanisms are yet to be understood with confidence. Generally supernovae and AGB stars are taken as the prime dust producers in the galaxy. In order to explain the large masses of dust at high redshift galaxies, Dwek and Cherchneff (2011) assumed the contribution from all possible dust producers in the galaxy and its balance with the galactic star formation rate. Massive stars evolve much faster compared to the low mass stars. Therefore, in the timescales of the lifetime of the early galaxies, supernovae are estimated to be the main source of dust. If 1 M of dust is produced by corecollapse supernovae before getting reprocessed in the remnant phase and a top-heavy initial Fig. 11 Evidence of large masses of dust reported from the far-IR spectrum of hyperluminous galaxy SDSS J1148+5251 at high redshift z = 6.4 (Bertoldi et al. 2003;Dwek et al. 2007) (Photo courtesy: NASA WMAP) mass function is assumed for the primeval stars, then the dust masses observed at high redshift galaxies can be explained. However, theoretical models infer a mass of 0.1-0.15 M of silicate dust formed per massive pair-instability supernova as an upper limit (Cherchneff and Dwek 2009;Cherchneff 2010). Moreover, the destruction of dust by the reverse shock in the ejecta or by UV photolysis in the ISM are somewhat underestimated. AGB stars can only contribute as a major dust producer if the galactic age is assumed to be much larger (400 Myrs). Therefore, in case of high z galaxies, SNe are expected to have obvious advantages (Maiolino et al. 2004;Calura et al. 2014;Gall et al. 2011) over AGB stars in terms of relative contributions owing to the large abundances of massive stars and their shorter evolutionary timescales (Dwek 1998;Dwek and Cherchneff 2011).

Presolar Stardust from Meteorites
Laboratory studies have revealed the presence of surviving dust grains of circumstellar origin within primitive meteorites (Bernatowicz and Zinner 1997). The deviation from solar isotopic composition confirms the presolar origin of these dust particles. Further, depending on the isotopic abundances, the stellar source of the grains as well as some limitations on the physical conditions of the environment can also be inferred (Hoppe and Ott 1997). Dust formed in the ejecta of a supernova carries signatures of the nuclear processes from the pre-explosion and post-explosion phases of the progenitor (Ott 2001).
Out of several types of chemical species found in the meteorites, silicon carbide grains have been profoundly studied and grouped into categories (A, B, Y, Z, X etc.) depending on their possible sources of origin (Hoppe and Ott 1997). The presolar SiC grains of type X (∼1% of all presolar SiC samples collected from meteorites) has been attributed to have SN origin (Amari et al. 1997), as shown in Fig. 12, based on the enrichment of 12 C and 15 N, high abundance of 28 Si, and importantly the presence of 44 Ca, a remnant of short-lived radioisotope 44 Ti (Ott 2001;Hoppe 2010). Similar to silicon carbide, a few of the silicon nitride (Si 3 N 4 ) samples found also have isotopic signatures of that of SNe (Amari et al. 1997). Apart from these, a small subsection of presolar aluminum oxide grains (Nittler et al. 1998; Even though, SN has been established as a source of cosmic dust in many previous analyses, the isotopic compositions of the grains are not exactly coherent with the existing theoretical models (Herant and Woosley 1994) and also observational studies of the SN ejecta (Hughes et al. 2000). In order to match the isotopic signatures, selective yet microscopic mixing between the ejected materials from different zones of the ejecta is necessary (Nittler 2003;Zinner et al. 2010). In order to understand the mixing scales compatible with the isotopic compositions in the grains, observations of the SN ejecta at high spatial resolutions are required. New theoretical studies of the post-explosion ejecta using 3D hydrodynamic models (Hammer et al. 2010b;Wongwathanarat et al. 2015) also hint at more intense mixing between the stratified zones.

The Chemical Nature of Dust Grains
Compared to dust in earth's atmosphere, dust grains in space have relatively simpler chemical forms as they exist in more hostile and pristine surroundings. The chemical types of cosmic dust and their precursor molecules are not well defined through observations due to lack of high-resolution mid-IR data (Cherchneff 2013). Composition of dust in space can be analysed by (a) the laboratory study of presolar grains from meteorites, (b) fitting of mid-IR spectral energy distribution using known opacities of probable dust species, (c) modelling of dust condensation from its gas phase precursors. Cosmic dust from circumstellar sources can be classified chemically as (a) oxides (silicates of quartz or metal oxides), (b) carbides,  (2013), gives a summary of known dust species and their mid-IR bands with the references. In case of SN ejecta, hydrogen mixing is understood to occur in macroscopic scales only (Joggerst et al. 2010;McCray 1993). Therefore, dust formation pathways are controlled by simplified H-free chemistry leading to a small set of possible dust species. From bottomup perspectives, the ejecta chemistry is modelled to predict the types and abundances of molecules and dust grains in the metal-rich ejecta (Bianchi and Schneider 2007;Cherchneff and Dwek 2009;Sarangi and Cherchneff 2013;Sluder et al. 2016). However, in order to conclusively comment on the type(s) of dust present in a certain environment, radiative transfer (RT) modelling of the ejecta and fitting of the observed spectra are required (Sugerman et al. 2006;Ercolano et al. 2007). For optically thin media, the spectral features of the dust grains are prominent whereas in case of large optical depths, the dust continuum behaves more like a black body spectrum making the detection of a particular dust type more ambiguous. In the following sections of this chapter, the theoretical aspects of the study of SNe have been addressed in more detail.

Theoretical Models
Observational studies have confirmed the presence of dust signatures in SN ejecta, which can either be contribution of echo from pre-existing dust or newly formed dust in the ejecta. In order to understand the origin of dust in such environments, we shall now discuss the mechanisms which aid the synthesis of dust grains in SN ejecta. Attempts to theoretically model an astrophysical environment or event focus on the temporal and spatial variance of the following: (a) the morphology/geometry, (b) relevant physical processes, (c) chemical pathways. To study these aspects with respect to a typical SN ejecta, several simplifications/ approximations are necessarily imposed on the models. The following subsections describe the existing models of SN ejecta and their pertinence compared to the actual environments.

Models Based on CNT and LTE
The earliest approach to model dust synthesis in SN ejecta involved a Classical Nucleation Theory (CNT) based formalism. The CNT based approach was developed using the 'liquiddrop model' introduced to study water droplets in earth's atmosphere (Feder et al. 1966). In this model, the nucleation of gas-phase metals is addressed by a free-energy controlled nucleation current that leads to the formation of critical clusters. The initial growth of these clusters by the attachment of new molecules is an energetically unfavorable process as the surface energy of these species (the missing binding energy of surface molecules) exceeds the volume binding energy. The work required to form a single cluster of k-monomers is a function of (a) the change in the energy of molecules due to its transition from gas to solid phase, (b) the binding energy of the involved species, and (c) the surface energy of the cluster. In macroscopic scale, the energy to form a cluster of k-monomers is equal to the change in the Gibbs free energy of a gas containing k-number of monomers that is undergoing a phase transition at a given pressure and temperature, which is given as, where W (k) is the work required, k B is the Boltzmann constant, S ≡ P /P v is the supersaturation ratio and u s is the surface energy term. These supersaturated compounds (the critical clusters), synthesized in this process, are assumed to be in thermodynamic equilibrium (TE) with the solid phase of dust. The accretion and detachment rates of gas monomers are therefore calculated considering a steady state environment.
Dust formation in SNe was discussed in an empirical approach by Clayton (1979). Thereafter, a pioneering attempt to model dust formation was made by Kozasa et al. (1989) based on the early results of SN1987A. This study considered a fully microscopically mixed ejecta where the gas and dust phase was connected by CNT and TE assumptions. Following that, CNT based models were applied to study dust formation in SNe (CCSN & PISN) by Todini and Ferrara (2001), , Bianchi and Schneider (2007), and more recently by Marassi et al. (2015) for various metallicities and mass of progenitors as well as for various levels of mixing.
Models based on kinetic nucleation theory (KNT) provide an alternative formalism to study dust formation in such environments. In this approach, the number densities of clusters are traced explicitly. The condensation of grains is computed using accretion, where as the destruction of clusters occur through evaporation. The rate coefficients are calculated using kinetic balance. Therefore, it does not necessarily consider a steady state environment for the nucleation processes. Several studies over the last decade have applied this approach to model dust formation in SN ejecta (Nozawa et al. 2003(Nozawa et al. , 2008Nozawa and Kozasa 2013).
However, the use of nucleation theory to study dust formation in environments such as SN ejecta was questioned first by Donn and Nuth (1985) and later elaborated by Cherchneff and Lilly (2008). CNT is applicable to steady-state conditions, which are usually not found in dynamic environments out of equilibrium (Cherchneff and Dwek 2009). All CNT-based models ignore the non-equilibrium chemistry related to the formation of molecules and dust clusters and the specific physics of ejecta where radioactivity greatly impacts the gas-phase chemistry through cascades of Compton electrons. Moreover, the CNT approach includes the quantities like specific surface energy, sticking coefficient, supersaturation ratio, partial vapour pressure, critical clusters, etc. in analysing the dust formation scenario. These all are essentially bulk properties to matter. However, the critical clusters sizes are of molecular dimensions and therefore, application of the bulk properties like surface energy or sticking coefficient on such small particles is not appropriate (Cherchneff 2013). Additionally, in CNT-models, the condensation sequence is based on assuming condensation temperatures at equilibrium as references. Estimating the condensation temperatures at equilibrium for the seed nuclei is also not an appropriate approach in most cases, as the temperatures of the newly formed grains heavily depend on the ambient gas temperatures and densities, controlled by the amount of radioactive 56 Ni produced at the explosion. Even though the KNT based models can take the non-steady state environment into account, they still ignore the actual chemical pathways that lead to the synthesis of the dust-seed nuclei. Further, important condensation mechanisms such as coagulation are also not included in such analyses.

Chemical Kinetics Model
A formalism based on chemical kinetics can overcome most of the limitations imposed by previous models. In this case, dynamics and chemistry control the synthesis of condensates (Cherchneff and Lilly 2008;Dwek 2009, 2010;Cherchneff 2013, 2015;Cherchneff 2014, 2016;Sluder et al. 2016). A chemical kinetic approach deals with the dynamics of chemical reactions and the conditions affecting the reaction rates. It is controlled by the concentration of the reactants, temperature of the medium and the enthalpy of the individual reactions. The SN ejecta has a very dynamic and unstable physical environment. The passage of the supernovae blast wave through the helium core and the deposition of radioactive energy in the ejecta creates a complex physical environment whose chemistry is neither at steady state nor can be assumed to be at thermodynamic equilibrium (TE). Following the explosion, the ejecta moves at a velocity ∼2000 km/s and the density and temperature of the ejecta drop rapidly as power-laws. The rich chemical composition of the ejecta is controlled by a large set of chemical reactions with diverse reaction rates. All these chemical reactions occur simultaneously and are coupled to each other. Each individual species is linked to the chemical network through a series of formation and destruction pathways, which comprise of both thermal and non-thermal processes. Fig. 13 A cartoon shows the chemical paths leading to the synthesis of molecules, clusters and dust grains in the ejecta through nucleation and condensation phases Therefore, with the alteration of temperature and density in the ejecta, the relative importance of specific reactions in the network may also change. These factors not only ensure the complexity of the expected chemical yields but also create a case which is far from the steady state condition.
The chemical pathway proceeds through simultaneous phases of nucleation and condensation (Fig. 13). The nucleation phase comprises of chemical reactions which lead to the synthesis of molecules and small clusters. The condensation phase couples the gas phase chemistry to the solid phase of dust. The small clusters formed in the gas phase condense through coagulation and coalescence to form large grains of dust, provided suitable conditions are prevalent. Therefore, we first discuss the physical environment of a typical SN ejecta before going through the detail chemical scheme.

Physical Conditions in the Ejecta
We choose a typical 15 M Type IIP SN as a standard case study considering Type IIP SN ejecta to be the most conducive to dust formation. The ejecta of a typical PISN has some morphological differences with a CCSNe owing to the metallicity of the Population III stars. The contrasts in the physical conditions of such objects will be addressed separately. For the sake of simplicity, a 1D stratified model of the ejecta is assumed and we shall ignore the details of the explosion mechanism in this discussion. The abundances of post-explosion elements are provided by nucleosynthesis studies and hydrodynamic simulation models of the ejecta (Woosley and Weaver 1995;Thielemann et al. 1996Thielemann et al. , 2011Thielemann 2017;Limongi et al. 2000;Rauscher et al. 2002). The final yields from these models correspond to a timescale of 6-7 hours from the explosion. These elemental yields are then used as initial inputs to the chemical models of the ejecta (Rauscher et al. 2002).
Following the expulsion of the helium core, the ejected materials are crossed by the blast wave, which deposits energy in the gas. The neutrino-driven convection during the collapse leads to the Rayleigh-Taylor instabilities and mixing in the helium core. The mixing scenario continues over a few hours after the explosion and leads to partial fragmentation of the core (Joggerst et al. 2010). Thereafter, the chemical species produced by nucleosynthesis flows ballistically as carried by the blast wave. Hence the abundances of the elements presented as functions of mass zones (Rauscher et al. 2002) in the ejecta (Fig. 14) can be safely assumed to remain constant over the next few days until chemical reactions becomes relevant. The explosion energy, E exp (∼10 51 ergs) is assessed using SN1987A as reference (Woosley 1988), which gives an estimate of the ejecta velocity. Based on the distribution of elements, the ejecta can be stratified into several regions or zones, assuming no further microscopic mixing takes place after a few hours from explosion.  Fig. 14 The post-explosion elemental abundances inside the He-core as a function of mass zones for a 15 M progenitor star which has exploded as a Type IIP supernova (Rauscher et al. 2002) and Cherchneff 2013). In this study, the species in individual zones are assumed to be microscopically mixed. The dominant chemical processes in each zone vary depending on the temperature, density and the abundances.
The SN light curve is powered by radioactive decays (Woosley et al. 1989). The energy deposited by the SN explosion-flash is radiated away in a couple of weeks (Woosley 1988;Blinnikov et al. 2000;Ensman and Burrows 1992;. After that, the bolometric luminosity, which includes the optical, IR and UV emissions, is controlled by interactions between high energy radiations with the supernova ejecta. There are several radioactive nuclei produced in the explosion. On the basis of the abundances and the halflives, the most important radioactive nuclei in the ejecta are 56 Ni, 57 Ni, 55 Co, 44 Ti and 22 Na (Kozasa et al. 2009). 56 Ni has a very short lifetime of 6 days and quickly decays to 56 Co through electron capture. 56 Co in turn decays to 56 Fe with a half-life is 77.12 days. The light curve of the first few years of SN1987A is dominated by 56 Co and a mass of 0.075 M reproduces the bolometric luminosity (Woosley et al. 1989). Hence an equal amount of 56 Ni is estimated to be produced in the explosion. The 56 Co decays through electron capture ( 56 Co −→ 56 Fe + γ + ν e ) and produces energetic γ -rays and electron neutrinos. The energetic γ -rays interact with the ejecta and get degraded to hard X-rays through Compton scattering. The fast Compton electrons thermalize the ejecta through heating, excitation and ionisation (Liu and Dalgarno 1995). The hard X-ray gets further degraded to soft X-rays and UV photons through cascades of inverse Compton scattering and ionization-recombination processes. The entire energy of the γ -rays is not deposited in the ejecta. Study of the SN1987A light curve indicates the presence of macroscopic mixing of 56 Ni in the ejecta, which is also confirmed by the early emergence of γ -rays and X-rays within a few hundred years from explosion. The ejecta is presumably clumpy and hence a small fraction of the high energy photons can possibly pass through the less dense interclumps unabsorbed (McCray 1993). However, the opacity of the ejecta is found to increase rapidly at early times. Hence, to a first order approximation, the fraction of escaping energy becomes negligible. For the sake of generality, it can be assumed that the entire radioactive energy is uniformly deposited in the ejecta (Cherchneff and Dwek 2009;Sarangi and Cherchneff 2013). The rate of energy deposition (L γ ) in the ejecta by 56 Co is calculated following the study by Woosley et al. (1989) in terms of the γ -ray optical depth (τ 0 ), γ -ray mass-absorption coefficient (κ γ ), the deposition function (f (τ )) and the energy emitted per 56 Co nucleus (E γ ) (Liu and Dalgarno 1994;Cherchneff and Dwek 2009;Sugerman et al. 2006). The rate varies linearly with the Fig. 15 The schematic diagram, adopted from Cherchneff and Sarangi (2011), represents the time evolution of a typical CCSN ejecta. The scale in the bottom of the figure illustrates the temperature/density variation with post-explosion time. On the top left corner, the stratification of the stellar core and envelope is shown which retains the signatures from the burning phases of the massive star. The mixing and the Rayleigh-Taylor instabilities in the ejecta that occur until a few hours after the explosion, is also presented in the picture, adopted from Kifonidis et al. (2006) total mass of radioactive 56 Ni produced by the explosion Cherchneff 2013, 2015).
The ejecta is assumed to undergo a quasi-adiabatic expansion (Cherchneff and Dwek 2009) over a span of few years. The temperature profile is obtained with reference to observations and models based on SN1987A (Fransson and Chevalier 1989;Kozasa et al. 1989). The temporal variation of gas temperature in the ejecta is represented by an 'adiabatic' index γ (∼1.42) for SN Type IIP, as derived by Nozawa et al. (2010). The propagation of the shock causes a velocity shear in the post-explosion ejecta and within one day of explosion, the ejecta shell expands homologously (Shigeyama and Nomoto 1990;Nozawa et al. 2003). The initial gas densities and temperatures at ∼ day 100 (Sarangi and Cherchneff 2013;Biscaro and Cherchneff 2014) used for the model are derived from Nozawa et al. (2010). Figure 15, adopted from Cherchneff and Sarangi (2011), presents a cartoon illustrating the post-explosion ejecta expanding outwards from the stellar remnant.

Nucleation Phase
Molecules and clusters act as the building block of dust in supernova ejecta. The nucleation scheme, which deals with the formation and destruction of such species in the gas phase, therefore acts the bottleneck to the entire mechanism (Cherchneff and Dwek 2010;Sarangi and Cherchneff 2013). The He-core, rich in metal concentration and poor in H abundance, is the prime site of dust synthesis in SN ejecta. The environment in concern being H-poor, simplifies the relevant chemical pathways to a great extent. Moreover it diminishes the variety of chemical compounds synthesized as intermediate or end products. Table 3 provides the list of all possible chemical species (∼120) in the gas phase that are important for this study (Cherchneff and Dwek 2009;Sarangi and Cherchneff 2013;Biscaro and Cherchneff 2014). A chemical network is developed for the study comprising of all feasible chemical pathways Table 3 The set of atoms, ions, molecules, and molecular clusters important for the study Cherchneff 2013, 2015;Cherchneff 2014, 2016)  relevant to the hot and dense environments. The network consists of about 600 reactions of different kinds in the gas phase. They can be further categorized into thermal (termolecular or 3-body processes, thermal fragmentation, bimolecular reactions, ion-molecule reactions, radiative association etc.) and non-thermal processes (photon induced, Compton-electron induced etc.) (Cherchneff and Dwek 2009). The reactions are characterized by the respective reactions rates. Hence the temporal variation of number density of a certain species i located at mass zone M r in the ejecta at a certain time t , is given by, where P i represents the sum over the all the production (formation) processes whereas L i represents the sum over all the destruction processes. The term k ij represents the tempera- Fig. 16 The nucleation scheme for silicates leading to the formation of the enstatite and forsterite dimers, which act as the seed in the condensation formalism. The reactant species are given on the side of the directional arrows in the schematic diagram (Goumans and Bromley 2012;Zachariah and Tsang 1993) ture dependent reaction rates between the ith and j th species (Cherchneff and Dwek 2009). T 0 is a reference temperature chosen as 300 K for the study. A ij is the Arrhenius rate coefficient expressed in s −1 , cm 3 s −1 or cm 6 s −1 depending on the order of reaction and E ij is the activation energy in the unit of K −1 . The factor ν represents the time dependence of the pre-exponential factor of the general Arrhenius equation. Typical values of the νfactor range between −1 to 1. The gas phase chemistry is overall controlled by N number of stiff coupled nonlinear ordinary differential equations similar to (2), where N is the number of species that has been included in the study. This entire set of differential equations is solved simultaneously and the time evolution of the abundances for each species is estimated. The reaction rates are either theoretically calculated using translation theory or estimated in the laboratory. Standard chemical databases such as NIST, KIDA, UMIST, etc., are used as sources for the compilation of the network. The rate coefficients for the non-thermal processes are estimated using the rate of deposition of radioactive energy by the thermalized γ -rays in the ejecta (Cherchneff and Dwek 2009;Dalgarno 1994, 1995). The branching ratio between the energy deposited by electron through excitation, ionisation (A → A + + e − ) and fragmentation (AB → A + B + + e − ) depends on the mean energy per ion pair (W d = E p X i /N d , where N d and X i are the number of ion pairs produced by the species in this process and the abundance of that species in the ejecta). The detailed account of all the non-thermal processes related to SN ejecta can be found in Cherchneff and Dwek (2009), Dalgarno (1995) &Clayton et al. (1999). A stochastic, kinetically driven approach has been adopted to study the nucleation of gas molecules to small clusters. The vaporisation experiments in the laboratory and the study of ceramic synthesis in the flames provide insights into similar nucleation processes under nonsteady state environment. The important species that demand detailed analysis are oxygen rich clusters of magnesium silicates, alumina and Si m O n molecules, carbon chains, rings and fullerenes, silicon carbide, metallic sulphides and oxides, and pure metallic clusters. The upper limit of the cluster size in the nucleation network is not defined very specifically. The geometry of the molecules is the backbone to determine the final product of nucleation and the dust precursor in the network (Table 3 in bold font). The clusters that attain a three dimensional structure like a cage or twisted rings, has a larger collision cross section for coagulation (Cherchneff and Dwek 2010;Sarangi and Cherchneff 2015). Further, they are energetically more stable, and hence likely to be resilient under the destruction processes in the gas phase. In some cases, the favoured chemical routes for nucleation network are derived from quantum chemical calculations. Figure 17 presents the structures of a few stable Fig. 17 (1) (SiO) 5 (Reber et al. 2008), (2) (MgSiO 3 ) 2 (Goumans and Bromley 2012), (3) C 28 (Makurin et al. 2001), (4) Mg 4 (Cherchneff and Dwek 2010), (5) (Al 2 O 3 ) 2 (Gobrecht et al. 2016, private communications) gas phase molecular clusters, which acts as precursor seeds to the condensation phase. The nucleation schemes vary vastly depending on the type of species. Figure 16 is a schematic presentation of the nucleation pathways that lead to stable (SiO) n compounds and magnesium silicates (Mg 2 SiO 4 ) 2 . In a hydrogen free environment, pure carbon chains are formed of sp-hybridization up to C 9 followed by monocyclic and fullerenes. Cluster beam experiments suggest, the smallest stable structure of fullerene present in carbon vapour is C 28 (Dunk et al. 2012). All the possible nucleation networks relevant for such environments are described in Cherchneff and Dwek (2009) & Sarangi and Cherchneff (2013).
The model predicts the formation of several diatomic molecules in the ejecta, such as CO, SiO, SO, O 2 , CS, SiS, N 2 etc., over a span of a few (between 1 and 4) years after the explosion. The epoch of synthesis and the abundances vary with the zones in the ejecta. The findings are congruent with the observations of CO (Danziger et al. 1988;Catchpole et al. 1988) and SiO (Roche et al. 1991;Kotak et al. 2005 molecules in local supernovae, mainly SN1987A. Figure 18 shows the time evolution of CO and SiO masses in a 15 M progenitor SN ejecta compared to the observations. CO molecules are essentially synthesized in the C/O rich zone 4, followed by the central region on the O-core (zone 2). CO formation follows very efficient chemical routes. Hence, the abundance increases as early as 150 days and becomes constant in mass by 700 days post-explosion. Due to the strong chemical bonds, once they are formed in the ejecta, the CO molecules are found to withstand the harsh environment conditions and possible destruction mechanisms. The chemical pathways important for CO formation are mainly: C + O 2 → CO + O, C + O → CO + hν etc. The final masses estimated by the model are in accord with the findings of Kamenetzky et al. (2013) and also it matches well at ∼300 days (4 × 10 −3 M ) with the IR observations of SN1987A. There is no evidence of depletion of CO molecules with time. Hence it can be safely claimed that CO does not act as a precursor to the carbon dust (Sarangi and Cherchneff 2013;Biscaro and Cherchneff 2014).

Condensation Phase
The dust precursors act as seeds in the condensation scheme, where they grow by coagulation and coalescence to form large grains. The gas phase reactions and grain condensation occur in tandem. The model is designed to maintain the balance between gas phase  (Liu et al. 1992), thermal assumption: circles (Liu and Dalgarno 1995), ALMA data: triangle (Kamenetzky et al. 2013). The ALMA data is not true to the time axis, but rather indicates the final mass at late times as compared to the model results. Right panel: Evolution of SiO mass (M ) with post-explosion time in days as a function of different zones in the ejecta. The black dotted black line corresponds to the sum over the SiO mass from all the zones. The solid square, circle, and triangle indicates the mass of SiO reported from observations in IR from several supernovae as specified in the legend box (Sarangi and Cherchneff 2013) processes, the solid phase interactions and the smooth transition between the two interactive phases at all times. The formalism has been developed by Jacobson (2005), where the variation of the grain number density of volume v is described by the integro-differential coagulation equation, where β is the rate coefficient of coagulation between particles. The rate coefficient is controlled by physical processes such as Brownian diffusion, convective Brownian motion enhancement, gravitational collection, turbulent inertial motion and Van der Waals forces. In free molecular regimes (where mean free-path of a particle in the gas is much larger compared to the mean radius of a grain), Brownian diffusion prevails as the most important phenomena along with Van der Waals interactions. Brownian diffusion accounts for the scattering, collision and coalescence of the grains through Brownian motion. On the other hand, the Van der Waals forces develop weak, local charge fluctuations that enhance the rate of coagulation for particles with sizes in the molecular range. Solving the appropriate equations, the concentrations and the grain size distribution function for individual dust types are estimated. This formalism has been used to model dust in SN ejecta and winds of AGB stars by Sarangi and Cherchneff (2015), Gobrecht et al. (2016). The dust components comprise of silicates, alumina, carbon, iron sulphide, silicon carbide and pure metallic clusters of magnesium, iron and silicon. The analysis traces the evolution mechanism of all these dust components individually. The formation of dust grains in the ejecta is found to be highly sensitive to the gas densities and the concentration of the specific clusters in consideration. Hence, the dust formation scenario varies to a great extent from one zone to the other even for a single dust type. The estimated mass of each dust component and the size distributions of the grains were calculated as a function of post-explosion time. There is no specific threshold in size to differentiate between the large From the observational point of view, the dust formation in the ejecta is characterised by an excess in IR spectral energy distribution and blue shift of optical emission lines. From the modeling perspective, a fiducial size of 10 Å has been chosen as the minimum size of grains in the ejecta, which can be categorised as dust. The grains smaller in size are considered as uncondensed molecular clusters, which remain in gas phase. Figure 19 illustrates the evolution of dust masses of all the dust components from day 100 to 2000 after explosion. The masses of dust components in the figure are presented as the sum over all the ejecta zones. Formation of dust in the ejecta dawns after 200 days post-explosion time with the formation of iron sulphide grains. Later on in the period of 300-400 days silicates (forsterite) starts to form mainly in the inner zones of the O-core. The mass of dust surges up at ∼550 days due to synthesis of silicates in the central zone 2. The process of gradual increase of dust mass continue with the advent of alumina at day 800 and lately by amorphous carbon at 1100 days. This way the dust mass increases from 10 −5 M at 400 days to 0.035 M at the end of four years. Table 5 presents all the evolution pattern of each dust type over a span of a 2000 days after explosion. Amorphous carbon, alumina and forsterite are the most abun- Fig. 20 The final size distribution f (a) of dust grains at 2000 days post-explosion is presented for all the dust components (Sarangi and Cherchneff 2015). The distribution function is expressed in terms of total number of particles per unit radius (N tot (a). Å −1 ) dant dust types in the ejecta. The range of dust masses closely agrees with the estimated masses of dust derived from observational data in the IR. At times between 200-600 days, the mid-IR observations predict the presence of dust in the ejecta with masses typically in the order of 10 −5 to 10 −3 M (Kotak et al. , 2005Elmhamdi et al. 2003b;Wooden et al. 1993;Inserra et al. 2011;Sugerman et al. 2006;Gallagher et al. 2012). These findings are consistent with the model results of the total dust mass over time (Fig. 19). Submm observations however report the presence of much larger masses of dust in the remnants. In Cas A, ∼0.08 M of dust is inferred from Herschel data (Barlow et al. 2010). In case of SN1987A, Matsuura et al. 2011 has derived masses of dust in the range of 0.4-0.7 M to be present in the ejecta after 28 years from explosion. Recent ALMA observation concluded the final dust mass to be about 0.5 M (Matsuura et al. 2015). Also for the filaments of Crab nebula an estimated mass of 0.1-0.2 M is required to fit the Herschel data as explained by Gomez et al. 2012. The final mass of dust estimated in this study ends up in 0.035 M , and we can suppose the dust mass will remain consistent till the advent of reverse shock in the remnant. These trends derived from the model match well with the results of submm observations (Sarangi and Cherchneff 2015;Biscaro and Cherchneff 2014). However the final masses reported from these late time observational data are 2-10 times larger than our assessment. Importantly though, all the masses reported from observations are derived by fitting the spectral energy distributions at IR and submm wavelengths where some composition of dust has been assumed. The dust masses from observations hence can be sensitive to the choice of dust components, their compositions and size at a λ. Finally, there are several options for the choice of optical constants related to each dust type, and the results vary significantly from one to the other. In most of the cases, the size distributions were assumed to either follow a power law profile like the standard Mathis-Rumpl-Nordsieck (MRN) distribution (Mathis et al. 1977) or the grains were chosen to be of identical sizes. All these factors add to some uncertainty in the estimation of dust masses fitting the observed spectra.
The size distribution function defines the number of grains per unit of particle radius, assuming a spherical geometry. The size distribution function f (a), for all the important dust species at day 2000 is presented in Fig. 20. The main three dust components carbon, alumina and forsterite are presented in solid colours in the figure. The distribution function at a given time depends on (a) the number of monomers present in the gas phase, (b) the rate of coagulation at the given temperature and (c) the concentration of grains from the previous time step. The prime dust components being more abundant, also grow to form larger grains in the order of 200-500 Å. The size distribution function of one dust species  (Sarangi and Cherchneff 2013) differs in behaviour from the other and hence cannot be defined by a single mathematical expression. Depending on the formation zones in the ejecta, the distribution functions can have a single peak or sometimes are found to have multiple peaks. The findings suggest that the dust masses reside mainly in the larger grains even though the population of large grains may be less than that of relatively smaller grains in the ejecta. Moreover, the condensation phase is an efficient process in the ejecta and once the large grains are formed, they are unlikely to segregate to smaller size fragments.
The size distributions of dust produced by SNe do not follow a MRN size variation both for single dust type and for the total dust distribution. The size distributions derived from the study have less small grains with radius a < 30 Å and are skewed towards large grains in the size range 40-1000 Å, compared with the MRN distribution. This result clearly indicates that the use of a MRN size distribution for modeling the fluxes from a SN ejecta is inappropriate. Further, using MRN for dust sputtering calculations in the SNe remnants may lead to erroneous results on the derived dust masses and the dust mass that survives shock processing in the remnants. Ideally if there at all exists a distribution function for supernovae dust grains, that should be defined by the dust size patterns of a single dust component in a given zone, and not over a sum of all the species. The analysis clearly indicates none of the dust components follow a power law pattern of dust sizes at any epoch.
The ions of the nobel gases, namely helium, neon and argon present in the ejecta, play an important role in moulding the chemistry (Lepp et al. 1990). The γ -rays in the ejecta produced by radioactivity undergo Compton scattering and degrade to X-rays. The energetic Compton electrons further lose energy by excitation and ionisation of the gas. By day 100, most of the atoms in the ejecta gets ionised to their corresponding cation through collision with the fast Compton electrons. The ground state energy level of their inert parents being very low, the cations of He, Ne and Ar actively participate in charge exchange processes with newly synthesised molecules, leading to their destruction (e.g. AB + He + → A + B + + He). The impact is clearly evident in the outermost zones of the He-core, where the detrimental effects of He + delay the appearance of stable molecules and clusters in the ejecta till ∼ day 1000. Figure 21 (left-panel) presents the correlation between the epoch of dust formation to the abundances of He + , Ne + or Ar + in the ejecta zones. With the passage of time, effect of radioactivity decreases as most of the 56 Co already gets converted to stable 56 Fe. The energy from the γ -rays and X-rays get deposited in the ejecta and a small fraction escapes outside the core. Thereafter, the abundance of Compton electrons also falls rapidly and all the cations eventually recombine back to their stable inert parents. The molecule and dust synthesis gain in efficiency when the ions of He, Ne and Ar are no longer present in the ejecta.
The amount of radioactive elements produced by the explosion also directly influences the chemical yields. The radioactive 56 Ni mass produced by a SN can be derived from the variation of the optical light curves and H α luminosities in the nebula phase (Elmhamdi et al. 2003a). Sarangi and Cherchneff (2013) provides a list of local SNe with their estimated progenitor masses as well as the masses of 56 Ni produced by the explosion. Observational findings support the trend derived by Hamuy (2003) that more massive supernovae produce more energetic explosions, and supernovae with greater energies produce larger 56 Ni masses. The energy deposited by the radioactive 56 Ni in the ejecta is linearly proportional to the its mass. Therefore, the rate coefficients for the destructions by energetic Compton electrons are directly connected to the mass of 56 Ni. Physically, a smaller 56 Ni mass in the ejecta reduces the number of Compton electrons resulting from the degradation of a lower amount of γ -rays. Therefore, fewer ions such as He + , Ne + or Ar + are produced, thereby enhancing the rate of recombination. The destruction rates for molecules and clusters (e.g., SiO, C 2 , AlO) are much less severe for a low 56 Ni mass case, because of the lower Ni + and He + concentration. So the chances of survival for the molecules and clusters increase, thereby also leading to early condensation of dust grains. Figure 21 (right-panel) shows the evolution of silicates, alumina and carbon dust masses corresponding to 56 Ni = 0.075 M and 56 Ni = 0.01 M cases and compares the respective yields.

Mixing and Clumpiness
A homogeneous, 1D model of the ejecta is capable to indicate trends, which are closely similar to the actual event. Nevertheless, a real supernova ejecta is neither homogeneous nor spherically symmetric. Observational findings suggest the possibility of macroscopic mixing of 56 Ni in the outer layers of the H-He core (McCray 1993;Hanuschik et al. 1988). Detection of hard X-rays and γ -ray lines in SN1987A, immediately after the explosion, indicates a presence of a high-velocity component of radioactive 56 Ni (Dotani et al. 1987;Syunyaev et al. 1987;Matz et al. 1988). The mass of the high velocity (>3000 km/s) 56 Ni clump is estimated to be about 10 −3 M , which is between 4-17 % of the total 56 Ni produced at the explosion (Haas et al. 1990;Utrobin et al. 1995), as explained in Fig. 22 (Wongwathanarat et al. 2015). Therefore, the mixing of 56 Ni in the earliest days post-explosion is intense yet not complete. As a result, a two-phase SN-ejecta structure is formed. Low-density blobs consisting of Fe-bubbles are created by the energy input from radioactive 56 Co and 56 Ni, which is surrounded by compressed high density metal-rich ejecta (Blondin et al. 2001).
The multi-D explosion models study the presence of Rayleigh-Taylor unstable filaments in the non-homogeneous ejecta (Kifonidis et al. 2003(Kifonidis et al. , 2006Joggerst et al. 2010;Hammer et al. 2010a;Wongwathanarat et al. 2015) of Type II supernovae. Figure 22 presents the isosurfaces of various metal-rich blobs after a few minutes post-explosion (Wongwathanarat et al. 2015). Moreover, the deviation from spherical symmetry is also evident from the structure of the ejecta.
The explosion energy is imparted in the ejecta in form of kinetic energy, which triggers an outward flow with velocities of around 2000-2500 km/s. If all the materials in the Hecore travel at a constant velocity, the density of the outflowing matter should decrease as Fig. 22 The isosurfaces of constant mass-fraction of Ni plumes and O/Mg zone at 9000 s post outburst is shown in the Figure adopted from Wongwathanarat et al. (2015). The 56 Ni blob is found to travel at velocities between 3800 and 4000 km/s t −2 law (Woosley 1988). However, the active radioactive processes (e.g. 56 Ni → 56 Co → 56 Fe) deposit the extra energy in the innermost Ni-Co core. Further, the neutrino-driven convection in the ejecta creates inhomogeneity between the layers (Joggerst et al. 2010). This creates a velocity shear across the zones of the ejecta, and the inner zone with extra energy tends to expand faster (Woosley 1988). The velocity shear induces the homologous expansion in the ejecta. Therefore, the fast moving low-density Ni-layers penetrate in the high density overlying zones of the ejecta. This density inversion occurs within a day of the explosion, thereby creating the Rayleigh-Taylor instabilities. The mixed Ni-bubbles further deposit energy in the O-core, which is in turn accelerated through the outer He/N/C zones. These instabilities lead to the fragmentation of the He-core, but the chemical stratification persists over time. Importantly though, this velocity shear results in clumping of the ejecta very early after the explosion.
Radiative transfer models trying to fit the IR contribution from the ejecta dust have repeatedly emphasized on the clumpy nature of the medium (Ercolano et al. 2007;Gallagher et al. 2012). Further, in order to account for the isotopic anomalies in the presolar grains in meteorites of SN origin, a small percentage of element mixing in microscopic scales is necessary. The scenario of mixing and inhomogeneity in the ejecta can only be modelled using hydrodynamic simulations, which is a very computationally rigorous process when all the important metals in the ejecta are being traced. Macroscopic mixing between the ejecta zones following the explosion is evident from observations and simulations. However, the possibility or the mechanism of microscopic diffusion of metals between the zone boundaries is not yet well understood.
A quasi-clumpy model based on SN1987A was formulated by Sarangi and Cherchneff (2015) using the volume filling factors for individual zones derived by Jerkstrand et al. (2011) through the study of the ultraviolet, optical, and near-IR emission lines. Dust formation is found to be faster and more efficient inside a high-density clump. However, depending on the optical depth of the clumpy region and rapid cooling of the newly formed grains, the dust can remain hidden and escape mid-IR detection. Microscopic mixing between inner (Si/S) and outer (H/He/C) zones would lead to formation of SiC grains (important constituent in presolar grains) in larger abundance and H-bearing species such as recently detected HCO + ).

Summary of Current State-of-the-Art
Dust formation scenario in SNe is a complex interplay between several physical and chemical processes, which may vary depending on minimal alteration of conditions. Therefore, formulation of an all-inclusive single model of dust in core-collapse SNe most of the time over-simplifies the actual scenario. On the contrary, a more acceptable approach is to understand the general trends and then apply the knowledge to study any specific SN with its given constraints. Here we shall summarise the aspects related to the parameter space, which characterises the physics and chemistry of dust formation in a typical SN-ejecta.
(a) Core-collapse supernovae, mainly Type II-P, are the most efficient dust producers among all SN-types. On the other hand, the amount of dust formed in Type Ia SN is negligible. (b) Dust formation in the ejecta is characterised by the chemistry prevailing in a non-steady state non-thermodynamic equilibrium environment. (c) In a typical SN ejecta, dust formation commences within a year from the explosion and the dust mass continues to increase for a couple of years. The span depends essentially on the density and clumpiness of the medium. (d) Dust synthesis is controlled by simultaneous processes of nucleation and condensation.
Nucleation of stable dust clusters act as the bottleneck to the entire chain of chemical network. These molecular clusters resulting from nucleation determines the chemical composition of the dust as well as the earliest epochs of dust formation. (e) Theory predicts of about 30-40% of the total mass in the ejecta. CO, SiO, O 2 , CS, AlO, SO etc. are important molecular species. Dust, on the other hand, is about 3-10% of the total ejecta mass. (f) The dust budget in the SN ejecta is comprises of silicates, amorphous carbon, alumina and pure metallic clusters of silicon, iron and magnesium. There is also the possibility of mixed species as the coexisting grains in a certain zone might coagulate heterogeneously. (g) Dust masses and grain sizes are inversely related to the amount of radioactive 56 Ni produced at the time of explosion. On the other hand, they are directly proportional to the degree of clumpiness. (h) Relative abundances of C-rich and O-rich dust depend on the mass of the progenitor.
The larger the progenitor, the larger is the O-core, and therefore O-rich dust dominates over other components. (i) The size distributions of dust grains in SN ejecta are functions of dust type and epochs of formation. They cannot be generalised as power laws of any exponent, contrary to the general assumption. (j) A clumpy ejecta of supernova forms larger grains at high densities. These grains are most likely to survive the passage of the reverse shock in the later phase of the evolution. (k) Theory predicts that a typical type II-P supernova can produce dust masses between 0.02 to 0.2 M , while measured dust masses range from 10 −6 to 1 M . (l) A 1-D model can precisely determine the trends and upper or lower limits for the scenario. However, in order to incorporate the impact of mixing and inhomogeneous deposition of radioactive energy in the ejecta, multi-D dynamical modeling of the ejecta is necessary.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.