Posts Tagged dust grains

Recent Postings from dust grains

New observations and models of circumstellar CO line emission of AGB stars in the Herschel SUCCESS programme

CONTEXT: Asymptotic giant branch (AGB) stars are in one of the latest evolutionary stages of low to intermediate-mass stars. Their vigorous mass loss has a significant effect on the stellar evolution, and is a significant source of heavy elements and dust grains for the interstellar medium. The mass-loss rate can be well traced by carbon monoxide (CO) line emission. AIMS: We present new \textit{Herschel} HIFI and IRAM 30m telescope CO line data for a sample of 53 galactic AGB stars. The lines cover a fairly large range of excitation energy from the $J=1\to0$ line to the $J=9\to8$ line, and even the $J=14\to13$ line in a few cases. We perform radiative transfer modelling for 38 of these sources to estimate their mass-loss rates. METHODS: We used a radiative transfer code based on the Monte Carlo method to model the CO line emission. We assume spherically symmetric circumstellar envelopes that are formed by a constant mass-loss rate through a smoothly accelerating wind. RESULTS: We find models that are consistent across a broad range of CO lines for most of the stars in our sample, i.e., a large number of the circumstellar envelopes can be described with a constant mass-loss rate. We also find that an accelerating wind is required to fit, in particular, the higher-J lines and that a velocity law will have a significant effect on the model line intensities. The results cover a wide range of mass-loss rates ($\sim 10^{-8}$ to $2\times 10^{-5}~\mathrm{M}_\odot~\mathrm{ yr}^{-1}$) and gas expansion velocities (2 to $21.5$~ km s$^{-1}$), and include M-, S-, and C-type AGB stars. Our results generally agree with those of earlier studies, although we tend to find slightly lower mass-loss rates by about 40\%, on average. We also present "bonus" lines detected during our CO observations.

Dust Cooling in Supernova Remnants in the Large Magellanic Cloud

The infrared-to-X-ray (IRX) flux ratio traces the relative importance of dust cooling to gas cooling in astrophysical plasma such as supernova remnants (SNRs). We derive IRX ratios of SNRs in the LMC using Spitzer and Chandra SNR survey data and compare them with those of Galactic SNRs. IRX ratios of all the SNRs in the sample are found to be moderately greater than unity, indicating that dust grains are a more efficient coolant than gas although gas cooling may not be negligible. The IRX ratios of the LMC SNRs are systematically lower than those of the Galactic SNRs. As both dust cooling and gas cooling pertain to the properties of the interstellar medium, the lower IRX ratios of the LMC SNRs may reflect the characteristics of the LMC, and the lower dust-to- gas ratio (a quarter of the Galactic value) is likely to be the most significant factor. The observed IRX ratios are compared with theoretical predictions that yield IRX ratios an order of magnitude larger. This discrepancy may originate from the dearth of dust in the remnants due to either the local variation of the dust abundance in the preshock medium with respect to the canonical abundance or the dust destruction in the postshock medium. The non-equilibrium ionization cooling of hot gas, in particular for young SNRs, may also cause the discrepancy. Finally, we discuss implications for the dominant cooling mechanism of SNRs in low-metallicity galaxies.

A SMACK Model of Colliding Planetesimals and Dust in the $\beta$ Pictoris Debris Disk: Thermal Radiation and Scattered Light

We present a new model of the $\beta$ Pictoris disk-and-planet system that simulates both the planetesimal collisions and the dynamics of the resulting dust grains, allowing us to model features and asymmetries in both thermal and scattered light images of the disk. Given the observed inclination and eccentricity of the $\beta$ Pictoris b planet, the model neatly ties together several features of the disk: the central hole in the submillimeter images, the two-disk "x"-pattern seen in scattered light, the "wing-tilt" asymmetry, and possibly even the clumpy gas seen by ALMA. We also find that most of the dust in the $\beta$ Pictoris system is likely produced outside the ring at 60-100 AU. Instead of a birth ring, this disk has a "stirring ring" at 60-100 AU where the high-velocity collisions produced by the secular wave launched by the planet are concentrated. The two-disk x-pattern arises because collisions occur more frequently at the peaks and troughs of the secular wave. The perturbations of the disk in this region create an azimuthally and vertically asymmetric spatial distribution of collisions, which could yield an azimuthal clump of gas without invoking resonances or an additional planet.

Benchmarking the Calculation of Stochastic Heating and Emissivity of Dust Grains in the Context of Radiative Transfer Simulations

We define an appropriate problem for benchmarking dust emissivity calculations in the context of radiative transfer (RT) simulations, specifically including the emission from stochastically heated dust grains. Our aim is to provide a self-contained guide for implementors of such functionality, and to offer insights in the effects of the various approximations and heuristics implemented by the participating codes to accelerate the calculations. The benchmark problem definition includes the optical and calorimetric material properties, and the grain size distributions, for a typical astronomical dust mixture with silicate, graphite and PAH components; a series of analytically defined radiation fields to which the dust population is to be exposed; and instructions for the desired output. We process this problem using six RT codes participating in this benchmark effort, and compare the results to a reference solution computed with the publicly available dust emission code DustEM. The participating codes implement different heuristics to keep the calculation time at an acceptable level. We study the effects of these mechanisms on the calculated solutions, and report on the level of (dis)agreement between the participating codes. For all but the most extreme input fields, we find agreement within 10% across the important wavelength range from 3 to 1000 micron. We conclude that the relevant modules in RT codes can and do produce fairly consistent results for the emissivity spectra of stochastically heated dust grains.

A stochastic model and Monte Carlo algorithm for fluctuation-induced H$_2$ formation on the surface of interstellar dust grains

A stochastic algorithm for simulation of fluctuation-induced kinetics of H$_2$ formation on grain surfaces is suggested as a generalization of the technique developed in our recent studies where this method was developed to describe the annihilation of spatially separate electrons and holes in a disordered semiconductor. The stochastic model is based on the spatially inhomogeneous, nonlinear integro-differential Smoluchowski equations with random source term. In this paper we derive the general system of Smoluchowski type equations for the formation of H$_2$ from two hydrogen atoms on the surface of interstellar dust grains with physisorption and chemisorption sites. We focus in this study on the spatial distribution, and numerically investigate the segregation in the case of a source with a continuous generation in time and randomly distributed in space. The stochastic particle method presented is based on a probabilistic interpretation of the underlying process as a stochastic Markov process of interacting particle system in discrete but randomly progressed time instances. The segregation is analyzed through the correlation analysis of the vector random field of concentrations which appears to be isotropic in space and stationary in time.

The accumulation and trapping of grains at planet gaps: effects of grain growth and fragmentation

We model the dust evolution in protoplanetary disks with full 3D, Smoothed Particle Hydrodynamics (SPH), two-phase (gas+dust) hydrodynamical simulations. The gas+dust dynamics, where aerodynamic drag leads to the vertical settling and radial migration of grains, is consistently treated. In a previous work, we characterized the spatial distribution of non-growing dust grains of different sizes in a disk containing a gap-opening planet and investigated the gap’s detectability with the Atacama Large Millimeter/submillimeter Array (ALMA). Here we take into account the effects of grain growth and fragmentation and study their impact on the distribution of solids in the disk. We show that rapid grain growth in the two accumulation zones around planet gaps is strongly affected by fragmentation. We discuss the consequences for ALMA observations.

Runaway growth of fractal dust grains

Fractal grains have large surface area, which leads to more efficient condensation. The special limit case where the volume-area ratio is constant (corresponding to, e.g., a very rough grain surface or non-compacts aggregates) is particularly interesting, as well as convenient, from a mathematical point of view. If dust grains from AGB stars have rough surfaces’, it may have important implications for our understanding of dust and wind formation in AGB stars.

Homogeneous dust emission and jet structure near active cometary nuclei: the case of 67P/Churyumov-Gerasimenko

We compute trajectories of dust grains starting from a homogeneous surface activity-profile. Despite the homogeneous initial distribution a collimation in jet-like structures becomes visible. The fine structure is caused by topographical features with similar bundles of normal vectors. Gravitational forces are accurately determined from the triangular surface mesh. For the comet 67P/Churyumov-Gerasimenko, we find several areas of good agreement between the homogeneous dust emission model and Rosetta observation of dust jets.

Tables of phase functions, opacities, albedos, equilibrium temperatures, and radiative accelerations of dust grains in exoplanets

There has been growing observational evidence for the presence of condensates in the atmospheres and/or comet-like tails of extrasolar planets. As a result, systematic and homogeneous tables of dust properties are useful in order to facilitate further observational and theoretical studies. In this paper we present calculations and analysis of non-isotropic phase functions, asymmetry parameter (mean cosine of the scattering angle), absorption and scattering opacities, single scattering albedos, equilibrium temperatures, and radiative accelerations of dust grains relevant for extrasolar planets. Our assumptions include spherical grain shape, Deirmendjian particle size distribution, and Mie theory. We consider several species: corundum/alumina, iron, olivines with 0% and 50% iron content, pyroxenes with 0%, 20% and 60% iron content, carbon at two different temperatures, water ice, liquid water, and ammonia. The presented tables cover the wavelength range of 0.2 to 500 micron and modal particle radii from 0.01 micron to 100 micron. Equilibrium temperatures and radiative accelerations assume irradiation by a non-black-body source of light with temperatures from 7000K to 700K seen at solid angles from 2$\pi$ to $10^{-6}$ sr. The tables are provided to the community together with a simple code which allows for an optional, finite, angular dimension of the source of light (star) in the phase function.

A compact concentration of large grains in the HD142527 protoplanetary dust trap

A pathway to the formation of planetesimals, and eventually giant planets, may occur in concentrations of dust grains trapped in pressure maxima. Dramatic crescent-shaped dust concentrations have been seen in recent radio images at sub-mm wavelengths. These disk asymmetries could represent the initial phases of planet formation in the dust trap scenario, provided that grain sizes are spatially segregated. A testable prediction of azimuthal dust trapping is that progressively larger grains should be more sharply confined and furthermore the trapped grains should follow a distribution that is markedly different from the gas. However, gas tracers such as CO and the infrared emission from small grains are both very optically thick where the submm continuum originates, so observations have been unable to test the trapping predictions or to identify compact concentrations of larger grains required for planet formation by core-accretion. Here we report multifrequency observations of HD142527, from 34GHz to 700GHz, that reveal a compact concentration of ~cm-sized grains, with a few Earth masses, embedded in a large-scale crescent of ~mm-sized particles. The emission peaks at wavelengths shorter than ~1mm are optically thick and trace the temperature structure resulting from shadows cast by the inner regions. Given this temperature structure, we infer that the largest dust grains are concentrated in the 34 GHz clump. We conclude that dust trapping is efficient for approximately cm-sized grains and leads to enhanced concentrations, while the smaller grains largely reflect the gas distribution.

Wind acceleration in AGB stars: Solid ground and loose ends

The winds of cool luminous AGB stars are commonly assumed to be driven by radiative acceleration of dust grains which form in the extended atmospheres produced by pulsation-induced shock waves. The dust particles gain momentum by absorption or scattering of stellar photons, and they drag along the surrounding gas particles through collisions, triggering an outflow. This scenario, here referred to as Pulsation-Enhanced Dust-DRiven Outflow (PEDDRO), has passed a range of critical observational tests as models have developed from empirical and qualitative to increasingly self-consistent and quantitative. A reliable theory of mass loss is an essential piece in the bigger picture of stellar and galactic chemical evolution, and central for determining the contribution of AGB stars to the dust budget of galaxies. In this review, I discuss the current understanding of wind acceleration and indicate areas where further efforts by theorists and observers are needed.

Compact dust concentration in the MWC 758 protoplanetary disk

The formation of planetesimals requires that primordial dust grains grow from micron- to km-sized bodies. Dust traps caused by gas pressure maxima have been proposed as regions where grains can concentrate and grow fast enough to form planetesimals, before radially migrating onto the star. We report new VLA Ka & Ku observations of the protoplanetary disk around the Herbig Ae/Be star MWC 758. The Ka image shows a compact emission region in the outer disk indicating a strong concentration of big dust grains. Tracing smaller grains, archival ALMA data in band 7 continuum shows extended disk emission with an intensity maximum to the north-west of the central star, which matches the VLA clump position. This segregation of grains sizes is expected in the context of dust trapping, where big grains are trapped more easily than smaller grains in gas pressure maxima. We develop a non-axisymmetric parametric model inspired by a steady state vortex solution which reproduces the observations, including the spectral energy distribution. Finally, we compare the radio continuum with SPHERE scattered light data. The ALMA continuum spatially coincides with a region devoid of scattered polarised emission and the VLA clump is offset to the north of the north-western spiral-like feature, indicating moderate or no flaring in the outer disk.

Far-ultraviolet study of the local supershell GSH 006-15+7

We have analyzed the archival data of FUV observations for the region of GSH 006-15+7, a large shell-like structure discovered by Moss et al. (2012) from the H I velocity maps. FUV emission is seen to be enhanced in the lower supershell region. The FUV emission is considered to come mainly from the scattering of interstellar photons by dust grains. A corresponding Monte Carlo simulation indicates that the distance to the supershell is 1300 +- 800 pc, which is similar to the previous estimation of 1500 +- 500 pc based on kinematic considerations. The spectrum at lower Galactic latitudes of the supershell exhibits molecular hydrogen fluorescence lines; a simulation model for this candidate photodissociation region (PDR) yields an H_2 column density of N(H_2) = 10^{18.0-20.0} cm^{-2} with a rather high total hydrogen density of n_H ~ 30 cm^{-3}.

Inhomogeneous cloud coverage through the Coulomb explosion of dust in substellar atmospheres

Recent observations of brown dwarf spectroscopic variability in the infrared infer the presence of patchy cloud cover. This paper proposes a mechanism for producing inhomogeneous cloud coverage due to the depletion of cloud particles through the Coulomb explosion of dust in atmospheric plasma regions. Charged dust grains Coulomb-explode when the electrostatic stress of the grain exceeds its mechanical tensile stress, which results in grains below a critical radius $a<a^{\rm Coul}_{\rm crit}$ being broken up. This work outlines the criteria required for the Coulomb explosion of dust clouds in substellar atmospheres, the effect on the dust particle size distribution function, and the resulting radiative properties of the atmospheric regions. Our results show that for an atmospheric plasma region with an electron temperature of $T_{e}=10$~eV ($\approx10^{5}$~K), the critical grain radius varies from $10^{-7}$ to $10^{-4}$~cm, depending on the grains’ tensile strength. Higher critical radii up to $10^{-3}$~cm are attainable for higher electron temperatures. We find that the process produces a bimodal particle size distribution composed of stable nanoscale seed particles and dust particles with $a\geq a^{\rm Coul}_{\rm crit}$, with the intervening particle sizes defining a region devoid of dust. As a result, the dust population is depleted, and the clouds become optically thin in the wavelength range $0.1-10~\mu$m, with a characteristic peak that shifts to higher wavelengths as more sub-micrometer particles are destroyed. In an atmosphere populated with a distribution of plasma volumes, this will yield regions of contrasting radiative properties, thereby giving a source of inhomogeneous cloud coverage. The results presented here may also be relevant for dust in supernova remnants and protoplanetary disks.

Dust Dynamics in Protoplanetary Disk Winds Driven by Magneto-Rotational Turbulence: A Mechanism for Floating Dust Grains with Characteristic Size

We investigate the dynamics of dust grains with various sizes in protoplanetary disk winds driven by magnetorotational turbulence, by simulating the time evolution of the dust grain distribution in the vertical direction. Small dust grains, which are well coupled to the gas, are dragged upward with the upflowing gas, while large grains remain near the midplane of a disk. Intermediate–size grains float at several scale heights from the midplane in time-averated force balance between the downward gravity and the upward gas drag. For the minimum mass solar nebula at 1 AU, dust grains with size of 20 — 40 $\mu m$ float at 5-10 scale heights from the midplane. Considering the dependence on the distance from the central star, smaller-size grains remain only in an outer region of the disk, while larger-size grains are distributed in a broader region. This implies that the dust depletion is expected to take place in small-to-large and inside-out manners. We also discuss the implication of our result to the observation of dusty material around young stellar objects.

A New Dust Budget In The Large Magellanic Cloud

The origin of dust in a galaxy is poorly understood. Recently, the surveys of the Large Magellanic Cloud (LMC) provide astrophysical laboratories for the dust studies. By a method of population synthesis, we investigate the contributions of dust produced by asymptotic giant branch (AGB) stars, common envelope (CE) ejecta and type II supernovae (SNe II) to the total dust budget in the LMC. Based on our models, the dust production rates (DPRs) of AGB stars in the LMC are between about $2.5\times10^{-5}$ and $4.0\times10^{-6}M_\odot{\rm yr^{-1}}$. The uncertainty mainly results from different models for the dust yields of AGB stars. The DPRs of CE ejecta are about $6.3\times10^{-6}$(The initial binary fraction is 50\%). These results are within the large scatter of several observational estimates. AGB stars mainly produce carbon grains, which is consistent with the observations. Most of dust grains manufactured by CE ejecta are silicate and iron grains. The contributions of SNe II are very uncertain. Compared with SNe II without reverse shock, the DPRs of AGB stars and CE ejecta are negligible. However, if only 2 \% of dust grains produced by SNe II can survive after reverse shock, the contributions of SNe II are very small. The total dust masses produced by AGB stars in the LMC are between $2.8\times10^4$ and $3.2\times10^5M_\odot$, and those produced by CE ejecta are about $6.3\times10^4$. They are much lower than the values estimated by observations. Therefore, there should be other dust sources in the LMC.

Survival and Structure of Dusty Vortices in Protoplanetary Discs

We have studied the impact of dust feedback on the survival and structure of vortices in protoplanetary discs using 2-D shearing box simulations with Lagrangian dust particles. We consider dust with a variety of sizes (stopping time $t_s = 10^{-2}\Omega^{-1} – 10^{2}\Omega^{-1}$, from fully coupled with the gas to the decoupling limit. We find that a vortex is destroyed by dust feedback when the total dust-to-gas mass ratio within the vortex is larger than 30-50%, independent of the dust size. The dust distribution can still be asymmetric in some cases after the vortex has been destroyed. With smaller amounts of dust, a vortex can survive for at least 100 orbits, and the maximum dust surface density within the vortex can be more than 100 times larger than the gas surface density, potentially facilitating planetesimal formation. On the other hand, in these stable vortices, small ($t_s < \Omega^{-1}$) and large ($t_s > \Omega^{-1}$) dust grains concentrate differently and affect the gas dynamics in different ways. The distribution of large dust is more elongated than that of small dust. Large dust ($t_s > \Omega^{-1}$) concentrates in the centre of the vortex and feedback leads to turn-over in vorticity towards the centre, forming a quiescent region within an anticyclonic vortex. Such a turn-over is absent if the vortex is loaded with small grains. We demonstrate that, in protoplanetary discs where both large and small dust grains are present and under the right condition, the concentration of large dust towards the vortex centre can lead to a quiescent centre, repelling the small dust and forming a small dust ring around the vortex centre. Such anticorrelations between small and large dust within vortices may explain the discrepancy between ALMA and near-IR scattered light observations in the asymmetric region of transitional discs.

Quantifying the gas inside dust cavities in transitional disks: implications for young planets

ALMA observations of a small sample of transitional disks with large dust cavities observed in Cycle 0 and 1 are summarized. The gas and dust surface density structures are inferred from the continuum and 12CO, 13CO and C18O line data using the DALI physical-chemical code. Thanks to its ability to self-shield, CO can survive inside dust cavities in spite of being exposed to intense UV radiation and can thus be used as a probe of the gas structure. Modeling of the existing data shows that gas is present inside the dust cavities in all cases, but at a reduced level compared with the gas surface density profile of the outer disk. The gas density decrease inside the dust cavity radius by factors of up to 10^4 suggests clearing by one or more planetary-mass companions. The accompanying pressure bumps naturally lead to trapping of the mm-sized dust grains observed in the ALMA images.

A fast and explicit algorithm for simulating the dynamics of small dust grains with smoothed particle hydrodynamics [Replacement]

We describe a simple method for simulating the dynamics of small grains in a dusty gas, relevant to micron-sized grains in the interstellar medium and grains of centimetre size and smaller in protoplanetary discs. The method involves solving one extra diffusion equation for the dust fraction in addition to the usual equations of hydrodynamics. This "diffusion approximation for dust" is valid when the dust stopping time is smaller than the computational timestep. We present a numerical implementation using Smoothed Particle Hydrodynamics (SPH) that is conservative, accurate and fast. It does not require any implicit timestepping and can be straightforwardly ported into existing 3D codes.

A fast and explicit algorithm for simulating the dynamics of small dust grains with smoothed particle hydrodynamics

We describe a simple method for simulating the dynamics of small grains in a dusty gas, relevant to micron-sized grains in the interstellar medium and grains of centimetre size and smaller in protoplanetary discs. The method involves solving one extra diffusion equation for the dust fraction in addition to the usual equations of hydrodynamics. This "diffusion approximation for dust" is valid when the dust stopping time is smaller than the computational timestep. We present a numerical implementation using Smoothed Particle Hydrodynamics (SPH) that is conservative, accurate and fast. It does not require any implicit timestepping and can be straightforwardly ported into existing 3D codes.

Dust trapping by spiral arms in gravitationally unstable protostellar discs

In this paper we discuss the influence of gravitational instabilities in massive protostellar discs on the dynamics of dust grains. Starting from a Smoothed Particle Hydrodynamics (SPH) simulation, we have computed the evolution of the dust in a quasi-static gas density structure typical of self-gravitating disc. For different grain size distributions we have investigated the capability of spiral arms to trap particles. We have run 3D radiative transfer simulations in order to construct maps of the expected emission at (sub-)millimetre and near-infrared wavelengths. Finally, we have simulated realistic observations of our disc models at (sub-)millimetre and near-infrared wavelengths as they may appear with the Atacama Large Millimetre/sub-millimetre Array (ALMA) and the High-Contrast Coronographic Imager for Adaptive Optics (HiCIAO) in order to investigate whether there are observational signatures of the spiral structure. We find that the pressure inhomogeites induced by gravitational instabilities produce a non-negligible dynamical effect on centimetre sized particles leading to significant overdensities in spiral arms. We also find that the spiral structure is readily detectable by ALMA over a wide range of (sub-)millimetre wavelengths and by HiCIAO in near-infrared scattered light for non-face-on discs located in the Ophiucus star-forming region. In addition, we find clear spatial spectral index variations across the disc, revealing that the dust trapping produces a migration of large grains that can be potentially investigated through multi-wavelenghts observations in the (sub-)millimetric. Therefore, the spiral arms observed to date in protoplanetary disc might be interpreted as density waves induced by the development of gravitational instabilities.

Millimeter-wave polarization of protoplanetary disks due to dust scattering

We present a new method to constrain the grain size in protoplanetary disks with polarization observations at millimeter wavelengths. If dust grains are grown to the size comparable to the wavelengths, the dust grains are expected to have a large scattering opacity and thus the continuum emission is expected to be polarized due to self-scattering. We perform 3D radiative transfer calculations to estimate the polarization degree for the protoplanetary disks having radial Gaussian-like dust surface density distributions, which have been recently discovered. The maximum grain size is set to be $100 {\rm~\mu m}$ and the observing wavelength to be 870 ${\rm \mu m}$. We find that the polarization degree is as high as 2.5% with a subarcsec spatial resolution, which is likely to be detected with near-future ALMA observations. The emission is polarized due to scattering of anisotropic continuum emission. The map of the polarization degree shows a double peaked distribution and the polarization vectors are in the radial direction in the inner ring and in the azimuthal direction in the outer ring. We also find the wavelength dependence of the polarization degree: the polarized emission is strongest if dust grains have a maximum size of $a_{\rm max}\sim\lambda/2\pi$, where $\lambda$ is the observing wavelength. Hence, multi-wave and spatially resolved polarization observations toward protoplanetary disks enable us to put a constraint on the grain size. The constraint on the grain size from polarization observations is independent of or may be even stronger than that from the opacity index.

Scattered Light from Dust in the Cavity of the V4046 Sgr Transition Disk

We report the presence of scattered light from dust grains located in the giant planet formation region of the circumbinary disk orbiting the ~20-Myr-old close (~0.045 AU separation) binary system V4046 Sgr AB based on observations with the new Gemini Planet Imager (GPI) instrument. These GPI images probe to within ~7 AU of the central binary with linear spatial resolution of ~3 AU, and are thereby capable of revealing dust disk structure within a region corresponding to the giant planets in our solar system. The GPI imaging reveals a relatively narrow (FWHM ~10 AU) ring of polarized near-infrared flux whose brightness peaks at ~14 AU. This ~14 AU radius ring is surrounded by a fainter outer halo of scattered light extending to ~45 AU, which coincides with previously detected mm-wave thermal dust emission. The presence of small grains that efficiently scatter starlight well inside the mm-wavelength disk cavity supports current models of planet formation that suggest planet-disk interactions can generate pressure traps that impose strong radial variations in the particle size distribution throughout the disk.

Impulsive Spot Heating and Thermal Explosion of Interstellar Grains Revisited

The problem of impulsive heating of dust grains in cold, dense interstellar clouds is revisited theoretically, with the aim to better understand leading mechanisms of the explosive desorption of icy mantles. It is rigorously shown that if the heating of a reactive medium occurs within a sufficiently localized spot (e.g., heating of mantles by cosmic rays), then the subsequent thermal evolution is characterized by a single dimensionless number $\lambda$. This number identifies a bifurcation between two distinct regimes: When $\lambda$ exceeds a critical value (threshold), the heat equation exhibits the explosive solution, i.e., the thermal (chemical) explosion is triggered. Otherwise, thermal diffusion causes the deposited heat to spread over the entire grain — this regime is commonly known as the whole-grain heating. The theory allows us to find a critical combination of the physical parameters that govern the explosion of icy mantles due to impulsive spot heating. In particular, the calculations suggest that heavy cosmic ray species (e.g., iron ions) colliding with dust are able to trigger the explosion. Based on the recently calculated local cosmic-ray spectra, the expected rate of the explosive desorption is estimated. The efficiency of the desorption, which affects all solid species independent of their binding energy, is shown to be comparable with other cosmic-ray desorption mechanisms typically considered in the literature. Also, the theory allows us to estimate maximum abundances of reactive species that may be stored in the mantles, which provides important constraints on available astrochemical models.

The silicate absorption profile in the ISM towards the heavily obscured nucleus of NGC 4418

The 9.7-micron silicate absorption profile in the interstellar medium provides important information on the physical and chemical composition of interstellar dust grains. Measurements in the Milky Way have shown that the profile in the diffuse interstellar medium is very similar to the amorphous silicate profiles found in circumstellar dust shells around late M stars, and narrower than the silicate profile in denser star-forming regions. Here, we investigate the silicate absorption profile towards the very heavily obscured nucleus of NGC 4418, the galaxy with the deepest known silicate absorption feature, and compare it to the profiles seen in the Milky Way. Comparison between the 8-13 micron spectrum obtained with TReCS on Gemini and the larger aperture spectrum obtained from the Spitzer archive indicates that the former isolates the nuclear emission, while Spitzer detects low surface brightness circumnuclear diffuse emission in addition. The silicate absorption profile towards the nucleus is very similar to that in the diffuse ISM in the Milky Way with no evidence of spectral structure from crystalline silicates or silicon carbide grains.

Variation of the ultraviolet extinction law across the Taurus-Auriga star forming complex. A GALEX based study

The Taurus-Auriga molecular complex (TMC) is the main laboratory for the study of low mass star formation. The density and properties of interstellar dust are expected to vary across the TMC. These variations trace important processes such as dust nucleation or the magnetic field coupling with the cloud. In this article, we show how the combination of near ultraviolet (NUV) and infrared (IR) photometry can be used to derive the strength of the 2175 \AA\ bump and thus any enhancement in the abundance of small dust grains and PAHs in the dust grains size distribution. This technique is applied to the envelope of the TMC, mapped by the GALEX All Sky Survey (AIS). UV and IR photometric data have been retrieved from the GALEX-AIS and the 2MASS catalogues. NUV and K-band star counts have been used to identify the areas in the cloud envelope where the 2175 \AA\ bump is weaker than in the diffuse ISM namely, the low column density extensions of L1495, L1498 and L1524 in Taurus, L1545, L1548, L1519, L1513 in Auriga and L1482-83 in the California region. This finding agrees with previous results on dust evolution derived from Spitzer data and suggests that dust grains begin to decouple from the environmental galactic magnetic field already in the envelope.

The Impact of Dust Evolution and Photoevaporation on Disk Dispersal

Protoplanetary disks are dispersed by viscous evolution and photoevaporation in a few million years; in the interim small, sub-micron sized dust grains must grow and form planets. The time-varying abundance of small grains in an evolving disk directly affects gas heating by far-ultraviolet photons, while dust evolution affects photoevaporation by changing the disk opacity and resulting penetration of FUV photons in the disk. Photoevaporative flows, in turn, selectively carry small dust grains leaving the larger particles—which decouple from the gas—behind in the disk. We study these effects by investigating the evolution of a disk subject to viscosity, photoevaporation by EUV, FUV and X-rays, dust evolution, and radial drift using a 1-D multi-fluid approach (gas + different dust grain sizes) to solve for the evolving surface density distributions. The 1-D evolution is augmented by 1+1D models constructed at each epoch to obtain the instantaneous disk structure and determine photoevaporation rates. The implementation of a dust coagulation/fragmentation model results in a marginal decrease in disk lifetimes when compared to models with no dust evolution; the disk lifetime is thus found to be relatively insensitive to the evolving dust opacity. We find that photoevaporation can cause significant reductions in the gas/dust mass ratio in the planet-forming regions of the disk as it evolves, and may result in a corresponding increase in heavy element abundances relative to hydrogen. We discuss implications for theories of planetesimal formation and giant planet formation, including the formation of gas-poor giants. After gas disk dispersal, $\sim 3\times 10^{-4}$ \ms\ of mass in solids typically remain, comparable to the solids inventory of our solar system.

Non-conservative evolution in Algols: where is the matter?

There is gathering indirect evidence suggesting non-conservative evolutions in Algols. However, the systemic mass-loss rate is poorly constrained by observations and generally set as a free parameter in binary-star evolution simulations. Moreover, systemic mass loss may lead to observational signatures that are still to be found. We investigate the impact of the outflowing gas and the possible presence of dust grains on the spectral energy distribution (SED). We used the 1D plasma code Cloudy and compared the results with the 3D Monte-Carlo radiative transfer code Skirt for dusty simulations. The circumbinary mass-distribution and binary parameters are computed with state-of-the-art binary calculations done with the Binstar evolution code. The outflowing material reduces the continuum flux-level of the stellar SED in the optical and UV. Due to the time-dependence of this effect, it may help to distinguish between different ejection mechanisms. Dust, if present, leads to observable infrared excesses even with low dust-to-gas ratios and traces the cold material at large distances from the star. By searching for such dust emission in the WISE catalogue, we found a small number of Algols showing infrared excesses, among which the two rather surprising objects SX Aur and CZ Vel. We find that some binary B[e] stars show the same strong Balmer continuum as we predict with our models. However, direct evidence of systemic mass loss is probably not observable in genuine Algols, since these systems no longer eject mass through the hotspot mechanism. Furthermore, owing to its high velocity, the outflowing material dissipates in a few hundred years. If hot enough, the hotspot may produce highly ionised species such as SiIV and observable characteristics that are typical of W Ser systems.

On the equivalent width of the Fe K$\alpha$ line produced by a dusty absorber in active galactic nuclei

Obscured AGNs provide an opportunity to study the material surrounding the central engine. Geometric and physical constraints on the absorber can be deduced from the reprocessed AGN emission. In particular, the obscuring gas may reprocess the nuclear X-ray emission producing a narrow Fe K$\alpha$ line and a Compton reflection hump. In recent years, models of the X-ray reflection from an obscuring torus have been computed; however, although the reflecting gas may be dusty, the models do not yet take into account the effects of dust on the predicted spectrum. We study this problem by analyzing two sets of models, with and without the presence of dust, using the one dimensional photo-ionization code Cloudy. The calculations are performed for a range of column densities ($22 <{\rm log}[N_H(\rm cm^{-2})]< 24.5$ ) and hydrogen densities ( $6 <{\rm log}[n_H(\rm cm^{-3})]< 8$). The calculations show the presence of dust can enhance the Fe K$\alpha$ equivalent width (EW) in the reflected spectrum by factors up to $\approx$ 8 for Compton thick (CT) gas and a typical ISM grain size distribution. The enhancement in EW with respect to the reflection continuum is due to the reduction in the reflected continuum intensity caused by the anisotropic scattering behaviour of dust grains. This effect will be most relevant for reflection from distant, predominately neutral gas, and is a possible explanation for AGNs which show a strong Fe K$\alpha$ EW and a relatively weak reflection continuum. Our results show it is an important to take into account dust while modeling the X-ray reflection spectrum, and that inferring a CT column density from an observed Fe K$\alpha$ EW may not always be valid. Multi-dimensional models are needed to fully explore the magnitude of the effect.

Numerical code for multi-component galaxies: from N-body to chemistry and magnetic fields [Cross-Listing]

We present a numerical code for multi-component simulation of the galactic evolution. Our code includes the following parts: $N$-body is used to evolve dark matter, stellar dynamics and dust grains, gas dynamics is based on TVD-MUSCL scheme with the extra modules for thermal processes, star formation, magnetic fields, chemical kinetics and multi-species advection. We describe our code in brief, but we give more details for the magneto-gas dynamics. We present several tests for our code and show that our code have passed the tests with a reasonable accuracy. Our code is parallelized using the MPI library. We apply our code to study the large scale dynamics of galactic discs.

The Pseudo-zodi Problem for Edge-on Planetary Systems

Future direct observations of extrasolar Earth-sized planets in the habitable zone could be hampered by a worrisome source of noise, starlight-reflecting exozodiacal dust. Mid-infrared surveys are currently underway to constrain the amount of exozodiacal dust in the habitable zones around nearby stars. However, at visible wavelengths another source of dust, invisible to these surveys, may dominate over exozodiacal dust. For systems observed near edge-on, a cloud of dust with face-on optical depth 10^-7 beyond ~5 AU can mimic the surface brightness of a cloud of exozodiacal dust with equal optical depth if the dust grains are sufficiently forward-scattering. We posit that dust migrating inward from cold debris belts via Poynting-Robertson drag could produce this "pseudo-zodiacal" effect, potentially making it ~50% as common as exozodiacal clouds. We place constraints on the disk radii and scattering phase function required to produce the effect.

Destruction of Interstellar Dust in Evolving Supernova Remnant Shock Waves

Supernova generated shock waves are responsible for most of the destruction of dust grains in the interstellar medium (ISM). Calculations of the dust destruction timescale have so far been carried out using plane parallel steady shocks, however that approximation breaks down when the destruction timescale becomes longer than that for the evolution of the supernova remnant (SNR) shock. In this paper we present new calculations of grain destruction in evolving, radiative SNRs. To facilitate comparison with the previous study by Jones et al. (1996), we adopt the same dust properties as in that paper. We find that the efficiencies of grain destruction are most divergent from those for a steady shock when the thermal history of a shocked gas parcel in the SNR differs significantly from that behind a steady shock. This occurs in shocks with velocities >~ 200 km/s for which the remnant is just beginning to go radiative. Assuming SNRs evolve in a warm phase dominated ISM, we find dust destruction timescales are increased by a factor of ~2 compared to those of Jones et al. (1996), who assumed a hot gas dominated ISM. Recent estimates of supernova rates and ISM mass lead to another factor of ~3 increase in the destruction timescales, resulting in a silicate grain destruction timescale of ~2-3 Gyr. These increases, while not able resolve the problem of the discrepant timescales for silicate grain destruction and creation, are an important step towards understanding the origin, and evolution of dust in the ISM.

Cold condensation of dust in the ISM

The condensation of complex silicates with pyroxene and olivine composition at conditions prevailing in molecular clouds has been experimentally studied. For this purpose, molecular species comprising refractory elements were forced to accrete on cold substrates representing the cold surfaces of surviving dust grains in the interstellar medium. The efficient formation of amorphous and homogeneous magnesium iron silicates at temperatures of about 12 K has been monitored by IR spectroscopy. The gaseous precursors of such condensation processes in the interstellar medium are formed by erosion of dust grains in supernova shock waves. In the laboratory, we have evaporated glassy silicate dust analogs and embedded the released species in neon ice matrices that have been studied spectroscopically to identify the molecular precursors of the condensing solid silicates. A sound coincidence between the 10 micron band of the interstellar silicates and the 10 micron band of the low-temperature siliceous condensates can be noted.

Study of grain alignment efficiency and distance estimate for CB4 cloud

We study the polarization efficiency (defined as ratio of polarization to extinction) of stars background to a small, nearly spherical and isolated Bok globule CB4 to understand the grain alignment process. A decrease in polarization efficiency with increase in visual extinction is noticed. This suggests that the observed polarization in lines of sight which intercept a Bok globule tends to show dominance of dust grains in the outer layers of the globule. This finding is consistent with the results obtained for other clouds in past. We determined the distance to the CB4 cloud using near-infrared photometry (2MASS $JHK_S$ colors) of moderately obscured stars located at the peripheries of the cloud. From the extinction-distance plot, the distance to this cloud is estimated to be ($459 \pm 85$) parsec.

Study of grain alignment efficiency and a distance estimate for small globule CB4 [Replacement]

We study the polarization efficiency (defined as the ratio of polarization to extinction) of stars in the background of the small, nearly spherical and isolated Bok globule CB4 to understand the grain alignment process. A decrease in polarization efficiency with an increase in visual extinction is noticed. This suggests that the observed polarization in lines of sight which intercept a Bok globule tends to show dominance of dust grains in the outer layers of the globule. This finding is consistent with the results obtained for other clouds in the past. We determined the distance to the cloud CB4 using near-infrared photometry (2MASS JHKS colors) of moderately obscured stars located at the periphery of the cloud. From the extinction-distance plot, the distance to this cloud is estimated to be ($459 \pm 85$) pc.

Dust Heating by Low-mass Stars in Massive Galaxies at z<1

Using the Hubble Space Telescope/Wide Field Camera 3 imaging data and multi-wavelength photometric catalog, we investigated the dust temperature of passively evolving and star-forming galaxies at 0.2<z<1.0 in the CANDELS fields. We estimated the stellar radiation field by low-mass stars from the stellar mass and surface brightness profile of these galaxies and then calculated their steady-state dust temperature. At first, we tested our method using nearby early-type galaxies with the deep FIR data by the Herschel Virgo cluster survey and confirmed that the estimated dust temperatures are consistent with the observed temperatures within the uncertainty. We then applied the method to galaxies at 0.2<z<1.0, and found that most of passively evolving galaxies with Mstar > 10^{10} Msun have a relatively high dust temperature of Tdust > 20 K, for which the formation efficiency of molecular hydrogen on the surface of dust grains in the diffuse ISM is expected to be very low from the laboratory experiments. The fraction of passively evolving galaxies strongly depends on the expected dust temperature at all redshifts and increases rapidly with increasing the temperature around Tdust ~ 20 K. These results suggest that the dust heating by low-mass stars in massive galaxies plays an important role for the continuation of their passive evolution, because the lack of the shielding effect of the molecular hydrogen on the UV radiation can prevent the gas cooling and formation of new stars.

Formation of hydroxylamine on dust grains via ammonia oxidation

The quest to detect prebiotic molecules in space, notably amino acids, requires an understanding of the chemistry involving nitrogen atoms. Hydroxylamine (NH$_2$OH) is considered a precursor to the amino acid glycine. Although not yet detected, NH$_2$OH is considered a likely target of detection with ALMA. We report on an experimental investigation of the formation of hydroxylamine on an amorphous silicate surface via the oxidation of ammonia. The experimental data are then fed into a simulation of the formation of NH$_2$OH in dense cloud conditions. On ices at 14 K and with a modest activation energy barrier, NH$_2$OH is found to be formed with an abundance that never falls below a factor 10 with respect to NH$_3$. Suggestions of conditions for future observations are provided.

Capture and evolution of dust in planetary mean-motion resonances: a fast, semi-analytic method for generating resonantly trapped disk images

Dust grains migrating under Poynting-Robertson drag may be trapped in mean-motion resonances with planets. Such resonantly trapped grains are observed in the solar system. In extrasolar systems, the exozodiacal light produced by dust grains is expected to be a major obstacle to future missions attempting to directly image terrestrial planets. The patterns made by resonantly trapped dust, however, can be used to infer the presence of planets, and the properties of those planets, if the capture and evolution of the grains can be modelled. This has been done with N-body methods, but such methods are computationally expensive, limiting their usefulness when considering large, slowly evolving grains, and for extrasolar systems with unknown planets and parent bodies, where the possible parameter space for investigation is large. In this work, we present a semi-analytic method for calculating the capture and evolution of dust grains in resonance, which can be orders of magnitude faster than N-body methods. We calibrate the model against N-body simulations, finding excellent agreement for Earth to Neptune mass planets, for a variety of grain sizes, initial eccentricities, and initial semimajor axes. We then apply the model to observations of dust resonantly trapped by the Earth. We find that resonantly trapped, asteroidally produced grains naturally produce the trailing blob’ structure in the zodiacal cloud, while to match the intensity of the blob, most of the cloud must be composed of cometary grains, which owing to their high eccentricity are not captured, but produce a smooth disk.

The relationship between polycyclic aromatic hydrocarbon emission and far-infrared dust emission from NGC 2403 and M83

We examine the relation between polycyclic aromatic hydrocarbon (PAH) emission at 8 microns and far-infrared emission from hot dust grains at 24 microns and from large dust grains at 160 and 250 microns in the nearby spiral galaxies NGC 2403 and M83 using data from the Spitzer Space Telescope and Herschel Space Observatory. We find that the PAH emission in NGC 2403 is better correlated with emission at 250 microns from dust heated by the diffuse interstellar radiation field (ISRF) and that the 8/250 micron surface brightness ratio is well-correlated with the stellar surface brightness as measured at 3.6 microns. This implies that the PAHs in NGC 2403 are intermixed with cold large dust grains in the diffuse interstellar medium (ISM) and that the PAHs are excited by the diffuse ISRF. In M83, the PAH emission appears more strongly correlated with 160 micron emission originating from large dust grains heated by star forming regions. However, the PAH emission in M83 is low where the 24 micron emission peaks within star forming regions, and enhancements in the 8/160 micron surface brightness ratios appear offset relative to the dust and the star forming regions within the spiral arms. This suggests that the PAHs observed in the 8 micron band are not excited locally within star forming regions but either by light escaping non-axisymmetrically from star forming regions or locally by young, non-photoionising stars that have migrated downstream from the spiral density waves. The results from just these two galaxies show that PAHs may be excited by different stellar populations in different spiral galaxies.

Five steps in the evolution from protoplanetary to debris disk

The protoplanetary disks of Herbig Ae stars eventually dissipate leaving a tenuous debris disk comprised of planetesimals and dust, as well as possibly gas and planets. This paper uses the properties of 10-20Myr A star debris disks to consider the protoplanetary to debris disk transition. The physical distinction between these two classes is argued to rest on the presence of primordial gas in sufficient quantities to dominate the motion of small dust grains (not the secondary nature of the dust or its level of stirring). This motivates an observational classification based on the dust spectrum, empirically defined so that A star debris disks require fractional excesses <3 at 12um and <2000 at 70um. We also propose a hypothesis to test, that the main sequence planet/planetesimal structures are already in place (but obscured) during the protoplanetary disk phase. This may be only weakly true if planetary architectures change until frozen during disk dispersal, or completely false if planets and planetesimals form during disk dispersal. Five steps in the transition are discussed: (i) carving an inner hole to form a transition disk; (ii) depletion of mm-sized dust in outer disk, noting the importance of determining whether this mass ends up in planetesimals or is collisionally depleted; (iii) final clearing of inner regions, noting that many mechanisms replenish moderate hot dust levels at later phases, and likely also operate in protoplanetary disks; (iv) disappearence of gas, noting recent discoveries of primordial and secondary gas in debris disks that highlight our ignorance and its impending enlightenment by ALMA; (v) formation of ring-like planetesimal structures, noting these are shaped by interactions with planets, and that the location of planetesimals in protoplanetary disks may be unrelated to the dust concentrations therein that are set by gas interactions.

Condensation of dust in the ejecta of type II-P supernovae

Aims: We study the production of dust in Type II-P supernova by coupling the gas-phase chemistry to the dust nucleation and condensation phases. We consider two supernova progenitor masses with homogeneous and clumpy ejecta to assess the chemical type and quantity of dust that forms. Grain size distributions are derived as a function of post-explosion time. Methods: The chemistry of the gas phase and the simultaneous formation of dust clusters are described by a chemical network. The formation of key species (CO, SiO) and dust clusters of silicates, alumina, silica, metal carbides and sulphides, pure metals, and amorphous carbon is considered. The master equations describing the chemistry of the nucleation phase are coupled to a dust condensation formalism based on Brownian coagulation. Results: Type II-P supernovae produce dust grains of various chemical compositions and size distributions as a function of time. The grain size distributions gain in complexity with time, are slewed towards large grains, and differ from the usual MRN power-law distribution used for interstellar dust. Gas density enhancements in the form of clumps strongly affect the dust chemical composition and the grain size distributions. Silicates and pure metallic grains are highly dependent on clumpiness. Specifically, clumpy ejecta produce grains over 0.1 micron, and the final dust mass reaches 0.14 Msun. Conversely, carbon and alumina dust masses are controlled by the mass yields of alumina and carbon in the zones where the dust is produced. Several dust components form in the ejecta and the total dust mass gradually builds up over a time span of 3 to 5 years post-outburst. This gradual growth provides a possible explanation for the discrepancy between the small dust masses formed at early post-explosion times and the high dust masses derived from recent observations of supernova remnants.

Physical Processes in the Interstellar Medium

Interstellar space is filled with a dilute mixture of charged particles, atoms, molecules and dust grains, called the interstellar medium (ISM). Understanding its physical properties and dynamical behavior is of pivotal importance to many areas of astronomy and astrophysics. Galaxy formation and evolution, the formation of stars, cosmic nucleosynthesis, the origin of large complex, prebiotic molecules and the abundance, structure and growth of dust grains which constitute the fundamental building blocks of planets, all these processes are intimately coupled to the physics of the interstellar medium. However, despite its importance, its structure and evolution is still not fully understood. Observations reveal that the interstellar medium is highly turbulent, consists of different chemical phases, and is characterized by complex structure on all resolvable spatial and temporal scales. Our current numerical and theoretical models describe it as a strongly coupled system that is far from equilibrium and where the different components are intricately linked together by complex feedback loops. Describing the interstellar medium is truly a multi-scale and multi-physics problem. In these lecture notes we introduce the microphysics necessary to better understand the interstellar medium. We review the relations between large-scale and small-scale dynamics, we consider turbulence as one of the key drivers of galactic evolution, and we review the physical processes that lead to the formation of dense molecular clouds and that govern stellar birth in their interior.

Studies of Anomalous Microwave Emission (AME) with the SKA

In this chapter, we will outline the scientific motivation for studying Anomalous Microwave Emission (AME) with the SKA. AME is thought to be due to electric dipole radiation from small spinning dust grains, although thermal fluctuations of magnetic dust grains may also contribute. Studies of this mysterious component would shed light on the emission mechanism, which then opens up a new window onto the interstellar medium (ISM). AME is emitted mostly in the frequency range $\sim 10$–100\,GHz, and thus the SKA has the potential of measuring the low frequency side of the AME spectrum, particularly in band 5. Science targets include dense molecular clouds in the Milky Way, as well as extragalactic sources. We also discuss the possibility of detecting rotational line emission from Poly-cyclic Aromatic Hydrocarbons (PAHs), which could be the main carriers of AME. Detecting PAH lines of a given spacing would allow for a definitive identification of specific PAH species.

AGB stars in the LMC: evolution of dust in circumstellar envelopes

We calculated theoretical evolutionary sequences of asymptotic giant branch (AGB) stars, including formation and evolution of dust grains in their circumstellar envelope. By considering stellar populations of the Large Magellanic Cloud (LMC), we calculate synthetic colour-colour and colour-magnitude diagrams, which are compared with those obtained by the Spitzer Space Telescope. The comparison between observations and theoretical predictions outlines that extremely obscured carbon-stars and oxygen-rich sources experiencing hot bottom burning (HBB) occupy well defined, distinct regions in the colour-colour ($[3.6]-[4.5]$, $[5.8]-[8.0]$) diagram. The C-rich stars are distributed along a diagonal strip that we interpret as an evolutionary sequence, becoming progressively more obscured as the stellar surface layers enrich in carbon. Their circumstellar envelopes host solid carbon dust grains with size in the range $0.05 < a < 0.2 \mu m$. The presence of SiC particles is expected only in the more metal-rich stars. The reddest sources, with $[3.6]-[4.5] > 2$, are the descendants of stars with initial mass $M_{in} \sim 2.5 – 3 M_{\odot}$ in the very latest phases of the AGB life. The oxygen-rich stars with the reddest colours ($[5.8]-[8.0] > 0.6$) are those experiencing HBB, the descendants of $\sim 5 M_{\odot}$ objects formed $10^{8}$ yr ago; alumina and silicates dust start forming at different distances from the central star. The overall dust production rate in the LMC is $\sim 4.5 \times 10^{-5} M_{\odot}/yr$, the relative percentages due to C- and M- star being respectively 85$%$ and 15 $%$.

On origin and destruction of relativistic dust and its implication for ultrahigh energy cosmic rays

Dust grains may be accelerated to relativistic speeds by radiation pressure of luminous sources, diffusive shocks, and other acceleration mechanisms. Such relativistic grains have been suggested as potential primary particles of ultrahigh energy cosmic rays (UHECRs). In this paper, we reexamine this idea by studying in detail different destruction mechanisms for relativistic grains moving with Lorentz factor $\gamma$ through a variety of environment conditions. For the solar radiation field, we find that sublimation/melting is a dominant destruction mechanism for silicate grains and large graphite grains. Using an improved treatment of photoelectric emission, we calculate the closest distance that relativistic grains can approach the Sun before destroyed by Coulomb explosions. A range of survival parameters for relativistic grains (size $a$ and $\gamma$) against both sublimation and Coulomb explosions by the solar radiation field is identified. We also study collisional destruction mechanisms, consisting of electronic sputtering by ions and grain-grain collisions. Electronic sputtering by light ions is found to be rather inefficient, whereas the evaporation induced by grain-grain collisions is shown to be an important mechanism for which the $a\le 1\mu$m grains can be completely destroyed after sweeping a column of gas and dust $N_{coll}\le 4\times 10^{20}$cm$^{-2}$. The destruction of relativistic dust by the interstellar radiation field (ISRF) and cosmic microwave background (CMB) in the intergalactic medium by melting is inefficient, while Coulomb explosions are only important for grains of very large $\gamma$. The obtained results indicate that relativistic dust grains from extragalactic sources would likely be destroyed in the interstellar medium, but the grains accelerated to relativistic speeds in our Galaxy are not completely ruled out as primary particles of UHECRs.

On origin and destruction of relativistic dust and its implication for ultrahigh energy cosmic rays [Replacement]

Dust grains may be accelerated to relativistic speeds by radiation pressure, diffusive shocks, and other acceleration mechanisms. Such relativistic grains have been suggested as primary particles of ultrahigh energy cosmic rays (UHECRs). In this paper, we first revisit the problem of acceleration by radiation pressure and calculate maximum grain velocities achieved. We find that grains can be accelerated to relativistic speeds with Lorentz factor $\gamma < 2$ by powerful radiation sources, which is lower than earlier estimates showing that $\gamma$ could reach $\sim 10$. We then investigate different destruction mechanisms for relativistic grains traversing a variety of environments. In solar radiation, we find that the destruction by thermal sublimation and Coulomb explosions is important. We also quantify grain destruction due to electronic sputtering by ions and grain-grain collisions. Electronic sputtering is found to be rather inefficient, whereas the evaporation following grain-grain collisions is shown to be an important mechanism for which the $a \sim 0.01- 1\mu$m grains would be destroyed after sweeping a gas column $N_{Coul}\sim 5\times 10^{19}-5\times 10^{20}cm^{-2}$. Relativistic dust in the interstellar medium and intergalactic medium (IGM) would be disrupted by Coulomb explosions due to collisional charging after traversing a gas column $N_{Coul} \sim 10^{17}cm^{-2}$ unless grain material is very strong. We show that photoelectric emission by optical and ultraviolet background radiation is also significant for the destruction of relativistic dust in the IGM. The obtained results indicate that relativistic dust from galaxies would be destroyed before reaching the Earth’s atmosphere and unlikely to account for UHECRs.

Depletion of chlorine into HCl ice in a protostellar core

The freezeout of gas-phase species onto cold dust grains can drastically alter the chemistry and the heating-cooling balance of protostellar material. In contrast to well-known species such as carbon monoxide (CO), the freezeout of various carriers of elements with abundances $<10^{-5}$ has not yet been well studied. Our aim here is to study the depletion of chlorine in the protostellar core, OMC-2 FIR 4. We observed transitions of HCl and H2Cl+ towards OMC-2 FIR 4 using the Herschel Space Observatory and Caltech Submillimeter Observatory facilities. Our analysis makes use of state of the art chlorine gas-grain chemical models and newly calculated HCl-H$_{2}$ hyperfine collisional excitation rate coefficients. A narrow emission component in the HCl lines traces the extended envelope, and a broad one traces a more compact central region. The gas-phase HCl abundance in FIR 4 is 9e-11, a factor of only 0.001 that of volatile elemental chlorine. The H2Cl+ lines are detected in absorption and trace a tenuous foreground cloud, where we find no depletion of volatile chlorine. Gas-phase HCl is the tip of the chlorine iceberg in protostellar cores. Using a gas-grain chemical model, we show that the hydrogenation of atomic chlorine on grain surfaces in the dark cloud stage sequesters at least 90% of the volatile chlorine into HCl ice, where it remains in the protostellar stage. About 10% of chlorine is in gaseous atomic form. Gas-phase HCl is a minor, but diagnostically key reservoir, with an abundance of <1e-10 in most of the protostellar core. We find the 35Cl/37Cl ratio in OMC-2 FIR 4 to be 3.2\pm0.1, consistent with the solar system value.

Outward Motion of Porous Dust Aggregates by Stellar Radiation Pressure in Protoplanetary Disks

We study the dust motion at the surface layer of protoplanetary disks. Dust grains in surface layer migrate outward due to angular momentum transport via gas-drag force induced by the stellar radiation pressure. In this study, we calculate mass flux of the outward motion of compact grains and porous dust aggregates by the radiation pressure. The radiation pressure force for porous dust aggregates is calculated using the T-Matrix Method for the Clusters of Spheres. First, we confirm that porous dust aggregates are forced by strong radiation pressure even if they grow to be larger aggregates in contrast to homogeneous and spherical compact grains to which efficiency of radiation pressure becomes lower when their sizes increase. In addition, we find that the outward mass flux of porous dust aggregates with monomer size of 0.1 $\mu$m is larger than that of compact grains by an order of magnitude at the disk radius of 1 AU, when their sizes are several microns. This implies that large compact grains like calcium-aluminum rich inclusions (CAIs) are hardly transported to outer region by stellar radiation pressure, whereas porous dust aggregates like chondritic-porous interplanetary dust particles (CP-IDPs) are efficiently transported to comet formation region. Crystalline silicates are possibly transported in porous dust aggregates by stellar radiation pressure from inner hot region to outer cold cometary region in the protosolar nebula.