The KOSMA-$\tau$ PDR Model -- I. Recent updates to the numerical model of photo-dissociated regions

Abstract

Numerical models of Photodissociation Regions (PDRs) are an essential tool to quantitatively understand observations of massive star forming regions through simulations. Few mature PDR models are available and the Cologne KOSMA-$\tau$ PDR model is the only sophisticated model that uses a spherical cloud geometry thereby allowing us to simulate clumpy PDRs. We present the current status of the code as reference for modelers and for observers that plan to apply KOSMA-$\tau$ to interpret their data. For the numerical solution of the chemical problem we present a superior Newton-Raphson stepping algorithm and discuss strategies to numerically stabilize the problem and speed up the iterations. The chemistry in KOSMA-$\tau$ is upgraded to include the full surface chemistry in an up-to-date formulation and we discuss a novel computation of branching ratios in chemical desorption reactions. The high dust temperature in PDRs leads to a selective freeze-out of oxygen-bearing ice species due to their higher condensation temperatures and we study changes in the ice mantle structures depending on the PDR parameters, in particular the impinging UV field. Selective freeze-out can produce enhanced C abundances and higher gas temperatures resulting in a fine-structure line emission of atomic carbon [C] enhanced by up to 50% if surface reactions are considered. We show how recent ALMA observations of HCO${}^{+}$ emission in the Orion Bar with high spatial resolution on the scale of individual clumps can be interpreted in the context of non-stationary, clumpy PDR ensembles. Additionally, we introduce WL-PDR, a simple plane-parallel PDR model written in Mathematica to act as numerical testing environment of PDR modeling aspects.

39 pages, accepted by A&A

References in corpus (79)

• The James Webb Space Telescope
• Complex Chemistry in Star-Forming Regions: An Expanded Gas-Grain Warm-up Chemical Model
• Formation of methyl formate and other organic species in the warm-up phase of hot molecular cores
• The Atomic to Molecular Transition in Galaxies. II: HI and H_2 Column Densities
• The 2014 KIDA network for interstellar chemistry
• Early Science with SOFIA, the Stratospheric Observatory for Infrared Astronomy
• Simulation of the Formation and Morphology of Ice Mantles on Interstellar Grains
• GREAT: the SOFIA high-frequency heterodyne instrument
• Binding energies: new values and impact on the efficiency of chemical desorption
• HI-to-H2 Transitions and H I Column Densities in Galaxy Star-Forming Regions
• The Atomic to Molecular Transition in Galaxies. I: An Analytic Approximation for Photodissociation Fronts in Finite Clouds
• [CII] 158$\mu$m Emission and Metallicity in PDRs
• Photodesorption of water ice: a molecular dynamics study
• Modeling Complex Organic Molecules in dense regions: Eley-Rideal and complex induced reaction
• Tracing the total molecular gas in galaxies: [CII] and the CO-dark gas
• Line Emission from Gas in Optically Thick Dust Disks around Young Stars
• Compression and ablation of the photo-irradiated cloud the Orion Bar
• Desorption From Interstellar Ices
• The Perils of Clumpfind: The Mass Spectrum of Sub-structures in Molecular Clouds
• Sensitivity analysis of grain surface chemistry to binding energies of ice species
• Properties of the dense core population in Orion B as seen by the Herschel Gould Belt survey
• Gas-grain chemistry in cold interstellar cloud cores with a microscopic Monte Carlo approach to surface chemistry
• The Herschel Dwarf Galaxy Survey: II. Physical conditions, origin of [CII] emission, and porosity of the multiphase low-metallicity ISM
• 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 upGREAT dual frequency heterodyne arrays for SOFIA
• Surface chemistry in the Interstellar Medium II. ${\mathrm{H}}_2$ formation on dust with random temperature fluctuations
• The penetration of FUV radiation into molecular clouds
• Consistent dust and gas models for protoplanetary disks: II. Chemical networks and rates
• Chemical Modelling of Young Stellar Objects, I. Method and Benchmarks
• From diffuse gas to dense molecular cloud cores
• Efficiency of non-thermal desorptions in cold-core conditions. Testing the sputtering of grain mantles induced by cosmic rays
• Multidimensional chemical modelling, II. Irradiated outflow walls
• A clumpy-cloud PDR model of the global far-infrared line emission of the Milky Way
• Dust dynamics and evolution in expanding HII regions. I. Radiative drift of neutral and charged grains
• H2 formation on PAHs in photodissociation regions: a high-temperature pathway to molecular hydrogen
• Photodissociation Region Diagnostics Across Galactic Environments
• The Astrochemical Impact of Cosmic Rays in Protoclusters I: Molecular Cloud Chemistry
• Full SED fitting with the KOSMA-τ PDR code - I. Dust modelling
• Prestellar Core Formation, Evolution, and Accretion from Gravitational Fragmentation in Turbulent Converging Flows
• The H+ Region Contribution to [C II] 158 Micron Emission
• Photoevaporation of Jeans-unstable molecular clumps
• A three-phase approach to grain surface chemistry in protoplanetary disks: Gas, ice surfaces and ice mantles of dust grains
• Bottlenecks to interstellar sulfur chemistry: Sulfur-bearing hydrides in UV-illuminated gas and grains
• A New and Simple Approach to Determine the Abundance of Hydrogen Molecules on Interstellar Ice Mantles
• Ion-molecule reactions involving HCO${}^{+}$ and N${}_2$H${}^{+}$: Isotopologue equilibria from new theoretical calculations and consequences for interstellar isotope fractionation
• Herschel / HIFI spectral line survey of the Orion Bar - Temperature and density differentiation near the PDR surface
• Submillimeter Line Emission from LMC N159W: a Dense, Clumpy PDR in a Low Metallicity Environment
• Self-absorption in [CII], ${}^{12}$CO, and HI in RCW120. Building up a geometrical and physical model of the region
• Molecular globules in the Veil bubble of Orion. IRAM 30m 12CO, 13CO, and C18O 2-1 expanded maps of Orion A
• The interstellar medium of the Antennae Galaxies
• Modelling the ArH${}^{+}$ emission from the Crab Nebula
• Assessing the Performance of Molecular Gas Clump Identification Algorithms
• SOFIA/FIFI-LS Full-disk [CII] Mapping and CO-dark Molecular Gas across the Nearby Spiral Galaxy NGC 6946
• The PDR structure and kinematics around the compact HII regions S235A and S235C with [CII], [13CII], [OI] and HCO+ line profiles
• Characterizing the Multi-Phase Origin of [CII] Emission in M101 and NGC 6946 with Velocity Resolved Spectroscopy
• A Revised Description of the Cosmic Ray-Induced Desorption of Interstellar Ices
• Study of Photon Dominated Regions in Cepheus B
• Molecular clumps photoevaporation in ionized regions
• The Interstellar Medium and Star Formation of Galactic Disks. I. ISM and GMC properties with Diffuse FUV and Cosmic Ray Backgrounds
• A Search for O_2 in CO-depleted Molecular Cloud Cores with Herschel
• Opening the Treasure Chest in Carina
• High Spectral and Spatial Resolution Observations of the PDR Emission in the NGC2023 Reflection Nebula with SOFIA and APEX
• HI-to-H${}_2$ Transitions in Dust-Free Interstellar Gas
• Dynamical effects of the radiative stellar feedback on the H I-to-H2 transition
• Effect of binding energies on the encounter desorption
• The Photon Dominated Region in the IC 348 molecular cloud
• First detection of [13CII] in the Large Magellanic Cloud
• Globules and pillars in Cygnus X III. Herschel and upGREAT/SOFIA far-infrared spectroscopy of the globule IRAS 20319+3958 inCygnus X
• Photoevaporating PDR models with the Hydra PDR Code
• Observations of the [CI] (${}^3{P}_1$-${}^3{P}_0$) emission toward the massive star-forming region RCW38: further evidence for highly-clumped density distribution of the molecular gas
• The CO universe: Modelling CO emission and H${}_2$ abundance in cosmological galaxy formation simulations
• Molecular dark matter in galaxies
• The distribution of ionized, atomic and PDR gas around S1 in Rho Ophiuchus
• Sustained oscillations in Interstellar chemistry models
• Astrochemical Bistability: Autocatalysis in Oxygen Chemistry
• Warm ISM in the Sgr A Complex. II. The [C/N] abundance ratio traced by [CII] 158 um and [NII] 205 um observations toward the Arched Filaments at the Galactic center
• Determining the H+ Region/PDR Equation of State in Star Forming Regions
• Species cycling and the enhancement of ammonia in prestellar cores
• Rate Constants for Fine-Structure Excitations in O-H Collisions with Error Bars Obtained by Machine Learning

Cited by in corpus (9)

• The PhotoDissociation Region Toolbox: Software and Models for Astrophysical Analysis
• PDFchem: A new fast method to determine ISM properties and infer environmental parameters using probability distributions
• High angular resolution near-IR view of the Orion Bar revealed by Keck/NIRC2
• Neural network-based emulation of interstellar medium models
• Impact of HAC evolution on the formation of small hydrocarbons in the Orion Bar and the Horsehead PDRs
• Ionized Carbon in Galaxies: The [C II] 158 $\mu$m Line as a Total Molecular Gas Mass Tracer Revisited
• Beetroots: spatially-regularized Bayesian inference of physical parameter maps -- Application to Orion
• Extended atomic carbon around molecular clouds
• First Submillimeter Lights from Dome A: Tracing the Carbon Cycle in the Feedback of Massive Stars