Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA model for atomic and molecular interstellar gas: The Meudon PDR code
108
0
0.0
(
0
)
Authors
Franck Le Petit
Ask ChatGPT about the research
No Arabic abstract
We present the revised ``Meudon model of Photon Dominated Region (PDR code), presently available on the web under the Gnu Public Licence at: http://aristote.obspm.fr/MIS. General organisation of the code is described down to a level that should allow most observers to use it as an interpretation tool with minimal help from our part. Two grids of models, one for low excitation diffuse clouds and one for dense highly illuminated clouds, are discussed, and some new results on PDR modelisation highlighted.
rate research
Read More
Our main purpose is to estimate the effect of assuming uniform density on the line-of-sight in PDR chemistry models, compared to a more realistic distribution for which total gas densities may well vary by several orders of magnitude. A secondary goal of this paper is to estimate the amount of molecular hydrogen which is not properly traced by the CO (J = 1 -> 0) line, the so-called dark molecular gas. We use results from a magnetohydrodynamical (MHD) simulation as a model for the density structures found in a turbulent diffuse ISM with no star-formation activity. The Meudon PDR code is then applied to a number of lines of sight through this model, to derive their chemical structures. It is found that, compared to the uniform density assumption, maximal chemical abundances for H2, CO, CH and CN are increased by a factor 2 to 4 when taking into account density fluctuations on the line of sight. The correlations between column densities of CO, CH and CN with respect to those of H2 are also found to be in better overall agreement with observations. For instance, at N(H2) > 2.10^{20} cm-2, while observations suggest that d[log N(CO)]=d[log N(H2)] = 3.07 +/- 0.73, we find d[log N(CO)]=d[log N(H2)] =14 when assuming uniform density, and d[log N(CO)]=d[log N(H2)] = 5.2 when including density fluctuations.
In this contribution, we test our previously published one-dimensional PDR model for deriving total hydrogen volume densities from HI column density measurements in extragalactic regions by applying it to the Taurus molecular cloud, where its predictions can be compared to available data. Also, we make the first direct detailed comparison of our model to CO(1-0) and far-infrared emission. Using an incident UV flux G0 of 4.25 ({chi} = 5) throughout the main body of the cloud, we derive total hydrogen volume densities of approx 430 cm-3, consistent with the extensive literature available on Taurus. The distribution of the volume densities shows a log-normal shape with a hint of a power-law shape on the high density end. We convert our volume densities to H2 column densities assuming a cloud depth of 5 parsec and compare these column densities to observed CO emission. We find a slope equivalent to a CO conversion factor relation that is on the low end of reported values for this factor in the literature (0.9 x 1020 cm-2 (K km s-1)-1), although this value is directly proportional to our assumed value of G0 as well as the cloud depth. We seem to under-predict the total hydrogen gas as compared to 100 {mu}m dust emission, which we speculate may be caused by a higher actual G0 incident on the Taurus cloud than is generally assumed.
E(A+M)PEC traces the ionization structure, cooling and emission spectra of plasmas. It is written in OpenCL, runs in NVIDIA Graphics Processor Units and can be coupled to any HD or MHD code to follow the dynamical and thermal evolution of any plasma in, e.g., the interstellar medium (ISM).
We present a comparison between independent computer codes, modeling the physics and chemistry of interstellar photon dominated regions (PDRs). Our goal was to understand the mutual differences in the PDR codes and their effects on the physical and chemical structure of the model clouds, and to converge the output of different codes to a common solution. A number of benchmark models have been created, covering low and high gas densities and far ultraviolet intensities. The benchmark models were computed in two ways: one set assuming constant temperatures, thus testing the consistency of the chemical network and photo-processes, and a second set determining the temperature selfconsistently. We investigated the impact of PDR geometry and agreed on the comparison of results from spherical and plane-parallel PDR models. We identified a number of key processes governing the chemical network which have been treated differently in the various codes, and defined a proper common treatment. We established a comprehensive set of reference models for ongoing and future PDR model bench-marking and were able to increase the agreement in model predictions for all benchmark models significantly.
Recent Herschel and ALMA observations of Photodissociation Regions (PDRs) have revealed the presence of a high thermal pressure (P ~ 10^7-10^8 K cm-3) thin compressed layer at the PDR surface where warm molecular tracer emission (e.g. CH+, SH+, high-J CO, H2,...) originate. These high pressures (unbalanced by the surrounding environment) and a correlation between pressure and incident FUV field (G0) seem to indicate a dynamical origin with the radiation field playing an important role in driving the dynamics. We investigate whether photoevaporation of the illuminated edge of a molecular cloud could explain these high pressures and pressure-UV field correlation. We developed a 1D hydrodynamical PDR code coupling hydrodynamics, EUV and FUV radiative transfer and time-dependent thermo-chemical evolution. We applied it to a 1D plane-parallel photoevaporation scenario where a UV-illuminated molecular cloud can freely evaporate in a surrounding low-pressure medium. We find that photoevaporation can produce high thermal pressures and the observed P-G0 correlation, almost independently from the initial gas density. In addition, we find that constant-pressure PDR models are a better approximation to the structure of photoevaporating PDRs than constant-density PDR models, although moderate pressure gradients are present. Strong density gradients from the molecular to the neutral atomic region are found, which naturally explain the large density contrasts (1-2 orders of magnitude) derived from observations of different tracers. The photoevaporating PDR is preceded by a low velocity shock (a few km/s) propagating into the molecular cloud. Photoevaporating PDR models offer a promising explanation to the recent observational evidence of dynamical effects in PDRs.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
