The interstellar chemistry of protostellar disks

A study has been undertaken of the gas-grain chemistry of protostellar disks which are sufficiently cool that in the outer regions, where the gas density is less than 10¹³ cm^–3 and the ionization rate highest, a bimolecular chemistry resembling that of dark clouds can occur. Since the gas-grain collision rate is so high, outgassing mantle molecules effectively determine the gas phase composition at any position in the disk. In contrast to previous work, a detailed gas phase chemistry is considered along with the accretion and desorption of mantle species which is controlled locally by the dust temperature.     

The first galaxies: assembly, cooling and the onset of turbulence

We investigate the properties of the first galaxies at   z ≳ 10  with highly resolved numerical simulations, starting from cosmological initial conditions and taking into account all relevant primordial chemistry and cooling. A first galaxy is characterized by the onset of atomic hydrogen cooling, once the virial temperature exceeds  ≃10⁴ K  , and its ability to retain photoheated gas. We follow the complex accretion and star formation history of a  ≃5 × 10⁷ M_⊙  system by means of a detailed merger tree and derive an upper limit on the number of Population III (Pop III) stars formed prior to its assembly. We investigate the thermal and chemical evolution of infalling gas and find that partial ionization at temperatures  ≳10⁴ K  catalyses the formation of  H₂  and hydrogen deuteride, allowing the gas to cool to the temperature of the cosmic microwave background. Depending on the strength of radiative and chemical feedback, primordial star formation might be dominated by intermediate-mass Pop III stars formed during the assembly of the first galaxies. Accretion on to the nascent galaxy begins with hot accretion, where gas is accreted directly from the intergalactic medium and shock heated to the virial temperature, but is quickly accompanied by a phase of cold accretion, where the gas cools in filaments before flowing into the parent halo with high velocities. The latter drives supersonic turbulence at the centre of the galaxy and could lead to very efficient chemical mixing. The onset of turbulence in the first galaxies thus likely marks the transition to Pop II star formation.     

Localized magnetic reconnection as a cause of extraplanar diffuse ionized gas in the Galactic halo

Many observations indicate the occurrence of ionized gas in the distant haloes of galaxies (including our own). Since photoionization by stars (mainly O stars, young stars or evolved low-mass stars depending on the kind of galaxy) does not seem to be exclusively responsible for the ionization of the hydrogen filaments that should otherwise cool fast and recombine quickly, the question arises which extra energy source can produce the quasi-stationary ionization. We show that stationary localized magnetic reconnection in current filaments may contribute to the ionization of the extraplanar halo gas. In these filaments magnetic energy is dissipated. Consequently, the ionized as well as the neutral component is heated and re-ionized on a time-scale significantly shorter than the recombination time-scale. The amount of energy required for efficient re-ionization can in principle easily be provided by the free magnetic energy. We present quasi-static models that are characterized by plasma temperatures and densities that agree well with the observed values for the diffuse ionized gas component of the interstellar medium. Plasma–neutral gas fluid simulations are made to show that the recombination-induced dynamical reconnection process indeed works in a self-regulatory way.     

How big were the first cosmological objects?

We calculate the cooling times at constant density for haloes with virial temperatures from 100 K to  1×10⁵ K  that originate from a 3 σ fluctuation of a CDM power spectrum in three different cosmologies. Our intention is to determine the first objects that can cool to low temperatures, but not to follow their dynamical evolution. We identify two generations of haloes: those with low virial temperatures,   T _vir≲9000 K  that remain largely neutral, and those with larger virial temperatures that become ionized. The lower temperature, lower mass haloes are the first to cool to 75 per cent of their virial temperature. The precise temperature and mass of the first objects are dependent upon the molecular hydrogen (H₂) cooling function and the cosmological model. The higher mass haloes collapse later but, in this paradigm, cool much more efficiently once they have done so, first via electronic transitions and then via molecular cooling: in fact, a greater residual ionization once the haloes cool below 9000 K results in an enhanced H₂ production and hence a higher cooling rate at low temperatures than for the lower mass haloes, so that within our constant-density model it is the former that are the first to cool to really low temperatures. We discuss the possible significance of this result in the context of CDM models in which the shallow slope of the initial fluctuation spectrum on small scales leads to a wide range of halo masses (of differing overdensities) collapsing over a small redshift interval. This 'crosstalk' is sufficiently important that both high- and low-mass haloes collapse during the lifetimes of the massive stars which may be formed at these epochs. Further investigation is thus required to determine which generation of haloes plays the dominant role in early structure formation.     

The temperature structure and pressure balance of magnetic loops in active regions

EUV observations show many active region loops in lines formed at temperatures between 10⁴K and 2×l0⁶K. The brightest loops are associated with flux tubes leading to the umbrae of sunspots. It is shown that the high visibility of certain loops in transition region lines is due principallly to a sharp radial decrease of temperature to chromospheric values toward the loop axis. The plasma density of these cool loops is not significantly greater than in the hot gas immediately surrounding it. Consequently, the internal gas pressure of the cool material is clearly lower. The hot material immediately surrounding the cool loops is generally denser than the external corona by a factor 3–4. When the active region is examined in coronal lines, this hot high pressure plasma shows up as loops that are generally parallel to the cool loops but significantly displaced laterally. In general the loop phenomenon in an active region is the result of temperature variations by two orders of magnitude and density variations of around a factor five between adjacent flux tubes in the corona.     

On the soft X‐ray emission of M82

We present a spatial analysis of the soft X‐ray and Hα emissions from the outflow of the starburst galaxy M82. We find that the two emissions are tightly correlated on various scales. The O VII triplet of M82, as resolved by X‐ray grating observations of XMM‐Newton, is dominated by the forbidden line, inconsistent with the thermal prediction. The O VII triplet also shows some spatial variations. We discuss three possible explanations for the observed O VII triplet, including the charge exchange at interfaces between the hot outflow and neutral cool gas, a collisional non‐equilibrium‐ionization recombining plasma, and resonance scattering (© 2012 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Quantitative Analysis of H2OComa Images Using a Multiscale MHD Model with Detailed Ion Chemistry

The observation of ions created by ionization of cometary gas, either by ground-based observations or byin situmeasurements can give us useful information about the gas production and composition of comets. However, due to the interaction of ions with the magnetized solar wind and their high chemical reactivity, it is not possible to relate measured ion densities (or column densities) directly to the parent gas densities. In order to quantitatively analyze measured ion abundances in cometary comae it is necessary to understand their dynamics and chemistry. We have developed a detailed ion–chemical network of cometary atmospheres. We include production of ions by photo- and electron impact-ionization of a background neutral atmosphere, charge exchange of solar wind ions with cometary atoms/molecules, reactions between ions and molecules, and dissociative recombination of molecular ions with thermal electrons. By combining the ion–chemical network with the three-dimensional plasma flow as computed by a new fully three-dimensional MHD model of cometary plasma environments (Gombosiet al.1996) we are able to compute the density of the major cometary ions everywhere in the coma. The input parameters for our model are the solar wind conditions (density, speed, temperature, magnetic field) and the composition and production rate of the gas. We applied our model to Comet P/Halley in early March 1986, for which the input parameters are reasonably well known. We compare the resulting column density of H₂O⁺with ground-based observations of H₂O⁺from DiSantiet al.(1990). The results of our model are in good agreement with both the spatial distribution and the absolute abundance of H₂O⁺and with their variations with the changing overall water production rate between two days. The results are encouraging that it will be possible to obtain production rates of neutral cometary constituents from observations of their ion products.     

Comments on some possible models of TMC-1

Chemical differences between cores in the dark ridge TMC-1 have been attributed to the cores being in different stages of chemical evolution with those having high NH₃ to cyanopolyyne abundance ratios being the most evolved. We suggest several alternative models including one in which the highest NH₃ to cyanopolyyne abundance ratio obtains in the youngest TMC-1 core; in this and in one of the other models the evolution of the chemistry as depletion increases is supposed to lead to a lower NH₃ to cyanopolyyne ratio. The possibility that the cyanopolyynes exist primarily in an interface between dark core material and the wind of a low mass star is considered; this wind interface model may account for the sharp cyanopolyyne emission gradient on the side of the ridge away from the star. Implosion of the cores by the ram pressure of the wind may have caused them to collapse more rapidly than gravity could and more rapidly than chemistry evolves so that the chemistry reflects a core's state at a lower density.     

The two basic components of the interstellar neutral hydrogen and its relation with the galactic structure

The two basic components of the neutral hydrogen, cool dense clouds merged in a hotter tenuous medium, are studied using 21 cm absorption data of the Parkes Survey. The mean parameters obtained for the typical clouds next to the galactic plane are τ_p = 1.7, velocity half-width=3.3 km s^?1. Their temperatures areT _sc ≥40 K with a meanT _sc =63±12 K and the obtained hot gas density isn _HH=(0.15±0.05) atom cm^?3. Theoretical analysis following Giovanelli and Brown (1973) reveals that the pressure equilibrium condition (n _HH+2n _e )·T _SH≈n _HC·T _sc is compatible with the quoted values if it is assumed that the cosmic abundances in the interstellar medium are below the adopted normal solar abundance. This lack of heavy elements suggests accretion to grains which is consistent with the observed narrow concentration of the dark matter on the galactic layer (≤100 pc halfwidth). The same pressure condition leads to a mean cool cloud density ofn _HC～30 atom cm^?3 and a hot gas temperature ofT _SH～10 500 K. Comparison with data from Hii regions suggests that the cool clouds are somewhat denser and less extensive than such regions. An explanation for it is the expansion that the Hii regions went through in their origin. Comparison with 21 cm emission data shows that the cloud galactic layer is only about a quarter as thick as the hot gas layer. All the present results suggest that only such clouds can be spatially related with the typical I population associated with the spiral structure.     

A multi‐step model for the origin of E3 (enstatite) chondrites

Abstract— It appears that the mineralogy and chemical properties of type 3 enstatite chondrites could have been established by fractionation processes (removal of a refractory component, and depletion of water) in the solar nebula, and by equilibration with nebular gas at low‐to‐intermediate temperatures (approximately 700–950 K). We describe a model for the origin of type 3 enstatite chondrites that for the first time can simultaneously account for the mineral abundances, bulk‐chemistry, and phase compositions of these chondrites by the operation of plausible processes in the solar nebula. This model, which assumes a representative nebular gas pressure of 10^?5 bar, entails three steps: (1) initial removal of 56% of the equilibrium condensed phases in a system of solar composition at 1270 K; (2) an average loss of 80–85% water vapor in the remaining gas; and (3) two different closure temperatures for the condensed phases. The first step involves a “refractory element fractionation” and is needed to account for the overall major element composition of enstatite chondrites, assuming an initial system with a solar composition. The second step, water‐vapor depletion, is needed to stabilize Si‐bearing metal, oldhamite, and niningerite, which are characteristic minerals of the enstatite chondrites. Variations in closure temperatures are suggested by the way in which the bulk chemistry and mineral assemblages of predicted condensates change with temperature, and how these parameters correlate with the observations of enstatite chondrites. In general, most phases in type 3 enstatite chondrites appear to have ceased equilibrating with nebular gas at approximately 900–950 K, except for Fe‐metal, which continued to partially react with nebular gas to temperatures as low as ～700 K.     

Gas-phase models for the evolved planetary nebulae NGC 6781, M4-9 and NGC 7293

We have studied the chemistry of the molecular gas in evolved planetary nebulae. Three pseudo-time-dependent gas-phase models have been constructed for dense (10⁴–10⁵ cm⁻³) and cool ( T ∼15 K) clumpy envelopes of the evolved nebulae NGC 6781, M4-9 and NGC 7293. The three nebulae are modelled as carbon-rich stars evolved from the asymptotic giant branch to the late planetary nebula phase. The clumpy neutral envelopes are subjected to ultraviolet radiation from the central star and X-rays that enhance the rate of ionization in the clumps. With the ionization rate enhanced by four orders of magnitude over that of the ISM, we find that resultant abundances of the species HCN, HNC, HC₃N and SiC₂ are in good agreement with observations, while those of CN, HCO⁺, CS and SiO are in rough agreement. The results indicate that molecular species such as CH, CH₂, CH₂⁺ , HCl, OH and H₂O are anticipated to be highly abundant in these objects.     

Intrinsic properties of carbon stars. II. Spectra,colours, and HR diagram of cool carbon stars

On the basis of the effective temperature scale proposed previously for cool carbon stars (Paper I), other intrinsic properties of them are examined in detail. It is shown that the major spectroscopic properties of cool carbon stars, including those of molecular bands due to polyatomic species (SiC₂, HCN, C₂H₂ etc.), can most consistently be understood on the basis of our new effective temperature scale and the theoretical prediction of chemical equilibrium. Various photometric indices of cool carbon stars also appear to be well correlated with the new effective temperatures. Furthermore, as effective temperatures of some 30 carbon stars are now obtained, the calibration of any photometric index is straightforward, and some examples of such a calibration are given. In general, colour index-effective temperature calibrations for carbon stars are quite different from those for K-M giant stars. It is found that the intrinsic (R —I)₀ colour is nearly the same for N-irregular variables in spite of a considerable spread in effective temperatures, and this fact is used to estimate the interstellar reddening of carbon stars. An observational HR diagram of red giant stars, including carbon stars as well as K-M giant stars, is obtained on the basis of our colour index-effective temperature calibrations and the best estimations of luminosities. It is shown that carbon stars and M giant stars are sharply divided in the HR diagram by a nearly vertical line at aboutT _eff = 3200 K (logT _eff = 3.50) and the carbon stars occupy the upper right region of M giant stars (except for some high luminosity, high temperature J-type stars in the Magellanic Clouds; also Mira variables are not considered). Such an observational HR diagram of red giant stars shows rather a poor agreement with the current stellar evolution models. Especially, a more efficient mixing process in red giant stars, as compared with those ever proposed, is required to explain the formation of carbon stars.     

The effect of particle precipitation events on the neutral and ion chemistry of the middle atmosphere: II. Odd hydrogen

A one dimensional time-dependent model of the neutral and ion chemistry of the middle atmosphere has been used to examine the production of odd hydrogen (H, OH, and HO₂) during charged particle precipitation. At altitudes above about 65 km, odd hydrogen production depends on the ionization rate, and the atomic oxygen and water vapor densities. Odd hydrogen production is shown to exhibit diurnal and other time dependent variations during such an event at these altitudes, and the assumption that two odd hydrogen particles are always produced per ionization is reexamined.     

Chemical evolution of the gas in C-type shocks in dark clouds

A magnetohydrodynamic model of a steady, transverse C-type shock in a dense molecular cloud is presented. A complete gas–grain chemical network is taken into account: the gas-phase chemistry, the adsorption of gas species on dust grains, various desorption mechanisms, the grain surface chemistry, the ion neutralization on dust grains, the sputtering of grain mantles. The population densities of energy levels of ions CI, CII and OI and molecules H₂, CO, H₂O are computed in parallel with the dynamical and chemical rate equations. The large velocity gradient approximation is used in the line radiative transfer calculations. The simulations consist of two steps: (i) modelling of the chemical and thermal evolution of a static molecular cloud and (ii) shock simulations. A comparison is made with the results of publicly available models of similar physical systems.The focus of the paper is on the chemical processing of gas material and ice mantles of dust grains by the shock. Sputtering of ice mantles takes place in the shock region close to the temperature peak of the neutral gas. At high shock speeds, molecules ejected from ice mantles are effectively destroyed in hot gas, and their survival time is low—of the order of dozens of years. After a passage of high-speed C-type shock, a zone of high abundance of atomic hydrogen appears in the cooling postshock gas that triggers formation of complex organic species such as methanol. It is shown that abundances of some complex organic molecules (COMs) in the postshock region can be much higher than in the preshock gas. These results are important for interpretation of observations of COMs in protostellar outflows.     

Spaceborne ultraviolet 251–384 nm spectroscopy of a meteor during the 1997 Leonid shower

Abstract— We used the ultraviolet to visible spectrometers onboard the midcourse space experiment to obtain the first ultraviolet spectral measurements of a bright meteor during the 1997 Leonid shower. The meteor was most likely a Leonid with a brightness of about‐2 magnitude at 100 km altitude. In the region between 251 and 310 nm, the two strongest emission lines are from neutral and ionized magnesium. Ionized Ca lines, indicative of a hot T ? 10 000 K plasma, are not detected. The Mg and Mg+ line intensity ratio alone does not yield the ionization temperature, which can be determined only by assuming the electron density. A typical air plasma temperature of T = 4400 K would imply a very high electron density: n_e = 2.2 times 10¹⁸ m‐³, but at chondritic abundances of Fe/Mg and Si/Mg ? 1. For a more reasonable local‐thermodynamic‐equilibrium (LTE) air plasma electron density, the Mg and Mg+ line ratio implies a less than chondritic Fe/Mg = 0.06 abundance ratio and a cool non‐LTE T = 2830 K ionization temperature for the ablation vapor plasma. The present observations do not permit a choice between these alternatives. The new data provide also the first spectral confirmation of the presence of molecular OH and NO emission in meteor spectra.     

The effect of displaced geomagnetic and geographic poles on the thermospheric neutral winds

A three-dimensional, semi-empirical dynamic model of the neutral thermosphere is used to examine the effect of the displaced geomagnetic and geographic poles on the daily variation of neutral gas motion. The global-scale pressure distribution to drive the neutral gas motion is derived from the empirical model of Jacchia (1965). The ionization distribution is obtained from the Pennsylvania State M.K 1 model ionosphere using the first few longitudinal Fourier coefficients of the ionization distribution. The calculations were made at various latitudes at equinox and solstice and for various values of solar activity. The results show that the calculated neutral winds for the case where the geomagnetic and geographic poles are coincident differ at most only a few per cent from the winds calculated assuming the poles displaced. With the poles coincident, longitude and local time are interchangeable, and one dimension in any dynamic model of the thermosphere may be eliminated.     

Photoionization of cool MHD disk winds

Recent ultraviolet observations point out that there is hot, dense plasma associated with the optical jet in some T Tauri stars. In this contribution, cool MHD disk wind physics is reviewed by means of a self-similar analytical model to analyze whether hot (Te ? 80,000 K) and dense (n_e ? 10⁹ cm^-3) plasma can be produced in disk winds. It is shown that these high densities can only be achieved at the base of the wind where the stellar X-rays radiation field is strong. The propagation of the X-rays radiation through the disk wind is analyzed: a cocoon of photoionized gas is generated around the star. However, it is difficult to foresee how temperatures as high as ～ 5 × 10⁴ can be reached unless a significant fraction of the X-rays radiation is produced by magnetic reconnection at the boundary between the stellar magnetosphere and the accretion disk.     

X-ray temperature spectroscopy of simulated cooling clusters

Results from a large sample of hydrodynamical/N-body simulations of galaxy clusters in a ΛCDM cosmology are used to simulate cluster X-ray observations. The physical modeling of the gas includes radiative cooling, star formation, energy feedback, and metal enrichment that follow from supernova explosions. Mock cluster samples are constructed grouping simulation data according to a number of constraints which would be satisfied by a data set of X-ray measurements of cluster temperatures as expected from Chandra observations. The X-ray spectra from simulated clusters are fitted into different energy bands using the XSPEC mekal model. The biasing of spectral temperatures with respect to mass-weighted temperatures is found to be influenced by two independent processes. The first scale dependency is absent in adiabatic runs and is due to cooling, whose efficiency to transform cold gas into stars is higher for cool clusters and this in turn implies a strong dependency of the spectral versus mass-weighted temperature relation on the cluster mass. The second dependency is due to photon emission because of cool gas which is accreted during merging events and biases the spectral fits. These events have been quantified according to the power ratio method and a robust correlation is found to exist between the spectral bias and the amount of cluster substructure.The shape of the simulated temperature profiles is not universal and it is steeper at the cluster center for cool clusters than for the massive ones. This follows owing to the scale dependency introduced by cooling which implies for cool clusters higher central temperatures, in scaled units, than for massive clusters. The profiles are in good agreement with data in the radial range between ∼0.1r_vir and ∼0.4r_vir; at small radii (r ≲ 0.1r_vir) the cooling runs fail to reproduce the shape of the observed profiles. The fit is improved if one considers a hierarchical merging scenario in which cluster cores can accrete cooler gas through merging with cluster subclumps, though the shape of the temperature profiles is modified in a significant way only in the regime where the mass of the substructure is a large fraction of the cluster mass.     

Structure and chemistry in the hot molecular core G34.3+0.15

New multifrequency spatial and spectral studies of the hot molecular core associated with the ultracompact HII region G34.3+0.15 have demonstrated an extremely rich chemistry in this archetypal hot core and revealed differing spatial structure between certain species which may be a dynamical effect of chemical evolution. The structure of the hot core has been studied with the JCMT in the high excitation J=19-18 and J=18-17 lines of CH₃CN and with the Nobeyama Millimetre Array at 4 arc resolution in the J=6-5 transition. Comparison with a VLA NH₃(3,3) map shows a displacement between peak emission in the two chemical species which is consistent with chemical processing on a time scale comparable to the dynamical time scale of 10⁵ yrs.A 330-360 GHz spectral survey of the hot core with the JCMT has detected 358 spectral lines from at least 46 distinct chemical species, including many typical of shocked chemistry while other species indicate abundances that reflect the chemistry of a previous cold phase. The first unambiguous detection of ethanol in hot gas has been made. Observations of 14 rotational transitions of this molecule yield a temperature of 125 K and column density 2×10¹⁵ cm^–2. This large abundance cannot be made by purely gas-phase processes and it is concluded that ethanol must have formed by grain surface chemistry.     

The collapse of self-gravitating clouds of pure hydrogen

Exploratory models of the collapse of spherical self-gravitating clouds are studied in relation to the problem of the formation of first generation star-systems. The masses which were considered are in the range of 83 to 5.2×10¹⁰ M _⊙. For simplicity, the assumed composition includes hydrogen only, which could be in the form of H, H₂, H⁺ or H^?. Since the physical conditions that might have prevailed in a primeval nebula are not well known, rather simple initial conditions were chosen: The gas starts from rest and has initially a uniform temperature. We consider the case of rather cool (T ₀～100 K) neutral clouds with different initial ionization degrees. Some of the initial density-distributions here considered are uniform while others are decreasing from the center outwards. The assumed initial values for the densities are ～10^?24 g cm^?3, except for one of the models, for which it is ～10^?26 g cm^?3. Several atomic processes within the gas, including physical-chemical reactions and the evaluation of radiative emission coefficients are considered. A system of differential equations is set up in order to evaluate the concentrationsn _H,n _{H 2},n _H ⁺,n _H ^? andn _e as a function of time. The treatment makes possible the study of the cooling and heating properties of the gas. Furthermore, the dynamical, thermal and chemical evolution of the cloud can be followed during the collapse. The computations apply only to the optically thin stages. The models show the importance of a correct evaluation of the chemical reactions and dissipative mechanisms, which cannot be ignored in a realistic treatment of the collapse of self-gravitating clouds. The influence of the initial conditions on the dynamical and thermal properties during evolution are also analysed.