Formation of chondrules in radiative shock waves I. First results, spherical dust particles, stationary shocks

The formation of chondrules in the protoplanetary nebulae causes many questions concerning the formation process, the source of energy for melting the rims, and the composition of the origin material. The aim of this work is to explore the heating of the chondrule in a single precursor as is typical for radiation hydrodynamical shock waves. We take into account the gas-particle friction for the duration of the shock transition and calculate the heat conduction into the chondrules. These processes are located in the protoplanetary nebulae at a region around 2.5 AU, which is considered to be the most likely place of chondrule formation. The present models are a first step towards computing radiative shock waves occurring in a particle-rich environment. We calculated the shock waves using one-dimensional, time-independent equations of radiation hydrodynamics involving realistic gas and dust opacities and gas-particle friction. The evolution of spherical chondrules was followed by solving the heat conduction equation on an adaptive grid. The results for the shock-heating event are consistent with the cosmochemical constraints of chondrule properties. The calculations yield a relative narrow range for density or temperature to meet the requested heating rates of $R > 10^4\,K\, h^{-1}$ as extracted from cosmochemical constraints. Molecular gas, opacities with dust, and a protoplanetary nebula with accretion are necessary requirements for a fast heating process. The thermal structure in the far post-shock region is not fully consistent with experimental constraints on chondrule formation since the models do not include additional molecular cooling processes.Comment: 8 pages,5 figure

A cosmic ray driven instability

The interaction between energetic charged particles and thermal plasma which forms the basis of diffusive shock acceleration leads also to interesting dynamical phenomena. For a compressional mode propagating in a system with homogeneous energetic particle pressure it is well known that friction with the energetic particles leads to damping. The linear theory of this effect has been analyzed in detail by Ptuskin. Not so obvious is that a non-uniform energetic particle pressure can addition amplify compressional disturbances. If the pressure gradient is sufficiently steep this growth can dominate the frictional damping and lead to an instability. It is important to not that this effect results from the collective nature of the interaction between the energetic particles and the gas and is not connected with the Parker instability, nor with the resonant amplification of Alfven waves

Three-component modeling of C-rich AGB star winds I. Method and first results

Radiative acceleration of newly-formed dust grains and transfer of momentum from the dust to the gas plays an important role for driving winds of AGB stars. Therefore a detailed description of the interaction of gas and dust is a prerequisite for realistic models of such winds. In this paper we present the method and first results of a three-component time-dependent model of dust-driven AGB star winds. With the model we plan to study the role and effects of the gas-dust interaction on the mass loss and wind formation. The wind model includes separate conservation laws for each of the three components of gas, dust and the radiation field and is developed from an existing model which assumes position coupling between the gas and the dust. As a new feature we introduce a separate equation of motion for the dust component in order to fully separate the dust phase from the gas phase. The transfer of mass, energy and momentum between the phases is treated by interaction terms. We also carry out a detailed study of the physical form and influence of the momentum transfer term (the drag force) and three approximations to it. In the present study we are interested mainly in the effect of the new treatment of the dust velocity on dust-induced instabilities in the wind. As we want to study the consequences of the additional freedom of the dust velocity on the model we calculate winds both with and without the separate dust equation of motion. The wind models are calculated for several sets of stellar parameters. We find that there is a higher threshold in the carbon/oxygen abundance ratio at which winds form in the new model. The winds of the new models, which include drift, differ from the previously stationary winds, and the winds with the lowest mass loss rates no longer form.Comment: 15 pages, 5 figures, accepted by A&

Approximate supernova remnant dynamics with cosmic ray production

Supernova explosions are the most violent and energetic events in the galaxy and have long been considered probably sources of Cosmic Rays. Recent shock acceleration models treating the Cosmic Rays (CR's) as test particles nb a prescribed Supernova Remnant (SNR) evolution, indeed indicate an approximate power law momentum distribution f sub source (p) approximation p(-a) for the particles ultimately injected into the Interstellar Medium (ISM). This spectrum extends almost to the momentum p = 1 million GeV/c, where the break in the observed spectrum occurs. The calculated power law index approximately less than 4.2 agrees with that inferred for the galactic CR sources. The absolute CR intensity can however not be well determined in such a test particle approximation

Hydrodynamic simulations of captured protoatmospheres around Earth-like planets

Young terrestrial planets, when they are still embedded in a circumstellar disk, accumulate an atmosphere of nebula gas. The evolution and eventual evaporation of the protoplanetary disk affect the structure and dynamics of the planetary atmosphere. These processes, combined with other mass loss mechanisms, such as thermal escape driven by extreme ultraviolet and soft X-ray radiation (XUV) from the young host star, determine how much of the primary atmosphere, if anything at all, survives into later stages of planetary evolution. Our aim is to explore the structure and the dynamic outflow processes of nebula-accreted atmospheres in dependency on changes in the planetary environment. We integrate stationary hydrostatic models and perform time-dependent dynamical simulations to investigate the effect of a changing nebula environment on the atmospheric structure and the timescales on which the protoatmosphere reacts to these changes. We find that the behavior of the atmospheres strongly depends on the mass of the planetary core. For planets of about Mars-mass the atmospheric structure, and in particular the atmospheric mass, changes drastically and on very short timescales whereas atmospheres around higher mass planets are much more robust and inert

Acceleration of cosmic rays in supernova-remnants

It is commonly accepted that supernova-explosions are the dominant source of cosmic rays up to an energy of 10 to the 14th power eV/nucleon. Moreover, these high energy particles provide a major contribution to the energy density of the interstellar medium (ISM) and should therefore be included in calculations of interstellar dynamic phenomena. For the following the first order Fermi mechanism in shock waves are considered to be the main acceleration mechanism. The influence of this process is twofold; first, if the process is efficient (and in fact this is the cas) it will modify the dynamics and evolution of a supernova-remnant (SNR), and secondly, the existence of a significant high energy component changes the overall picture of the ISM. The complexity of the underlying physics prevented detailed investigations of the full non-linear selfconsistent problem. For example, in the context of the energy balance of the ISM it has not been investigated how much energy of a SN-explosion can be transfered to cosmic rays in a time-dependent selfconsistent model. Nevertheless, a lot of progress was made on many aspects of the acceleration mechanism

Solar-Like Cycle in Asymptotic Giant Branch Stars

I propose that the mechanism behind the formation of concentric semi-periodic shells found in several planetary nebulae (PNs) and proto-PNs, and around one asymptotic giant branch (AGB) star, is a solar-like magnetic activity cycle in the progenitor AGB stars. The time intervals between consecutive ejection events is about 200-1,000 years, which is assumed to be the cycle period (the full magnetic cycle can be twice as long, as is the 22-year period in the sun). The magnetic field has no dynamical effects; it regulates the mass loss rate by the formation of magnetic cool spots. The enhanced magnetic activity at the cycle maximum results in more magnetic cool spots, which facilitate the formation of dust, hence increasing the mass loss rate. The strong magnetic activity implies that the AGB star is spun up by a companion, via a tidal or common envelope interaction. The strong interaction with a stellar companion explains the observations that the concentric semi-periodic shells are found mainly in bipolar PNs.Comment: 10 pages, submitted to Ap

Hydrodynamic Simulation of Supernova Remnants Including Efficient Particle Acceleration

A number of supernova remnants (SNRs) show nonthermal X-rays assumed to be synchrotron emission from shock accelerated TeV electrons. The existence of these TeV electrons strongly suggests that the shocks in SNRs are sources of galactic cosmic rays (CRs). In addition, there is convincing evidence from broad-band studies of individual SNRs and elsewhere that the particle acceleration process in SNRs can be efficient and nonlinear. If SNR shocks are efficient particle accelerators, the production of CRs impacts the thermal properties of the shock heated, X-ray emitting gas and the SNR evolution. We report on a technique that couples nonlinear diffusive shock acceleration, including the backreaction of the accelerated particles on the structure of the forward and reverse shocks, with a hydrodynamic simulation of SNR evolution. Compared to models which ignore CRs, the most important hydrodynamical effects of placing a significant fraction of shock energy into CRs are larger shock compression ratios and lower temperatures in the shocked gas. We compare our results, which use an approximate description of the acceleration process, with a more complete model where the full CR transport equations are solved (i.e., Berezhko et al., 2002), and find excellent agreement for the CR spectrum summed over the SNR lifetime and the evolving shock compression ratio. The importance of the coupling between particle acceleration and SNR dynamics for the interpretation of broad-band continuum and thermal X-ray observations is discussed.Comment: Accepted for publication in A & A; 14 pages including 11 figure

TeV Gamma Rays Expected from Supernova Remnants in Different Uniform Interstellar Media

Calculations of the expected TeV $\gamma$-ray emission, produced by accelerated cosmic rays (CRs) in nuclear collisions, from supernova remnants evolving in a uniform interstellar medium (ISM) are presented. The aim is to study the sensitivity of $\gamma$-ray production to a physical parameter set. Apart from its general proportionality to N_H, it is shown that the $\gamma$-ray production essentially depends upon the ratio of the CR diffusion coefficient $\kappa$ to a critical value $\kappa_{crit}=10(B_0/5 \mu{G})(N_H/0.3 {cm}^{-3})^{-1/3}\kappa_B$, where B_0 and N_H are the magnetic field and hydrogen number density of the ISM, and $\kappa_B$ denotes the Bohm diffusion coefficient. If $\kappa$ is of the same order or lower than $\kappa_{crit}$, then the peak TeV $\gamma$-ray flux in the Sedov evolutionary phase, normalized to a distance of 1 kpc, is about 10^{-10}(N_H/0.3 {cm}^{-3}) photons cm^{-2} s^{-1}. For a CR diffusion coefficient that is significantly larger than $\kappa_{crit}$, the CR cutoff energy is less than 10 TeV and the expected $\gamma$-ray flux at 1 TeV is considerably below the presently detectable level of 10^{-12} photons cm^{-2} s^{-1}. The same is of course true for a SNR in the rarified, so-called hot ISM.Comment: 9 pages, 2 figures, to appear in Astroparticle Physic