Observable effects of dust formation in dynamic atmospheres of M-type Mira variables

The formation of dust with temperature-dependent non-grey opacity is considered in a series of self-consistent model atmospheres at different phases of an O-rich Mira variable of mass  1.2 M_⊙  . Photometric and interferometric properties of these models are predicted under different physical assumptions regarding the dust formation. The iron content of the initial silicate that forms and the availability of grain nuclei are found to be critical parameters that affect the observable properties. For certain plausible parameter values, dust would form at 2–3 times the average continuum photospheric radius. This work provides a consistent physical explanation for the larger apparent size of Mira variables at wavelengths shorter than 1 μm than that predicted by dust-free fundamental-mode pulsation models.     

Modelling the Extinction Properties of Galaxies

Recently (Granato, Lacey, Silva et al., 2000, astro-ph/0001308)we have combined our spectrophotometric galaxy evolution code which includes dust reprocessing (GRASIL, Silva et al., 1998) with semi-analytical galaxy formation models (GALFORM, Cole et al. 1999). One of the most characteristic features of the former is that the dust is divided in two main phases: molecular cloud complexes, where stars are assumed to be born, and the diffuse interstellar medium. As a consequence, stellar populations of different ages have different geometrical relationships with the two phases, which is essential in understanding several observed properties of galaxies, in particular those undergoing major episodes of star formation at any redshift. Indeed, our merged GRASIL+GALFORM model reproduces fairly well the SEDs of normal spirals and starbursts from the far-UV to the sub-mm and their internal extinction properties. In particular in the model the observed starburst attenuation law (Calzetti, 1999) is accounted for as an effect of geometry of stars and dust, and has nothing to do with the optical properties of dust grains.     

Isophot Mapping Observations of Carbon Stars

Mapping observations have been made toward five carbon stars in the far-infrared using ISOPHOT, an imaging photo-polarimeter on board the Infrared Space Observatory. Cold, very extended dust shells are clearly revealed in two of them. Y CVn is surrounded by a very extended, detached dust shell, which indicates a sudden decline of the mass-loss by two orders of magnitude in the last (1-2) × 104 years on a short time scale. The Hipparcos parallax reinforces our previous conclusion that Y CVn is a J-type carbon star on the asymptotic giant branch. U Ant shows a double shell structure, a compact dust shell surrounded by a very extended one. The outer shell has a brightness comparable to the dust shell of Y CVn. The structure indicates that there were two different high mass-loss phases separated by a period with a much lower mass-loss rate in between the two. The structure is consistent with the double dust shell proposed for this star by Izumiura et al. (1997) based on a detailed investigation of IRAS survey data. This revised version was published online in August 2006 with corrections to the Cover Date.     

Gravitational Instability in the Dust Layer of a Protoplanetary Disk with Interaction between the Layer and the Surrounding Gas

We consider gravitational instability of the dust layer in the midplane of a protoplanetary disk with turbulence and shear stresses between the gas in the disk and that in the dust layer. We solve a linearized system of hydrodynamic equations for perturbations of dust (monodisperse) and gas phases in the incompressible gas approximation. We take into account the gas drag of solid particles (dust aggregates), turbulent diffusion and the velocity dispersion of particles, and the perturbation of the azimuthal velocity of gas in the layer upon the transfer of angular momentum from solid particles to it and from this gas to the surrounding gas in the disk. We obtain and solve the dispersion equation for the layer with the ratio of surface densities of the dust phase and gas being well above unity. The following parameters of gravitational instability in the dust layer are calculated: the critical surface density of solid matter and the Stokes number of particles corresponding to the onset of instability, the wavelength range in which instability occurs, and the rate of its growth as a function of the perturbation wavelength in the circumsolar disk at radial distances of 1 and 10 AU. We show that at 10 AU, the maximum instability growth rate increases due to the transfer of angular momentum of gas in the layer to gas outside it, a new maximum emerges at a longer wavelength, a long-wavelength instability “tail” forms, and the critical surface density initiating instability decreases relative to that determined without the transfer of angular momentum to gas outside the layer. None of these effects are observed at 1 AU, since instability in this region probably develops faster than the transfer of angular momentum to the surrounding gаs of a protoplanetary disk occurs.     

Assessing uncertainties in dust emission in the Aral Sea region caused by meteorological fields predicted with a mesoscale model

This study investigates how the choice of the planetary boundary layer (PBL) parameterization and dust emission scheme affects the prediction of dust entrainment simulated with a regional mesoscale model. For this analysis we consider a representative dust episode which occurred on April 2001 in the Aral Sea region. The meteorological fields were simulated using the PSU/NCAR MM5 modeling system considering two different boundary layer parameterizations. In each case, emitted dust fluxes were computed off-line by incorporating MM5 meteorological fields into the dust module DuMo. Several dust emission schemes with a prescribed erodible fraction and fixed threshold wind speed were the subject of our analysis. Implications to assessment of the anthropogenic fraction of dust emitted in the Aral region were investigated by conducting the full, half, and no lake modeling experiments.Our results show that the discrepancies in dust fluxes between the two different PBLs are much higher compared to the discrepancy associated with the use of considered dust production schemes. Furthermore, the choice of the PBL affects the timing and duration of a modeled dust event. We demonstrate that different combinations of the PBL parameterization and wind- or friction velocity-driven dust emission schemes can result in up to about a 50% difference in predicted dust mass caused by the Aral Sea desiccation. We found that the drying-up of the Aral cannot only affect the dust emission by expanding the source area, but also by affecting atmospheric characteristics, especially winds. These competitive factors add further complexity to quantification of the anthropogenic dust fraction in the region.     

Dust ring formation due to sublimation of dust grains drifting radially inward by the Poynting-Robertson drag: An analytical model

Dust particles exposed to the stellar radiation and wind drift radially inward by the Poynting-Robertson (P-R) drag and pile up at the zone where they begin to sublime substantially. The reason they pile up or form a ring is that their inward drifts due to the P-R drag are suppressed by stellar radiation pressure when the ratio of radiation pressure to stellar gravity on them increases during their sublimation phases. We present analytic solutions to the orbital and mass evolution of such subliming dust particles, and find their drift velocities at the pileup zone are almost independent of their initial semimajor axes and masses. We derive analytically an enhancement factor of the number density of the particles at the outer edge of the sublimation zone from the solutions. We show that the formula of the enhancement factor reproduces well numerical simulations in the previous studies. The enhancement factor for spherical dust particles of silicate and carbon extends from 3 to more than 20 at stellar luminosities L_?=0.8-500L_⊙, where L_⊙ is solar luminosity. Although the enhancement factor for fluffy dust particles is smaller than that for spherical particles, sublimating particles inevitably form a dust ring as long as their masses decrease faster than their surface areas during sublimation. The formulation is applicable to dust ring formation for arbitrary shape and material of dust in dust-debris disks as well as in the Solar System.     

Gap formation in the dust layer of 3D protoplanetary disks

We numerically model the evolution of dust in a protoplanetary disk using a two-phase (gas+dust) Smoothed Particle Hydrodynamics (SPH) code, which is non-self-gravitating and locally isothermal. The code follows the three dimensional distribution of dust in a protoplanetary disk as it interacts with the gas via aerodynamic drag. In this work, we present the evolution of a disk comprising 1% dust by mass in the presence of an embedded planet for two different disk configurations: a small, minimum mass solar nebular (MMSN) disk and a larger, more massive Classical T Tauri star (CTTS) disk. We then vary the grain size and planetary mass to see how they effect the resulting disk structure. We find that gap formation is much more rapid and striking in the dust layer than in the gaseous disk and that a system with a given stellar, disk and planetary mass will have a different appearance depending on the grain size and that such differences will be detectable in the millimetre domain with ALMA. For low mass planets in our MMSN models, a gap can open in the dust disk while not in the gas disk. We also note that dust accumulates at the external edge of the planetary gap and speculate that the presence of a planet in the disk may facilitate the growth of planetesimals in this high density region.     

Dust coma morphology in the Deep Impact images of Comet 9P/Tempel 1

We present an overview of the dust coma observations of Comet Tempel 1 that were obtained during the approach and encounter phases of the Deep Impact mission. We use these observations to set constraints on the pre-impact activity of the comet and discuss some preliminary results. The temporal and spatial changes that were observed during approach reveal three distinct jets rotating with a 1.7-day periodicity. The brightest jet produces an arcuate feature that expands outward with a projected velocity of about 12 m s⁻¹, suggesting that the ambient dust coma is dominated by millimeter-sized dust grains. As the spatial resolution improves, more jets and fans are revealed. We use stereo pairs of high-resolution images to put some crude constraints on the source locations of some of the brightest features. We also present a number of interesting coma features that were observed, including surface jets detected at the limb of the nucleus when the exposed ice patches are passing over the horizon, and features that appear to be jets emanating from unilluminated sources near the negative pole. We also provide a list of 10 outbursts of various sizes that were observed in the near-continuous monitoring during the approach phase.     

Constraints on the nucleus and dust properties from mid-infrared imaging of comet Hyakutake

We use Mie scattering theory to determine the expected thermal emission from dust grains in cometary comae and apply these results to mid-infrared images of comet Hyakutake (C/1996 B2) obtained preperihelion in 1996 March. Calculations were performed for dust grains in the size range from 0.1 to 10 micrometers for two different compositions: amorphous olivine (a silicate glass) and an organic residue mixture. The resulting emission efficiencies are complicated functions of wavelength and particle size and are significantly different for the two materials in question. The Hyakutake data set consists of three nights of high-resolution imaging (100-150 km pixel-1 at the comet) of the inner coma at 8.7, 11.7, 12.5, and 19.7 micrometers. Attempts to fit the observed colors (ratios of fluxes at different wavelengths) using a single grain composition failed. However, fits to the data were achieved for all three nights using a mixture of approximately 1 micrometer olivine grains and approximately 7 micrometers organic grains. The resulting olivine mass fraction was between 8% and 16% of the total dust mass-loss rate. We also estimate the radius of the nucleus to be r = 2.1 +/- 0.4 km.     

Dust models confronted with observations

We critically discuss the two major models for interstellar dust grains developed by Greenberg and Mathis and their co-workers. So far, knowledge about the dust has been obtained mainly from extinction and polarization measurements in the UV, visual and near IR. Analysis of dust emission at longer wavelengths generally requires an evaluation of the radiative transfer. Such a procedure is necessary to assess optical grain properties at long wavelengths and also to interpret the emission from the different kinds of IR objects. The 10 m silicate resonance and the submm dust absorptivity are presented as two examples. We then briefly summarize the implications which follow from the existence of extremely small dust grains.Paper presented at a Workshop on The Role of Dust in Dense Regions of Interstellar Matter, held at Georgenthal, G.D.R., in March 1986.     

ISO-SWS Observations of the Time Variability of Oxygen-righ Mira Variables

We present the first observation of variation of the infrared spectrum of circumstellar dust emission around an oxygen-rich Mira variable star, Z Cyg, over an entire light variation cycle, based on the periodic SWS01 observations. The 10 μm and 20 μm ‘silicate features’ become stronger relative to the photospheric emission at maximum that at minimum. In addition, the relative intensity of the 10 μm to the 20 μm features increases at maximum, indicating an increase of the temperature of circumstellar dust grain. A simple model analysis suggests enhanced dust formation near the photosphere around maximum, leading to a scenario that dust nucleation may have occurred near minimum and the dust grains may have subsequently grown till maximum. The strong observed variation suggests that the variability must be taken into account in the interpretation of the infrared spectra of oxygen-rich Mira variables. This revised version was published online in August 2006 with corrections to the Cover Date.     

Axisymmetry in protoplanetary nebulae – II. A near-infrared imaging polarimetric survey

In an imaging polarimetry survey of candidate post-AGB stars, scattering envelopes are detected around 20 objects. These envelopes represent the final mass-loss phases at the end of the AGB. In all cases, evidence for axisymmetry in the dust density is seen, suggesting that the presence of an axisymmetric outflow may be a ubiquitous phenomenon of the AGB to post-AGB transition. We use the polarized flux images to classify the objects into detached shell, bipolar and unresolved targets. Modelling based on a simple axisymmetric shell geometry supports the contention that post-AGB objects fall into one of two classes that differ in the amount of dust in the circumstellar environment: the detached shells correspond to stars with an optically thin expanding circumstellar envelope (CSE) whereas the bipolar and unresolved targets have optically thick dust structures, probably in the form of discs, which remain bound to the star, rather than partaking in the expansion of the AGB CSE. It is suggested that this bifurcation in morphology is rooted in the presence or absence of a binary companion, which determines whether or not a disc forms. Because the detached shell objects also appear axisymmetric, an additional mechanism for generating the axisymmetry, such as a magnetically shaped outflow, is needed if they do indeed have single star progenitors.     

An explanation for time dependent variability of the solar dust ring

There are several observations showing an enhancement of infrared emission and optical polarization at a distance of 4R , (R is the angular radius of the Sun) implying a ring of dust in near-ecliptic orbit about the Sun; but there is an almost equal number of observations, which do not show any such enhancement. We plotted the observational results for the detection and the non-detection of the circum-solar dust on a diagram for the variation of the sun-spot number with time, and found that its detection and non-detection occurred near the solar minimum and maximum phases, respectively. We present the possibility that this phenomena is caused by an additional process, i.e., the Lorentz force acting on a charged dust particle. Since, at the maximum phase a dust particle in a near-solar region acquires a higher positive potential, and the solar magnetic field is very strong, the Lorentz force becomes dominant and can affect the orbit of the dust particle.     

Optical polarimetric study of open clusters: distribution of interstellar matter towards NGC 654

We present new B ,  V and R linear polarimetric observations for 61 stars towards the region of the young open cluster NGC 654. In this study we found evidence for the presence of at least two layers of dust along the line of sight to the cluster. The distances to the two dust layers are estimated to be ∼200 pc and ∼1 kpc which are located much closer to the Sun than the cluster (∼2.4 kpc). Both the dust layers have their local magnetic field orientation nearly parallel to the direction of the Galactic plane. The foreground dust layer is found to have a ring morphology with the central hole coinciding with the centre of the cluster. The foreground dust grains are suggested to be mainly responsible for both the observed differential reddening and the polarization towards the cluster.     

Assessing the power law hypothesis for the size-frequency distribution of terrestrial and martian dust devils

Competing hypotheses for the diameter dependence of terrestrial and martian dust devil frequency are assessed using new field observations from two sites in the southwestern United States. We show that at diameters less than 12 m, our observed dust devil size-frequency distributions are better fit by an exponential function than by a power law formulation, and discuss the implications for larger dust devils on Earth and Mars.     

Settling and growth of dust particles in a laminar phase of a low-mass solar nebula

It is considered that some vertical convection as well as possible turbulence in an early phase of solar nebula soon terminates owing to diminution of the temperature dependence of dust opacity due to rapid growth of dust particles. We reexamine settling and growth of dust particles in the subsequent laminar phase of the solar nebula in detail, treating a dust layer as a two-component fluid composed of the dust and the gas. We obtain analytic expressions for the settling path, the growing size, and the settling time. The settling process is divided into two phases, i.e., an early gas-dominant phase and a later dust-dominant phase. So far, only the former phase, where the particle path finally turns from vertical to radial, has been investigated. In the latter phase, dust particles drag the gas, rather than the gas does dust particles. Consequently, the particle path turns from radial to vertical. Dust particles grow most appreciably and rapidly in a radially sweeping phase. The final radii of dust particles at the onset of gravitational instability of the dust layer are 20, 5.9, and 0.60 cm in the Earth's, Jupiter's, and Neptune's zones, respectively. These values are much smaller than those obtained previously by S.J. Weidenschilling [1980, Icarus44, 172–189] and Y. Nakagawa et al. [1981, Icarus45, 517–528]. The total settling times are 1.9 × 10³, 4.6 × 10³, and 2.8 × 10⁴ years in the above-mentioned three zones, respectively. These are somewhat smaller than those obtained by the previous studies. Most of the settling time is spent in the early vertically settling phase.     

Preliminary analysis of the MER magnetic properties experiment using a computational fluid dynamics model

Motivated by questions raised by the magnetic properties experiments on the NASA Mars Pathfinder and Mars Exploration Rover (MER) missions, we have studied in detail the capture of airborne magnetic dust by permanent magnets using a computational fluid dynamics (CFD) model supported by laboratory simulations. The magnets studied are identical to the capture magnet and filter magnet on MER, though results are more generally applicable.The dust capture process is found to be dependent upon wind speed, dust magnetization, dust grain size and dust grain mass density. Here we develop an understanding of how these parameters affect dust capture rates and patterns on the magnets and set bounds for these parameters based on MER data and results from the numerical model. This results in a consistent picture of the dust as containing varying amounts of at least two separate components with different physical properties.     

High angular resolution millimeter observations of circumstellar disks

In this lecture, we review the properties of protoplanetary disks as derived from high angular resolution observations at millimeter wavelengths. We discuss how the combination of several different high angular resolution techniques allow us to probe different regions of the disk around young stellar objects and to derive the properties of the dust when combined with sophisticated disk models. The picture that emerges is that the dust in circumstellar disks surrounding pre-main sequence stars is in many cases significantly evolved compared to the dust in molecular clouds and the interstellar medium. It is however still difficult to derive a consistent picture and timeline for dust evolution in disks as the observations are still limited to small samples of objects.We also review the evidence for and properties of disks around high-mass young stellar objects and the implications on their formation mechanisms. The study of massive YSOs is complicated by their short lifetimes and larger average distances. In most cases high angular resolution data at millimeter wavelengths are the only method to probe the structure of disks in these objects.We provide a summary of the characteristics of available high angular resolution millimeter and submillimeter observatories. We also describe the characteristics of the ALMA observatory being constructed in the Chilean Andes. ALMA is going to be the world leading observatory at millimeter wavelengths in the coming decades, the project is now in its main construction phase with early science activities envisaged for 2010 and full science operations for 2012.     

Modelling the circumnuclear coma of comets: objectives,methods and recent results

We investigate whether the modelling of the immediate vicinity of an active nucleus—currently unobservable—can, as the modelling of the outer, observable coma, be based on unrealistic simple assumptions such as those of nucleus and dust grains sphericity. We point out the inconsistency of models based on such assumptions, which, to manage compatibility with the observations, have to introduce additional assumptions that conflict with the previous ones, such as the existence of active areas of the nucleus. We argue that, while the outer coma models being phenomenological in nature, can perhaps tolerate such inconsistencies, the circumnuclear coma models must be predictive, having to obviate the lack of observational data, and therefore must exclude implausible and ad hoc assumptions, and advocate only well-understood physical processes and duly validated modelling methods. We describe the first steps of development of a predictive circumnuclear coma model, and present a set of results obtained with parameters fitted to comet P/Wirtanen, the target of the Rosetta mission, but of a quite general significance. Considering, first an inhomogeneous spherical nucleus with spherical dust grains, and then an aspherical homogeneous nucleus with spherical dust grains, we show that, in both cases (1) the surface temperature and initial gas parameters differ considerably from the Hertz–Knudsen values; (2) the near-surface gas and dust flows are not in general vertical, (3) the gas and dust density do not always monotonically decrease outwards, (4) the gas and dust velocity vary strongly from point to point, (5) shock structures are formed, which result in the formation of pseudo-jets of dust grains originating from various points of the surface. No simple method to distinguish between dust structures created by the surface inhomogeneity and by the surface orography is found. We show, for the first time, the deformation of the near-nucleus dust coma during a full rotation of an homogeneous, aspherical nucleus. We also show that identical active regions located at different points of an inhomogeneous spherical nucleus produce very different dust distributions, suggesting that the dust distribution is also strongly deformed during the rotation of such a nucleus. Finally, we consider, for the first time, a spherical homogeneous nucleus emitting aspherical dust grains. We show that, in such a case, the terminal grain velocity depends upon the shape, initial position, and even possibly upon the initial orientation of the grain at the surface, so that there cannot exist a precise relation between terminal velocity and dust grain mass. We conclude that, far from giving an approximate or average representation of the circumnuclear coma, the classical modelling approach yields in this region predictions that are in total conflict with the real behaviour of the gas and dust. As a most dramatic consequence, the use of this classical approach may have obscured completely the significance of the few direct and of the many indirect informations acquired hitherto on the nucleus activity.     

Effects of dust abundance on the far-infrared colours of blue compact dwarf galaxies

We investigate the far-infrared (FIR) properties of a sample of blue compact dwarf galaxies (BCDs) observed by AKARI . By utilizing the data at wavelengths of  λ= 65  , 90 and 140 μm, we find that the FIR colours of the BCDs are located at the natural high-temperature extension of those of the Milky Way and the Magellanic Clouds. This implies that the optical properties of dust in BCDs are similar to those in the Milky Way. Indeed, we explain the FIR colours by assuming the same grain optical properties, which may be appropriate for amorphous dust grains, and the same size distribution as those adopted for the Milky Way dust. Since both interstellar radiation field and dust optical depth affect the dust temperature, it is difficult to distinguish which of these two physical properties is responsible for the change of FIR colours. Then, in order to examine if the dust optical depth plays an important role in determining the dust temperature, we investigate the correlation between FIR colour (dust temperature) and dust-to-gas ratio. We find that the dust temperature tends to be high as the dust-to-gas ratio decreases but that this trend cannot be explained by the effect of dust optical depth. Rather, it indicates a correlation between dust-to-gas ratio and interstellar radiation field. Although the metallicity may also play a role in this correlation, we suggest that the dust optical depth could regulate the star formation activities, which govern the interstellar radiation field. We also mention the importance of submillimetre data in tracing the emission from highly shielded low-temperature dust.