Phenomenological Theory of the Photoevaporation Front Instability

The dynamics of photoevaporated molecular clouds is determined by the ablative pressure acting on the ionization front. An important step in the understanding of the ensuing motion is to develop the linear stability theory for an initially flat front. Despite the simplifications introduced by linearization, the problem remains quite complex and still draws a lot of attention. The complexity is related to the large number of effects that have to be included in the analysis: acceleration of the front, possible temporal variation of the intensity of the ionizing radiation, the tilt of the radiation flux with respect to the normal to the surface, and partial absorption of the incident radiation in the ablated material. In this paper, we describe a model where all these effects can be taken into account simultaneously, and a relatively simple and universal dispersion relation can be obtained. The proposed phenomenological model may prove to be a helpful tool in assessing the feasibility of the laboratory experiments directed towards scaled modeling of astrophysical phenomena. PACS Numbers: 98.38.Dq, 98.38.Hv, 52.38.Mf, 5257.FG, 52.72.+v     

Blowing in the wind. II. Creation and redistribution of refractory inclusions in a turbulent protoplanetary nebula

Ca-Al rich refractory mineral inclusions (CAIs) found at 1-6% mass fraction in primitive chondrites appear to be 1-3 million years older than the dominant (chondrule) components which were accreted into the same parent bodies. A prevalent concern is that it is difficult to retain CAIs for this long against gas-drag-induced radial drift into the Sun. We reassess the situation in terms of a hot inner (turbulent) nebula context for CAI formation, using analytical models of nebula evolution and particle diffusion. We show that outward radial diffusion in a weakly turbulent nebula can overcome inward drift, and prevent significant numbers of CAI-size particles from being lost into the Sun for times on the order of 10⁶ years. CAIs can form early, when the inner nebula was hot, and persist in sufficient abundance to be incorporated into primitive planetesimals at a much later time. Small (?0.1 mm diameter) CAIs persist for longer times than large (?5 mm diameter) ones. To obtain a quantitative match to the observed volume fractions of CAIs in chondrites, another process must be allowed for: a substantial enhancement of the inner hot nebula in silicate-forming material, which we suggest was caused by rapid inward drift of meter-sized objects. This early in nebula history, the drifting rubble would have a carbon content probably an order of magnitude larger than even the most primitive (CI) carbonaceous chondrites. Abundant carbon in the evaporating material would help keep the nebula oxygen fugacity low, plausibly solar, as inferred for the formation environment of CAIs. The associated production of a larger than canonical amount of CO₂ might also play a role in mass-independent fractionation of oxygen isotopes, leaving the gas rich in ¹⁶O as inferred from CAIs and other high temperature condensates.     

Formation of terrestrial planets in a dissipating gas disk with Jupiter and Saturn

We have performed N-body simulations on final accretion stage of terrestrial planets, including the eccentricity and inclination damping effect due to tidal interaction with a gas disk. We investigated the dependence on a depletion time scale of the disk, and the effect of secular perturbations by Jupiter and Saturn. In the final stage, terrestrial planets are formed through coagulation of protoplanets of about the size of Mars. They would collide and grow in a decaying gas disk. Kominami and Ida [Icarus 157 (2002) 43-56] showed that it is plausible that Earth-sized, low-eccentricity planets are formed in a mostly depleted gas disk. In this paper, we investigate the formation of planets in a decaying gas disk with various depletion time scales, assuming disk surface density of gas component decays exponentially with time scale of τ_gas. Fifteen protoplanets with are initially distributed in the terrestrial planet regions. We found that Earth-sized planets with low eccentricities are formed, independent of initial gas surface density, when the condition (τ_cross+τ_growth)/2?τ_gas?τ_cross is satisfied, where τ_cross is the time scale for initial protoplanets to start orbit crossing in a gas-free case and τ_growth is the time scale for Earth-sized planets to accrete during the orbit crossing stage. In the cases satisfying the above condition, the final masses and eccentricities of the largest planets are consistent with those of Earth and Venus. However, four or five protoplanets with the initial mass remain. In the final stage of terrestrial planetary formation, it is likely that Jupiter and Saturn have already been formed. When Jupiter and Saturn are included, their secular perturbations pump up eccentricities of protoplanets and tend to reduce the number of final planets in the terrestrial planet regions. However, we found that the reduction is not significant. The perturbations also shorten τ_cross. If the eccentricities of Jupiter and Saturn are comparable to or larger than present values (∼0.05), τ_cross become too short to satisfy the above condition. As a result, eccentricities of the planets cannot be damped to the observed value of Earth and Venus. Hence, for the formation of terrestrial planets, it is preferable that the secular perturbations from Jupiter and Saturn do not have significant effect upon the evolution. Such situation may be reproduced by Jupiter and Saturn not being fully grown, or their eccentricities being smaller than the present values during the terrestrial planets' formation. However, in such cases, we need some other mechanism to eliminate the problem that numerous Mars-sized planets remain uncollided.     

Simulations of self-magnetization in dusty space plasmas comparing different charge polarities and charge numbers of the dust

In partially ionized dusty space plasmas collisional momentum transfer between neutral and charged components of different inertia yields significant self-induced magnetic fields on remarkable short time scales of several dust Alfven times. Considering the proto-solar accretion disk previous first self-consistent plasma-neutral gas-dust simulations have shown, that this process yields a self-magnetization of 10^-5- on time scales of about 100 days. Thus, this mechanism is able to explain the remanent magnetization of chondrite type meteorite matter. New simulations show the quantitative dependence of the self-magnetization on polarity and charge numbers of the dust grains.     

Kinematics of the young stellar objects associated with the cometary globules in the Gum Nebula

An analysis of the proper motion measurements of the young stellar objects (YSOs) associated with the cometary globules (CGs) in the Gum Nebula is presented. While earlier studies based on the radial-velocity measurements of the CGs suggested expansion of the system of the CGs, the observed proper motion of the YSOs shows no evidence for expansion. In particular, the kinematics of two YSOs embedded in CGs are inconsistent with the supernova explosion of the companion of ζ Pup, about 1.5 Myr ago being the cause of the expansion of the system of the CGs. YSOs associated with the CGs share the average proper motion of the member stars of the Vela OB2 association. A few YSOs that have relatively large proper motions are found to show relatively low infrared excesses.     

Structure and kinematics of the interstellar medium around WR 134 and WR 135

Based on our H _α interferometry and 21-cm and CO observations, we analyze the structure and kinematics of the interstellar medium around the stars WR 134 and WR 135. We conclude that the HI bubble found here previously is associated with WR 135, not with WR 134. High-velocity motions of ionized gas that can be interpreted as expansion of the gas swept up by the stellar wind with a velocity up to 50–80 km s^?1 are observed around both stars. The line-of-sight velocity field of the ionized hydrogen in the Cygnus arm is shown to agree with the large-scale line-of-sight velocity distribution of the CO emission.     

Origin of the ocean on the Earth: Early evolution of water D/H in a hydrogen-rich atmosphere

The origin of the Earth's ocean has been discussed on the basis of deuterium/hydrogen ratios (D/H) of several sources of water in the Solar System. The average D/H of carbonaceous chondrites (CC's) is known to be close to the current D/H of the Earth's ocean, while those of comets and the solar nebula are larger by about a factor of two and smaller by about a factor of seven, respectively, than that of the Earth's ocean. Thus, the main source of the Earth's ocean has been thought to be CC's or adequate mixing of comets and the solar nebula. However, those conclusions are correct only if D/H of water on the Earth has remained unchanged for the past 4.5 Gyr. In this paper, we investigate evolution of D/H in the ocean in the case that the early Earth had a hydrogen-rich atmosphere, the existence of which is predicted by recent theories of planet formation no matter whether the nebula remains or not. Then we show that D/H in the ocean increases by a factor of 2-9, which is caused by the mass fractionation during atmospheric hydrogen loss, followed by deuterium exchange between hydrogen gas and water vapor during ocean formation. This result suggests that the apparent similarity in D/H of water between CC's and the current Earth's ocean does not necessarily support the CC's origin of water and that the apparent discrepancy in D/H is not a good reason for excluding the nebular origin of water.     

Interpretation of the carbon abundance in Saturn measured by Cassini

Spectral observations of Saturn from the far infrared spectrometer aboard the Cassini spacecraft [Flasar, F.M., et al., 2005. Temperatures, winds, and composition in the Saturnian system. Science 307, 1247-1251] have revealed that the C/H ratio in the planet is in fact about twice higher than previously derived from ground based observations and in agreement with the C/H value derived from Voyager IRIS by Courtin et al. [1984. The composition of Saturn's atmosphere at northern temperate latitudes from Voyager IRIS spectra - NH₃, PH₃, C₂H₂, C₂H₆, CH₃D, CH₄, and the Saturnian D/H isotopic ratio. Astrophys. J. 287, 899-916]. The implications of this measurement are reanalyzed in the present report on the basis that volatiles observed in cometary atmospheres, namely CO₂, CH₄, NH₃ and H₂S may have been trapped as solids in the feeding zone of the planet. CH₄ and H₂S may have been in the form of clathrate hydrates while CO₂ presumably condensed in the cooling solar nebula. Carbon may also have been incorporated in organics. Conditions of temperature and pressure ease the hydratation of NH₃. Such icy grains were included in planetesimals which subsequently collapsed into the hydrogen envelope of the planet, then resulting in C, N and S enrichments with respect to the solar abundance. Our calculations are consistent, within error bars, with observed elemental abundances on Saturn provided that the carbon trapped in planetesimals was mainly in the form of CH₄ clathrate and CO₂ ice (and maybe as organics) while nitrogen was in the form of NH₃ hydrate. Our approach has implications on the possible pattern of noble gases in Saturn, since we predict that contrary to what is observed in Jupiter, Ar and Kr should be in solar abundance while Xe might be strongly oversolar. The only way to verify this scenario is to send a probe making in situ mass spectrometer measurements. Our scenario also predicts that the ¹⁴N/¹⁵N ratio should be somewhat smaller in Saturn than measured in Jupiter by Galileo.     

Mechanism of generation of the emission bands in the dynamic spectrum of the Crab pulsar

We show that the proportionately spaced emission bands in the dynamic spectrum of the Crab pulsar fit the oscillations of the square of a Bessel function whose argument exceeds its order. This function has already been encountered in the analysis of the emission from a polarization current with a superluminal distribution pattern: a current whose distribution pattern rotates (with an angular frequency ω) and oscillates (with a frequency  Ω > ω  differing from an integral multiple of ω) at the same time. Using the results of our earlier analysis, we find that the dependence on frequency of the spacing and width of the observed emission bands can be quantitatively accounted for by an appropriate choice of the value of the single free parameter  Ω/ω  . In addition, the value of this parameter, thus implied by Hankins & Eilek's data, places the last peak in the amplitude of the oscillating Bessel function in question at a frequency  (∼Ω³/ω²)  that agrees with the position of the observed ultraviolet peak in the spectrum of the Crab pulsar. We also show how the suppression of the emission bands by the interference of the contributions from differing polarizations can account for the differences in the time and frequency signatures of the interpulse and the main pulse in the Crab pulsar. Finally, we put the emission bands in the context of the observed continuum spectrum of the Crab pulsar by fitting this broad-band spectrum (over 16 orders of magnitude of frequency) with that generated by an electric current with a superluminally rotating distribution pattern.