Dust associated with various stellar sources in galaxies at all cosmic epochs remains a controversial topic, particularly whether supernovae (SNe) play an important role in dust production. We report evidence of dust formation in the cold, dense shell behind the ejecta-circumstellar medium (CSM) interaction in the Type Ia-CSM SN 2018evt three years after the explosion, characterized by a rise in the mid-infrared (MIR) emission accompanied by an accelerated decline in the optical radiation of the SN. Such a dust-formation picture is also corroborated by the concurrent evolution of the profiles of the Ha emission line. Our model suggests enhanced CSM dust concentration at increasing distances from the SN as compared to what can be expected from the density profile of the mass loss from a steady stellar wind. By the time of the last MIR observations at day +1041, a total amount of 1.2+-0.2x10^{-2} Msun of new dust has been formed by SN 2018evt, making SN 2018evt one of the most prolific dust factories among SNe with evidence of dust formation. The unprecedented witness of the intense production procedure of dust may shed light on the perceptions of dust formation in cosmic history.

We present initial results from a JWST survey of the youngest Galactic core-collapse supernova remnant Cassiopeia A (Cas A), made up of NIRCam and MIRI imaging mosaics that map emission from the main shell, interior, and surrounding circumstellar/interstellar material (CSM/ISM). We also present four exploratory positions of MIRI/MRS IFU spectroscopy that sample ejecta, CSM, and associated dust from representative shocked and unshocked regions. Surprising discoveries include: 1) a web-like network of unshocked ejecta filaments resolved to 0.01 pc scales exhibiting an overall morphology consistent with turbulent mixing of cool, low-entropy matter from the progenitor's oxygen layer with hot, high-entropy matter heated by neutrino interactions and radioactivity, 2) a thick sheet of dust-dominated emission from shocked CSM seen in projection toward the remnant's interior pockmarked with small (approximately one arcsecond) round holes formed by knots of high-velocity ejecta that have pierced through the CSM and driven expanding tangential shocks, 3) dozens of light echoes with angular sizes between 0.1 arcsecond to 1 arcminute reflecting previously unseen fine-scale structure in the ISM. NIRCam observations place new upper limits on infrared emission from the neutron star in Cas A's center and tightly constrain scenarios involving a possible fallback disk. These JWST survey data and initial findings help address unresolved questions about massive star explosions that have broad implications for the formation and evolution of stellar populations, the metal and dust enrichment of galaxies, and the origin of compact remnant objects.

ABSTRACT JWST/NIRCam obtained high angular resolution (0.05–0.1 arcsec), deep near-infrared 1–5 $\mu$m imaging of Supernova (SN) 1987A taken 35 yr after the explosion. In the NIRCam images, we identify: (1) faint H2 crescents, which are emissions located between the ejecta and the equatorial ring, (2) a bar, which is a substructure of the ejecta, and (3) the bright 3–5 $\mu$m continuum emission exterior to the equatorial ring. The emission of the remnant in the NIRCam 1–2.3 $\mu$m images is mostly due to line emission, which is mostly emitted in the ejecta and in the hotspots within the equatorial ring. In contrast, the NIRCam 3–5 $\mu$m images are dominated by continuum emission. In the ejecta, the continuum is due to dust, obscuring the centre of the ejecta. In contrast, in the ring and exterior to the ring, synchrotron emission contributes a substantial fraction to the continuum. Dust emission contributes to the continuum at outer spots and diffuse emission exterior to the ring, but little within the ring. This shows that dust cooling and destruction time-scales are shorter than the synchrotron cooling time-scale, and the time-scale of hydrogen recombination in the ring is even longer than the synchrotron cooling time-scale. With the advent of high sensitivity and high angular resolution images provided by JWST/NIRCam, our observations of SN 1987A demonstrate that NIRCam opens up a window to study particle-acceleration and shock physics in unprecedented details, probed by near-infrared synchrotron emission, building a precise picture of how an SN evolves.

We present JWST NIRCam (F356W and F444W filters) and MIRI (F770W) images and NIRSpec- IFU spectroscopy of the young supernova remnant Cassiopeia A (Cas A). We obtained the data as part of a JWST survey of Cas A. The NIRCam and MIRI images map the spatial distributions of synchrotron radiation, Ar-rich ejecta, and CO on both large and small scales, revealing remarkably complex structures. The CO emission is stronger at the outer layers than the Ar ejecta, which indicates the reformation of CO molecules behind the reverse shock. NIRSpec-IFU spectra (3 - 5.5 microns) were obtained toward two representative knots in the NE and S fields. Both regions are dominated by the bright fundamental rovibrational band of CO in the two R and P branches, with strong [Ar VI] and relatively weaker, variable strength ejecta lines of [Si IX], [Ca IV], [Ca V] and [Mg IV]. The NIRSpec-IFU data resolve individual ejecta knots and filaments spatially and in velocity space. The fundamental CO band in the JWST spectra reveals unique shapes of CO, showing a few tens of sinusoidal patterns of rovibrational lines with pseudo-continuum underneath, which is attributed to the high-velocity widths of CO lines. The CO also shows high J lines at different vibrational transitions. Our results with LTE modeling of CO emission indicate a temperature of 1080 K and provide unique insight into the correlations between dust, molecules, and highly ionized ejecta in supernovae, and have strong ramifications for modeling dust formation that is led by CO cooling in the early Universe.

ABSTRACT We have estimated the mass of metals in the molecular gas in 13 dusty star-forming galaxies at $z \sim 4$ in which the gas, based on previous observations, lies in a cold rotating disc. We estimated the metal masses using either the submillimetre line or continuum emission from three tracers of the overall metal content – carbon atoms, carbon monoxide molecules, and dust grains – using the first simultaneous calibration of all three tracers. We obtain very similar mass estimates from the different tracers, which are similar to the entire metal content of a present-day massive early-type galaxy. We used the dynamical masses of these galaxies to estimate an upper limit on the mass of the molecular gas in each galaxy, allowing us to estimate a lower limit on the metal abundance of the gas, finding values for many of the galaxies well above the solar value. We show that the high metal masses and metal abundances are what is expected shortly after the formation of a galaxy for a top-heavy IMF. We suggest a scenario for galaxy evolution in which massive galaxies reach a high metal abundance during their formation phase, which is then gradually reduced by dry mergers with lower mass galaxies. We show that the metals in the outflows from high-redshift dusty star-forming galaxies can quantitatively explain the long-standing puzzle that such a large fraction of the metals in galaxy clusters ($\simeq$0.75) is in the intracluster gas rather than in the galaxies themselves.

The amount of evolution in the dust content of galaxies over the past 5 billion years of cosmic history is contested in the literature. Here, we present a far-infrared (FIR) census of dust based on a sample of 29 241 galaxies with redshifts ranging from 0 < z < 0.5 using data from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). We use the spectral energy distribution fitting tool magphys and a stacking analysis to investigate the evolution of dust mass and temperature of FIR-selected galaxies as a function of both luminosity and redshift. At low redshifts, we find that the mass-weighted and luminosity-weighted dust temperatures from the stacking analysis both exhibit a trend for brighter galaxies to have warmer dust. In higher redshift bins, we see some evolution in both mass-weighted and luminosity-weighted dust temperatures with redshift, but the effect is strongest for luminosity-weighted temperature. The measure of dust content in galaxies at z < 0.1 (the dust mass function) has a different shape to that derived using optically selected galaxies from the same region of sky. We revise the local dust mass density (z < 0.1) to rho(d )=(1.37 +/- 0.08) x 10(5) M(circle dot)Mpc(-3)h(70)(-1); corresponding to an overall fraction of baryons (by mass) stored in dust of f(mb)(dust) = (2.22 +/- 0.13) x 10(-5). We confirm evolution in both the luminosity density and dust mass density over the past few billion years (rho d proportional to (1+z)(2.6 +/- 0.6)), with a flatter evolution than observed in previous FIR-selected studies. We attribute the evolution in rho(L) and rho(m) to an evolution in the dust mass.

JWST Near Infrared Camera (NIRCam) observations at 1.5–4.5 μ m have provided broadband and narrowband imaging of the evolving remnant of SN 1987A with unparalleled sensitivity and spatial resolution. Comparing with previous marginally spatially resolved Spitzer Infrared Array Camera (IRAC) observations from 2004 to 2019 confirms that the emission arises from the circumstellar equatorial ring (ER), and the current brightness at 3.6 and 4.5 μ m was accurately predicted by extrapolation of the declining brightness tracked by IRAC. Despite the regular light curve, the NIRCam observations clearly reveal that much of this emission is from a newly developing outer portion of the ER. Spots in the outer ER tend to lie at position angles in between the well-known ER hotspots. We show that the bulk of the emission in the field can be represented by five standard spectral energy distributions, each with a distinct origin and spatial distribution. This spectral decomposition provides a powerful technique for distinguishing overlapping emission from the circumstellar medium and the supernova ejecta, excited by the forward and reverse shocks, respectively.

We present polarization observations of the young supernova remnant (SNR) Cas A using the High-resolution Airborne Wideband Camera-Plus (HAWC+) instrument onboard the Stratospheric Observatory for Infrared Astronomy (SOFIA). The polarization map at 154 microns reveals dust grains with strong polarization fractions (5 - 30 percent), supporting previous measurements made over a smaller region of the remnant at 850 microns. The 154 microns emission and the polarization signal is coincident with a region of cold dust observed in the southeastern shell and in the unshocked central ejecta. The highly polarized far-IR emission implies the grains are large (greater than 0.14 microns) and silicate-dominated. The polarization level varies across the SNR, with an inverse correlation between the polarization degree and the intensity and smaller polarization angle dispersion for brighter SNR emission. Stronger polarization is detected between the bright structures. This may result from a higher collision rate between the gas and dust producing a lower grain alignment efficiency where the gas density is higher. We use the dust emission to provide an estimate of the magnetic field strength in Cas A using the Davis-Chandrasekhar-Fermi method. The high polarization level is direct evidence that grains are highly elongated and strongly aligned with the magnetic field of the SNR. The dust mass from the polarized region is 0.14+-0.04 Msun, a lower limit of the amount of dust present within the ejecta of Cas A. This result strengthens the hypothesis that core-collapse SNe are an important contributor to the dust mass in high redshift galaxies.

Abstract The amount of cosmic dust contributed by stellar sources in galaxies at all cosmic epochs remains a controversial topic, particularly whether or not supernovae (SNe) have an important role to play given the dust-hostile environments provided by SNe. To date, freshly-formed dust has been observed in a handful of core collapse (CC) SNe, both in the ejecta in-situ and in the interactions between the ejecta and circumstellar medium (CSM). As yet, there exists no clear observational evidence for dust formation in Type Ia SNe despite predictions of $3×10^{-4}$- $0.2$ Msun of dust forming per Ia. Here we report evidence of dust formation in the ejecta-CSM interaction in the Type Ia (SNIa) SN2018evt just three years after the explosion, characterized by a staggering rise in the mid-infrared (MIR) flux accompanied by an accelerated decline in the optical. This hypothesis is strengthened by the concurrent evolution of the profiles of the Hα emission lines. A preexisting hydrogen-rich torus which itself may be dusty before the SN explosion, provides a natural explanation of the observed data. SN 2018evt is the first SNIa with clear evidence of both dust destruction and formation in its ejecta and surroundings. The amount of the newly formed dust follows a steep power law rise of index 4 with time after the explosion. This steep rise indicates that the newly formed dust is formed in the SN ejecta, which is compressed by the shock interaction with the CSM.

We have carried out the first spatially-resolved investigation of the multi-phase interstellar medium (ISM) at high redshift, using the z=4.24 strongly-lensed sub-millimetre galaxy H-ATLASJ142413.9+022303 (ID141). We present high-resolution (down to ~350 pc) ALMA observations in dust continuum emission and in the CO(7-6), H_2O (2_{1,1} - 2_{0,2}), CI(1-0) and CI(2-1) lines, the latter two allowing us to spatially resolve the cool phase of the ISM for the first time. Our modelling of the kinematics reveals that the system appears to be dominated by a rotationally-supported gas disk with evidence of a nearby perturber. We find that the CI(1-0) line has a very different distribution to the other lines, showing the existence of a reservoir of cool gas that might have been missed in studies of other galaxies. We have estimated the mass of the ISM using four different tracers, always obtaining an estimate in the range (3.2-3.8) x 10^{11} M_sol, significantly higher than our dynamical mass estimate of (0.8-1.3) x 10^{11} M_sol. We suggest that this conflict and other similar conflicts reported in the literature is because the gas-to-tracer ratios are ~4 times lower than the Galactic values used to calibrate the ISM in high-redshift galaxies. We demonstrate that this could result from a top-heavy initial mass function and strong chemical evolution. Using a variety of quantitative indicators, we show that, extreme though it is at z=4.24, ID141 will likely join the population of quiescent galaxies that appears in the Universe at z~3.

The Nearby Evolved Stars Survey (NESS) is a volume-complete sample of $\sim$850 Galactic evolved stars within 3\,kpc at (sub-)mm wavelengths, observed in the CO $J = $ (2$-$1) and (3$-$2) rotational lines, and the sub-mm continuum, using the James Clark Maxwell Telescope and Atacama Pathfinder Experiment. NESS consists of five tiers, based on distances and dust-production rate (DPR). We define a new metric for estimating the distances to evolved stars and compare its results to \emph{Gaia} EDR3. Replicating other studies, the most-evolved, highly enshrouded objects in the Galactic Plane dominate the dust returned by our sources, and we initially estimate a total DPR of $4.7\times 10^{-5}$ M$_\odot$ yr$^{-1}$ from our sample. Our sub-mm fluxes are systematically higher and spectral indices are typically shallower than dust models typically predict. The 450/850 $\mu$m spectral indices are consistent with the blackbody Rayleigh--Jeans regime, suggesting a large fraction of evolved stars have unexpectedly large envelopes of cold dust.

Supernova remnants (SNRs) are well-recognised dust producers, but their net dust production rate remains elusive due to uncertainties in grain properties that propagate into observed dust mass uncertainties, and determine how efficiently these grains are processed by reverse shocks. In this paper, we present a detection of polarised dust emission in the Crab pulsar wind nebula, the second SNR with confirmed polarised dust emission after Cassiopeia A. We constrain the bulk composition of the dust with new SOFIA/HAWC+ polarimetric data in band C 89 um and band D 154 um. After correcting for synchrotron polarisation, we report dust polarisation fractions ranging between 3.7-9.6 per cent and 2.7-7.6 per cent in three individual dusty filaments at 89 and 154 um, respectively. The detected polarised signal suggests the presence of large (> 0.05-0.1 um) grains in the Crab Nebula. With the observed polarisation, and polarised and total fluxes, we constrain the temperatures and masses of carbonaceous and silicate grains. We find that the carbon-rich grain mass fraction varies between 12 and 70 per cent, demonstrating that carbonaceous and silicate grains co-exist in this SNR. Temperatures range from 40 K to 70 K and from 30 K to 50 K for carbonaceous and silicate grains, respectively. Dust masses range from 10^{-4} Msol to 10^{-2} Msol for carbonaceous grains and to 10^{-1} Msol for silicate grains, in three individual regions.

ABSTRACT Shockwaves driven by supernovae both destroy dust and reprocess the surviving grains, greatly affecting the resulting dust properties of the interstellar medium (ISM). While these processes have been extensively studied theoretically, observational constraints are limited. We use physically motivated models of dust emission to fit the infrared (IR) spectral energy distributions of seven Galactic supernova remnants, allowing us to determine the distribution of dust mass between diffuse and dense gas phases, and between large and small grain sizes. We find that the dense ($\sim \! 10^3\ {\rm cm}^{-3}$), relatively cool ($\sim \! 10^3\ {\rm K}$) gas phase contains $\gt 90{{\ \rm per\ cent}}$ of the dust mass, making the warm dust located in the X-ray emitting plasma ($\sim \! 1\ {\rm cm}^{-3}$/$10^6\ {\rm K}$) a negligible fraction of the total, despite dominating the mid-IR emission. The ratio of small ($\lesssim\!{10}\ {\rm nm}$) to large ($\gtrsim \! 0.1\ {\rm \mu m}$) grains in the cold component is consistent with that in the ISM, and possibly even higher, whereas the hot phase is almost entirely devoid of small grains. This suggests that grain shattering, which processes large grains into smaller ones, is ineffective in the low-density gas, contrary to model predictions. Single-phase models of dust destruction in the ISM, which do not account for the existence of the cold swept-up material containing most of the dust mass, are likely to greatly overestimate the rate of dust destruction by supernovae.

At a distance of 50 kpc, Supernova 1987A is an ideal target to study how a young supernova (SN) evolves in time. Its equatorial ring, filled with material expelled from the progenitor star about 20,000 years ago, has been engulfed with SN blast waves. Shocks heat dust grains in the ring, emitting their energy at mid-infrared (IR) wavelengths We present ground-based 10--18$\mu$m monitoring of the ring of SN 1987A from day 6067 to 12814 at a resolution of 0.5", together with SOFIA photometry at 10-30 $\mu$m. The IR images in the 2000's (day 6067-7242) showed that the shocks first began brightening the east side of the ring. Later, our mid-IR images from 2017 to 2022 (day 10952-12714) show that dust emission is now fading in the east, while it has brightened on the west side of the ring. Because dust grains are heated in the shocked plasma, which can emit X-rays, the IR and X-ray brightness ratio represent shock diagnostics. Until 2007 the IR to X-ray brightness ratio remained constant over time, and during this time shocks seemed to be largely influencing the east side of the ring. However, since then, the IR to X-ray ratio has been declining, due to increased X-ray brightness. Whether the declining IR brightness is because of dust grains being destroyed or being cooled in the post-shock regions will require more detailed modelling.

ABSTRACT We present a self-consistent cross-calibration of the three main molecular gas mass tracers in galaxies, namely the 12CO(1–0), [C i](3P1–3P0) lines, and the submm dust continuum emission, using a sample of 407 galaxies, ranging from local discs to submillimetre-selected galaxies (SMGs) up to z ≈ 6. A Bayesian statistical method is used to produce galaxy-scale universal calibrations of these molecular gas indicators, that hold over 3–4 orders of magnitude in infrared luminosity, LIR. Regarding the dust continuum, we use a mass-weighted dust temperature, Tmw, determined using new empirical relations between temperature and luminosity. We find the average L/Mmol gas mass conversion factors (including He) to be ${\alpha _{850}}=6.9\times 10^{12}\, \rm W\, Hz^{-1}\, M_{\odot }^{-1}$, αCO = 4.0 M⊙ (κ km s−1 pc2)−1, and $\alpha _{\rm C\, I}$ = 17.0 M⊙ (κ km s−1 pc2)−1, based on the assumption that the mean dust properties of the sample (κH = gas-to-dust ratio/dust emissivity) will be similar to those of local metal rich galaxies and the Milky Way. The tracer with the least intrinsic scatter is [C i](1–0), while CO(1–0) has the highest. The conversion factors show a weak but significant correlation with LIR which is not apparent when Tmw is held constant. Assuming dust properties typical of metal-rich galaxies, we infer a neutral carbon abundance ${X_{\rm C\, I}}=\rm [C^0/\rm{H_2}]=1.6\times 10^{-5}$, similar to that in the Milky Way. We find no evidence for bi-modality of αCO between main-sequence (MS) galaxies and those with extreme star-formation intensity, i.e. ultraluminous infrared galaxies (ULIRGs) and SMGs. The means of the three conversion factors are found to be similar between MS galaxies and ULIRGs/SMGs, to within 10–20 per cent. The overarching conclusion of our work is that, for metal-rich galaxies, near-universal average values for αCO, $X_{\rm C\, I}$, and κH are adequate for global molecular gas estimates within the expected uncertainties. The 1σ scatter in our optimized values for αCO, $X_{\rm C\, I}$, and κH are 0.14, 0.11, and 0.15 dex, respectively.

Dust destruction by supernovae is one of the main processes removing dust from the interstellar medium (ISM). Estimates of the efficiency of this process, both theoretical and observational, typically assume a shock propagating into a homogeneous medium, whereas the ISM possesses significant substructure in reality. We self-consistently model the dust and gas properties of the shocked ISM in three supernova remnants (SNRs), using X-ray and infrared (IR) data combined with corresponding emission models. Collisional heating by gas with properties derived from X-ray observations produces dust temperatures too high to fit the far-IR fluxes from each SNR. An additional colder dust component is required, which has a minimum mass several orders of magnitude larger than that of the warm dust heated by the X-ray emitting gas. Dust-to-gas mass ratios indicate that the majority of the dust in the X-ray emitting material has been destroyed, while the fraction of surviving dust in the cold component is plausibly close to unity. As the cold component makes up virtually all the total dust mass, destruction time-scales based on homogeneous models, which cannot account for multiple phases of shocked gas and dust, may be significantly overestimating actual dust destruction efficiencies, and subsequently underestimating grain lifetimes.

ABSTRACT We investigate the evolution in galactic dust mass over cosmic time through (i) empirically derived dust masses using stacked submillimetre fluxes at 850 μm in the COSMOS field and (ii) dust masses derived using a robust post-processing method on the results from the cosmological hydrodynamical simulation IllustrisTNG. We effectively perform a ‘self-calibration’ of the dust mass absorption coefficient by forcing the model and observations to agree at low redshift and then compare the evolution shown by the observations with that predicted by the model. We create dust mass functions (DMFs) based on the IllustrisTNG simulations from 0 < z < 0.5 and compare these with previously observed DMFs. We find a lack of evolution in the DMFs derived from the simulations, in conflict with the rapid evolution seen in empirically derived estimates of the low-redshift DMF. Furthermore, we observe a strong evolution in the observed mean ratio of dust mass to stellar mass of galaxies over the redshift range 0 < z < 5, whereas the corresponding dust masses from IllustrisTNG show relatively little evolution, even after splitting the sample into satellites and centrals. The large discrepancy between the strong observed evolution and the weak evolution predicted by IllustrisTNG plus post-processing may be explained by either strong cosmic evolution in the properties of the dust grains or limitations in the model. In the latter case, the limitation may be connected to previous claims that the neutral gas content of galaxies does not evolve fast enough in IllustrisTNG.

Observing nearby galaxies with submillimeter telescopes on the ground has two major challenges. First, the brightness is significantly reduced at long submillimeter wavelengths compared to the brightness at the peak of the dust emission. Second, it is necessary to use a high-pass spatial filter to remove atmospheric noise on large angular scales, which has the unwelcome by-product of also removing the galaxy's large-scale structure. We have developed a technique for producing high-resolution submillimeter images of galaxies of large angular size by using the telescope on the ground to determine the small-scale structure (the large Fourier components) and a space telescope (Herschel or Planck) to determine the large-scale structure (the small Fourier components). Using this technique, we are carrying out the HARP and SCUBA-2 High Resolution Terahertz Andromeda Galaxy Survey (HASHTAG), an international Large Program on the James Clerk Maxwell Telescope, with one aim being to produce the first high-fidelity high-resolution submillimeter images of Andromeda. In this paper, we describe the survey, the method we have developed for combining the space-based and ground-based data, and present the first HASHTAG images of Andromeda at 450 and 850um. We also have created a method to predict the CO(J=3-2) line flux across M31, which contaminates the 850um band. We find that while normally the contamination is below our sensitivity limit, the contamination can be significant (up to 28%) in a few of the brightest regions of the 10 kpc ring. We therefore also provide images with the predicted line emission removed.

We present ALMA observations of a small but statistically complete sample of twelve 250 micron selected galaxies at $z=0.35$ designed to measure their dust submillimeter continuum emission as well as their CO(1-0) and atomic carbon [CI](3P1-3P0) spectral lines. This is the first sample of galaxies with global measures of all three $H_2$-mass tracers and which show star formation rates (4-26 Msun yr$^{-1}$) and infra-red luminosities ($1-6\times10^{11}$ Lsun) typical of star forming galaxies in their era. We find a surprising diversity of morphology and kinematic structure; one-third of the sample have evidence for interaction with nearby smaller galaxies, several sources have disjoint dust and gas morphology. Moreover two galaxies have very high $L_{CI}/L_{CO}$ ratios for their global molecular gas reservoirs; if confirmed, such extreme intensity ratios in a sample of dust selected, massive star forming galaxies presents a challenge to our understanding of ISM. Finally, we use the emission of the three molecular gas tracers, to determine the carbon abundance, $X_{ci}$, and CO-$\rm{H_2}$ conversion $\alpha_{co}$ in our sample, using a weak prior that the gas-to-dust ratio is similar to that of the Milky Way for these massive and metal rich galaxies. Using a likelihood method which simultaneously uses all three gas tracer measurements, we find mean values and errors on the mean of $\alpha_{co}=3.0\pm0.5\,\rm{Msun\,(K\,kms^{-1}\,pc^2)^{-1}}$ and $X_{ci}=1.6\pm0.1\times 10^{-5}$ (or $\alpha_{ci}=18.8\,K kms^{-1}\,pc^2 (Msun)^{-1}$) and $\delta_{GDR}=128\pm16$ (or $\alpha_{850}=5.9\times10^{12}\,\rm{W\,Hz^{-1}\, Msun^{-1}}$), where our starting assumption is that these metal rich galaxies have an average gas-to-dust ratio similar to that of the Milky Way centered on $\delta_{GDR}=135$.

We build a rigorous statistical framework to provide constraints on the chemical and dust evolution parameters for nearby late-type galaxies with a wide range of gas fractions (3 per cent < f(g) < 94 per cent). A Bayesian Monte Carlo Markov Chain framework provides statistical constraints on the parameters used in chemical evolution models. Nearly a million one-zone chemical and dust evolution models were compared to 340 galaxies. Relative probabilities were calculated from the chi(2) between data and models, marginalized over the different time-steps, galaxy masses, and star formation histories. We applied this method to find 'best-fitting' model parameters related to metallicity, and subsequently fix these metal parameters to study the dust parameters. For the metal parameters, a degeneracy was found between the choice of initial mass function, supernova metal yield tables, and outflow prescription. For the dust parameters, the uncertainties on the best-fitting values are often large except for the fraction of metals available for grain growth, which is well constrained. We find a number of degeneracies between the dust parameters, limiting our ability to discriminate between chemical models using observations only. For example, we show that the low dust content of low-metallicity galaxies can be resolved by either reducing the supernova dust yields and/or including photofragmentation. We also show that supernova dust dominates the dust mass for low-metallicity galaxies and grain growth dominates for high-metallicity galaxies. The transition occurs around 12 + log(O/H) = 7.75, which is lower than found in most studies in the literature.