Source author record

J. Le Bourlot

J. Le Bourlot appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2025arXiv

JWST observations of photodissociation regions: II. Warm molecular Hydrogen spectroscopy in the Horsehead nebula

H2 is the most abundant molecule in the interstellar medium and is a useful tool to study photodissociation regions, where radiative feedback from massive stars on molecular clouds is dominant. The James Webb Space Telescope, with its high spatial resolution, sensitivity, and wavelength coverage provides unique access to the detection of most of H2 lines and the analysis of its spatial morphology. Our goal is to use H2 line emission detected with the JWST in the Horsehead nebula to constrain the physical parameters (e.g., extinction, gas temperature, thermal pressure) throughout the PDR and its geometry. The study of H2 morphology reveals that FUV-pumped lines peak closer to the edge of the PDR than thermalized lines. From H2 lines, we estimate the value of extinction throughout the PDR. We find that AV is increasing from the edge of the PDR to the second and third H2 filaments. Then, we study the H2 excitation in different regions across the PDR. The temperature profile shows that the observed gas temperature is quite constant throughout the PDR, with a slight decline in each of the dissociation fronts. This study also reveals that the OPR is far from equilibrium. We observe a spatial separation of para and ortho rovibrational levels, indicating that efficient ortho-para conversion and preferential ortho self-shielding are driving the spatial variations of the OPR. Finally, we derive a thermal pressure in the first filament around P > 6x10$^6$ K cm$^{-3}$, about ten times higher than that of the ionized gas. We highlight that template stationary 1D PDR models cannot account for the intrinsic 2D structure and the very high temperature observed in the Horsehead nebula. We argue the highly excited, over-pressurized H2 gas at the edge of the PDR interface could originate from the mixing between the cold and hot phase induced by the photo-evaporation of the cloud.

preprint2016arXiv

A milestone toward understanding PDR properties in the extreme environment of LMC-30Dor

More complete knowledge of galaxy evolution requires understanding the process of star formation and interaction between the interstellar radiation field and the interstellar medium in galactic environments traversing a wide range of physical parameter space. Here we focus on the impact of massive star formation on the surrounding low metallicity ISM in 30 Doradus in the Large Magellanic Cloud. A low metal abundance, as is the case of some galaxies of the early universe, results in less ultra-violet shielding for the formation of the molecular gas necessary for star formation to proceed. The half-solar metallicity gas in this region is strongly irradiated by the super star cluster R136, making it an ideal laboratory to study the structure of the ISM in an extreme environment. Our spatially resolved study investigates the gas heating and cooling mechanisms, particularly in the photo-dissociation regions where the chemistry and thermal balance are regulated by far-ultraviolet photons (6 eV< hν<13.6 eV). We present Herschel observations of far-infrared fine-structure lines obtained with PACS and SPIRE/FTS. We have combined atomic fine-structure lines from Herschel and Spitzer observations with ground-based CO data to provide diagnostics on the properties and the structure of the gas by modeling it with the Meudon PDR code. We derive the spatial distribution of the radiation field, the pressure, the size, and the filling factor of the photodissociated gas and molecular clouds. We find a range of pressure of ~ 10^5 - 1.7x10^6 cm^{-3} K and a range of incident radiation field G_UV ~ 10^2 - 2.5x10^4 through PDR modeling. Assuming a plane-parallel geometry and a uniform medium, we find a total extinction of 1-3 mag , which correspond to a PDR cloud size of 0.2 to 3pc, with small CO depth scale of 0.06 to 0.5pc. We also determine the three dimensional structure of the gas. (Abridged)

preprint2015arXiv

$Herschel$ SPIRE-FTS observations of RCW 120

The expansion of Galactic HII regions can trigger the formation of a new generation of stars. However, little is know about the physical conditions that prevail in these regions. We study the physical conditions that prevail in specific zones towards expanding HII regions that trace representative media such as the photodissociation region, the ionized region, and condensations with and without ongoing star formation. We use the SPIRE Fourier Transform Spectrometer (FTS) on board $Herschel$ to observe the HII region RCW 120. Continuum and lines are observed in the $190-670\,μ$m range. Line intensities and line ratios are obtained and used as physical diagnostics of the gas. We used the Meudon PDR code and the RADEX code to derive the gas density and the radiation field at nine distinct positions including the PDR surface and regions with and without star-formation activity. For the different regions we detect the atomic lines [NII] at $205\,μ$m and [CI] at $370$ and $609\,μ$m, the $^{12}{\rm CO}$ ladder between the $J=4$ and $J=13$ levels and the $^{13}{\rm CO}$ ladder between the $J=5$ and $J=14$ levels, as well as CH$ ^{+} $ in absorption. We find gas temperatures in the range $45-250\,$K for densities of $10^4-10^6\,{\rm cm}^{-3}$, and a high column density on the order of $N_{\rm H}\sim10^{22}\,{\rm cm}^{-2}$ that is in agreement with dust analysis. The ubiquitousness of the atomic and CH$ ^{+} $ emission suggests the presence of a low-density PDR throughout RCW 120. High-excitation lines of CO indicate the presence of irradiated dense structures or small dense clumps containing young stellar objects, while we also find a less dense medium ($N_{\rm H}\sim10^{20}\,{\rm cm}^{-2}$) with high temperatures ($80-200\,$K).

preprint2014arXiv

Photodissociation of interstellar ArH+

Aims. Following the recent detection of 36ArH+ in the Crab nebula spectrum, we have computed the photodissociation rate of ArH+ in order to constrain the physical processes at work in this environment. Methods. Photodissociation cross sections of ArH+ are computed in an ab initio approach including explicit account of spin-orbit coupling. Results. We report the photodissociation cross section of ArH+ as a function of wavelength. Photodissociation probabilities are derived for different impinging radiation fields.The photodissociation probability of for a very small unshielded cloud surrounded on all sides by the unshielded InterStellar Radiation Field (ISRF) model described by Draine (1978) is equal to 9.9e-12 s-1 and 1.9e-9 s-1 in the Crab nebula conditions. The dependence on the visual extinction is obtained by using the Meudon Photon Dominated Region (PDR) code and corresponding analytical fits are provided. Conclusions. These data will help to produce a realistic chemical network to interpret the observations. Photodissociation of ArH+ is found to be moderate and the presence of this molecular ion is mainly dependent on the molecular fraction

preprint2013arXiv

Spatial distribution of small hydrocarbons in the neighborhood of the Ultra Compact HII region Monoceros R2

We study the chemistry of small hydrocarbons in the photon-dominated regions (PDRs) associated with the ultra-compact HII region Mon R2. Our goal is to determine the variations of the abundance of small hydrocarbons in a high-UV irradiated PDR and investigate their chemistry. We present an observational study of CH, CCH and c-C$_3$H$_2$ in Mon R2 combining data obtained with the IRAM 30m telescope and Herschel. We determine the column densities of these species, and compare their spatial distributions with that of polycyclic aromatic hydrocarbon (PAH). We compare the observational results with different chemical models to explore the relative importance of gas-phase, grain-surface and time-dependent chemistry in these environments. The emission of the small hydrocarbons show different patterns. The CCH emission is extended while CH and c-C$_3$H$_2$ are concentrated towards the more illuminated layers of the PDR. The ratio of the column densities of c-C$_3$H$_2$ and CCH shows spatial variations up to a factor of a few, increasing from $N(c-C$_3$H$_2$)/N(CCH)\approx0.004$ in the envelope to a maximum of $\sim0.015-0.029$ towards the 8$μ$m emission peak. Comparing these results with other galactic PDRs, we find that the abundance of CCH is quite constant over a wide range of G$_0$, whereas the abundance of c-C$_3$H$_2$ is higher in low-UV PDRs. In Mon R2, the gas-phase steady-state chemistry can account relatively well for the abundances of CH and CCH in the most exposed layers of the PDR, but falls short by a factor of 10 to reproduce c-C$_3$H$_2$. In the molecular envelope, time-dependent effects and grain surface chemistry play a dominant role in determining the hydrocarbons abundances. Our study shows that CCH and c-C$_3$H$_2$ present a complex chemistry in which UV photons, grain-surface chemistry and time dependent effects contribute to determine their abundances.

preprint2013arXiv

The chemistry of ions in the Orion Bar I. - CH+, SH+, and CF+: The effect of high electron density and vibrationally excited H2 in a warm PDR surface

The abundances of interstellar CH+ and SH+ are not well understood as their most likely formation channels are highly endothermic. Using data from Herschel, we study the formation of CH+ and SH+ in a typical high UV-illumination photon-dominated region (PDR), the Orion Bar. Herschel/HIFI provides velocity-resolved data of CH+ 1-0 and 2-1 and three hyperfine transitions of SH+. Herschel/PACS provides information on the excitation and spatial distribution of CH+ (up to J=6-5). The widths of the CH+ 2-1 and 1-0 transitions are of ~5 km/s, significantly broader than the typical width of dense gas tracers in the Orion Bar (2-3 km/s) and are comparable to the width of tracers of the interclump medium such as C+ and HF. The detected SH+ transitions are narrower compared to CH+ and have line widths of 3 km/s, indicating that SH+ emission mainly originates in denser condensations. Non-LTE radiative transfer models show that electron collisions affect the excitation of CH+ and SH+, and that reactive collisions need to be taken into account to calculate the excitation of CH+. Comparison to PDR models shows that CH+ and SH+ are tracers of the warm surface region (AV<1.5) of the PDR with temperatures between 500-1000 K. We have also detected the 5-4 transition of CF+ (FWHM=1.9 km/s) with an intensity that is consistent with previous observations of lower-J CF+ transitions toward the Orion Bar. A comparison to PDR models indicate that the internal vibrational energy of H2 can explain the formation of CH+ for typical physical conditions in the Orion Bar near the ionization front. H2 vibrational excitation is the most likely explanation of SH+ formation as well. The abundance ratios of CH+ and SH+ trace the destruction paths of these ions, and through that, indirectly, the ratios of H, H2 and electron abundances as a function of depth into the cloud.

preprint2013arXiv

The IRAM-30m line survey of the Horsehead PDR: IV. Comparative chemistry of H2CO and CH3OH

Aims. We investigate the dominant formation mechanism of H2CO and CH3OH in the Horsehead PDR and its associated dense core. Methods. We performed deep integrations of several H2CO and CH3OH lines at two positions in the Horsehead, namely the PDR and dense core, with the IRAM-30m telescope. In addition, we observed one H2CO higher frequency line with the CSO telescope at both positions. We determine the H2CO and CH3OH column densities and abundances from the single-dish observations complemented with IRAM-PdBI high-angular resolution maps (6") of both species. We compare the observed abundances with PDR models including either pure gas-phase chemistry or both gas-phase and grain surface chemistry. Results. We derive CH3OH abundances relative to total number of hydrogen atoms of ~1.2e-10 and ~2.3e-10 in the PDR and dense core positions, respectively. These abundances are similar to the inferred H2CO abundance in both positions (~2e-10). We find an abundance ratio H2CO/CH3OH of ~2 in the PDR and ~1 in the dense core. Pure gas-phase models cannot reproduce the observed abundances of either H2CO or CH3OH at the PDR position. Both species are therefore formed on the surface of dust grains and are subsequently photodesorbed into the gas-phase at this position. At the dense core, on the other hand, photodesorption of ices is needed to explain the observed abundance of CH3OH, while a pure gas-phase model can reproduce the observed H2CO abundance. The high-resolution observations show that CH3OH is depleted onto grains at the dense core. CH3OH is thus present in an envelope around this position, while H2CO is present in both the envelope and the dense core itself. Conclusions. Photodesorption is an efficient mechanism to release complex molecules in low FUV-illuminated PDRs, where thermal desorption of ice mantles is ineffective.

preprint2012arXiv

Herschel / HIFI observations of CO, H2O and NH3 in Mon R2

Context. Mon R2 is the only ultracompact HII region (UCHII) where the associated photon-dominated region (PDR) can be resolved with Herschel. Due to its brightness and proximity, it is the best source to investigate the chemistry and physics of highly UV-irradiated PDRs. Aims. Our goal is to estimate the abundance of H2O and NH3 in this region and investigate their origin. Methods. We present new observations obtained with HIFI and the IRAM-30m telescope. Using a large velocity gradient approach, we model the line intensities and derive an average abundance of H2O and NH3 across the region. Finally, we model the line profiles with a non-local radiative transfer model and compare these results with the abundance predicted by the Meudon PDR code. Results. The variations of the line profiles and intensities indicate complex geometrical and kinematical patterns. The H2O lines present a strong absorption at the ambient velocity and emission in high velocity wings towards the HII region. The spatial distribution of the o-H2^18O line shows that the its emission arises in the PDR surrounding the HII region. By modeling the o-H2^18O emission we derive a mean abundance of o-H2O of ~10^-8 relative to H2. The ortho-H2O abundance is however larger, ~1x10^-7, in the high velocity wings. Possible explanations for this larger abundance include an expanding hot PDR and/or an outflow. Ammonia seems to be present only in the envelope with an average abundance of ~2x10^-9 relative to H2. Conclusions. The Meudon PDR code can account for the measured water abundance in the high velocity gas as long as we assume that it originates from a <1 mag hot expanding layer of the PDR, i.e. that the outflow has only a minor contribution to this emission. To explain the abundances in the rest of the cloud the molecular freeze out and grain surface chemistry would need to be included.

preprint2012arXiv

Multi-line detection of O2 toward rho Oph A

Models of pure gas-phase chemistry in well-shielded regions of molecular clouds predict relatively high levels of molecular oxygen, O2, and water, H2O. Contrary to expectation, the space missions SWAS and Odin found only very small amounts of water vapour and essentially no O2 in the dense star-forming interstellar medium. Only toward rho Oph A did Odin detect a weak line of O2 at 119 GHz in a beam size of 10 arcmin. A larger telescope aperture such as that of the Herschel Space Observatory is required to resolve the O2 emission and to pinpoint its origin. We use the Heterodyne Instrument for the Far Infrared aboard Herschel to obtain high resolution O2 spectra toward selected positions in rho Oph A. These data are analysed using standard techniques for O2 excitation and compared to recent PDR-like chemical cloud models. The 487.2GHz line was clearly detected toward all three observed positions in rho Oph A. In addition, an oversampled map of the 773.8GHz transition revealed the detection of the line in only half of the observed area. Based on their ratios, the temperature of the O2 emitting gas appears to vary quite substantially, with warm gas (> 50 K) adjacent to a much colder region, where temperatures are below 30 K. The exploited models predict O2 column densities to be sensitive to the prevailing dust temperatures, but rather insensitive to the temperatures of the gas. In agreement with these model, the observationally determined O2 column densities seem not to depend strongly on the derived gas temperatures, but fall into the range N(O2) = (3 to >6)e15/cm^2. Beam averaged O2 abundances are about 5e-8 relative to H2. Combining the HIFI data with earlier Odin observations yields a source size at 119 GHz of about 4 - 5 arcmin, encompassing the entire rho Oph A core.

preprint2012arXiv

The IRAM-30m line survey of the Horsehead PDR: II. First detection of the l-C3H+ hydrocarbon cation

We present the first detection of the l-C3H+ hydrocarbon in the interstellar medium. The Horsehead WHISPER project, a millimeter unbiased line survey at two positions, namely the photo-dissociation region (PDR) and the nearby shielded core, revealed a consistent set of eight unidentified lines toward the PDR position. Six of them are detected with a signal-to-noise ratio from 6 to 19, while the two last ones are tentatively detected. Mostly noise appears at the same frequency toward the dense core, located less than 40" away. We simultaneously fit 1) the rotational and centrifugal distortion constants of a linear rotor, and 2) the Gaussian line shapes located at the eight predicted frequencies. The observed lines can be accurately fitted with a linear rotor model, implying a 1Sigma ground electronic state. The deduced rotational constant value is Be= 11244.9512 +/- 0.0015 MHz, close to that of l-C3H. We thus associate the lines to the l-C3H+ hydrocarbon cation, which enables us to constrain the chemistry of small hydrocarbons. A rotational diagram is then used to infer the excitation temperature and the column density. We finally compare the abundance to the results of the Meudon PDR photochemical model.

preprint2011arXiv

OH emission from warm and dense gas in the Orion Bar PDR

As part of a far-infrared (FIR) spectral scan with Herschel/PACS, we present the first detection of the hydroxyl radical (OH) towards the Orion Bar photodissociation region (PDR). Five OH rotational Lambda-doublets involving energy levels out to E_u/k~511 K have been detected (at ~65, ~79, ~84, ~119 and ~163um). The total intensity of the OH lines is I(OH)~5x10^-4 erg s^-1 cm^-2 sr^-1. The observed emission of rotationally excited OH lines is extended and correlates well with the high-J CO and CH^+ J=3-2 line emission (but apparently not with water vapour), pointing towards a common origin. Nonlocal, non-LTE radiative transfer models including excitation by the ambient FIR radiation field suggest that OH arises in a small filling factor component of warm (Tk~160-220 K) and dense (n_H~10^{6-7} cm^-3) gas with source-averaged OH column densities of ~10^15 cm^-2. High density and temperature photochemical models predict such enhanced OH columns at low depths (A_V<1) and small spatial scales (~10^15 cm), where OH formation is driven by gas-phase endothermic reactions of atomic oxygen with molecular hydrogen. We interpret the extended OH emission as coming from unresolved structures exposed to far-ultraviolet (FUV) radiation near the Bar edge (photoevaporating clumps or filaments) and not from the lower density "interclump" medium. Photodissociation leads to OH/H2O abundance ratios (>1) much higher than those expected in equally warm regions without enhanced FUV radiation fields.

preprint2010arXiv

Gas morphology and energetics at the surface of PDRs: new insights with Herschel observations of NGC 7023

We investigate the physics and chemistry of the gas and dust in dense photon-dominated regions (PDRs), along with their dependence on the illuminating UV field. Using Herschel-HIFI observations, we study the gas energetics in NGC 7023 in relation to the morphology of this nebula. NGC 7023 is the prototype of a PDR illuminated by a B2V star and is one of the key targets of Herschel. Our approach consists in determining the energetics of the region by combining the information carried by the mid-IR spectrum (extinction by classical grains, emission from very small dust particles) with that of the main gas coolant lines. In this letter, we discuss more specifically the intensity and line profile of the 158 micron (1901 GHz) [CII] line measured by HIFI and provide information on the emitting gas. We show that both the [CII] emission and the mid-IR emission from polycyclic aromatic hydrocarbons (PAHs) arise from the regions located in the transition zone between atomic and molecular gas. Using the Meudon PDR code and a simple transfer model, we find good agreement between the calculated and observed [CII] intensities. HIFI observations of NGC 7023 provide the opportunity to constrain the energetics at the surface of PDRs. Future work will include analysis of the main coolant line [OI] and use of a new PDR model that includes PAH-related species.

preprint2010arXiv

Herschel/HIFI discovery of interstellar chloronium (H$_2$Cl$^+$)

We report the first detection of chloronium, H$_2$Cl$^+$, in the interstellar medium, using the HIFI instrument aboard the \emph{Herschel} Space Observatory. The $2_{12}-1_{01}$ lines of ortho-H$_2^{35}$Cl$^+$ and ortho-H$_2^{37}$Cl$^+$ are detected in absorption towards NGC~6334I, and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ is detected in absorption towards NGC~6334I and Sgr~B2(S). The H$_2$Cl$^+$ column densities are compared to those of the chemically-related species HCl. The derived HCl/H$_2$Cl$^+$ column density ratios, $\sim$1--10, are within the range predicted by models of diffuse and dense Photon Dominated Regions (PDRs). However, the observed H$_2$Cl$^+$ column densities, in excess of $10^{13}$~cm$^{-2}$, are significantly higher than the model predictions. Our observations demonstrate the outstanding spectroscopic capabilities of HIFI for detecting new interstellar molecules and providing key constraints for astrochemical models.

preprint2010arXiv

Sensitivity analyses of dense cloud chemical models

Because of new telescopes that will dramatically improve our knowledge of the interstellar medium, chemical models will have to be used to simulate the chemistry of many regions with diverse properties. To make these models more robust, it is important to understand their sensitivity to a variety of parameters. In this article, we report a study of the sensitivity of a chemical model of a cold dense core, with homogeneous and time-independent physical conditions, to variations in the following parameters: initial chemical inventory, gas temperature and density, cosmic-ray ionization rate, chemical reaction rate coefficients, and elemental abundances. From the results of the parameter variations, we can quantify the sensitivity of the model to each parameter as a function of time. Our results can be used in principle with observations to constrain some parameters for different cold clouds. We also attempted to use the Monte Carlo approach with all parameters varied collectively. Within the parameter ranges studied, the most critical parameters turn out to be the reaction rate coefficients at times up to 4e5 yr and elemental abundances at later times. At typical times of best agreement with observation, models are sensitive to both of these parameters. The models are less sensitive to other parameters such as the gas density and temperature. The improvement of models will require that the uncertainties in rate coefficients of important reactions be reduced. As the chemistry becomes better understood and more robust, it should be possible to use model sensitivities concerning other parameters, such as the elemental abundances and the cosmic ray ionization rate, to yield detailed information on cloud properties and history. Nevertheless, at the current stage, we cannot determine the best values of all the parameters simultaneously based on purely observational constraints.

preprint2010arXiv

The global dust SED: Tracing the nature and evolution of dust with DustEM

The Planck and Herschel missions are currently measuring the farIR-mm emission of dust, which combined with existing IR data, will for the first time provide the full SED of the galactic ISM dust emission with an unprecedented sensitivity and angular resolution. It will allow a systematic study of the dust evolution processes that affect the SED. Here we present a versatile numerical tool, DustEM, that predicts the emission and extinction of dust given their size distribution and their optical and thermal properties. In order to model dust evolution, DustEM has been designed to deal with a variety of grain types, structures and size distributions and to be able to easily include new dust physics. We use DustEM to model the dust SED and extinction in the diffuse interstellar medium at high-galactic latitude (DHGL), a natural reference SED. We present a coherent set of observations for the DHGL SED. The dust components in our DHGL model are (i) PAHs, (ii) amorphous carbon and (iii) amorphous silicates. We use amorphous carbon dust, rather than graphite, because it better explains the observed high abundances of gas-phase carbon in shocked regions of the interstellar medium. Using the DustEM model, we illustrate how, in the optically thin limit, the IRAS/Planck HFI (and likewise Spitzer/Herschel for smaller spatial scales) photometric band ratios of the dust SED can disentangle the influence of the exciting radiation field intensity and constrain the abundance of small grains relative to the larger grains. We also discuss the contributions of the different grain populations to the IRAS, Planck and Herschel channels. Such information is required to enable a study of the evolution of dust as well as to systematically extract the dust thermal emission from CMB data and to analyze the emission in the Planck polarized channels. The DustEM code described in this paper is publically available.

preprint2009arXiv

The ionization fraction gradient across the Horsehead edge: An archetype for molecular clouds

The ionization fraction plays a key role in the chemistry and dynamics of molecular clouds. We study the H13CO+, DCO+ and HOC+ line emission towards the Horsehead, from the shielded core to the UV irradiated cloud edge, i.e., the Photodissociation Region (PDR), as a template to investigate the ionization fraction gradient in molecular clouds. We analyze a PdBI map of the H13CO+ J=1-0 line, complemented with IRAM-30m H13CO+ and DCO+ higher-J line maps and new HOC+ and CO+ observations. We compare self-consistently the observed spatial distribution and line intensities with detailed depth-dependent predictions of a PDR model coupled with a nonlocal radiative transfer calculation. The chemical network includes deuterated species, 13C fractionation reactions and HCO+/HOC+ isomerization reactions. The role of neutral and charged PAHs in the cloud chemistry and ionization balance is investigated. The detection of HOC+ reactive ion towards the Horsehead PDR proves the high ionization fraction of the outer UV irradiated regions, where we derive a low [HCO+]/[HOC+]~75-200 abundance ratio. In the absence of PAHs, we reproduce the observations with gas-phase metal abundances, [Fe+Mg+...], lower than 4x10(-9) (with respect to H) and a cosmic-rays ionization rate of zeta=(5+/-3)x10(-17) s(-1). The inclusion of PAHs modifies the ionization fraction gradient and increases the required metal abundance. The ionization fraction in the Horsehead edge follows a steep gradient, with a scale length of ~0.05 pc (or ~25''), from [e-]~10(-4) (or n_e ~ 1-5 cm(-3)) in the PDR to a few times ~10(-9) in the core. PAH^- anions play a role in the charge balance of the cold and neutral gas if substantial amounts of free PAHs are present ([PAH] >10(-8)).

preprint2007arXiv

A PDR-Code Comparison Study

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.