There is currently a severe discrepancy between theoretical models of dust formation in core-collapse supernovae (CCSNe), which predict $\gtrsim 0.01$ M$_\odot$ of ejecta dust forming within $\sim 1000$ days, and observations at these epochs, which infer much lower masses. We demonstrate that, in the optically thin case, these low dust masses are robust despite significant observational and model uncertainties. For a sample of 11 well-observed CCSNe, no plausible model reaches carbon dust masses above $10^{-4}$ M$_\odot$, or silicate masses above $\sim 10^{-3}$ M$_\odot$. Optically thick models can accommodate larger dust masses, but the dust must be clumped and have a low ($<0.1$) covering fraction to avoid conflict with data at optical wavelengths. These values are insufficient to reproduce the observed infrared fluxes, and the required covering fraction varies not only between SNe but between epochs for the same object. The difficulty in reconciling large dust masses with early-time observations of CCSNe, combined with well-established detections of comparably large dust masses in supernova remnants, suggests that a mechanism for late-time dust formation is necessary.

The luminous Type IIn SN 2010jl shows strong signs of interaction between the SN ejecta and dense circumstellar material. Dust may be present in the unshocked ejecta, the cool, dense shell between the shocks in the interaction region, or in the circumstellar medium. We present and model new optical and infrared photometry and spectroscopy of SN 2010jl from 82 to 1367 days since explosion. We evaluate the photometric and spectroscopic evolution using the radiative transfer codes MOCASSIN and DAMOCLES, respectively. We propose an interaction scenario and investigate the resulting dust formation scenarios and dust masses. We find that SN 2010jl has been continuously forming dust based on the evolution of its infrared emission and optical spectra. There is evidence for pre-existing dust in the circumstellar medium as well as new dust formation in the cool, dense shell and/or ejecta. We estimate that 0.005--0.01M$_{\odot}$ of predominantly carbon dust grains has formed in SN 2010jl by ~1400d post-outburst.

Observations of the ejecta of core-collapse supernovae have shown that dust grains form in over-dense gas clumps in the expanding ejecta. The clumps are later subject to the passage of the reverse shock and a significant amount of the newly formed dust material can be destroyed due to the high temperatures and high velocities in the post-shock gas. To determine dust survival rates, we have performed a set of hydrodynamic simulations using the grid-based code AstroBEAR in order to model a shock wave interacting with a clump of gas and dust. Afterwards, dust motions and dust destruction rates are computed using our newly developed external, post-processing code Paperboats, which includes gas and plasma drag, grain charging, kinematic and thermal sputtering as well as grain-grain collisions. We have determined dust survival rates for the oxygen-rich supernova remnant Cassiopeia A as a function of initial grain sizes, dust materials and clump gas densities.

The reverse shock in the ejecta of core-collapse supernovae is potentially able to destroy newly formed dust material. In order to determine dust survival rates, we have performed a set of hydrodynamic simulations using the grid-based code AstroBEAR in order to model a shock wave interacting with clumpy supernova ejecta. Dust motions and destruction rates were computed using our newly developed external, post-processing code Paperboats, which includes gas drag, grain charging, sputtering and grain-grain collisions. We have determined dust destruction rates for the oxygen-rich supernova remnant Cassiopeia A as a function of initial grain sizes and clump gas density. We found that up to 30 % of the carbon dust mass is able to survive the passage of the reverse shock if the initial grain size distribution is narrow with radii around ~10 - 50 nm for high gas densities, or with radii around ~0.5 - 1.5 ${\mu}$m for low and medium gas densities. Silicate grains with initial radii around 10 - 30 nm show survival rates of up to 40 % for medium and high density contrasts, while silicate material with micron sized distributions is mostly destroyed. For both materials, the surviving dust mass is rearranged into a new size distribution that can be approximated by two components: a power-law distribution of small grains and a log-normal distribution of grains having the same size range as the initial distribution. Our results show that grain-grain collisions and sputtering are synergistic and that grain-grain collisions can play a crucial role in determining the surviving dust budget in supernova remnants.

We have modelled the near-infrared to radio images of the Crab Nebula with a Bayesian SED model to simultaneously fit its synchrotron, interstellar and supernova dust emission. We infer an interstellar dust extinction map with an average $A_{\text{V}}$=1.08$\pm$0.38 mag, consistent with a small contribution (<22%) to the Crab's overall infrared emission. The Crab's supernova dust mass is estimated to be between 0.032 and 0.049 M$_{\odot}$ (for amorphous carbon grains) with an average dust temperature $T_{\text{dust}}$=41$\pm$3K, corresponding to a dust condensation efficiency of 8-12%. This revised dust mass is up to an order of magnitude lower than some previous estimates, which can be attributed to our different interstellar dust corrections, lower SPIRE flux densities, and higher dust temperature than were used in previous studies. The dust within the Crab is predominantly found in dense filaments south of the pulsar, with an average V-band dust extinction of $A_{\text{V}}$=0.20-0.39 mag, consistent with recent optical dust extinction studies. The modelled synchrotron power-law spectrum is consistent with a radio spectral index $\alpha_{\text{radio}}$=0.297$\pm$0.009 and an infrared spectral index $\alpha_{\text{IR}}$=0.429$\pm$0.021. We have identified a millimetre excess emission in the Crab's central regions, and argue that it most likely results from two distinct populations of synchrotron emitting particles. We conclude that the Crab's efficient dust condensation (8-12%) provides further evidence for a scenario where supernovae can provide substantial contributions to the interstellar dust budgets in galaxies.

In order to understand the contribution of core-collapse supernovae to the dust budget of the early universe, it is important to understand not only the mass of dust that can form in core-collapse supernovae but also the location and rate of dust formation. SN 2005ip is of particular interest since dust has been inferred to have formed in both the ejecta and the post-shock region behind the radiative reverse shock. We have collated eight optical archival spectra that span the lifetime of SN 2005ip and we additionally present a new X-shooter optical-near-IR spectrum of SN 2005ip at 4075d post-discovery. Using the Monte Carlo line transfer code DAMOCLES, we have modelled the blueshifted broad and intermediate width H$\alpha$, H$\beta$ and He I lines from 48d to 4075d post-discovery using an ejecta dust model. We find that dust in the ejecta can account for the asymmetries observed in the broad and intermediate width H$\alpha$, H$\beta$ and He I line profiles at all epochs and that it is not necessary to invoke post-shock dust formation to explain the blueshifting observed in the intermediate width post-shock lines. Using a Bayesian approach, we have determined the evolution of the ejecta dust mass in SN 2005ip over 10 years presuming an ejecta dust model, with an increasing dust mass from ~10$^{-8}$ M$_{\odot}$ at 48d to a current dust mass of $\sim$0.1 M$_{\odot}$.

Optical and near-IR (NIR) line profiles of many ageing core-collapse supernovae (CCSNe) exhibit an apparently asymmetric bluewards shift often attributed to greater extinction by internal dust of redshifted radiation emitted from the receding regions of the SN ejecta. The DAMOCLES Monte Carlo line radiative transfer code models the extent and shape of these dust-affected line profiles to determine the dust mass that has condensed, in addition to other properties of the dusty ejecta. I present here the application of an affine invariant Markov Chain Monte Carlo (MCMC) ensemble sampler (emcee) to the DAMOCLES code in order to investigate the multi-dimensional parameter space rigorously and characterise the posterior probability distribution. A likelihood function is formulated that handles both Monte Carlo and observational uncertainties. This Bayesian approach is applied to four simulated line profiles in order to test the method and investigate its efficacy. The majority of parameters can be tightly constrained using this method, and a strong (predictable) dependence between the grain size and the dust mass is quantified. The new approach is also applied to the H$\alpha$ line and [O I]6300,6363 \AA\ doublet of SN 1987A at 714d post-outburst, re-examining a previous 5-dimensional smooth model and also investigating a new, more complex, 10-dimensional model that treats both features simultaneously. The dust mass, dust grain size and a range of other parameters can be well constrained using this technique, representing a significant improvement over the previous manual approach.

The equatorial ring of Supernova (SN) 1987A has been exposed to forward shocks from the SN blast wave, and it has been suggested that these forward shocks have been causing on-going destruction of dust in the ring. We obtained SOFIA FORCAST 11.1, 19.7 and 31.5 micron photometry of SN\,1987A in 2016. Compared with Spitzer measurements 10 years earlier, the 31.5 micron flux has significantly increased. The excess at 31.5 micron appears to be related to the Herschel 70 micron excess, which was detected 5 years earlier. The dust mass needed to account for the the 31.5--70 micron excess is 3--7x10^-4 Msun, more than ten times larger than the ring dust mass (1x10^-5 Msun) estimate from the data 10-years earlier. We argue that dust grains are re-formed or grown in the post-shock regions in the ring after forward shocks have destroyed pre-existing dust grains in the ring and released refractory elements into gas. In the post-shock region, atoms can stick to surviving dust grains, and the dust mass may have increased (grain growth), or dust grains might have condensed directly from the gas. An alternative possibility is that the outer part of the expanding ejecta dust might have been heated by X-ray emission from the circumstellar ring. The future development of this excess could reveal whether grains are reformed in the post-shocked region of the ring or eject dust is heated by X-ray.

We present Monte Carlo line transfer models that investigate the effects of dust on the very late time emission line spectra of the core collapse supernovae SN 1980K and SN 1993J and the young core collapse supernova remnant Cassiopeia A. Their blue-shifted emission peaks, resulting from the removal by dust of redshifted photons emitted from the far sides of the remnants, and the presence of extended red emission wings are used to constrain dust compositions and radii and to determine the masses of dust in the remnants. We estimate dust masses of between 0.08 - 0.15 M$_\odot$ for SN 1993J at year 16, 0.12 - 0.30 M$_\odot$ for SN 1980K at year 30 and $\sim$1.1 M$_\odot$ for Cas A at year $\sim$330. Our models for the strong oxygen forbidden lines of Cas A require the overall modelled profiles to be shifted to the red by between 700 - 1000 km s$^{-1}$, consistent with previous estimates for the shift of the dynamical centroid of this remnant.

The late time optical and near-IR line profiles of many core-collapse supernovae exhibit a red-blue asymmetry as a result of greater extinction by internal dust of radiation emitted from the receding parts of the supernova ejecta. We present here a new code, DAMOCLES, that models the effects of dust on the line profiles of core-collapse supernovae in order to determine the masses of newly formed dust. As noted by Lucy et al. (1989), the presence of an extended red scattering wing in late-time line profiles can also indicate dust formation. We find that dust-affected line profiles need not necessarily be flux-biased towards to the blue, although the profile peak will always be blue-shifted. We have collated optical spectra of SN 1987A from a variety of archival sources and have modelled the evolution of the H$\alpha$ line from days 714 to 3604, as well as that of the [OI] 6300,6363A doublet between days 714 and 1478. A variety of evidence points to the presence of clumping and we find that our clumped dust models require significantly higher dust masses than smoothly distributed dust models. Our line profile fits imply day 714 dust masses of <3 $\times$ 10$^{-3}$ M$_{\odot}$ for all grain types apart from very high albedo pure magnesium silicates, for which up to 0.07M$_{\odot}$ can be accommodated. Large grain radii (>0.6$\mu$m) are generally required to fit the line profiles even at the earlier epochs. We find that a large dust mass (>0.1M$_{\odot}$) had formed by day 3604 and infer that the majority of the present dust mass must have formed after this epoch. Our findings agree with recent estimates from SED fits for the dust mass evolution of SN 1987A and support the inference that the majority of SN 1987A's dust formed many years after the initial explosion.

We investigate the effect of ionising radiation from Massive Young Stellar Objects impinging on their emerging spectral energy distribution. By means of detailed radiative transfer calculations including both the gaseous and dust phase of their surrounding protoplanetary discs we highlight the importance of modelling both phases simultaneously when interpreting observations from such objects. In particular we find that models that only include dust may lead to incorrect conclusions about the inner disc evolution. Furthermore the omission of gas from models overproduces far-infrared and sub-millimiter fluxes with the result that derived dust masses may be underestimated by a factor of two in some cases. Finally free-free emission from the ionised component of gaseous discs causes the slope of the dust emission in the sub-mm and mm regime to appear flatter, resulting in incorrectly modelling the dust properties, with consequences on the derived disc masses, power law index of the surface density profile and other disc properties.