Self-gravity driven instabilities of interfaces in the interstellar medium

In order to understand star formation it is important to understand the dynamics of atomic and molecular clouds in the interstellar medium (ISM). Non-linear hydrodynamic flows are a key component to the ISM. One route by which non-linear flows arise is the onset and evolution of interfacial instabilities. Interfacial instabilities act to modify the interface between gas components at different densities and temperatures. Such an interface may be subject to a host of instabilities, including the Rayleigh–Taylor, Kelvin–Helmholtz, and Richtmyer–Meshkov instabilities. Recently, a new density interface instability was identified. This self-gravity interfacial instability (SGI) causes any displacement of the interface to grow on roughly a free-fall time-scale, even when the perturbation wavelength is much less than the Jeans length. In previous work, we used numerical simulations to confirm the expectations of linear theory and examine the non-linear evolution of the SGI. We now continue our study by generalizing our initial conditions to allow the acceleration due to self-gravity to be non-zero across the interface. We also consider the behaviour of the SGI for perturbation wavelengths near the Jeans wavelength. We conclude that the action of self-gravity across a density interface may play a significant role in the ISM either by fuelling the growth of new instabilities or modifying the evolution of existing instabilities.     

The chemistry of transient microstructure in the diffuse interstellar medium

Transient microstructure in the diffuse interstellar medium (ISM) has been observed towards Galactic and extragalactic sources for decades, usually in lines of atoms and ions, and, more recently, in molecular lines. Evidently, there is a molecular component to the transient microstructure. In this paper, we explore the chemistry that may arise in such microstructure. We use a photodissociation region (PDR) code to model the conditions of relatively high density, low temperature, very low visual extinction and very short elapsed time that are appropriate for these objects. We find that there is a well-defined region of parameter space where detectable abundances of molecular species might be found. The best matching models are those where the interstellar microstructure is young (<100 yr), small (∼100 au) and dense  (>10⁴ cm⁻³)  .     

Deuterium enrichment of the interstellar medium

Despite the low elemental abundance of atomic deuterium in the interstellar medium (ISM), observational evidence suggests that several species, both in the gas phase and in ices, could be heavily fractionated. We explore various aspects of deuterium enrichment by constructing a chemical evolution model in both gaseous and granular phases. Depending on various physical parameters, gases and grains are allowed to interact with each other through the exchange of their chemical species. It is known that HCO⁺ and N₂H⁺ are two abundant gas phase ions in the ISM and, their deuterium fractionation is generally used to predict the degree of ionization in the various regions of a molecular cloud. For a more accurate estimation, we consider the density profile of a collapsing cloud. The radial distributions of important interstellar molecules, along with their deuterated isotopomers, are presented. Quantum chemical simulations are computed to study the effects of isotopic substitution on the spectral properties of these interstellar species. We calculate the vibrational (harmonic) frequencies of the most important deuterated species (neutral and ions). The rotational and distortional constants of these molecules are also computed in order to predict the rotational transitions of these species. We compare vibrational (harmonic) and rotational transitions as computed by us with existing experimental and theoretical results. It is hope that our results will assist observers in detecting several hitherto unobserved deuterated species.     

ISO observations of 3–200 μm emission by three dust populations in an isolated local translucent cloud

We present isophot spectrophotometry of three positions within the isolated high-latitude cirrus cloud G 300.2−16.8, spanning from the near- to far-infrared (NIR to FIR). The positions exhibit contrasting emission spectrum contributions from the unidentified infrared bands (UIBs), very small grains (VSGs) and large classical grains, and both semi-empirical and numerical models are presented. At all three positions, the UIB spectrum shapes are found to be similar and the large grain emission may be fitted by an equilibrium temperature of  ∼17.5 K  . The energy requirements of both the observed emission spectrum and optical scattered light are shown to be satisfied by the incident local interstellar radiation field (ISRF). The FIR emissivity of dust in G 300.2−16.8 is found to be lower than in globules or dense clouds and is even lower than model predictions for dust in the diffuse interstellar medium (ISM). The results suggest physical differences in the ISM mixtures between positions within the cloud, possibly arising from grain coagulation processes.     

Charge transfer reactions at interfaces between neutral gas and plasma: Dynamical effects and X‐ray emission

Charge‐transfer is the main process linking neutrals and charged particles in the interaction regions of neutral (or partly ionized) gas with a plasma. In this paper we illustrate the importance of charge‐transfer with respect to the dynamics and the structure of neutral gas‐plasma interfaces. We consider the following phenomena: (1) the heliospheric interface ‐ region where the solar wind plasma interacts with the partly‐ionized local interstellar medium (LISM) and (2) neutral interstellar clouds embedded in a hot, tenuous plasma such as the million degree gas that fills the so‐called “Local Bubble”. In (1), we discuss several effects in the outer heliosphere caused by charge exchange of interstellar neutral atoms and plasma protons. In (2) we describe the role of charge exchange in the formation of a transition region between the cloud and the surrounding plasma based on a two‐component model of the cloud‐plasma interaction. In the model the cloud consists of relatively cold and dense atomic hydrogen gas, surrounded by hot, low density, fully ionized plasma. We discuss the structure of the cloud‐plasma interface and the effect of charge exchange on the lifetime of interstellar clouds. Charge transfer between neutral atoms and minor ions in the plasma produces X‐ray emission. Assuming standard abundances of minor ions in the hot gas surrounding the cold interstellar cloud, we estimate the X‐ray emissivity consecutive to the charge transfer reactions. Our model shows that the charge‐transfer X‐ray emission from the neutral cloud‐plasma interface may be comparable to the diffuse thermal X‐ray emission from the million degree gas cavity itself (© 2012 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the influence of cooling and heating processes on the Parker instability – II. Numerical simulations

We investigate the Parker instability (PI) influenced by thermal processes in a non-adiabatic, gravitationally stratified interstellar medium and discuss a model including the photoionization heating together with the supplemental heating mechanisms postulated by Reynolds, Haffner and Tufte. A cooling rate due to radiative losses is described by an approximation to the realistic cooling function of Dalgarno and McCray for ionized interstellar gas. An unperturbed initial state of the system simultaneously represents both a magnetohydrostatic and thermal equilibrium, and is thermally stable. We perform a set of 3D numerical magnetohydrodynamic simulations using the zeusmp code. We find that PI developing in the presence of non-adiabatic effects promotes a transition of gas in magnetic valleys to a thermally unstable regime. We find that the region of initially enhanced density due to PI starts to condense more as the result of thermal instability action. The density in this region rises above the classical isothermal limit of two times the equilibrium value at the mid-plane. The maximum density in an evolved system reaches 10–40 times the equilibrium value at the mid-plane, and the structures so formed attain oval shapes. These results lead to the conclusion that PI, operating in the presence of realistic cooling and heating processes, can trigger the formation of dense clouds, which may give rise to giant molecular complexes.     

A simple model for the relationship between star formation and surface density

We investigate the relationship between the star formation rate per unit area and the surface density of the interstellar medium (ISM; the local Kennicutt–Schmitt law) using a simplified model of the ISM and a simple estimate of the star formation rate based on the mass of gas in bound clumps, the local dynamical time-scales of the clumps and an efficiency parameter of around  ε≈ 5  per cent. Despite the simplicity of the approach, we are able to reproduce the observed linear relation between star formation rate and surface density of dense (molecular) gas. We use a simple model for the dependence of H₂ fraction on total surface density to argue why neither total surface density nor the H  i surface density is a good local indicator of star formation rate. We also investigate the dependence of the star formation rate on the depth of the spiral potential. Our model indicates that the mean star formation rate does not depend significantly on the strength of the spiral potential, but that a stronger spiral potential, for a given mean surface density, does result in more of the star formation occurring close to the spiral arms. This agrees with the observation that grand design galaxies do not appear to show a larger degree of star formation compared to their flocculent counterparts.     

A model of supernova feedback in galaxy formation

A model of supernova feedback during disc galaxy formation is developed. The model incorporates infall of cooling gas from a halo, and outflow of hot gas from a multiphase interstellar medium (ISM). The star formation rate is determined by balancing the energy dissipated in collisions between cold gas clouds with that supplied by supernovae in a disc marginally unstable to axisymmetric instabilities. Hot gas is created by thermal evaporation of cold gas clouds in supernova remnants, and criteria are derived to estimate the characteristic temperature and density of the hot component and hence the net mass outflow rate. A number of refinements of the model are investigated, including a simple model of a galactic fountain, the response of the cold component to the pressure of the hot gas, pressure-induced star formation and chemical evolution. The main conclusion of this paper is that low rates of star formation can expel a large fraction of the gas from a dwarf galaxy. For example, a galaxy with circular speed 50 km s¹ can expel 6080 per cent of its gas over a time-scale of 1 Gyr, with a star formation rate that never exceeds 0.1 M yr¹. Effective feedback can therefore take place in a quiescent mode and does not require strong bursts of star formation. Even a large galaxy, such as the Milky Way, might have lost as much as 20 per cent of its mass in a supernova-driven wind. The models developed here suggest that dwarf galaxies at high redshifts will have low average star formation rates and may contain extended gaseous discs of largely unprocessed gas. Such extended gaseous discs might explain the numbers, metallicities and metallicity dispersions of damped Lyman systems.     

Continuum-driven shocks in active galactic nuclei

In this paper, we extend the study of instabilities in flows driven by the radiation pressure of an ionizing continuum to flows that are not plane parallel. It is well known that the plane-parallel instability leads eventually to the formation of continuum-driven shocks backed by a sonic transition. If these structures are thin, we find that they are unstable to a corrugation mode, and evolve to form sharp-peaked triangular profiles. Once this has occurred, the thin-shock approximation is no longer valid.
We study the further development of the shocks by numerical hydrodynamic simulations. The flow tends to break up into numerous discrete bow-shaped components. The speed of these components through the upstream material is almost constant. As a result, the maximal velocity of radiatively driven shocks through the upstream gas may be determined by instabilities rather than by other physical effects. Interactions between gas in the wings of neighbouring bowshocks can, however, form subsequent generations of bowshocks that are faster and more acute than their predecessors.
One likely location where continuum-driven shocks may occur is in the broad-line regions of active nuclei. We discuss the application of our results to such flows.     

Parker instability in a self-gravitating magnetized gas disk: II. Competition between gravitational and convective instabilities

Under influence of external gravity generated by Galactic all components excluding ISM, a magnetized gas disk may experience both Parker and convective instabilities. Growth rate of the convective instability increases with decreasing perturbation wavelength, and the convective motion makes sheet-like structures all over before the Parker instability forms structures of any meaningful size in the disk. Yet the Parker instability is thought to be an ideal route to form large-scale condensations in the Galaxy. In search of a means to curb convective activities in the Galactic ISM disk, the external gravity is replaced by self-gravity as a driving force of the Parker instability and the gravitational instability is invoked to reinforce the Parker instability. Perturbation of interchange mode is known to trigger convective instability in such disk and the one of undular mode to activate the Parker instability, while the gravitational instability can be triggered by both modes. Therefore, the resulting Jeans instability would help the Parker instability to overcome disrupting behavior of the convection. Dynamical properties of the disk can be characterized by ratio α of magnetic to gas pressure, adiabatic exponent γ, scale height H of the ISM, and disk thickness z_a. A linear stability analysis has been done to the disk, and the maximum growth rate of the Parker–Jeans instability is compared with that of the convective instability. The latter may or may not be higher than the former, depending on the disk parameters. The Parker–Jeans instability has chances to override the convective instability, when the disk is thicker than a certain value. In the disk thinner than the critical one, the Jeans instability can always suppress the convection. Since the growth rate of the convective instability is proportional to local gravitational acceleration, thereby in the general Galactic gravity, the convective instability works actively only in upper regions, we expect chaotic features to appear in regions of low density far from Galactic mid-plane.     

Formation of interstellar clouds: Parker instability with phase transitions

We numerically follow the nonlinear evolution of the Parker instability in the presence of phase transitions from a warm to a cold H  i interstellar medium in two spatial dimensions. The nonlinear evolution of the system favours modes that allow the magnetic field lines to cross the galactic plane. Cold H  i clouds form with typical masses  ≃10⁵ M_⊙  , mean densities  ≃20 cm⁻³  , mean magnetic-field strengths  ≃4.3 μG  (rms field strengths  ≃6.4 μG  ), mass-to-flux ratios  ≃0.1–0.3  relative to critical, temperatures  ≃50 K  , (two-dimensional) turbulent velocity dispersions  ≃1.6 km s⁻¹  and separations  ≃500 pc  , in agreement with observations. The maximum density and magnetic-field strength are  ≃10³ cm⁻³  and  ≃20 μG  , respectively. Approximately 60 per cent of all H  i mass is in the warm neutral medium. The cold neutral medium is arranged into sheet-like structures both perpendicular and parallel to the galactic plane, but it is also found almost everywhere in the galactic plane, with the density being highest in valleys of the magnetic field lines. 'Cloudlets' also form whose physical properties are in quantitative agreement with those observed for such objects by Heiles. The nonlinear phase of the evolution takes ≲30 Myr, so that, if the instability is triggered by a nonlinear perturbation such as a spiral density shock wave, interstellar clouds can form within a time suggested by observations.     

银河系化学元素的丰度特征和合成图像(I):星际介质中轻元素丰度的观测及其特征

给出并解释了星际介质中轻元素D，3He，4He和Li的最新观测数据．星际介质中轻元素的丰度观测结果可以用来检验标准大爆炸核合成理论，因此对这些元素的丰度研究具有重要的天体物理意义．到目前为止，轻元素丰度的观测结果基本上支持开放宇宙的观点．根据最新的观测结果，在本地星际介质中D丰度可能存在小尺度不均匀性，而对类星体吸收云的观测表明不同观测者所获得的原初D丰度结果最大差别可达一个量级．如果观测是可靠的，那么在目前的标准大爆炸核合成理论和星系化学演化模型框架下还不能解释这种结果．另外种种迹象表明太阳系丰度可能不代表45亿年前本地星际介质的丰度．     

AMRVAC and relativistic hydrodynamic simulations for gamma-ray burst afterglow phases

We apply a novel adaptive mesh refinement (AMR) code, AMRVAC (Adaptive Mesh Refinement version of the Versatile Advection Code), to numerically investigate the various evolutionary phases in the interaction of a relativistic shell with its surrounding cold interstellar medium (ISM). We do this for both 1D isotropic and full 2D jet-like fireball models. This is relevant for gamma-ray bursts (GRBs), and we demonstrate that, thanks to the AMR strategy, we resolve the internal structure of the shocked shell–ISM matter, which will leave its imprint on the GRB afterglow. We determine the deceleration from an initial Lorentz factor  γ= 100  up to the almost Newtonian     phase of the flow. We present axisymmetric 2D shell evolutions, with the 2D extent characterized by their initial opening angle. In such jet-like GRB models, we discuss the differences with the 1D isotropic GRB equivalents. These are mainly due to thermally induced sideways expansions of both the shocked shell and shocked ISM regions. We found that the propagating 2D ultrarelativistic shell does not accrete all the surrounding medium located within its initial opening angle. Part of this ISM matter gets pushed away laterally and forms a wide bow-shock configuration with swirling flow patterns trailing the thin shell. The resulting shell deceleration is quite different from that found in isotropic GRB models. As long as the lateral shell expansion is merely due to ballistic spreading of the shell, isotropic and 2D models agree perfectly. As thermally induced expansions eventually lead to significantly higher lateral speeds, the 2D shell interacts with comparably more ISM matter and decelerates earlier than its isotropic counterpart.     

Interstellar deuteroammonia

Close to 30 deuterated molecules have now been detected in the ISM, including doubly-deuterated species D₂H⁺, ND₂H, D₂CO, CHD₂OH, D₂S, and D₂CS, as well as triply-deuterated ammonia and methanol. We review the current understanding of depletion and deuteration processes in cold, dense interstellar medium (ISM) and discuss the utility of deuteroammonia as a tracer of the physical conditions and kinematics of cold, dense gas.     

Thermal instability in a weakly ionized plasma

We revisit the problem of clump formation due to thermal instabilities in a weakly ionized plasma with the help of a linear perturbation analysis, as discussed by Nejad-Asghar & Ghanbari. In the absence of a magnetic field and ambipolar diffusion the characteristic equation reduces to the thermal instability described by Field. We derive the critical wavelengths, which separate the spatial ranges of stability and instability. Contrary to the original analysis of Nejad-Asghar & Ghanbari, perturbations with a wavelength larger than the critical wavelength destabilize the cloud. Moreover, the instability regime of isentropic perturbations is drastically reduced. Isobaric modes with real values of the critical wavelength appear only if the density dependence of the cooling rate is more pronounced than the temperature dependence. Isentropic modes arise only if the power of the density in the cooling rate is smaller than 1/2, which is not fulfilled for CO cooling. We find that ambipolar diffusion is not a dominating heating process in molecular gas.     

Millimetre/submillimetre-wave emission-line searches for high-redshift galaxies

The redshifted spectral line radiation emitted from both atomic fine-structure and molecular rotational transitions in the interstellar medium (ISM) of high-redshift galaxies can be detected in the centimetre, millimetre and submillimetre wavebands. Here we predict the counts of galaxies detectable in an array of molecular and atomic lines. This calculation requires a reasonable knowledge of both the surface density of these galaxies on the sky, and the physical conditions in their ISM. The surface density is constrained using the results of submillimetre-wave continuum surveys. Follow-up OVRO Millimeter Array observations of two of the galaxies detected in the dust continuum have provided direct measurements of CO rotational line emission at redshifts of 2.56 and 2.81. Based on these direct high-redshift observations and on models of the ISM that are constrained by observations of low-redshift ultraluminous infrared galaxies, we predict the surface density of line-emitting galaxies as a function of line flux density and observing frequency. We incorporate the sensitivities and mapping speeds of existing and future millimetre/submillimetre-wave telescopes and spectrographs, and so assess the prospects for blank-field surveys to detect this line emission from gas-rich high-redshift galaxies.     

Star formation in a multi-phase ISM

We present a 3d code for the dynamical evolution of a multi-phase interstellar medium (ISM) coupled to stars via star formation (SF) and feedback processes. The multi-phase ISM consists of clouds (sticky particles) and diffuse gas (SPH): exchange of matter, energy and momentum is achieved by drag (due to ram pressure) and condensation or evaporation processes. The cycle of matter is completed by SF and feedback by SNe and PNe. A SF scheme based on a variable SF efficiency as proposed by Elmegreen and Efremov (1997) is presented. For a Milky Way type galaxy we get a SF rate of ∼1 M_⊙ yr^-1 with an average SF efficiency of ∼5%. This revised version was published online in August 2006 with corrections to the Cover Date.     

The orientations of molecular clouds in the outer Galaxy: evidence for the scale of the turbulence driver?

Supernova (SN) explosions inject a considerable amount of energy into the interstellar medium (ISM) in regions with high-to-moderate star formation rates. In order to assess whether the driving of turbulence by supernovae is also important in the outer Galactic disc, where the star formation rates are lower, we study the spatial distribution of molecular cloud (MC) inclinations with respect to the Galactic plane. The latter contains important information on the nature of the mechanism of energy injection into the ISM. We analyse the spatial correlations between the position angles (PAs) of a selected sample of MCs (the largest clouds in the catalogue of the outer Galaxy published by Heyer et al). Our results show that when the PAs of the clouds are all mapped to values into the  [0°, 90°]  interval, there is a significant degree of spatial correlation between the PAs on spatial scales in the range of 100–800 pc. These scales are of the order of the sizes of individual SN shells in low-density environments such as those prevailing in the outer Galaxy and where the metallicity of the ambient gas is of the order of the solar value or smaller. These findings suggest that individual SN explosions, occurring in the outer regions of the Galaxy and in likewise spiral galaxies, albeit at lower rates, continue to play an important role in shaping the structure and dynamics of the ISM in those regions. The SN explosions we postulate here are likely associated with the existence of young stellar clusters in the far outer regions of the Galaxy and the ultraviolet emission and low levels of star formation observed with the Galaxy Evolution Explorer (GALEX) satellite in the outer regions of local galaxies.     

Evaporation and condensation of giant interstellar clouds in a hot-gas environment

Gas phases of the interstellar medium (ISM) coexist locally, penetrate each other and mix by means of dynamical and plasmaphysical processes. E.g. heat conduction from the hot to the cooler gas leads to energy and mass exchange between the gas phases. Analytical solutions exist under which evaporation of cloudy material or condensation of hot gas onto the clouds' surface dominate. Since these results are derived for stationary and static conditions and under ideal assumptions, they do not necessarily hold for a dynamical ISM. On the other hand, the mass and energy exchange between the gas phases is of great importance for the energy budget of the ISM and by this influences the evolution of galaxies. This led us to investigate the evolution of interstellar clouds in a hot gas by means of numerical simulations. At first, we compare static models with the analytical results and found that interstellar clouds with parameters requiring analytically evaporation are, in contrast, accreting surrounding material if self-gravitation and cooling are implied. For the more realistic case, where clouds are embedded in a streaming hot gas, the models show that Kelvin-Helmholtz instability which leads to the disruption of the clouds is suppressed by heat conduction so that the clouds are stabilized to survive. This revised version was published online in September 2006 with corrections to the Cover Date.