Carbonates in space - The challenge of low temperature data

Carbonates have repeatedly been discussed as possible carriers of stardust emission bands. However, the band assignments proposed so far were mainly based on room temperature powder transmission spectra of the respective minerals. Since very cold calcite grains have been claimed to be present in protostars and in Planetary Nebulae such as NGC 6302, the changes of their dielectric functions at low temperatures are relevant from an astronomical point of view. We have derived the IR optical constants of calcite and dolomite from reflectance spectra - measured at 300, 200, 100 and 10K - and calculated small particle spectra for different grain shapes, with the following results: i) The absorption efficiency factors both of calcite and dolomite are extremely dependent on the particle shapes. This is due to the high peak values of the optical constants of CaCO3 and CaMg[CO3]2. ii) The far infrared properties of calcite and dolomite depend also very significantly on the temperature. Below 200K, a pronounced sharpening and increase in the band strengths of the FIR resonances occurs. iii) In view of the intrinsic strength and sharpening of the 44 mum band of calcite at 200-100K, the absence of this band -- inferred from Infrared Space Observatory data -- in PNe requires dust temperatures below 45K. iv) Calcite grains at such low temperatures can account for the '92' mum band, while our data rule out dolomite as the carrier of the 60-65 mum band. The optical constants here presented are publicly available in the electronic database http://www.astro.uni-jena.de/Laboratory/OCDBComment: 20 pages, 10 figures, accepted by ApJ, corrected typo

From interstellar abundances to grain composition: the major dust constituents Mg, Si and Fe

We analyse observational correlations for three elements entering into the composition of interstellar silicate and oxide grains. Using current solar abundances (Asplund et al. 2009), we convert the gas-phase abundances into dust-phase abundances for 196 sightlines. We deduce a sharp difference in abundances for sightlines located at low (|b|<30\degr) and high (|b|>30\degr) galactic latitudes. For high-latitude stars the ratios Mg/Si and Fe/Si in dust are close to 1.5. For disk stars they are reduced to ${\rm Mg/Si} \sim 1.2$ and ${\rm Fe/Si} \sim 1.05$. The derived numbers indicate that 1) the dust grains cannot be the mixture of silicates with olivine and pyroxene composition only and some amount of magnesium or iron (or both) should be in another population and 2) the destruction of Mg-rich grains in the warm medium is more effective than of Fe-rich grains. We reveal a decrease of dust-phase abundances and correspondingly an increase of gas-phase abundances with distance $D$ for stars with D\ga 400\,pc. We attribute this fact to an observational selection effect: a systematic trend toward smaller observed hydrogen column density for distant stars. We find differences in abundances for disk stars with low (E({\rm B-V}) \la 0.2) and high (E({\rm B-V}) \ga 0.2) reddenings which reflect the distinction between the sightlines passing through diffuse and translucent interstellar clouds. For Scorpius-Ophiuchus we detect an uniform increase of dust-phase abundances of Mg and Si with an increase of the ratio of total to selective extinction $R_{\rm V}$ and a decrease of the strength of the far-UV extinction. This is the first evidence for a growth of Mg-Si grains due to accretion in the interstellar medium.Comment: 16 pages, 16 figures, accepted for publication in Astronomy and Astrophysic

Migration and Accretion of Protoplanets in 2D and 3D Global Hydrodynamical Simulations

Planet evolution is tightly connected to the dynamics of both distant and close disk material. Hence, an appropriate description of disk-planet interaction requires global and high resolution computations, which we accomplish by applying a Nested-Grid method. Through simulations in two and three dimensions, we investigate how migration and accretion are affected by long and short range interactions. For small mass objects, 3D models provide longer growth and migration time scales than 2D ones do, whereas time lengths are comparable for large mass planets.Comment: 4 pages, 4 figures; to appear in the Conference Proceedings of "Scientific Frontiers in Research on Extrasolar Planets

Chemo-dynamical deuterium fractionation in the early solar nebula: The origin of water on Earth and in asteroids and comets

Formation and evolution of water in the Solar System and the origin of water on Earth constitute one of the most interesting questions in astronomy. The prevailing hypothesis for the origin of water on Earth is by delivery through water-rich small Solar system bodies. In this paper, the isotopic and chemical evolution of water during the early history of the solar nebula, before the onset of planetesimal formation, is studied. A gas-grain chemical model that includes multiply-deuterated species and nuclear spin-states is combined with a steady-state solar nebula model. To calculate initial abundances, we simulated 1 Myr of evolution of a cold and dark TMC1-like prestellar core. Two time-dependent chemical models of the solar nebula are calculated over 1 Myr: (1) a laminar model and (2) a model with 2D turbulent mixing. We find that the radial outward increase of the H2O D/H ratio is shallower in the chemo-dynamical nebular model compared to the laminar model. This is related to more efficient de-fractionation of HDO via rapid gas-phase processes, as the 2D mixing model allows the water ice to be transported either inward and thermally evaporated or upward and photodesorbed. The laminar model shows the Earth water D/H ratio at r ~<2.5 AU, while for the 2D chemo-dynamical model this zone is larger, r ~<9 AU. Similarly, the water D/H ratios representative of the Oort-family comets, ~2.5-10 x 10-4, are achieved within ~2-6 AU and ~2-20 AU in the laminar and the 2D model, respectively. We find that with regards to the water isotopic composition and the origin of the comets, the mixing model seems to be favored over the laminar model.Comment: 28 pages, 5 figures, 2 tables (accepted for publication by ApJ

Do we need to know the temperature in prestellar cores?

Molecular line observations of starless (prestellar) cores combined with a chemical evolution modeling and radiative transfer calculations are a powerful tool to study the earliest stages of star formation. However, conclusions drawn from such a modeling may noticeably depend on the assumed thermal structure of the cores. The assumption of isothermality, which may work well in chemo-dynamical studies, becomes a critical factor in molecular line formation simulations. We argue that even small temperature variations, which are likely to exist in starless cores, can have a non-negligible effect on the interpretation of molecular line data and derived core properties. In particular, ``chemically pristine'' isothermal cores (low depletion) can have centrally peaked C$^{18}$O and C$^{34}$S radial intensity profiles, while having ring-like intensity distributions in models with a colder center and/or warmer envelope assuming the same underlying chemical structure. Therefore, derived molecular abundances based on oversimplified thermal models may lead to a mis-interpretation of the line data.Comment: ApJL, accepte

Fragmentation, infall, and outflow around the showcase massive protostar NGC7538 IRS1 at 500 AU resolution

Aims: Revealing the fragmentation, infall, and outflow processes in the immediate environment around massive young stellar objects is crucial for understanding the formation of the most massive stars. Methods: With this goal in mind we present the so far highest spatial-resolution thermal submm line and continuum observations toward the young high-mass protostar NGC7538 IRS1. Using the Plateau de Bure Interferometer in its most extended configuration at 843mum wavelength, we achieved a spatial resolution of 0.2"x0.17", corresponding to ~500AU at a distance of 2.7\,kpc. Results: For the first time, we have observed the fragmentation of the dense inner core of this region with at least three subsources within the inner 3000 AU. The outflow exhibits blue- and red-shifted emission on both sides of the central source indicating that the current orientation has to be close to the line-of-sight, which differs from other recent models. We observe rotational signatures in northeast-southwest direction; however, even on scales of 500 AU, we do not identify any Keplerian rotation signatures. This implies that during the early evolutionary stages any stable Keplerian inner disk has to be very small (<=500 AU). The high-energy line HCN(4-3)v2=1 (E_u/k=1050K) is detected over an extent of approximately 3000 AU. In addition to this, the detection of red-shifted absorption from this line toward the central dust continuum peak position allows us to estimate infall rates of ~1.8x10^(-3)Msun/yr on the smallest spatial scales. Although all that gas will not necessarily be accreted onto the central protostar, nevertheless, such inner core infall rates are among the best proxies of the actual accretion rates one can derive during the early embedded star formation phase. These data are consistent with collapse simulations and the observed high multiplicity of massive stars.Comment: Accepted for Astronomy & Astrophysics, 8 pages, also available at http://www.mpia.de/homes/beuther/papers.htm

Effective medium theories for irregular fluffy structures: aggregation of small particles

We study the extinction efficiencies as well as scattering properties of particles of different porosity. Calculations are performed for porous pseudospheres with small size (Rayleigh) inclusions using the discrete dipole approximation. Five refractive indices of materials covering the range from $1.20+0.00i$ to $1.75+0.58i$ were selected. They correspond to biological particles, dirty ice, silicate, amorphous carbon and soot in the visual part of spectrum. We attempt to describe the optical properties of such particles using Lorenz-Mie theory and a refractive index found from some effective medium theory (EMT) assuming the particle is homogeneous. We refer to this as the effective model. It is found that the deviations are minimal when utilizing the EMT based on the Bruggeman mixing rule. Usually the deviations in extinction factor do not exceed $\sim 5$ for particle porosity ${\cal P}=0 - 0.9$ and size parameters x_{\rm porous} = 2 \pi r_{\rm s, porous}/\lambda \la 25. The deviations are larger for scattering and absorption efficiencies and smaller for particle albedo and asymmetry parameter. Our calculations made for spheroids confirm these conclusions. Preliminary consideration shows that the effective model represents the intensity and polarization of radiation scattered by fluffy aggregates quite well. Thus, the effective models of spherical and non-spherical particles can be used to significantly simplify computations of the optical properties of aggregates containing only Rayleigh inclusions.Comment: 24 pages, 9 figures, accepted for publication in Applied Optic

Size distribution of dust grains: A problem of self-similarity

Distribution functions describing the results of natural processes frequently show the shape of power laws, e.g., mass functions of stars and molecular clouds, velocity spectrum of turbulence, size distributions of asteroids, micrometeorites and also interstellar dust grains. It is an open question whether this behavior is a result simply coming about by the chosen mathematical representation of the observational data or reflects a deep-seated principle of nature. The authors suppose the latter being the case. Using a dust model consisting of silicate and graphite grains Mathis et al. (1977) showed that the interstellar extinction curve can be represented by taking a grain radii distribution of power law type n(a) varies as a(exp -p) with 3.3 less than or equal to p less than or equal to 3.6 (example 1) as a basis. A different approach to understanding power laws like that in example 1 becomes possible by the theory of self-similar processes (scale invariance). The beta model of turbulence (Frisch et al., 1978) leads in an elementary way to the concept of the self-similarity dimension D, a special case of Mandelbrot's (1977) fractal dimension. In the frame of this beta model, it is supposed that on each stage of a cascade the system decays to N clumps and that only the portion beta N remains active further on. An important feature of this model is that the active eddies become less and less space-filling. In the following, the authors assume that grain-grain collisions are such a scale-invarient process and that the remaining grains are the inactive (frozen) clumps of the cascade. In this way, a size distribution n(a) da varies as a(exp -(D+1))da (example 2) results. It seems to be highly probable that the power law character of the size distribution of interstellar dust grains is the result of a self-similarity process. We can, however, not exclude that the process leading to the interstellar grain size distribution is not fragmentation at all. It could be, e.g., diffusion-limited growth discussed by Sander (1986), who applied the theory of fractal geometry to the classification of non-equilibrium growth processes. He received D=2.4 for diffusion-limited aggregation in 3d-space