# Posts Tagged dust grains

## Recent Postings from dust grains

### Chemical Enrichment of the Pre-Solar Cloud by Supernova Dust Grains

The presence of short-lived radioisotopes (SLRs) in solar system meteorites has been interpreted as evidence that the solar system was exposed to a supernova shortly before or during its formation. Yet results from hydrodynamical models of SLR injection into the proto-solar cloud or disc suggest that gas-phase mixing may not be efficient enough to reproduce the observed abundances. As an alternative, we explore the injection of SLRs via dust grains as a way to overcome the mixing barrier. We numerically model the interaction of a supernova remnant containing SLR-rich dust grains with a nearby molecular cloud. The dust grains are subject to drag forces and both thermal and non-thermal sputtering. We confirm that the expanding gas shell stalls upon impact with the dense cloud and that gas-phase SLR injection occurs slowly due to hydrodynamical instabilities at the cloud surface. In contrast, dust grains of sufficient size (> 1 micron) decouple from the gas and penetrate into the cloud within 0.1 Myr. Once inside the cloud, the dust grains are destroyed by sputtering, releasing SLRs and rapidly enriching the dense (potentially star-forming) regions. Our results suggest that SLR transport on dust grains is a viable mechanism to explain SLR enrichment.

### Modelling the nebular emission from primeval to present-day star-forming galaxies

We present a new model of the nebular emission from star-forming galaxies in a wide range of chemical compositions, appropriate to interpret observations of galaxies at all cosmic epochs. The model relies on the combination of state-of-the-art stellar population synthesis and photoionization codes to describe the ensemble of HII regions and the diffuse gas ionized by young stars in a galaxy. A main feature of this model is the self-consistent yet versatile treatment of element abundances and depletion onto dust grains, which allows one to relate the observed nebular emission from a galaxy to both gas-phase and dust-phase metal enrichment. We show that this model can account for the rest-frame ultraviolet and optical emission-line properties of galaxies at different redshifts and find that ultraviolet emission lines are more sensitive than optical ones to parameters such as C/O abundance ratio, hydrogen gas density, dust-to-metal mass ratio and upper cutoff of the stellar initial mass function. We also find that, for gas-phase metallicities around solar to slightly sub-solar, widely used formulae to constrain oxygen ionic fractions and the C/O ratio from ultraviolet and optical emission-line luminosities are reasonable faithful. However, the recipes break down at non-solar metallicities, making them inappropriate to study chemically young galaxies. In such cases, a fully self-consistent model of the kind presented in this paper is required to interpret the observed nebular emission.

### Modelling extragalactic extinction through gamma-ray burst afterglows

We analyze extragalactic extinction pro?les derived through gamma-ray burst afterglows, using a dust model speci?cally constructed on the assumption that dust grains are not immutable but respond time-dependently to the local physics. Such a model includes core-mantle spherical particles of mixed chemical composition (silicate core, sp2 and sp3 carbonaceous layers), and an additional molecular component, in the form of free-flying polycyclic aromatic hydrocarbons. We fit most of the observed extinction pro?les. Failures occur for lines of sight presenting remarkable rises blueward the bump. We find a tendency in the carbon chemical structure to become more aliphatic with the galactic activity, and to some extent with increasing redshifts. Moreover, the contribution of the moleclar component to the total extinction is more important in younger objects. The results of the ?tting procedure (either successes and failures) may be naturally interpreted through an evolutionary prescription based on the carbon cycle in the interstellar medium of galaxies.

### Chemical differentiation in a prestellar core traces non-uniform illumination

Dense cloud cores present chemical differentiation due to the different distribution of C-bearing and N-bearing molecules, the latter being less affected by freeze-out onto dust grains. In this letter we show that two C-bearing molecules, CH$_3$OH and $c$-C$_3$H$_2$, present a strikingly different (complementary) morphology while showing the same kinematics toward the prestellar core L1544. After comparing their distribution with large scale H$_2$ column density N(H$_2$) map from the Herschel satellite, we find that these two molecules trace different environmental conditions in the surrounding of L1544: the $c$-C$_3$H$_2$ distribution peaks close to the southern part of the core, where the surrounding molecular cloud has a N(H$_2$) sharp edge, while CH$_3$OH mainly traces the northern part of the core, where N(H$_2$) presents a shallower tail. We conclude that this is evidence of chemical differentiation driven by different amount of illumination from the interstellar radiation field: in the South, photochemistry maintains more C atoms in the gas phase allowing carbon chain (such as $c$-C$_3$H$_2$) production; in the North, C is mainly locked in CO and methanol traces the zone where CO starts to freeze out significantly. During the process of cloud contraction, different gas and ice compositions are thus expected to mix toward the central regions of the core, where a potential Solar-type system will form. An alternative view on carbon-chain chemistry in star-forming regions is also provided.

### Investigating dust trapping in transition disks with millimeter-wave polarization

Spatially resolved polarized (sub-)mm emission has been observed for example in the protoplanetary disk around HL Tau. Magnetically aligned grains are commonly interpreted as the source of polarization. However, self-scattering by large dust grains with a high enough albedo is another polarization mechanism, becoming a compelling method independent of the spectral index to constrain the dust grain size in protoplanetary disks. We study the dust polarization at mm wavelength in the dust trapping scenario proposed for transition disks, when a giant planet opens a gap in the disk. We investigate the characteristic polarization patterns and their dependence on disk inclination, dust size evolution, planet position, and observing wavelength. We combine two-dimensional hydrodynamical simulations of planet-disk interactions with self-consistent dust growth models. These size-dependent dust density distributions are used for follow-up three-dimensional radiative transfer calculations to predict the polarization degree at ALMA bands due to scattered thermal emission. We find that the polarization pattern of a disk with a planetary gap after 1 Myr of dust evolution shows a distinctive three ring structure. Two narrow inner rings are located at the planet gap edges. For increasing observing wavelengths all three rings slightly change their position, where the innermost and outermost rings move inward. This distance is detectable comparing the results at ALMA bands 3, 6 and 7. Within the highest polarized intensity regions the polarization vectors are oriented in the azimuthal direction. For an inclined disk there is an interplay between polarization originating from a flux gradient and inclination-induced quadrupole polarization. For intermediate inclined transition disks the polarization degree is as high as ~ 2% at band 3, which is well above the detection limit of future ALMA observations.

### Variations of the Interstellar Extinction Law within the Nearest Kiloparsec

Multicolor photometry from the Tycho-2 and 2MASS catalogues for 11 990 OB and 30 671 K-type red giant branch stars is used to detect systematic large-scale variations of the interstellar extinction law within the nearest kiloparsec. The characteristic of the extinction law, the total-to-selective extinction ratio $R_V$, which also characterizes the size and other properties of interstellar dust grains, has been calculated for various regions of space by the extinction law extrapolation method. The results for the two classes of stars agree: the standard deviation of the "red giants minus OB" $R_V$ differences within 500 pc of the Sun is 0.2. The detected $R_V$ variations between 2.2 and 4.4 not only manifest themselves in individual clouds but also span the entire space near the Sun, following Galactic structures. In the Local Bubble within about 100 pc of the Sun, $R_V$ has a minimum. In the inner part of the Gould Belt and at high Galactic latitudes, at a distance of about 150 pc from the Sun, $R_V$ reaches a maximum and then decreases to its minimum in the outer part of the Belt and other directions at a distance of about 500 pc from the Sun, returning to its mean values far from the Sun. The detected maximum of $R_V$ at high Galactic latitudes is important when allowance is made for the interstellar extinction toward extragalactic objects. In addition, a monotonic increase in $R_V$ by 0.3 per kpc toward the Galactic center has been found near the Galactic equator. It is consistent with the result obtained by Zasowski et al. (2009) for much of the Galaxy. Ignoring the $R_V$ variations and traditionally using a single value for the entire space must lead to systematic errors in the calculated distances reaching 10\%.

### THz Time-Domain Spectroscopy of Mixed CO2-CH3OH Interstellar Ice Analogs

The icy mantles of interstellar dust grains are the birthplaces of the primordial prebiotic molecular inventory that may eventually seed nascent solar systems and the planets and planetesimals that form therein. Here, we present a study of two of the most abundant species in these ices after water: carbon dioxide (CO2) and methanol (CH3OH) using TeraHertz (THz) time-domain spectroscopy and mid-infrared spectroscopy. We study pure and mixed-ices of these species, and demonstrate the power of the THz region of the spectrum to elucidate the long-range structure (i.e. crystalline versus amorphous) of the ice, the degree of segregation of these species within the ice, and the thermal history of the species within the ice. Finally, we comment on the utility of the THz transitions arising from these ices for use in astronomical observations of interstellar ices.

### Temperature spectra of interstellar dust grains heated by cosmic-rays I: translucent clouds

Heating of whole interstellar dust grains by cosmic-ray (CR) particles affects the gas-grain chemistry in molecular clouds by promoting molecule desorption, diffusion, and chemical reactions on grain surfaces. The frequency of such heating $f_T$, s$^{-1}$, determines how often a certain temperature $T_{\rm CR}$, K, is reached for grains hit by CR particles. This study aims to provide astrochemists with comprehensive and updated dataset on the CR-induced whole-grain heating. We present calculations of $f_T$ and $T_{\rm CR}$ spectra for bare olivine grains with radius $a$ of 0.05; 0.1; 0.2 $\mu$m, and such grains covered with ice mantles of thickness 0.1$a$ and 0.3$a$. Grain shape and structure effects are considered, as well as 30 CR elemental constituents with an updated energy spectrum corresponding to a translucent cloud with $A_V=2$ mag. Energy deposition by CRs in grain material was calculated with the SRIM program. We report full $T_{\rm CR}$ spectra for all nine grain types and consider initial grain temperatures of 10 K and 20 K. We also provide frequencies for a range of minimum $T_{\rm CR}$ values. The calculated dataset can be simply and flexibly implemented in astrochemical models. The results show that, in the case of translucent clouds, the currently adopted rate for heating of whole grains to temperatures in excess of 70 K is underestimated by approximately two orders of magnitude in astrochemical numerical simulations. Additionally, grains are heated by CRs to modest temperatures (20--30 K) with intervals of a few years, which reduces the possibility of ice chemical explosions.

### On the chemistry of hydrides of N atoms and O$^+$ ions [Replacement]

Previous work by various authors has suggested that the detection by Herschel/HIFI of nitrogen hydrides along the low density lines of sight towards G10.6-0.4 (W31C) cannot be accounted for by gas-phase chemical models. In this paper we investigate the role of surface reactions on dust grains in diffuse regions, and we find that formation of the hydrides by surface reactions on dust grains with efficiency comparable to that for H$_2$ formation reconciles models with observations of nitrogen hydrides. However, similar surface reactions do not contribute significantly to the hydrides of O$^+$ ions detected by Herschel/HIFI present along many sight lines in the Galaxy. The O$^+$ hydrides can be accounted for by conventional gas-phase chemistry either in diffuse clouds of very low density with normal cosmic ray fluxes or in somewhat denser diffuse clouds with high cosmic ray fluxes. Hydride chemistry in dense dark clouds appears to be dominated by gas-phase ion-molecule reactions.

### On the chemistry of hydrides of N atoms and O$^+$ ions

Previous work by various authors has suggested that the detection by Herschel/HIFI of nitrogen hydrides along the low density lines of sight towards G10.6-0.4 (W31C) cannot be accounted for by gas-phase chemical models. In this paper we investigate the role of surface reactions on dust grains in diffuse regions, and we find that formation of the hydrides by surface reactions on dust grains with efficiency comparable to that for H$_2$ formation reconciles models with observations of nitrogen hydrides. However, similar surface reactions do not contribute significantly to the hydrides of O$^+$ ions detected by Herschel/HIFI present along many sight lines in the Galaxy. The O$^+$ hydrides can be accounted for by conventional gas-phase chemistry either in diffuse clouds of very low density with normal cosmic ray fluxes or in somewhat denser diffuse clouds with high cosmic ray fluxes. Hydride chemistry in dense dark clouds appears to be dominated by gas-phase ion-molecule reactions.

### The role of ices in star-forming clouds

Ices play a critical role during the evolution of interstellar clouds. Their presence is ubiquitous in the dense molecular medium and their impact is not only limited to chemistry. Species adsorbed onto dust grains also affect cloud thermodynamics. It all depends on the interstellar conditions, the chemical parameters, and the composition of ice layers. In this work, I study the formation of ices by focusing on the interplay between gas and solid phase to determine their role on cloud evolution and star formation. I show that while the formation of ices greatly impacts the cloud chemistry, their role on the thermodynamics is more conservative, and their influence on star formation is only marginal.

### Evidence for a CO desorption front in the outer AS 209 disk

Millimeter observations of CO isotopologues are often used to make inferences about protoplanetary disk gas density and temperature structures. The accuracy of these estimates depends on our understanding of CO freezeout and desorption from dust grains. Most models of these processes indicate that CO column density decreases monotonically with distance from the central star due to a decrease in gas density and freezeout beyond the CO snowline. We present ALMA Cycle 2 observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O $J=2-1$ emission that instead suggest CO enhancement in the outer disk of T Tauri star AS 209. Most notably, the C$^{18}$O emission consists of a central peak and a ring at a radius of $\sim1''$ (120 AU), well outside the expected CO snowline. We propose that the ring arises from the onset of CO desorption near the edge of the millimeter dust disk. CO desorption exterior to a CO snowline may occur via non-thermal processes involving cosmic rays or high-energy photons, or via a radial thermal inversion arising from dust migration.

### The JCMT Gould Belt Survey: Evidence for Dust Grain Evolution in Perseus Star-forming Clumps

The dust emissivity spectral index, $\beta$, is a critical parameter for deriving the mass and temperature of star-forming structures, and consequently their gravitational stability. The $\beta$ value is dependent on various dust grain properties, such as size, porosity, and surface composition, and is expected to vary as dust grains evolve. Here we present $\beta$, dust temperature, and optical depth maps of the star-forming clumps in the Perseus Molecular Cloud determined from fitting SEDs to combined Herschel and JCMT observations in the 160 $\mu$m, 250 $\mu$m, 350 $\mu$m, 500 $\mu$m, and 850 $\mu$m bands. Most of the derived $\beta$, and dust temperature values fall within the ranges of 1.0 - 2.7 and 8 - 20 K, respectively. In Perseus, we find the $\beta$ distribution differs significantly from clump to clump, indicative of grain growth. Furthermore, we also see significant, localized $\beta$ variations within individual clumps and find low $\beta$ regions correlate with local temperature peaks, hinting at the possible origins of low $\beta$ grains. Throughout Perseus, we also see indications of heating from B stars and embedded protostars, as well evidence of outflows shaping the local landscape.

### A Primer on Unifying Debris Disk Morphologies [Replacement]

A "minimum model" for debris disks consists of a narrow ring of parent bodies, secularly forced by a single planet on a possibly eccentric orbit, colliding to produce dust grains that are perturbed by stellar radiation pressure. We demonstrate how this minimum model can reproduce a wide variety of disk morphologies imaged in scattered starlight. Five broad categories of disk shape can be captured: "rings," "needles," "ships-and-wakes," "bars," and "moths (a.k.a. fans)," depending on the viewing geometry. Moths can also sport "double wings." We explain the origin of morphological features from first principles, exploring the dependence on planet eccentricity, disk inclination dispersion, and the parent body orbital phases at which dust grains are born. A key determinant in disk appearance is the degree to which dust grain orbits are apsidally aligned. Our study of a simple steady-state (secularly relaxed) disk should serve as a reference for more detailed models tailored to individual systems. We use the intuition gained from our guidebook of disk morphologies to interpret, informally, the images of a number of real-world debris disks.

### A Primer on Unifying Debris Disk Morphologies [Replacement]

A "minimum model" for debris disks consists of a narrow ring of parent bodies, secularly forced by a single planet on a possibly eccentric orbit, colliding to produce dust grains that are perturbed by stellar radiation pressure. We demonstrate how this minimum model can reproduce a wide variety of disk morphologies imaged in scattered starlight. Five broad categories of disk shape can be captured: "rings," "needles," "ships-and-wakes," "bars," and "moths (a.k.a. fans)," depending on the viewing geometry. Moths can also sport "double wings." We explain the origin of morphological features from first principles, exploring the dependence on planet eccentricity, disk inclination dispersion, and the parent body orbital phases at which dust grains are born. A key determinant in disk appearance is the degree to which dust grain orbits are apsidally aligned. Our study of a simple steady-state (secularly relaxed) disk should serve as a reference for more detailed models tailored to individual systems. We use the intuition gained from our guidebook of disk morphologies to interpret, informally, the images of a number of real-world debris disks. These interpretations suggest that the farthest reaches of planetary systems are perturbed by eccentric planets, possibly just a few Earth masses each.

### A Primer on Unifying Debris Disk Morphologies

A "minimum model" for debris disks consists of a narrow ring of parent bodies, secularly forced by a single planet on a possibly eccentric orbit, colliding to produce dust grains that are perturbed by stellar radiation pressure. We demonstrate how this minimum model can reproduce a wide variety of disk morphologies imaged in scattered starlight. Five broad categories of disk shape can be captured: "rings," "needles," "ships-and-wakes," "bars," and "moths (a.k.a. fans)," depending on the viewing geometry. Moths can also sport "double wings." We explain the origin of morphological features from first principles, exploring the dependence on planet eccentricity, disk inclination dispersion, and the parent body orbital phases at which dust grains are born. A key determinant in disk appearance is the degree to which dust grain orbits are apsidally aligned. Our study of a simple steady-state (secularly relaxed) disk should serve as a reference for more detailed models tailored to individual systems. We use the intuition gained from our guidebook of disk morphologies to interpret, informally, the images of a number of real-world debris disks.

### Detection of Phosphorus, Sulphur, and Zinc in the Carbon-Enhanced Metal-Poor Star BD+44 493

The carbon-enhanced metal-poor star BD+44 493 ([Fe/H]=-3.9) has been proposed as a candidate second-generation star enriched by metals from a single Pop III star. We report the first detections of P and S and the second detection of Zn in any extremely metal-poor carbon-enhanced star, using new spectra of BD+44 493 collected by the Cosmic Origins Spectrograph on the Hubble Space Telescope. We derive [P/Fe]=-0.34 +/- 0.21, [S/Fe]=+0.07 +/- 0.41, and [Zn/Fe]=-0.10 +/- 0.24. We increase by ten-fold the number of Si I lines detected in BD+44 493, yielding [Si/Fe]=+0.15 +/- 0.22. The solar [S/Fe] and [Zn/Fe] ratios exclude the hypothesis that the abundance pattern in BD+44 493 results from depletion of refractory elements onto dust grains. Comparison with zero-metallicity supernova models suggests that the stellar progenitor that enriched BD+44 493 was massive and ejected much less than 0.07 Msun of Ni-56, characteristic of a faint supernova.

### Significant enhancement of ${\rm H_2}$ formation in disk galaxies under strong ram pressure

We show, for the first time, that ${\rm H_2}$ formation on dust grains can be enhanced in disk galaxies under strong ram-pressure (RP). We numerically investigate how the time evolution, of ${\rm H}$ {\sc i} and ${\rm H_2}$ components in disk galaxies orbiting a group/cluster of galaxies, can be influenced by hydrodynamical interaction between the gaseous components of the galaxies and the hot intra-cluster medium (ICM). We find that compression of ${\rm H}$ {\sc i} caused by RP increases ${\rm H_2}$ formation in disk galaxies, before RP rapidly strips ${\rm H}$ {\sc i}, cutting off the fuel supply and causing a drop in ${\rm H_2}$ density. We also find that the level of this ${\rm H_2}$ formation enhancement in a disk galaxy under RP depends on the mass of its host cluster dark matter (DM) halo, initial positions and velocities of the disk galaxy, and disk inclination angle with respect to the orbital plane. We demonstrate that dust growth is a key factor in the evolution of the ${\rm H}$ {\sc i} and ${\rm H_2}$ mass in disk galaxies under strong RP. We discuss how the correlation between ${\rm H_2}$ fractions and surface gas densities of disk galaxies evolves with time in the galaxies under RP. We also discuss whether or not galaxy-wide star formation rates (SFRs) in cluster disk galaxies can be enhanced by RP if the SFRs depend on ${\rm H_2}$ densities.

### Observations and analysis of a curved jet in the coma of comet 67P/Churyumov-Gerasimenko

We analyze the physical properties and dynamical origin of a curved jet of comet 67P/Churyumov-Gerasimenko that was observed repeatedly in several nucleus rotations starting on May 30 and persisting until early August, 2015. We simulated the motion of dust grains ejected from the nucleus surface under the influence of the gravity and viscous drag effect of the expanding gas flow from the rotating nucleus. The formation of the curved jet is a combination of the size of the dust particles (~0.1-1 mm) and the location of the source region near the nucleus equator. This enhances the spiral feature of the collimated dust stream after the dust is accelerated to a terminal speed on the order of m/s.

### Polarimetric Detection of Exoplanets Transiting T- and L- Brown Dwarfs

While scattering of light by atoms and molecules yields large amount of polarization at the B-band of both T- and L-dwarfs, scattering by dust grains in cloudy atmosphere of L-dwarfs gives rise to significant polarization at the far-optical and infra-red wavelengths where these objects are much brighter. However, the observable disk averaged polarization should be zero if the clouds are uniformly distributed and the object is spherically symmetric. Therefore, in order to explain the observed large polarization of several L-dwarfs, rotation-induced oblateness or horizontally inhomogeneous cloud distribution in the atmosphere is invoked. On the other hand, when an extra-solar planet of Earth-size or larger transits the brown dwarf along the line of sight, the asymmetry induced during the transit gives rise to a net non-zero, time dependent polarization. Employing atmospheric models for a range of effective temperature and surface gravity appropriate for T- and L-dwarfs, I derive the time dependent polarization profiles of these objects during transit phase and estimate the peak amplitude of polarization that occurs during the inner contact points of the transit ingress/egress phase. It is found that peak polarization in the range of 0.2-1.0 % at I- and J-band may arise of cloudy L dwarfs occulted by Earth-size or larger exoplanets. Such an amount of polarization is higher than that can be produced by rotation-induced oblateness of even the rapidly rotating L-dwarfs. Hence, I suggest that time resolved imaging polarization should be a potential technique to detect transiting exoplanets around L-dwarfs.

### Polarimetric Detection of Exoplanets Transiting T- and L- Brown Dwarfs [Replacement]

While scattering of light by atoms and molecules yields large amount of polarization at the B-band of both T- and L-dwarfs, scattering by dust grains in cloudy atmosphere of L-dwarfs gives rise to significant polarization at the far-optical and infra-red wavelengths where these objects are much brighter. However, the observable disk averaged polarization should be zero if the clouds are uniformly distributed and the object is spherically symmetric. Therefore, in order to explain the observed large polarization of several L-dwarfs, rotation-induced oblateness or horizontally inhomogeneous cloud distribution in the atmosphere is invoked. On the other hand, when an extra-solar planet of Earth-size or larger transits the brown dwarf along the line of sight, the asymmetry induced during the transit gives rise to a net non-zero, time dependent polarization. Employing atmospheric models for a range of effective temperature and surface gravity appropriate for T- and L-dwarfs, I derive the time dependent polarization profiles of these objects during transit phase and estimate the peak amplitude of polarization that occurs during the inner contact points of the transit ingress/egress phase. It is found that peak polarization in the range of 0.2-1.0 % at I- and J-band may arise of cloudy L dwarfs occulted by Earth-size or larger exoplanets. Such an amount of polarization is higher than that can be produced by rotation-induced oblateness of even the rapidly rotating L-dwarfs. Hence, I suggest that time resolved imaging polarization should be a potential technique to detect transiting exoplanets around L-dwarfs.

### Surface chemistry in photodissociation regions

The presence of dust can strongly affect the chemical composition of the interstellar medium. We model the chemistry in photodissociation regions (PDRs) using both gas-phase and dust-phase chemical reactions. Our aim is to determine the chemical compositions of the interstellar medium (gas/dust/ice) in regions with distinct (molecular) gas densities that are exposed to radiation fields with different intensities. We have significantly improved the Meijerink PDR code by including 3050 new gas-phase chemical reactions and also by implementing surface chemistry. In particular, we have included 117 chemical reactions occurring on grain surfaces covering different processes, such as adsorption, thermal desorption, chemical desorption, two-body reactions, photo processes, and cosmic-ray processes on dust grains. We obtain abundances for different gas and solid species as a function of visual extinction, depending on the density and radiation field. We also analyse the rates of the formation of CO2 and H2O ices in different environments. In addition, we study how chemistry is affected by the presence/absence of ice mantles (bare dust or icy dust) and the impact of considering different desorption probabilities. The type of substrate (bare dust or icy dust) and the probability of desorption can significantly alter the chemistry occurring on grain surfaces, leading to differences of several orders of magnitude in the abundances of gas-phase species, such as CO, H2CO, and CH3OH. The type of substrate, together with the density and intensity of the radiation field, also determine the threshold extinction to form ices of CO2 and H2O. We also conclude that H2CO and CH3OH are mainly released into the gas phase of low, far-ultraviolet illuminated PDRs through chemical desorption upon two-body surface reactions, rather than through photodesorption.

### A chemical solver to compute molecule and grain abundances and non-ideal MHD resistivities in prestellar core collapse calculations [Replacement]

We develop a detailed chemical network relevant to the conditions characteristic of prestellar core collapse. We solve the system of time-dependent differential equations to calculate the equilibrium abundances of molecules and dust grains, with a size distribution given by size-bins for these latter. These abundances are used to compute the different non-ideal magneto-hydrodynamics resistivities (ambipolar, Ohmic and Hall), needed to carry out simulations of protostellar collapse. For the first time in this context, we take into account the evaporation of the grains, the thermal ionisation of Potassium, Sodium and Hydrogen at high temperature, and the thermionic emission of grains in the chemical network, and we explore the impact of various cosmic ray ionisation rates. All these processes significantly affect the non-ideal magneto-hydrodynamics resistivities, which will modify the dynamics of the collapse. Ambipolar diffusion and Hall effect dominate at low densities, up to n_H = 10^12 cm^-3, after which Ohmic diffusion takes over. We find that the time-scale needed to reach chemical equilibrium is always shorter than the typical dynamical (free fall) one. This allows us to build a large, multi-dimensional multi-species equilibrium abundance table over a large temperature, density and ionisation rate ranges. This table, which we make accessible to the community, is used during first and second prestellar core collapse calculations to compute the non-ideal magneto-hydrodynamics resistivities, yielding a consistent dynamical-chemical description of this process.

### A chemical solver to compute molecule and grain abundances and non-ideal MHD resistivities in prestellar core collapse calculations

We develop a detailed chemical network relevant to the conditions characteristic of prestellar core collapse. We solve the system of time-dependent differential equations to calculate the equilibrium abundances of molecules and dust grains, with a size distribution given by size-bins for these latter. These abundances are used to compute the different non-ideal magneto-hydrodynamics resistivities (ambipolar, Ohmic and Hall), needed to carry out simulations of protostellar collapse. For the first time in this context, we take into account the evaporation of the grains, the thermal ionisation of Potassium, Sodium and Hydrogen at high temperature, and the thermionic emission of grains in the chemical network, and we explore the impact of various cosmic ray ionisation rates. All these processes significantly affect the non-ideal magneto-hydrodynamics resistivities, which will modify the dynamics of the collapse. Ambipolar diffusion and Hall effect dominate at low densities, up to n_H = 10^12 cm^-3, after which Ohmic diffusion takes over. We find that the time-scale needed to reach chemical equilibrium is always shorter than the typical dynamical (free fall) one. This allows us to build a large, multi-dimensional multi-species equilibrium abundance table over a large temperature, density and ionisation rate ranges. This table, which we make accessible to the community, is used during first and second prestellar core collapse calculations to compute the non-ideal magneto-hydrodynamics resistivities, yielding a consistent dynamical-chemical description of this process.

### A chemical solver to compute molecule and grain abundances and non-ideal MHD resistivities in prestellar core collapse calculations [Replacement]

We develop a detailed chemical network relevant to the conditions characteristic of prestellar core collapse. We solve the system of time-dependent differential equations to calculate the equilibrium abundances of molecules and dust grains, with a size distribution given by size-bins for these latter. These abundances are used to compute the different non-ideal magneto-hydrodynamics resistivities (ambipolar, Ohmic and Hall), needed to carry out simulations of protostellar collapse. For the first time in this context, we take into account the evaporation of the grains, the thermal ionisation of Potassium, Sodium and Hydrogen at high temperature, and the thermionic emission of grains in the chemical network, and we explore the impact of various cosmic ray ionisation rates. All these processes significantly affect the non-ideal magneto-hydrodynamics resistivities, which will modify the dynamics of the collapse. Ambipolar diffusion and Hall effect dominate at low densities, up to n_H = 10^12 cm^-3, after which Ohmic diffusion takes over. We find that the time-scale needed to reach chemical equilibrium is always shorter than the typical dynamical (free fall) one. This allows us to build a large, multi-dimensional multi-species equilibrium abundance table over a large temperature, density and ionisation rate ranges. This table, which we make accessible to the community, is used during first and second prestellar core collapse calculations to compute the non-ideal magneto-hydrodynamics resistivities, yielding a consistent dynamical-chemical description of this process.

### A chemical solver to compute molecule and grain abundances and non-ideal MHD resistivities in prestellar core collapse calculations [Replacement]

We develop a detailed chemical network relevant to the conditions characteristic of prestellar core collapse. We solve the system of time-dependent differential equations to calculate the equilibrium abundances of molecules and dust grains, with a size distribution given by size-bins for these latter. These abundances are used to compute the different non-ideal magneto-hydrodynamics resistivities (ambipolar, Ohmic and Hall), needed to carry out simulations of protostellar collapse. For the first time in this context, we take into account the evaporation of the grains, the thermal ionisation of Potassium, Sodium and Hydrogen at high temperature, and the thermionic emission of grains in the chemical network, and we explore the impact of various cosmic ray ionisation rates. All these processes significantly affect the non-ideal magneto-hydrodynamics resistivities, which will modify the dynamics of the collapse. Ambipolar diffusion and Hall effect dominate at low densities, up to n_H = 10^12 cm^-3, after which Ohmic diffusion takes over. We find that the time-scale needed to reach chemical equilibrium is always shorter than the typical dynamical (free fall) one. This allows us to build a large, multi-dimensional multi-species equilibrium abundance table over a large temperature, density and ionisation rate ranges. This table, which we make accessible to the community, is used during first and second prestellar core collapse calculations to compute the non-ideal magneto-hydrodynamics resistivities, yielding a consistent dynamical-chemical description of this process.

### A chemical solver to compute molecule and grain abundances and non-ideal MHD resistivities in prestellar core collapse calculations [Replacement]

We develop a detailed chemical network relevant to the conditions characteristic of prestellar core collapse. We solve the system of time-dependent differential equations to calculate the equilibrium abundances of molecules and dust grains, with a size distribution given by size-bins for these latter. These abundances are used to compute the different non-ideal magneto-hydrodynamics resistivities (ambipolar, Ohmic and Hall), needed to carry out simulations of protostellar collapse. For the first time in this context, we take into account the evaporation of the grains, the thermal ionisation of Potassium, Sodium and Hydrogen at high temperature, and the thermionic emission of grains in the chemical network, and we explore the impact of various cosmic ray ionisation rates. All these processes significantly affect the non-ideal magneto-hydrodynamics resistivities, which will modify the dynamics of the collapse. Ambipolar diffusion and Hall effect dominate at low densities, up to n_H = 10^12 cm^-3, after which Ohmic diffusion takes over. We find that the time-scale needed to reach chemical equilibrium is always shorter than the typical dynamical (free fall) one. This allows us to build a large, multi-dimensional multi-species equilibrium abundance table over a large temperature, density and ionisation rate ranges. This table, which we make accessible to the community, is used during first and second prestellar core collapse calculations to compute the non-ideal magneto-hydrodynamics resistivities, yielding a consistent dynamical-chemical description of this process.

### A chemical solver to compute molecule and grain abundances and non-ideal MHD resistivities in prestellar core collapse calculations [Replacement]

We develop a detailed chemical network relevant to the conditions characteristic of prestellar core collapse. We solve the system of time-dependent differential equations to calculate the equilibrium abundances of molecules and dust grains, with a size distribution given by size-bins for these latter. These abundances are used to compute the different non-ideal magneto-hydrodynamics resistivities (ambipolar, Ohmic and Hall), needed to carry out simulations of protostellar collapse. For the first time in this context, we take into account the evaporation of the grains, the thermal ionisation of Potassium, Sodium and Hydrogen at high temperature, and the thermionic emission of grains in the chemical network, and we explore the impact of various cosmic ray ionisation rates. All these processes significantly affect the non-ideal magneto-hydrodynamics resistivities, which will modify the dynamics of the collapse. Ambipolar diffusion and Hall effect dominate at low densities, up to n_H = 10^12 cm^-3, after which Ohmic diffusion takes over. We find that the time-scale needed to reach chemical equilibrium is always shorter than the typical dynamical (free fall) one. This allows us to build a large, multi-dimensional multi-species equilibrium abundance table over a large temperature, density and ionisation rate ranges. This table, which we make accessible to the community, is used during first and second prestellar core collapse calculations to compute the non-ideal magneto-hydrodynamics resistivities, yielding a consistent dynamical-chemical description of this process.

### Probing the CO and methanol snow lines in young protostars. Results from the CALYPSO IRAM-PdBI survey

Context. "Snow lines", marking regions where abundant volatiles freeze out onto the surface of dust grains, play an important role for planet growth and bulk composition in protoplanetary disks. They can already be observed in the envelopes of the much younger, low-mass Class 0 protostars that are still in their early phase of heavy accretion. Aims. We aim at using the information on the sublimation regions of different kinds of ices to understand the chemistry of the envelope, its temperature and density structure, and the history of the accretion process. Methods. As part of the CALYPSO IRAM Large Program, we have obtained observations of C$^{18}$O, N$_2$H$^+$ and CH$_3$OH towards nearby Class 0 protostars with the IRAM Plateau de Bure interferometer at sub-arcsecond resolution. For four of these sources we have modeled the emission using a chemical code coupled with a radiative transfer module. Results. We observe an anti-correlation of C$^{18}$O and N$_2$H$^+$ in NGC 1333-IRAS4A, NGC 1333-IRAS4B, L1157, and L1448C, with N$_2$H$^+$ forming a ring around the centrally peaked C$^{18}$O emission due to N$_2$H$^+$ being chemically destroyed by CO. The emission regions of models and observations match for a CO binding energy of 1200 K, which is higher than the binding energy of pure CO ices ($\sim$855 K). Furthermore, we find very low CO abundances inside the snow lines in our sources, about an order of magnitude lower than the total CO abundance observed in the gas on large scales in molecular clouds before depletion sets in. Conclusions. The high CO binding energy may hint at CO being frozen out in a polar ice environment like amorphous water ice or in non-polar CO$_2$-rich ice. The low CO abundances are comparable to values found in protoplanetary disks, which may indicate an evolutionary scenario where these low values are already established in the protostellar phase. (Abbr. Version)

### Detecting Exomoons Around Self-luminous Giant Exoplanets Through Polarization

Many of the directly imaged self-luminous gas giant exoplanets have been found to have cloudy atmospheres. Scattering of the emergent thermal radiation from these planets by the dust grains in their atmospheres should locally give rise to significant linear polarization of the emitted radiation. However, the observable disk averaged polarization should be zero if the planet is spherically symmetric. Rotation-induced oblateness may yield a net non-zero disk averaged polarization if the planets have sufficiently high spin rotation velocity. On the other hand, when a large natural satellite or exomoon transits a planet with cloudy atmosphere along the line of sight, the asymmetry induced during the transit should give rise to a net non-zero, time resolved linear polarization signal. The peak amplitude of such time dependent polarization may be detectable even for slowly rotating exoplanets. Therefore, we suggest that large exomoons around directly imaged self-luminous exoplanets may be detectable through time resolved imaging polarimetry. Adopting detailed atmospheric models for several values of effective temperature and surface gravity which are appropriate for self-luminous exoplanets, we present the polarization profiles of these objects in the infrared during transit phase and estimate the peak amplitude of polarization that occurs during the inner contacts of the transit ingress/egress phase. The peak polarization is predicted to range between 0.1 and 0.3 % in the infrared.

### Extrasolar comets : the origin of dust in exozodiacal disks?

Comets have been invoked in numerous studies as a potentially important source of dust and gas around stars, but none has studied the thermo-physical evolution, out-gassing rate, and dust ejection of these objects in such stellar systems. We investigate the thermo-physical evolution of comets in exo-planetary systems in order to provide valuable theoretical data required to interpret observations of gas and dust. We use a quasi 3D model of cometary nucleus to study the thermo-physical evolution of comets evolving around a single star from 0.1 to 50 AU, whose homogeneous luminosity varies from 0.1 to 70 solar luminosities. This paper provides mass ejection, lifetimes, and the rate of dust and water gas mass productions for comets as a function of the distance to the star and stellar luminosity. Results show significant physical changes to comets at high stellar luminosities. The models are presented in such a manner that they can be readily applied to any planetary system. By considering the examples of the Solar System, Vega and HD 69830, we show that dust grains released from sublimating comets have the potential to create the observed (exo)zodiacal emission. We show that observations can be reproduced by 1 to 2 massive comets or by a large number of comets whose orbits approach close to the star. Our conclusions depend on the stellar luminosity and the uncertain lifetime of the dust grains. We find, as in previous studies, that exozodiacal dust disks can only survive if replenished by a population of typically sized comets renewed from a large and cold reservoir of cometary bodies beyond the water ice line. These comets could reach the inner regions of the planetary system following scattering by a (giant) planet.

### Dust Evolution and the Formation of Planetesimals [Replacement]

The solid content of circumstellar disks is inherited from the interstellar medium: dust particles of at most a micrometer in size. Protoplanetary disks are the environment where these dust grains need to grow at least 13 orders of magnitude in size. Our understanding of this growth process is far from complete, with different physics seemingly posing obstacles to this growth at various stages. Yet, the ubiquity of planets in our galaxy suggests that planet formation is a robust mechanism. This chapter focuses on the earliest stages of planet formation, the growth of small dust grains towards the gravitationally bound "planetesimals", the building blocks of planets. We will introduce some of the key physics involved in the growth processes and discuss how they are expected to shape the global behavior of the solid content of disks. We will consider possible pathways towards the formation of larger bodies and conclude by reviewing some of the recent observational advances in the field.

### Dust Evolution and the Formation of Planetesimals

The solid content of circumstellar disks is inherited from the interstellar medium: dust particles of at most a micrometer in size. Protoplanetary disks are the environment where these dust grains need to grow at least 13 orders of magnitude in size. Our understanding of this growth process is far from complete, with different physics seemingly posing obstacles to this growth at various stages. Yet, the ubiquity of planets in our galaxy suggests that planet formation is a robust mechanism. This chapter focuses on the earliest stages of planet formation, the growth of small dust grains towards the gravitationally bound "planetesimals", the building blocks of planets. We will introduce some of the key physics involved in the growth processes and discuss how they are expected to shape the global behavior of the solid content of disks. We will consider possible pathways towards the formation of larger bodies and conclude by reviewing some of the recent observational advances in the field.

### Circumstellar Debris Disks: Diagnosing the Unseen Perturber

The first indication of the presence of a circumstellar debris disk is usually the detection of excess infrared emission from the population of small dust grains orbiting the star. This dust is short-lived, requiring continual replenishment, and indicating that the disk must be excited by an unseen perturber. Previous theoretical studies have demonstrated that an eccentric planet orbiting interior to the disk will stir the larger bodies in the belt and produce dust via interparticle collisions. However, motivated by recent observations, we explore another possible mechanism for heating a debris disk: a stellar-mass perturber orbiting exterior to and inclined to the disk and exciting the disk particles' eccentricities and inclinations via the Kozai-Lidov mechanism. We explore the consequences of an exterior perturber on the evolution of a debris disk using secular analysis and collisional N-body simulations. We demonstrate that a Kozai-Lidov excited disk can generate a dust disk via collisions and we compare the results of the Kozai-Lidov excited disk with a simulated disk perturbed by an interior eccentric planet. Finally, we propose two observational tests of a dust disk that can distinguish whether the dust was produced by an exterior brown dwarf or stellar companion or an interior eccentric planet.

### Extinction Laws toward Stellar Sources within a Dusty Circimstellar Medium and Implications for Type Ia Supernovae

Many astronomical objects are surrounded by dusty environments. In such dusty objects, multiple scattering processes of photons by circumstellar (CS) dust grains can effectively alter extinction properties. In this paper, we systematically investigate effects of multiple scattering on extinction laws for steady-emission sources surrounded by the dusty CS medium, using a radiation transfer simulation based on the Monte Carlo technique. In particular, we focus on whether and how the extinction properties are affected by properties of CS dust grains, adopting various dust grain models. We {\bf confirm} that behaviors of the (effective) extinction laws are highly dependent on the properties of CS grains. Especially, the total-to-selective extinction ratio $R_{V}$, which characterizes the extinction law, can be either increased or decreased, compared to the case without multiple scattering. We find that the criterion for this behavior is given by a ratio of albedos in the $B$ and $V$ bands. We also find that either small silicate grains or polycyclic aromatic hydrocarbons (PAHs) are necessary for realizing a low value of $R_{V}$ as often measured toward Type Ia supernovae, if the multiple scattering by CS dust is responsible for their non-standard extinction laws. Using the derived relations between the properties of dust grains and the resulting effective extinction laws, we propose that the extinction laws toward dusty objects could be used to constrain the properties of dust grains in CS environments.

### First Detection of Galactic Latitude Dependence of Near-Infrared Diffuse Galactic Light from DIRBE Reanalysis

Observational study on near-infrared (IR) scattering properties of interstellar dust grains has been limited due to its faintness. Using all-sky maps obtained from Diffuse Infrared Background Experiment (DIRBE), we investigate the scattering property from diffuse Galactic light (DGL) measurements at 1.25, 2.2, and 3.5 {\mu}m in addition to our recent analyses of diffuse near-IR emission (Sano et al. 2015; Sano et al. 2016). As a result, we first find that the intensity ratios of near-IR DGL to 100 {\mu}m emission increase toward low Galactic latitudes at 1.25 and 2.2 {\mu}m. The derived latitude dependence can be reproduced by a scattered light model of interstellar dust with a large scattering asymmetry factor g = <cos{\theta}> of $0.8^{+0.2}_{-0.3}$ at 1.25 and 2.2 {\mu}m, assuming an infinite Galaxy disk as an illuminating source. The derived asymmetry factor is comparable to the values obtained in the optical, but several times larger than that expected from a recent dust model. Since possible latitude dependence of ultraviolet-excited dust emission at 1.25 and 2.2 {\mu}m would reduce the large asymmetry factor to the reasonable value, our result may indicate the first detection of such an additional emission component in the diffuse interstellar medium.

### Constraints on Planetesimal Collision Models in Debris Disks

Observations of debris disks offer a window into the physical and dynamical properties of planetesimals in extrasolar systems through the size distribution of dust grains. In particular, the millimeter spectral index of thermal dust emission encodes information on the grain size distribution. We have made new VLA observations of a sample of seven nearby debris disks at 9 mm, with 3" resolution and $\sim5$ $\mu$Jy/beam rms. We combine these with archival ATCA observations of eight additional debris disks observed at 7 mm, together with up-to-date observations of all disks at (sub)millimeter wavelengths from the literature to place tight constraints on the millimeter spectral indices and thus grain size distributions. The analysis gives a weighted mean for the slope of the power law grain size distribution, $n(a)\propto a^{-q}$, of $\langle q \rangle = 3.36\pm0.02$, with a possible trend of decreasing $q$ for later spectral type stars. We compare our results to a range of theoretical models of collisional cascades, from the standard self-similar, steady-state size distribution ($q=3.5$) to solutions that incorporate more realistic physics such as alternative velocity distributions and material strengths, the possibility of a cutoff at small dust sizes from radiation pressure, as well as results from detailed dynamical calculations of specific disks. Such effects can lead to size distributions consistent with the data, and plausibly the observed scatter in spectral indices. For the AU Mic system, the VLA observations show clear evidence of a highly variable stellar emission component; this stellar activity obviates the need to invoke the presence of an asteroid belt to explain the previously reported compact millimeter source in this system.

### Gaps, rings, and non-axisymmetric structures in protoplanetary disks - Emission from large grains

Dust grains with sizes around (sub)mm are expected to couple only weakly to the gas motion in regions beyond 10 au of circumstellar disks. In this work, we investigate the influence of the spatial distribution of such grains on the (sub)mm appearance of magnetized protoplanetary disks. We perform non-ideal global 3D magneto-hydrodynamic stratified disk simulations including particles of different sizes (50 micron to 1 cm), using a Lagrangian particle solver. We calculate the spatial dust temperature distribution, including the dynamically coupled submicron-sized dust grains, and derive ideal continuum re-emission maps of the disk through radiative transfer simulations. Finally, we investigate the feasibility to observe specific structures in the thermal re-emission maps with the Atacama Large Millimeter/submillimeter Array (ALMA). The pressure bump close to the outer edge of the dead-zone leads to particle trapping in ring structures. More specifically, vortices in the disk concentrate the dust and create an inhomogeneous distribution of the solid material in the azimuthal direction. The large-scale disk perturbations are preserved in the (sub)mm re-emission maps. The observable structures are very similar to those expected to result from planet-disk interaction. The larger dust particles increase the brightness contrast between gap and ring structures. We find that rings, gaps and the dust accumulation in the vortex could be traced with ALMA down to a scale of a few astronomical units in circumstellar disks located in nearby star-forming regions. Finally, we present a brief comparison of these structures with those recently found with ALMA in the young circumstellar disks of HL Tau and Oph IRS 48.

### Resolving the Planetesimal Belt of HR 8799 with ALMA

The star HR 8799 hosts one of the largest known debris discs and at least four giant planets. Previous observations have found evidence for a warm belt within the orbits of the planets, a cold planetesimal belt beyond their orbits and a halo of small grains. With the infrared data, it is hard to distinguish the planetesimal belt emission from that of the grains in the halo. With this in mind, the system has been observed with ALMA in band 6 (1.34 mm) using a compact array format. These observations allow the inner edge of the planetesimal belt to be resolved for the first time. A radial distribution of dust grains is fitted to the data using an MCMC method. The disc is best fit by a broad ring between $145^{+12}_{-12}$ AU and $429^{+37}_{-32}$ AU at an inclination of $40^{+5}_{-6}${\deg} and a position angle of $51^{+8}_{-8}${\deg}. A disc edge at ~145 AU is too far out to be explained simply by interactions with planet b, requiring either a more complicated dynamical history or an extra planet beyond the orbit of planet b.

### Magnetic field geometry of an unusual cometary cloud Gal 110-13

We carried out optical polarimetry of an isolated cloud, Gal 110-13, to map the plane-of-the-sky magnetic field geometry. The main aim of the study is to understand the most plausible mechanism responsible for the unusual cometary shape of the cloud in the context of its magnetic field geometry. When unpolarized starlight passes through the intervening interstellar dust grains that are aligned with their short axes parallel to the local magnetic field, it gets linearly polarized. The plane-of-the-sky magnetic field component can therefore be traced by doing polarization measurements of background stars projected on clouds. Because the light in the optical wavelength range is most efficiently polarized by the dust grains typically found in the outer layers of the molecular clouds, optical polarimetry enables us to trace the magnetic field geometry of the outer layers of the clouds. We made R-band polarization measurements of 207 stars in the direction of Gal 110-13. The distance of Gal 110-13 was determined as $\sim450\pm80$ pc using our polarization and 2MASS near-infrared data. The foreground interstellar contribution was removed from the observed polarization values by observing a number of stars located in the vicinity of Gal 110-13 which has Hipparcos parallax measurements. The plane-of-the-sky magnetic field lines are found to be well ordered and aligned with the elongated structure of Gal 110-13. Using structure function analysis, we estimated the strength of the plane-of-the-sky component of the magnetic field as $\sim25\mu$G. Based on our results and comparing them with those from simulations, we conclude that compression by the ionization fronts from 10 Lac is the most plausible cause of the comet-like morphology of Gal 110-13 and of the initiation of subsequent star formation.

### Ferromagnetism and particle collisions: applications to protoplanetary disks and the meteoritical record

The meteoritical record shows both iron partitioning and tungsten isotopic partitioning between matrix and chondrules. Tungsten is not abundant enough to have driven its own isotopic partitioning, but if tungsten were correlated with iron, then ferromagnetic interactions grains could help explain both observations. We derive a practical parameterization for the increase in particle-particle collision rates caused by mutually attracting particle magnetic dipole moments. While the appropriate magnetic parameters remain uncertain, we show that ambient magnetic fields in protoplanetary disks are expected to be strong enough to magnetize iron metal bearing dust grains sufficiently to drive large increases in their collision rates. Such increased collision rates between iron metal rich grains could help preserve primordial iron and W isotopic inhomogeneities; and would help explain why the meteoritical record shows their partitioning in the solar nebula. The importance of magnetic interactions for larger grains whose growth is balanced by fragmentation is less clear, and will require future laboratory or numerical studies.

### On dust entrainment in photoevaporative winds [Replacement]

We investigate dust entrainment by photoevaporative winds in protoplanetary discs using dusty smoothed particle hydrodrodynamics (SPH). We use unequal-mass particles to resolve more than five orders of magnitude in disc/outflow density and a one-fluid formulation to efficiently simulate an equivalent magnitude range in drag stopping time. We find that only micron sized dust grains and smaller can be entrained in EUV driven winds. The maximum grain size is set by dust settling in the disc rather than aerodynamic drag in the wind. More generally, there is a linear relationship between the base flow density and the maximum entrainable grain size in the wind. A pileup of micron sized dust grains can occur in the upper atmosphere at critical radii in the disc as grains decouple from the low-density wind. Entrainment is a strong function of location in the disc, resulting in a size sorting of grains in the outflow---the largest grain being carried out between $10$--$20\,$AU. The peak dust density for each grain size occurs at the inner edge of its own entrainment region.

### How dusty are photoevaporative winds?

Not very. We perform dusty smoothed particle hydrodrodynamic (SPH) calculations of photoevaporation in protoplanetary discs. We use unequal-mass particles to resolve more than five orders of magnitude in disc/outflow density and a one-fluid formulation to efficiently simulate an equivalent magnitude range in drag stopping time. We find that only micron sized dust grains and smaller can be entrained in photoevaporative winds. This makes it difficult to explain the dust holes seen in transition discs using photoevaporation and implies that only small grains can be transported to the outer disc by this mechanism. A pileup of micron sized dust grains can occur in the upper atmosphere at critical radii in the disc as grains decouple from the low-density wind. Entrainment is a strong function of location in the disc, resulting in a size sorting of grains in the outflow---the largest grain being carried out between $10$--$20\,$AU. The peak dust density for each grain occurs at the inner edge of its entrainment region.

### Direct measurement of desorption and diffusion energies of O and N atoms physisorbed on amorphous surfaces

Physisorbed atoms on the surface of interstellar dust grains play a central role in solid state astrochemistry. Their surface reactivity is one source of the observed molecular complexity in space. In experimental astrophysics, the high reactivity of atoms also constitutes an obstacle to measuring two of the fundamental properties in surface physics, namely desorption and diffusion energies, and so far direct measurements are non-existent for O and N atoms. We investigated the diffusion and desorption processes of O and N atoms on cold surfaces in order to give boundary conditions to astrochemical models. Here we propose a new technique for directly measuring the N- and O-atom mass signals. Including the experimental results in a simple model allows us to almost directly derive the desorption and diffusion barriers of N atoms on amorphous solid water ice (ASW) and O atoms on ASW and oxidized graphite. We find a strong constraint on the values of desorption and thermal diffusion energy barriers. The measured barriers for O atoms are consistent with recent independent estimations and prove to be much higher than previously believed (E$_{des}=1410_{-160}^{+290}$; E$_{dif}=990_{-360}^{+530}$ K on ASW). As for oxygen atoms, we propose that the combination E$_{des}$-E$_{dif}$=1320-750 K is a sensible choice among the possible pairs of solutions. Also, we managed to measure the desorption and diffusion energy of N atoms for the first time (E$_{des}=720_{-80}^{+160}$; E$_{dif}=525_{-200}^{+260}$ K on ASW) in the thermal hopping regime and propose that the combination E$_{des}$-E$_{dif}$=720-400 K can be reasonably adopted in models. The value of E$_{dif}$ for N atoms is slightly lower than previously suggested, which implies that the N chemistry on dust grains might be richer.

### Dust diffusion and settling in the presence of collisions: Trapping (sub)micron grains in the midplane

In protoplanetary disks, the distribution and abundance of small (sub)micron grains are important for a range of physical and chemical processes. For example, they dominate the optical depth at short wavelengths and their surfaces are the sites of many important chemical reactions such as the formation of water. Based on their aerodynamical properties (i.e., their strong dynamical coupling with the surrounding gas) it is often assumed that these small grains are well-mixed with the gas. Our goal is to study the vertical (re)distribution of grains taking into account settling, turbulent diffusion, as well as collisions with other dust grains. Assuming a fragmentation-limited background dust population, we developed a Monte Carlo approach that follows single monomers as they move through a vertical column of gas and become incorporated in different aggregates as they undergo sticking and fragmenting collisions. We find that (sub)micron grains are not necessarily well-mixed vertically, but can become trapped in a thin layer with a scale-height that is significantly smaller than that of the gas. This collisional trapping occurs when the timescale for diffusion is comparable to or longer than the collision timescale in the midplane and its effect is strongest when the most massive particles in the size-distribution show significant settling. Based on simulations and analytical considerations we conclude that for typical dust-to-gas ratios and turbulence levels, the collisional trapping of small grains should be a relatively common phenomenon. The absence of trapping could then indicate a low dust-to-gas ratio, possibly because a large portion of the dust mass has been removed through radial drift or is locked up in planetesimals.

### Formation, Evolution, and Revolution of Galaxies by SKA: Activities of SKA-Japan Galaxy Evolution Sub-SWG

Formation and evolution of galaxies have been a central driving force in the studies of galaxies and cosmology. Recent studies provided a global picture of cosmic star formation history. However, what drives the evolution of star formation activities in galaxies has long been a matter of debate. The key factor of the star formation is the transition of hydrogen from atomic to molecular state, since the star formation is associated with the molecular phase. This transition is also strongly coupled with chemical evolution, because dust grains, i.e., tiny solid particles of heavy elements, play a critical role in molecular formation. Therefore, a comprehensive understanding of neutral-molecular gas transition, star formation and chemical enrichment is necessary to clarify the galaxy formation and evolution. Here we present the activity of SKA-JP galaxy evolution sub-science working group (subSWG) Our activity is focused on three epochs: z \sim 0, 1, and z > 3. At z \sim 0, we try to construct a unified picture of atomic and molecular hydrogen through nearby galaxies in terms of metallicity and other various ISM properties. Up to intermediate redshifts z \sim 1, we explore scaling relations including gas and star formation properties, like the main sequence and the Kennicutt-Schmidt law of star forming galaxies. To connect the global studies with spatially-resolved investigations, such relations will be plausibly a viable way. For high redshift objects, the absorption lines of HI 21-cm line will be a very promising observable to explore the properties of gas in galaxies. By these studies, we will surely witness a real revolution in the studies of galaxies by SKA.

### Mid-infrared imaging- and spectro-polarimetric subarcsecond observations of NGC 1068

We present sub-arcsecond 7.5$-$13 $\mu$m imaging- and spectro-polarimetric observations of NGC 1068 using CanariCam on the 10.4-m Gran Telescopio CANARIAS. At all wavelengths, we find: (1) A 90 $\times$ 60 pc extended polarized feature in the northern ionization cone, with a uniform $\sim$44$^{\circ}$ polarization angle. Its polarization arises from dust and gas emission in the ionization cone, heated by the active nucleus and jet, and further extinguished by aligned dust grains in the host galaxy. The polarization spectrum of the jet-molecular cloud interaction at $\sim$24 pc from the core is highly polarized, and does not show a silicate feature, suggesting that the dust grains are different from those in the interstellar medium. (2) A southern polarized feature at $\sim$9.6 pc from the core. Its polarization arises from a dust emission component extinguished by a large concentration of dust in the galaxy disc. We cannot distinguish between dust emission from magnetically aligned dust grains directly heated by the jet close to the core, and aligned dust grains in the dusty obscuring material surrounding the central engine. Silicate-like grains reproduce the polarized dust emission in this feature, suggesting different dust compositions in both ionization cones. (3) An upper limit of polarization degree of 0.3 per cent in the core. Based on our polarization model, the expected polarization of the obscuring dusty material is $\lesssim$0.1 per cent in the 8$-$13 $\mu$m wavelength range. This low polarization may be arising from the passage of radiation through aligned dust grains in the shielded edges of the clumps.

### Influence of the water content in protoplanetary discs on planet migration and formation

The temperature and density profiles of protoplanetary discs depend crucially on the mass fraction of micrometre-sized dust grains and on their chemical composition. A larger abundance of micrometre-sized grains leads to an overall heating of the disc, so that the water ice line moves further away from the star. An increase in the water fraction inside the disc, maintaining a fixed dust abundance, increases the temperature in the icy regions of the disc and lowers the temperature in the inner regions. Discs with a larger silicate fraction have the opposite effect. Here we explore the consequence of the dust composition and abundance for the formation and migration of planets. We find that discs with low water content can only sustain outwards migration for planets up to 4 Earth masses, while outwards migration in discs with a larger water content persists up to 8 Earth masses in the late stages of the disc evolution. Icy planetary cores that do not reach run-away gas accretion can thus migrate to orbits close to the host star if the water abundance is low. Our results imply that hot and warm super-Earths found in exoplanet surveys could have formed beyond the ice line and thus contain a significant fraction in water. These water-rich super-Earths should orbit primarily around stars with a low oxygen abundance, where a low oxygen abundance is caused by either a low water-to-silicate ratio or by overall low metallicity.

### On the origins of polarization holes in Bok globules

Context. Polarimetric observations of Bok globules frequently show a decrease in the degree of polarization towards their central dense regions (polarization holes). This behaviour is usually explained with increased disalignment owing to high density and temperature, or insufficient angular resolution of a possibly complex magnetic field structure. Aims. We investigate whether a significant decrease in polarized emission of dense regions in Bok globules is possible under certain physical conditions. For instance, we evaluate the impact of optical depth effects and various properties of the dust phase. Methods. We use radiative transfer modelling to calculate the temperature structure of an analytical Bok globule model and simulate the polarized thermal emission of elongated dust grains. For the alignment of the dust grains, we consider a magnetic field and include radiative torque and internal alignment. Results. Besides the usual explanations, selected conditions of the temperature and density distribution, the dust phase and the magnetic field are also able to significantly decrease the polarized emission of dense regions in Bok globules. Taking submm/mm grains and typical column densities of existing Bok globules into consideration, the optical depth is high enough to decrease the degree of polarization by up to {\Delta}P~10%. If limited to the densest regions, dust grain growth to submm/mm size and accumulated graphite grains decrease the degree of polarization by up to {\Delta}P~10% and {\Delta}P~5%, respectively. However, the effect of the graphite grains occurs only if they do not align with the magnetic field.

### On the origins of polarization holes in Bok globules [Replacement]

Context. Polarimetric observations of Bok globules frequently show a decrease in the degree of polarization towards their central dense regions (polarization holes). This behaviour is usually explained with increased disalignment owing to high density and temperature, or insufficient angular resolution of a possibly complex magnetic field structure. Aims. We investigate whether a significant decrease in polarized emission of dense regions in Bok globules is possible under certain physical conditions. For instance, we evaluate the impact of optical depth effects and various properties of the dust phase. Methods. We use radiative transfer modelling to calculate the temperature structure of an analytical Bok globule model and simulate the polarized thermal emission of elongated dust grains. For the alignment of the dust grains, we consider a magnetic field and include radiative torque and internal alignment. Results. Besides the usual explanations, selected conditions of the temperature and density distribution, the dust phase and the magnetic field are also able to significantly decrease the polarized emission of dense regions in Bok globules. Taking submm/mm grains and typical column densities of existing Bok globules into consideration, the optical depth is high enough to decrease the degree of polarization by up to {\Delta}P~10%. If limited to the densest regions, dust grain growth to submm/mm size and accumulated graphite grains decrease the degree of polarization by up to {\Delta}P~10% and {\Delta}P~5%, respectively. However, the effect of the graphite grains occurs only if they do not align with the magnetic field.

### You need to log in to vote

The blog owner requires users to be logged in to be able to vote for this post.

Alternatively, if you do not have an account yet you can create one here.