KOSMA-τ Photodissociation Model

• KOSMA‑τ Photodissociation Model is an advanced numerical tool that uses spherical cloud geometry to simulate fragmented, UV-irradiated gas clouds.
• It integrates dust radiative transfer, gas-phase and grain-surface chemistry, and thermal equilibrium to predict line and continuum observables.
• The model supports both single-clump and ensemble formulations, enabling detailed analysis of star-forming regions and cloud fragmentation.

Searching arXiv for recent and foundational KOSMA-$\tau$ papers to ground the article. KOSMA-$\tau$ is a mature numerical photodissociation region (PDR) model for computing the coupled chemical, thermal, and radiative structure of UV-irradiated neutral gas clouds. Its distinguishing feature within the mature PDR-model landscape is its spherical cloud or clump geometry, which makes it possible to treat finite-mass clumps and, by superposition, ensembles of clumps as an approximation to fragmented, porous, and fractal molecular material. In its modern form, KOSMA-$\tau$ combines dust radiative transfer, stationary gas-phase and grain-surface chemistry, thermal balance, and post-processed line transfer to predict line and continuum observables for massive star-forming environments (Röllig et al., 2022).

1. Position within PDR modeling

KOSMA-$\tau$ is designed to simulate regions in which far-ultraviolet photons control the chemistry and temperature of neutral gas. In that regime, the physics is strongly coupled: FUV photons photodissociate and photoionize species, the resulting abundances regulate heating and cooling, and the temperature in turn modifies chemistry and excitation. The code addresses that stationary problem in one spatial dimension, but does so in spherical rather than plane-parallel geometry. The 2022 reference description states that few mature PDR models are available and identifies KOSMA-$\tau$ as the only sophisticated model built around spherical cloud geometry, a choice motivated by the observational fact that molecular clouds are fragmented and clumpy (Röllig et al., 2022).

Historically, the code evolved from the Sternberg plane-parallel PDR code into a spherical Cologne/Tel Aviv model. That evolution is methodologically important because the spherical, finite-mass setup changes both the surface-to-volume ratio and the angular structure of attenuation. In practice, KOSMA-$\tau$ is used to compute depth-dependent gas and dust temperatures, abundances of gas-phase and surface species, local photorates and shielding, level populations, emergent line intensities, line spectral energy distributions, and, in ensemble mode, cloud-scale observables obtained by summing over a clump mass spectrum (Röllig et al., 2022).

2. Geometry, illumination, and radiative-transfer architecture

In the single-clump formulation, KOSMA-$\tau$ is a 1D spherical model with depth coordinate

$z = R_\mathrm{tot} - r,$

where $R_\mathrm{tot}$ is the total clump radius and $r$ is the radial distance from the center. The standard density structure is centrally condensed,

$\tau$0

with standard parameters $\tau$1 and $\tau$2. Illumination is taken to be isotropic FUV irradiation. In the 2022 code description the field strength is denoted $\tau$3 in units of the Draine field,

$\tau$4

integrated from 91.2 to 200 nm, with $\tau$5 in Habing units (Röllig et al., 2022).

Radiative transfer is split into three layers. First, a dust-continuum pre-processing stage uses MCDRT to compute the internal spectrally resolved FUV field, the dust temperature distribution, and emitted dust continuum. Second, during the PDR iterations, the model treats line-related UV transfer and local cooling transfer, including H$\tau$6 shielding and pumping, C photoionization and dissociation shielding as implemented, CO and isotopologue shielding, and line cooling via spherical escape probability. Third, final emergent line emission is obtained with the ONION post-processing code. The angular structure of attenuation enters explicitly through the directional average

$\tau$7

so local photorates depend on angle-dependent columns rather than a single slab sightline (Röllig et al., 2022).

KOSMA-$\tau$8 is also used in an ensemble formulation. In the S140 application, the PDR is represented by an ensemble of spherical clumps with a power-law clump-mass spectrum,

$\tau$9

and a mass-size relation

$\tau$0

The fixed structural parameters in that formulation are $\tau$1, $\tau$2, $\tau$3, and $\tau$4, with a clump density profile

$\tau$5

This ensemble formalism underlies the model’s use for clumpy, multi-component PDRs rather than uniform slabs (Dedes et al., 2010).

3. Chemistry, surface processes, and thermal balance

The physical foundation of KOSMA-$\tau$6 is the general PDR result that photodissociation rates depend on both molecular physics and radiative transfer, including dust extinction, scattering, self-shielding, and mutual shielding. In that sense, a PDR code cannot rely only on a scalar UV scaling factor unless the photorates have already been calibrated for a specific field shape and attenuation law (Dishoeck et al., 2011).

The modern gas-phase chemistry in KOSMA-$\tau$7 is modular and based on the UMIST 2012 network, with updates including isotopologues such as $\tau$8C and $\tau$9O, revised branching ratios, $\tau$0-type isomers, updated fractionation reactions, and low-temperature refits. For each gas-phase species $\tau$1, the steady-state chemistry is written as

$\tau$2

Thermal equilibrium is obtained from

$\tau$3

with implemented heating and cooling terms including photoelectric heating, H$\tau$4 vibrational de-excitation heating, H$\tau$5 photodissociation heating, H$\tau$6 formation heating, cosmic-ray heating, C photoionization heating, gas-grain collisional exchange, and cooling from [O I], [C II], [C I], CO, $\tau$7CO, H$\tau$8O, OH, Ly$\tau$9, and O I 6300 Å (Röllig et al., 2022).

A major 2022 update is the inclusion of full grain-surface chemistry in a quasi-three-phase model consisting of gas, chemically active surface, and inert bulk ice mantle. Only the top surface layer is mobile, and desorption is limited to the top two monolayers. The surface rate equation is

$\tau$0

The model uses a surface site density $\tau$1, a tunneling barrier width $\tau$2, Cazaux–Tielens H$\tau$3 formation, encounter desorption for $\tau$4 and $\tau$5, and updated cosmic-ray and chemical desorption prescriptions (Röllig et al., 2022).

One of the principal physical consequences of this upgrade is selective freeze-out. Because different ice species have different condensation temperatures—$\tau$6 at $\tau$7 K, $\tau$8 at $\tau$9 K, $\tau$0 at $\tau$1 K, $\tau$2 at $\tau$3 K, $\tau$4 at $\tau$5 K, and $\tau$6 at $\tau$7 K—warm PDR dust can preferentially lock oxygen-bearing species into ice while leaving much carbon in the gas. The resulting reduction of CO cooling and suppression of gas-phase destruction routes for atomic carbon can enhance [C I] fine-structure emission by up to $\tau$8 when surface reactions are included (Röllig et al., 2022).

4. Dust-consistent extension and full SED modeling

A major extension of KOSMA-$\tau$9 revised the treatment of interstellar dust so that dust-related physics is described consistently throughout the model. In that formulation, the code is coupled to MCDRT to solve frequency-dependent radiative transfer and the dust thermal-balance equation in a dusty clump under spherical symmetry. The treatment assumes thermal equilibrium dust temperatures, includes isotropic scattering, computes both line and dust-continuum emission, and deliberately neglects non-equilibrium or stochastic heating of very small grains and PAHs (Röllig et al., 2012).

Dust is represented as a mixture of grain sorts with explicit size distributions. The paper considers MRN dust and three Weingartner & Draine (2001) models: WD01-7 ($\tau$0), WD01-21 ($\tau$1), and WD01-25 ($\tau$2). Optical properties are computed with Mie theory from $\tau$3 to $\tau$4 on 333 wavelengths. The radiative-transfer output provides the local mean intensity

$\tau$5

and the local dust temperatures $\tau$6. These feed directly into photo-reaction rates,

$\tau$7

and into the thermal and chemical coupling between gas and grains (Röllig et al., 2012).

The dust extension also revised photoelectric heating and H$\tau$8 formation. The control parameter for photoelectric heating is

$\tau$9

and the updated fitting range extends to $z = R_\mathrm{tot} - r,$0, beyond the original WD01 domain. H$z = R_\mathrm{tot} - r,$1 formation on grains was reworked to include physisorption, chemisorption, and what the paper describes as the Eley-Rideal effect, allowing efficient high-temperature formation. The total formation rate is

$z = R_\mathrm{tot} - r,$2

and the authors explicitly set $z = R_\mathrm{tot} - r,$3 in the adopted efficiency expression in order to avoid unrealistically trapping newly formed H$z = R_\mathrm{tot} - r,$4 on cold grains (Röllig et al., 2012).

The paper identifies the H$z = R_\mathrm{tot} - r,$5 formation revision as the most influential modification. H$z = R_\mathrm{tot} - r,$6 formation heating, with roughly $z = R_\mathrm{tot} - r,$7 eV per formed H$z = R_\mathrm{tot} - r,$8 assumed to heat the gas, can dominate the thermal budget of outer cloud layers, drive temperatures above $z = R_\mathrm{tot} - r,$9 K in some models, and strongly boost high-$R_\mathrm{tot}$0 CO emission. PAH surfaces further enhance both photoelectric heating and H$R_\mathrm{tot}$1 formation. Increasing the abundance of small grains produces hotter outer layers because heating becomes more efficient, but cooler cloud centers because FUV extinction becomes more efficient. High-$R_\mathrm{tot}$2 CO transitions are correspondingly emphasized as the clearest line diagnostic of these dust-driven changes (Röllig et al., 2012).

5. Observational configurations and representative applications

KOSMA-$R_\mathrm{tot}$3 has been applied in at least two distinct interpretive modes: as a compact single-clump diagnostic based on a small set of line ratios, and as a clumpy ensemble model fit to a broad set of absolute line intensities.

Application	KOSMA-$R_\mathrm{tot}$4 configuration	Main inference
IC1396A	Single externally illuminated clump; ratios of $R_\mathrm{tot}$5 and $R_\mathrm{tot}$6	Densities of $R_\mathrm{tot}$7–$R_\mathrm{tot}$8, typically $R_\mathrm{tot}$9, with local shielding at one rim position
S140	Clumpy ensemble of spherical clumps; absolute intensities; two components at IRS1	A small, hot, highly irradiated component plus a cooler, more massive shielded component; the hot component dominates $r$0

In IC1396A, the model is used deliberately as a compact interpretive tool rather than a global fit. The observed quantities are the integrated-intensity ratios of $r$1, $r$2, and CO(4–3), compared to a single-clump KOSMA-$r$3 model with free parameters mean density $r$4, FUV field $r$5, and clump mass $r$6. Because only two independent ratios are available, the analysis explores two model families: Model 1 fixes $r$7 from the luminosity of HD 206267 and fits $r$8 and $r$9; Model 2 fixes $\tau$00 and fits $\tau$01 and $\tau$02. The [C II] profile is not integrated wholesale, but is fitted with the CO(4–3) profile at each position so that only the kinematically associated [C II] component enters the ratio analysis. After convolution to a common $\tau$03 beam and adopting $\tau$04 systematic uncertainties for each line, the principal result is that the inferred density is typically in the range $\tau$05–$\tau$06 at position B and about $\tau$07 at the other positions. The two model strategies give densities consistent within the errors. At position B, however, the fitted $\tau$08 is much lower than the geometric estimate, and the authors interpret that discrepancy as local shielding of the UV field at the rim of the globule (Okada et al., 2012).

In S140, the model is used in its clump-ensemble formulation. Each ensemble has five free parameters: average density $\tau$09, ensemble mass $\tau$10, UV field strength $\tau$11, and the minimum and maximum clump masses $\tau$12. At IRS1 the data require a two-component fit. The hot component has $\tau$13, $\tau$14, $\tau$15, $\tau$16, $\tau$17, $\tau$18–$\tau$19, and $\tau$20. The cool component has $\tau$21, $\tau$22–$\tau$23, $\tau$24, $\tau$25, $\tau$26–$\tau$27, $\tau$28–$\tau$29, and $\tau$30–$\tau$31. The hot component is interpreted as possibly associated with irradiated outflow cavity walls, and about $\tau$32 of the $\tau$33 emission around IRS1 is attributed to it. At the ionization front, no full fit is presented, but the observations imply a hot dense component with $\tau$34, $\tau$35, and a beam mass of order $\tau$36, together with an additional UV source beyond HD 211880 (Dedes et al., 2010).

6. Limitations, numerical behavior, and interpretive scope

KOSMA-$\tau$37 is fundamentally a stationary 1D spherical model. The 2022 reference paper lists several explicit scope limitations: isotropic FUV illumination, dust continuum radiative transfer precomputed once rather than iterated with evolving gas-line UV absorption, isotropic scattering, surface chemistry using a single average dust temperature rather than size-resolved temperatures, explicit line shielding only for selected species such as $\tau$38, C, and CO isotopologues, final line emission computed without line overlap or pumping between different molecules, and quasi-three-phase rather than full dynamic multilayer ice chemistry. The same paper also notes that multiple chemical and thermal solutions can exist, including hot atomic and warm molecular branches of the thermal equation, and that some models exhibit oscillatory convergence. WL-PDR, a simple plane-parallel PDR model written in Mathematica, is introduced specifically as a numerical testing environment for such issues (Röllig et al., 2022).

Application papers expose additional limitations of particular KOSMA-$\tau$39 use modes. In the S140 ensemble analysis, optical-depth effects are included only within individual clumps, not between clumps or between ensembles. The authors state explicitly that this approximation is acceptable for most species but not for $\tau$40, because the line is very optically thick, with $\tau$41, so the model is not able to give reliable estimates for the $\tau$42 intensities. In IC1396A, the inverse problem is underdetermined from the outset: three model parameters are confronted with only two independent line ratios, absolute intensities are not fitted, beam filling factors are not solved for independently, the source geometry is represented by a single unresolved clump, and the [C II] component used in the ratios must first be isolated by profile matching to CO(4–3), a physically sensible but not unique procedure (Dedes et al., 2010, Okada et al., 2012).

The dust-consistent extension improves internal consistency but retains approximations of its own. It assumes thermal-equilibrium grain temperatures, neglects stochastic heating of PAHs and very small grains, and uses isotropic rather than anisotropic scattering. A plausible implication is that KOSMA-$\tau$43 is best understood not as a universal geometric description of real PDR structure, but as a physically structured and observationally productive clump framework whose strengths are greatest when line selection, beam treatment, and source morphology are commensurate with its spherical, stationary assumptions (Röllig et al., 2012).