Life Cycle of Dust in the Magellanic Clouds and the Milky Way

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.Comment: 18 pages, 6 figures, published recently online in the Proceedings of Science for the conference: Life Cycle of Dust in the Universe, Observations, Theory and Laboratory Experiment

Dust input from AGB stars in the Large Magellanic Cloud

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.Comment: Accepted to A&

Collapse of massive magnetized dense cores using radiation-magneto-hydrodynamics: early fragmentation inhibition

We report the results of radiation-magneto-hydrodynamics calculations in the context of high mass star formation, using for the first time a self-consistent model for photon emission (i.e. via thermal emission and in radiative shocks) and with the high resolution necessary to resolve properly magnetic braking effects and radiative shocks on scales <100 AU. We investigate the combined effects of magnetic field, turbulence, and radiative transfer on the early phases of the collapse and the fragmentation of massive dense cores. We identify a new mechanism that inhibits initial fragmentation of massive dense cores, where magnetic field and radiative transfer interplay. We show that this interplay becomes stronger as the magnetic field strength increases. Magnetic braking is transporting angular momentum outwards and is lowering the rotational support and is thus increasing the infall velocity. This enhances the radiative feedback owing to the accretion shock on the first core. We speculate that highly magnetized massive dense cores are good candidates for isolated massive star formation, while moderately magnetized massive dense cores are more appropriate to form OB associations or small star clusters.Comment: Accepted for publication in ApJL, 19 pages, 4 figure

Hidden Q-structure and Lie 3-algebra for non-abelian superconformal models in six dimensions

We disclose the mathematical structure underlying the gauge field sector of the recently constructed non-abelian superconformal models in six spacetime dimensions. This is a coupled system of 1-form, 2-form, and 3-form gauge fields. We show that the algebraic consistency constraints governing this system permit to define a Lie 3-algebra, generalizing the structural Lie algebra of a standard Yang-Mills theory to the setting of a higher bundle. Reformulating the Lie 3-algebra in terms of a nilpotent degree 1 BRST-type operator Q, this higher bundle can be compactly described by means of a Q-bundle; its fiber is the shifted tangent of the Q-manifold corresponding to the Lie 3-algebra and its base the odd tangent bundle of spacetime equipped with the de Rham differential. The generalized Bianchi identities can then be retrieved concisely from Q^2=0, which encode all the essence of the structural identities. Gauge transformations are identified as vertical inner automorphisms of such a bundle, their algebra being determined from a Q-derived bracket.Comment: 51 pages, 3 figure

High-resolution simulations of planetesimal formation in turbulent protoplanetary discs

We present high-resolution computer simulations of dust dynamics and planetesimal formation in turbulence generated by the magnetorotational instability. We show that the turbulent viscosity associated with magnetorotational turbulence in a non-stratified shearing box increases when going from 256^3 to 512^3 grid points in the presence of a weak imposed magnetic field, yielding a turbulent viscosity of $\alpha\approx0.003$ at high resolution. Particles representing approximately meter-sized boulders concentrate in large-scale high-pressure regions in the simulation box. The appearance of zonal flows and particle concentration in pressure bumps is relatively similar at moderate (256^3) and high (512^3) resolution. In the moderate-resolution simulation we activate particle self-gravity at a time when there is little particle concentration, in contrast with previous simulations where particle self-gravity was activated during a concentration event. We observe that bound clumps form over the next ten orbits, with initial birth masses of a few times the dwarf planet Ceres. At high resolution we activate self-gravity during a particle concentration event, leading to a burst of planetesimal formation, with clump masses ranging from a significant fraction of to several times the mass of Ceres. We present a new domain decomposition algorithm for particle-mesh schemes. Particles are spread evenly among the processors and the local gas velocity field and assigned drag forces are exchanged between a domain-decomposed mesh and discrete blocks of particles. We obtain good load balancing on up to 4096 cores even in simulations where particles sediment to the mid-plane and concentrate in pressure bumps.Comment: Accepted for publication in Astronomy & Astrophysics, with some changes in response to referee repor

Iron and silicate dust growth in the Galactic interstellar medium: clues from element depletions

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.Comment: Accepted for publication in the ApJ, 15 pages, 9 figure