We study the destruction of interstellar dust via sputtering in supernova (SN) shocks using three-dimensional hydrodynamical simulations. With a novel numerical framework, we follow both sputtering and dust dynamics governed by direct collisions, plasma drag and betatron acceleration. Grain-grain collisions are not included and the grain-size distribution is assumed to be fixed. The amount of dust destroyed per SN is quantified for a broad range of ambient densities and fitting formulae are provided. Integrated over the grain-size distribution, nonthermal (inertial) sputtering dominates over thermal sputtering for typical ambient densities. We present the first simulations that explicitly follow dust sputtering within a turbulent multiphase interstellar medium. We find that the dust destruction timescales $\tau$ are 0.35 Gyr for silicate dust and 0.44 Gyr for carbon dust in solar neighborhood conditions. The SN environment has an important impact on $\tau$. SNe that occur in preexisting bubbles destroy less dust as the destruction is limited by the amount of dust in the shocked gas. This makes $\tau$ about 2.5 times longer than the estimate based on results from a single SN explosion. We investigate the evolution of the dust-to-gas mass ratio (DGR), and find that a spatial inhomogeneity of $\sim$ 14\% develops for scales below 10 pc. It locally correlates positively with gas density but negatively with gas temperature even in the exterior of the bubbles due to incomplete gas mixing. This leads to a $\sim$ 30\% lower DGR in the volume filling warm gas compared to that in the dense clouds.

The traditional picture of post-starburst galaxies as dust- and gas-poor merger remnants, rapidly transitioning to quiescence, has been recently challenged. Unexpected detections of a significant ISM in many post-starbursts raise important questions. Are they truly quiescent and, if so, what mechanisms inhibit further star formation? What processes dominate their ISM energetics? We present an infrared spectroscopic and photometric survey of 33 SDSS-selected E+A post-starbursts, aimed at resolving these questions. We find compact, warm dust reservoirs with high PAH abundances, and total gas and dust masses significantly higher than expected from stellar recycling alone. Both PAH/TIR and dust-to-burst stellar mass ratios are seen to decrease with post-burst age, indicative of the accumulating effects of dust destruction and an incipient transition to hot, early-type ISM properties. Their infrared spectral properties are unique, with dominant PAH emission, very weak nebular lines, unusually strong H$_{2}$ rotational emission, and deep ${\rm [C\, II]}$ deficits. There is substantial scatter among SFR indicators, and both PAH and TIR luminosities provide overestimates. Even as potential upper limits, all tracers show that the SFR has typically experienced a more than two order-of-magnitude decline since the starburst, and that the SFR is considerably lower than expected given both their stellar masses and molecular gas densities. These results paint a coherent picture of systems in which star formation was, indeed, rapidly truncated, but in which the ISM was $\textit{not}$ completely expelled, and is instead supported against collapse by latent or continued injection of turbulent or mechanical heating. The resulting aging burst populations provide a "high-soft" radiation field which seemingly dominates the E+As' unusual ISM energetics.

The interstellar abundances of refractory elements indicate a substantial depletion from the gas phase, that increases with gas density. Our recent model of dust evolution, based on hydrodynamic simulations of the lifecycle of giant molecular clouds (GMCs) proves that the observed trend for [Si$_{gas}$/H] is driven by a combination of dust growth by accretion in the cold diffuse interstellar medium (ISM) and efficient destruction by supernova (SN) shocks (Zhukovska et al. 2016). With an analytic model of dust evolution, we demonstrate that even with optimistic assumptions for the dust input from stars and without destruction of grains by SNe it is impossible to match the observed [Si$_{gas}$/H]$-n_H$ relation without growth in the ISM. We extend the framework developed in our previous work for silicates to include the evolution of iron grains and address a long-standing conundrum: ``Where is the interstellar iron?'. Much higher depletion of Fe in the warm neutral medium compared to Si is reproduced by the models, in which a large fraction of interstellar iron (70%) is locked as inclusions in silicate grains, where it is protected from sputtering by SN shocks. The slope of the observed [Fe$_{gas}$/H]$-n_H$ relation is reproduced if the remaining depleted iron resides in a population of metallic iron nanoparticles with sizes in the range of 1-10nm. Enhanced collision rates due to the Coulomb focusing are important for both silicate and iron dust models to match the observed slopes of the relations between depletion and density and the magnitudes of depletion at high density.

Dust grains are an important component of the interstellar medium (ISM) of galaxies. We present the first direct measurement of the residence times of interstellar dust in the different ISM phases, and of the transition rates between these phases, in realistic hydrodynamical simulations of the multi-phase ISM. Our simulations include a time-dependent chemical network that follows the abundances of H^+, H, H_2, C^+ and CO and take into account self-shielding by gas and dust using a tree-based radiation transfer method. Supernova explosions are injected either at random locations, at density peaks, or as a mixture of the two. For each simulation, we investigate how matter circulates between the ISM phases and find more sizeable transitions than considered in simple mass exchange schemes in the literature. The derived residence times in the ISM phases are characterised by broad distributions, in particular for the molecular, warm and hot medium. The most realistic simulations with random and mixed driving have median residence times in the molecular, cold, warm and hot phase around 17, 7, 44 and 1 Myr, respectively. The transition rates measured in the random driving run are in good agreement with observations of Ti gas-phase depletion in the warm and cold phases in a simple depletion model, although the depletion in the molecular phase is under-predicted. ISM phase definitions based on chemical abundance rather than temperature cuts are physically more meaningful, but lead to significantly different transition rates and residence times because there is no direct correspondence between the two definitions.

We develop a model of dust evolution in a multiphase, inhomogeneous ISM including dust growth and destruction processes. The physical conditions for grain evolution are taken from hydrodynamical simulations of giant molecular clouds in a Milky Way-like spiral galaxy. We improve the treatment of dust growth by accretion in the ISM to investigate the role of the temperature-dependent sticking coefficient and ion-grain interactions. From detailed observational data on the gas-phase Si abundances [Si/H]_{gas} measured in the local Galaxy, we derive a relation between the average [Si/H]_{gas} and the local gas density n(H) which we use as a critical constraint for the models. This relation requires a sticking coefficient that decreases with the gas temperature. The synthetic relation constructed from the spatial dust distribution reproduces the slope of -0.5 of the observed relation in cold clouds. This slope is steeper than that for the warm medium and is explained by the dust growth. We find that it occurs for all adopted values of the minimum grain size a_{min} from 1 to 5nm. For the classical cut-off of a_{min}=5 nm, the ion-grain interactions result in longer growth timescales and higher [Si/H]_{gas} than the observed values. For a_{min} below 3 nm, the ion-grain interactions enhance the growth rates, steepen the slope of [Si/H]_{gas}-n(H) relation and provide a better match to observations. The rates of dust re-formation in the ISM by far exceed the rates of dust production by stellar sources as expected from simple evolution models. After the cycle of matter in and out of dust reaches a steady state, the dust growth balances the destruction operating on similar timescales of 350 Myr.

We combine Herschel/SPIRE sub-millimeter (submm) observations with existing multi-wavelength data to investigate the characteristics of low redshift, optically red galaxies detected in submm bands. We select a sample of galaxies in the redshift range 0.01$\leq$z$\leq$0.2, having >5$\sigma$ detections in the SPIRE 250 micron submm waveband. Sources are then divided into two sub-samples of $red$ and $blue$ galaxies, based on their UV-optical colours. Galaxies in the $red$ sample account for $\approx$4.2 per cent of the total number of sources with stellar masses M$_{*}\gtrsim$10$^{10}$ Solar-mass. Following visual classification of the $red$ galaxies, we find that $\gtrsim$30 per cent of them are early-type galaxies and $\gtrsim$40 per cent are spirals. The colour of the $red$-spiral galaxies could be the result of their highly inclined orientation and/or a strong contribution of the old stellar population. It is found that irrespective of their morphological types, $red$ and $blue$ sources occupy environments with more or less similar densities (i.e., the $\Sigma_5$ parameter). From the analysis of the spectral energy distributions (SEDs) of galaxies in our samples based on MAGPHYS, we find that galaxies in the $red$ sample (of any morphological type) have dust masses similar to those in the $blue$ sample (i.e. normal spiral/star-forming systems). However, in comparison to the $red$-spirals and in particular $blue$ systems, $red$-ellipticals have lower mean dust-to-stellar mass ratios. Besides galaxies in the $red$-elliptical sample have much lower mean star-formation/specific-star-formation rates in contrast to their counterparts in the $blue$ sample. Our results support a scenario where dust in early-type systems is likely to be of an external origin.

We examine the fraction of massive asymptotic giant branch (AGB) stars remaining bound in their parent star clusters and the effect of irradiation of these stars by intracluster ultraviolet (UV) field. We employ a set of N-body models of dynamical evolution of star clusters rotating in a galactic potential at the solar galactocentric radius. The cluster models are combined with stellar evolution formulae, a library of stellar spectra, and simple models for SiO photodissociation in circumstellar environment (CSE). The initial stellar masses of clusters are varied from $50\rm M_\odot$ to $10^{5}\rm M_\odot$. Results derived for individual clusters are combined using a mass distribution function for young star clusters. We find that about 30% of massive AGB stars initially born in clusters become members of the field population, while the rest evolves in star clusters. They are irradiated by strong intracluster UV radiation resulting in the decrease of the photodissociation radius of SiO molecules, in many stars down to the dust formation zone. In absence of dust shielding, the UV photons penetrate in the CSE deeper than $10R_*$ in 64% and deeper than $2 R_*$ in 42% of all massive AGB stars. If this suppresses following dust formation, the current injection rate of silicate dust from AGB stars in the local Galaxy decreases from $2.2 \times 10^{-4}\rm M_\odot\,kpc^{-2}\,Gyr^{-1}$ to $1.8 \times 10^{-4}\rm M_\odot\,kpc^{-2}\,Gyr^{-1}$ at most. A lower revised value of 40% for the expected fraction of presolar silicate grains from massive AGB stars is still high to explain the non-detection of these grains in meteorites.

The spatial variations of the gas-to-dust ratio (GDR) provide constraints on the chemical evolution and lifecycle of dust in galaxies. We examine the relation between dust and gas at 10-50 pc resolution in the Large and Small Magellanic Clouds (LMC and SMC) based on Herschel far-infrared (FIR), H I 21 cm, CO, and Halpha observations. In the diffuse atomic ISM, we derive the gas-to-dust ratio as the slope of the dust-gas relation and find gas-to-dust ratios of 380+250-130 in the LMC, and 1200+1600-420 in the SMC, not including helium. The atomic-to-molecular transition is located at dust surface densities of 0.05 Mo pc-2 in the LMC and 0.03 Mo pc-2 in the SMC, corresponding to AV ~ 0.4 and 0.2, respectively. We investigate the range of CO-to-H2 conversion factor to best account for all the molecular gas in the beam of the observations, and find upper limits on XCO to be 6x1020 cm-2 K-1 km-1 s in the LMC (Z=0.5Zo) at 15 pc resolution, and 4x 1021 cm-2 K-1 km-1 s in the SMC (Z=0.2Zo) at 45 pc resolution. In the LMC, the slope of the dust-gas relation in the dense ISM is lower than in the diffuse ISM by a factor ~2, even after accounting for the effects of CO-dark H2 in the translucent envelopes of molecular clouds. Coagulation of dust grains and the subsequent dust emissivity increase in molecular clouds, and/or accretion of gas-phase metals onto dust grains, and the subsequent dust abundance (dust-to-gas ratio) increase in molecular clouds could explain the observations. In the SMC, variations in the dust-gas slope caused by coagulation or accretion are degenerate with the effects of CO-dark H2. Within the expected 5--20 times Galactic XCO range, the dust-gas slope can be either constant or decrease by a factor of several across ISM phases. Further modeling and observations are required to break the degeneracy between dust grain coagulation, accretion, and CO-dark H2.

To a great extent, our understanding of the life cycle of dust is based on the observational and theoretical studies of the Milky Way and the Magellanic Clouds, which will be the topic of this contribution. Over past years, a large volume of observations with unprecedented spatial resolution has been accumulated for the Milky Way. It permits investigations of different stages of the life cycle of dust, from its formation in stellar sources to destruction in star-forming regions and supernovae shocks. Observations of dust emission, extinction, polarisation of light, and interstellar element depletions in the solar neighbourhood provide the most accurate constraints for the reference dust models applied to study extragalactic systems. However, global spatial studies of the circumstellar and interstellar dust are complicated in the Milky Way disk because of high extinction, confusion along the line of sight and large uncertainties in distances. In contrast, the favourable location in the sky and the proximity of the Magellanic Clouds allow detailed multi-wavelength studies of the dust-forming stellar populations and the investigation of variations of the interstellar grain properties for the entire galaxies. They enable the first comparison between the global stardust production rates from theoretical calculations and those from observations, which confirm discrepancy between accumulated stardust mass and observed interstellar dust mass - "the missing dust-source" problem. Modelling of the life cycle of dust in the Large Magellanic Cloud showed that dust growth by mantle accretion in the ISM, a major dust source in the Milky Way, can be responsible for the existing dust mass in the LMC. We will present comparison of the dust input from different sources to the dust budgets of the Milky Way and the Magellanic Clouds, which reveals how the role of these dust sources depends on metallicity.

We re-evaluate the roles of different dust sources in dust production as a function of metallicity in late-type dwarf galaxies, with the goal of understanding the relation between dust content and metallicity. The dust content of late-type dwarf galaxies with episodic star formation is studied with a multicomponent model of dust evolution, which includes dust input from AGB stars, type II SNe and dust mass growth by accretion of gas species in the ISM. Dust growth in the ISM becomes an important dust source in dwarf galaxies, on the timescale of 0.1 - few Gyrs. It increases the dust-to-gas ratio (DGR) during post-burst evolution, unlike type II SNe, which eject grains into the ISM only during starbursts. Before the dust growth in the ISM overtakes the dust production, AGB stars can be major sources of dust in metal-poor dwarf galaxies. Our models reproduce the relation between the DGR and oxygen abundance, derived from observations of a large sample of dwarf galaxies. The steep decrease in the DGR at low O values is explained by the relatively low efficiency of dust condensation in stars. The scatter observed at higher O values is determined mainly by different critical metallicities for the transition from stardust- to ISM-growth dominated dust production, depending on the star formation history. In galaxies with episodic star formation, additional dispersion in the DGR is introduced by grain destruction during starbursts, followed by an increase of the dust mass due to dust growth in the ISM during post-burst evolution. We find that the carbon-to-silicate ratio changes dramatically, when the ISM growth becomes the dominant dust source, therefore this ratio can be used as an indicator of the transition. The observed DGR-O relation in dwarf galaxies favours low condensation efficiencies in type II SNe, together with an increase in the total dust mass by means of dust growth in the ISM.

The formation of grains in the interstellar medium, i.e., at low temperature, has been proposed as a possibility to solve the lifetime problem of cosmic dust. This process lacks a firm experimental basis, which is the goal of this study. We have investigated the condensation of SiO molecules at low temperature using neon matrix and helium droplet isolation techniques. The energies of SiO polymerization reactions have been determined experimentally with a calorimetric method and theoretically with calculations based on the density functional theory. The combined experimental and theoretical values have revealed the formation of cyclic (SiO)$_k$ ($k$ = 2--3) clusters inside helium droplets at $T$ = 0.37 K. Therefore, the oligomerization of SiO molecules is found to be barrierless and is expected to be fast in the low-temperature environment of the interstellar medium on the surface of dust grains. The incorporation of numerous SiO molecules in helium droplets leads to the formation of nanoscale amorphous SiO grains. Similarly, the annealing and evaporation of SiO-doped Ne matrices lead to the formation of solid amorphous SiO on the substrate. The structure and composition of the grains were determined by infrared absorption spectroscopy, transmission electron microscopy, and energy-dispersive X-ray spectroscopy. Our results support the hypothesis that interstellar silicates \textbf{can be formed} in the low temperature regions of the interstellar medium by accretion through barrierless reactions.

This paper analyses the behaviour of the gas-to-dust mass ratio (G/D) of local Universe galaxies over a large metallicity range. We combine three samples: the Dwarf Galaxy Survey, the KINGFISH survey and a subsample from Galametz et al. (2011) totalling 126 galaxies, covering a 2 dex metallicity range, with 30% of the sample with 12+log(O/H) < 8.0. The dust masses are homogeneously determined with a semi-empirical dust model, including submm constraints. The atomic and molecular gas masses are compiled from the literature. Two XCO are used to estimate molecular gas masses: the Galactic XCO, and a XCO depending on the metallicity (as Z^{-2}). Correlations with morphological types, stellar masses, star formation rates and specific star formation rates are discussed. The trend between G/D and metallicity is empirically modelled using power-laws (slope of -1 and free) and a broken power-law. We compare the evolution of the G/D with predictions from chemical evolution models. We find that out of the five tested galactic parameters, metallicity is the galactic property driving the observed G/D. The G/D versus metallicity relation cannot be represented by a power-law with a slope of -1 over the whole metallicity range. The observed trend is steeper for metallicities lower than ~ 8.0. A large scatter is observed in the G/D for a given metallicity, with a dispersion of 0.37 dex in metallicity bins of ~0.1 dex. The broken power-law reproduces best the observed G/D and provides estimates of the G/D that are accurate to a factor of 1.6. The good agreement of the G/D and its scatter with the three tested chemical evolution models shows that the scatter is intrinsic to galactic properties, reflecting the different star formation histories, dust destruction efficiencies, dust grain size distributions and chemical compositions across the sample. (abriged)

The dust-forming population of AGB stars and their input to the interstellar dust budget of the Large Magellanic Cloud (LMC) are studied with evolutionary dust models with the main goals (1) to investigate how the amount and composition of dust from AGB stars vary over galactic history; (2) to characterise the mass and metallicity distribution of the present population of AGB stars; (3) to quantify the contribution of AGB stars of different mass and metallicity to the present stardust population in the interstellar medium (ISM). We use models of the stardust lifecycle in the ISM developed and tested for the Solar neighbourhood. The first global spatially resolved reconstruction of the star formation history of the LMC from the Magellanic Clouds Photometric Survey is employed to calculate the stellar populations in the LMC. The dust input from AGB stars is dominated by carbon grains from stars with masses < 4 Msun almost over the entire history of the LMC. The production of silicate, silicon carbide and iron dust is delayed until the ISM is enriched to about half the present metallicity in the LMC. For the first time, theoretically calculated dust production rates of AGB stars are compared to those derived from IR observations of AGB stars for the entire galaxy. We find good agreement within scatter of various observational estimates. We show that the majority of silicate and iron grains in the present stardust population originate from a small population of intermediate-mass stars consisting of only about 4% of the total number of stars, whereas in the Solar neighbourhood they originate from low-mass stars. With models of the lifecycle of stardust grains in the ISM we confirm a large discrepancy between dust input from stars and the existing interstellar dust mass in the LMC reported in Matsuura et al. 2009.

We present a self-consistent model of a protoplanetary disk: 'ANDES' ('AccretioN disk with Dust Evolution and Sedimentation'). ANDES is based on a flexible and extendable modular structure that includes 1) a 1+1D frequency-dependent continuum radiative transfer module, 2) a module to calculate the chemical evolution using an extended gas-grain network with UV/X-ray-driven processes surface reactions, 3) a module to calculate the gas thermal energy balance, and 4) a 1+1D module that simulates dust grain evolution. For the first time, grain evolution and time-dependent molecular chemistry are included in a protoplanetary disk model. We find that grain growth and sedimentation of large grains to the disk midplane lead to a dust-depleted atmosphere. Consequently, dust and gas temperatures become higher in the inner disk (R < 50 AU) and lower in the outer disk (R > 50 AU), in comparison with the disk model with pristine dust. The response of disk chemical structure to the dust growth and sedimentation is twofold. First, due to higher transparency a partly UV-shielded molecular layer is shifted closer to the dense midplane. Second, the presence of big grains in the disk midplane delays the freeze-out of volatile gas-phase species such as CO there, while in adjacent upper layers the depletion is still effective. Molecular concentrations and thus column densities of many species are enhanced in the disk model with dust evolution, e.g., CO2, NH2CN, HNO, H2O, HCOOH, HCN, CO. We also show that time-dependent chemistry is important for a proper description of gas thermal balance.

We study the impact of dust evolution in a protoplanetary disk around a T Tauri star on the disk chemical composition. For the first time we utilize a comprehensive model of dust evolution which includes growth, fragmentation and sedimentation. Specific attention is paid to the influence of grain evolution on the penetration of the UV field in the disk. A chemical model that includes a comprehensive set of gas phase and grain surface chemical reactions is used to simulate the chemical structure of the disk. The main effect of the grain evolution on the disk chemical composition comes from sedimentation, and, to a lesser degree, from the reduction of the total grain surface area. The net effect of grain growth is suppressed by the fragmentation process which maintains a population of small grains, dominating the total grain surface area. We consider three models of dust properties. In model GS both growth and sedimentation are taken into account. In models A5 and A4 all grains are assumed to have the same size (10(-5) cm and 10(-4) cm, respectively) with constant gas-to-dust mass ratio of 100. Like in previous studies, the "three-layer" pattern (midplane, molecular layer, hot atmosphere) in the disk chemical structure is preserved in all models, but shifted closer to the midplane in models with increased grain size (GS and A4). Unlike other similar studies, we find that in models GS and A4 column densities of most gas-phase species are enhanced by 1-3 orders of magnitude relative to those in a model with pristine dust (A5), while column densities of their surface counterparts are decreased. We show that column densities of certain species, like C2H, HC(2n+1)N (n=0-3), H2O and some other molecules, as well as the C2H2/HCN abundance ratio which are accessible with Herschel and ALMA can be used as observational tracers of early stages of the grain evolution process in protoplanetary disks.

The basic mechanism responsible for the widespread condensation of MgS in the outflows from carbon rich stars on the tip of the AGB is discussed with the aim of developing a condensation model that can be applied in model calculations of dust formation in stellar winds. The different possibilities how MgS may be formed in the chemical environment of outflows from carbon stars are explored by some thermochemical calculations and by a detailed analysis of the growth kinetics of grains in stellar winds. The optical properties of core-mantle grains with a MgS mantle are calculated to demonstrate that such grains reproduce the structure of the observed 30 $\mu$m feature. These considerations are complemented by model calculations of circumstellar dust shells around carbon stars. It is argued that MgS is formed via precipitation on silicon carbide grains. This formation mechanism explains some of the basic observed features of MgS condensation in dust shells around carbon stars. A weak secondary peak at about 33 ... 36 $\mu$m is shown to exist in certain cases if MgS forms a coating on SiC.

The abundance evolution of interstellar dust species originating from stellar sources and from condensation in molecular clouds in the local interstellar medium of the Milky Way is studied and the input of dust material to the Solar System is determined. A one-zone chemical evolution model of the Milky Way for the elemental composition of the disk combined with an evolution model for its interstellar dust component similar to that of Dwek (1998) is developed. The dust model considers dust-mass return from AGB stars as calculated from synthetic AGB models combined with models for dust condensation in stellar outflows. Supernova dust formation is included in a simple parameterized form which is gauged by observed abundances of presolar dust grains with supernova origin. For dust growth in the ISM a simple method is developed for coupling this with disk and dust evolution models. The time evolution of the abundance of the following dust species is followed in the model: silicate, carbon, silicon carbide, and iron dust from AGB stars and from SNe as well as silicate, carbon, and iron dust grown in molecular clouds. It is shown that the interstellar dust population is dominated by dust accreted in molecular clouds; most of the dust material entering the Solar System at its formation does not show isotopic abundance anomalies of the refractory elements, i.e., inconspicuous isotopic abundances do not point to a Solar System origin of dust grains. The observed abundance ratios of presolar dust grains formed in SN ejecta and in AGB star outflows requires that for the ejecta from SNe the fraction of refractory elements condensed into dust is 0.15 for carbon dust and is quite small ($\sim10^{-4}$) for other dust species.