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.

Recent analyses of ${\it Herschel}$ observations suggest that in nearby disc galaxies the dust mass opacity at $500 \, {\rm \mu m}$, $\kappa_{500}$, decreases with increasing gas surface density, $\Sigma_{\rm ISM}$ (Clark et al. 2019). This apparent anti-correlation between $\kappa_{500}$ and $\Sigma_{\rm ISM}$ is opposite to the behaviour expected from theoretical dust evolution models; in such models, dust in denser, cooler regions (i.e. regions of increased $\Sigma_{\rm ISM}$) tends to grow and therefore to have increased $\kappa_{500}$. We show, using a toy model, that the presence of a range of dust temperatures along the line of sight can lead to spuriously low estimated values of $\kappa_{500}$. If in regions of higher $\Sigma_{\rm ISM}$ the range of dust temperatures extends to lower values (as seems likely), the magnitude of this effect may be sufficient to explain the apparent anti-correlation between $\kappa_{500}$ and $\Sigma_{\rm ISM}$. Therefore there may not be any need for spatial variation in the intrinsic dust properties that run counter to theoretical expectations.

We calculate dust spectral energy distributions (SEDs) for a range of grain sizes and compositions, using physical properties appropriate for five pulsar wind nebulae (PWNe) from which dust emission associated with the ejecta has been detected. By fitting the observed dust SED with our models, with the number of grains of different sizes as the free parameters, we are able to determine the grain size distribution and total dust mass in each PWN. We find that all five PWNe require large ($\ge 0.1 \, {\rm \mu m}$) grains to make up the majority of the dust mass, with strong evidence for the presence of micron-sized or larger grains. Only two PWNe contain non-negligible quantities of small ($<0.01 \, {\rm \mu m}$) grains. The size distributions are generally well-represented by broken power laws, although our uncertainties are too large to rule out alternative shapes. We find a total dust mass of $0.02-0.28 \, {\rm M}_\odot$ for the Crab Nebula, depending on the composition and distance from the synchrotron source, in agreement with recent estimates. For three objects in our sample, the PWN synchrotron luminosity is insufficient to power the observed dust emission, and additional collisional heating is required, either from warm, dense gas as found in the Crab Nebula, or higher temperature shocked material. For G$54.1$+$0.3$, the dust is heated by nearby OB stars rather than the PWN. Inferred dust masses vary significantly depending on the details of the assumed heating mechanism, but in all cases large mass fractions of micron-sized grains are required.

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 investigate differences in the molecular abundances between magnetically super- and sub-critical prestellar cores, performing three-dimensional non-ideal magnetohydrodynamical (MHD) simulations with varying densities and magnetic field strengths, and post-processing the results with a time-dependent gas-grain chemical code. Most molecular species show significantly more central depletion in subcritical models, due to the longer duration of collapse. However, the directly observable quantities - the molecule to hydrogen column density ratios - are generally too similar for observational data to discriminate between models. The profiles of N$_2$H$^+$ and HCO$^+$ show qualitative differences between supercritical and subcritical models on scales of $0.01 \, {\rm pc}$, which may allow the two cases to be distinguished. However, this requires knowledge of the hydrogen column density, which is not directly measureable, and predicted line intensity profiles from radiative transfer modelling are similar for these molecules. Other commonly observed species, such as HCN and CH$_3$OH, have line intensity profiles which differ more strongly between models, and so are more promising as tracers of the mechanism of cloud collapse.

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.

We model the thermal dust emission from dust grains heated by synchrotron radiation and by particle collisions, under conditions appropriate for four different shocked and unshocked gas components of the Cassiopeia A (Cas A) supernova remnant (SNR). By fitting the resulting spectral energy distributions (SEDs) to the observed SNR dust fluxes, we determine the required mass of dust in each component. We find the observed SED can be reproduced by $\sim 0.6 \, {\rm M_\odot}$ of silicate grains, the majority of which is in the unshocked ejecta and heated by the synchrotron radiation field. Warmer dust, located in the X-ray emitting reverse shock and blastwave regions, contribute to the shorter wavelength infrared emission but make only a small fraction of the total dust mass. Carbon grains can at most make up $\sim 25 \%$ of the total dust mass. Combined with estimates for the gas masses, we obtain dust-to-gas mass ratios for each component, which suggest that the condensation efficiency in the ejecta is high, and that dust in the shocked ejecta clumps is well protected from destruction by sputtering in the reverse shock.

We develop analytic approximations to the density evolution of prestellar cores, based on the results of hydrodynamical simulations. We use these approximations as input for a time-dependent gas-grain chemical code to investigate the effects of differing modes of collapse on the molecular abundances in the core. We confirm that our method can provide reasonable agreement with an exact numerical solution of both the hydrodynamics and chemistry while being significantly less computationally expensive, allowing a large grid of models varying multiple input parameters to be run. We present results using this method to illustrate how the chemistry is affected not only by the collapse model adopted, but also by the large number of unknown physical and chemical parameters. Models which are initially gravitationally unstable predict similar abundances despite differing densities and collapse timescales, while ambipolar diffusion produces more extended inner depleted regions which are not seen in observations of prestellar cores. Molecular observations are capable of discriminating between modes of collapse despite the unknown values of various input parameters. We also investigate the evolution of the ambipolar diffusion timescale for a range of collapse modes, metallicities and cosmic ray ionization rates, finding that it remains comparable to or larger than the collapse timescale during the initial stages for all models we consider, but becomes smaller at later evolutionary stages. This confirms that ambipolar diffusion is an important process for diffuse gas, but becomes less significant as cores collapse to higher densities.

We model the molecular emission from cometary knots in planetary nebulae (PNe) using a combination of photoionization and photodissociation region (PDR) codes, for a range of central star properties and gas densities. Without the inclusion of ionizing extreme ultraviolet (EUV) radiation, our models require central star temperatures $T_*$ to be near the upper limit of the range investigated in order to match observed H$_2$ and OH$^+$ surface brightnesses consistent with observations - with the addition of EUV flux, our models reproduce observed OH$^+$ surface brightnesses for $T_* \ge 100 \, {\rm kK}$. For $T_* < 80 \, {\rm kK}$, the predicted OH$^+$ surface brightness is much lower, consistent with the non-detection of this molecule in PNe with such central star temperatures. Our predicted level of H$_2$ emission is somewhat weaker than commonly observed in PNe, which may be resolved by the inclusion of shock heating or fluorescence due to UV photons. Some of our models also predict ArH$^+$ and HeH$^+$ rotational line emission above detection thresholds, despite neither molecule having been detected in PNe, although the inclusion of photodissociation by EUV photons, which is neglected by our models, would be expected to reduce their detectability.

We have performed combined photoionization and photodissociation region (PDR) modelling of a Crab Nebula filament subjected to the synchrotron radiation from the central pulsar wind nebula, and to a high flux of charged particles; a greatly enhanced cosmic ray ionization rate over the standard interstellar value, $\zeta_0$, is required to account for the lack of detected [C I] emission in published Herschel SPIRE FTS observations of the Crab Nebula. The observed line surface brightness ratios of the OH$^+$ and ArH$^+$ transitions seen in the SPIRE FTS frequency range can only be explained with both a high cosmic ray ionization rate and a reduced ArH$^+$ dissociative recombination rate compared to that used by previous authors, although consistent with experimental upper limits. We find that the ArH$^+$/OH$^+$ line strengths and the observed H$_2$ vibration-rotation emission can be reproduced by model filaments with $n_{\rm{H}} = 2 \times 10^4$ cm$^{-3}$, $\zeta = 10^7 \zeta_0$ and visual extinctions within the range found for dusty globules in the Crab Nebula, although far-infrared emission from [O I] and [C II] is higher than the observational constraints. Models with $n_{\rm{H}} = 1900$ cm$^{-3}$ underpredict the H$_2$ surface brightness, but agree with the ArH$^+$ and OH$^+$ surface brightnesses and predict [O I] and [C II] line ratios consistent with observations. These models predict HeH$^+$ rotational emission above detection thresholds, but consideration of the formation timescale suggests that the abundance of this molecule in the Crab Nebula should be lower than the equilibrium values obtained in our analysis.

We present a publicly available, open source version of the time dependent gas-grain chemical code UCLCHEM. UCLCHEM propagates the abundances of chemical species through a large network of chemical reactions in a variety of physical conditions. The model is described in detail along with its applications. As an example of possible uses, UCLCHEM is used to explore the effect of protostellar collapse on commonly observed molecules and to study the behaviour of molecules in C-type shocks. We find the collapse of a simple Bonnor-Ebert sphere successfully reproduces most of the behaviour of CO,CS and NH$_3$ from cores observed by \citet{Tafalla2004} but cannot predict the behaviour of N$_2$H$^+$. In the C-shock application, we find that molecules can be categorized so that they can be useful observational tracers of shocks and their physical properties. Whilst many molecules are enhanced in shocked gas, we identify two groups of molecules in particular. A small number of molecules are enhanced by the sputtering of the ices as the shock propagates and then remain high in abundance throughout the shock. A second, larger set are also enhanced by sputtering but then are destroyed as the gas temperature rises. Through these applications the general applicability of UCLCHEM is demonstrated.