Hot and cool water in Herbig Ae protoplanetary disks. A challenge for Herschel

The spatial origin and detectability of rotational H2O emission lines from Herbig Ae type protoplanetary disks beyond 70 micron is discussed. We use the recently developed disk code ProDiMo to calculate the thermo-chemical structure of a Herbig Ae type disk and apply the non-LTE line radiative transfer code Ratran to predict water line profiles and intensity maps. The model shows three spatially distinct regions in the disk where water concentrations are high, related to different chemical pathways to form the water: (1) a big water reservoir in the deep midplane behind the inner rim, (2) a belt of cold water around the distant icy midplane beyond the snow-line r>20AU, and (3) a layer of irradiated hot water at high altitudes z/r=0.1…0.3, extending from about 1AU to 30AU, where the kinetic gas temperature ranges from 200K to 1500K. Although region 3 contains only little amounts of water vapour (~3x10^-5 M_Earth), we find this warm layer to be almost entirely responsible for the rotational water emission lines, execpt for the 3 lowest excitation lines. Thus, Herschel will probe first and foremost the conditions and radial extension of region 3, where water is predominantly formed via neutral-neutral reactions and the gas is thermally decoupled from the dust T_gas>T_dust. The observations do not allow for a determination of the snow-line, because the snow-line truncates the radial extension of region 1, whereas the lines originate from region 3. Different line transfer approximations (LTE, escape probability, Monte Carlo) are discussed. A non-LTE treatment is required in most cases, and the results obtained with the escape probability method are found to underestimate the Monte Carlo results by 2%…45%.

💡 Research Summary

This paper investigates the origin, spatial distribution, and observability of rotational water (H₂O) emission lines in Herbig Ae protoplanetary disks at wavelengths longer than 70 µm, with a focus on what the Herschel Space Observatory can detect. The authors combine two state‑of‑the‑art modeling tools: ProDiMo, a thermo‑chemical disk code that solves for the 2‑D (radius‑height) structure of gas temperature, dust temperature, density, radiation field, and chemical abundances; and Ratran, a non‑LTE line radiative‑transfer code that uses Monte‑Carlo photon propagation to compute level populations and emergent line profiles.

ProDiMo is fed with a typical Herbig Ae stellar spectrum (effective temperature ≈ 9500 K, luminosity ≈ 30 L⊙), a disk mass of ~10⁻² M⊙, an inner radius of 0.1 AU, an outer radius of 200 AU, and a realistic dust size distribution. The chemical network includes thousands of reactions, allowing the model to self‑consistently determine where water is present in the gas phase versus frozen onto grains. Three distinct water‑rich zones emerge:

1. Deep mid‑plane reservoir (Region 1) – located just behind the puffed‑up inner rim, this region is cold and dense; water is almost entirely in ice, contributing the bulk of the total water mass but little to the gas‑phase emission.

2. Cold icy belt beyond the snow line (Region 2) – at radii larger than ~20 AU, temperatures drop to 30–100 K. Again, water is predominantly solid, with only trace gas‑phase abundances.

3. Warm, tenuous surface layer (Region 3) – extending from ~1 AU to ~30 AU at heights z/r ≈ 0.1–0.3, the gas temperature ranges from 200 K up to 1500 K. Here water is formed mainly through neutral‑neutral reactions (O + H₂ → OH, OH + H₂ → H₂O). Although the total water mass in this layer is modest (~3 × 10⁻⁵ M⊕), it dominates the rotational line emission, accounting for > 90 % of the flux for all but the three lowest‑excitation lines.

Using Ratran, the authors calculate the emergent spectra for ~30 key rotational transitions (e.g., the 179.5 µm and 538 µm lines). The Monte‑Carlo non‑LTE results show that LTE assumptions severely underestimate line strengths, especially for higher‑energy transitions where the critical densities exceed the local gas density. An escape‑probability approach improves on LTE but still falls short of the full Monte‑Carlo solution by 2 % to 45 %, depending on line opacity and excitation. Consequently, a proper non‑LTE treatment is mandatory for quantitative interpretation of Herschel data.

From the synthetic line profiles the authors infer that Herschel will primarily probe Region 3. The warm surface layer is thermally decoupled from the dust (T_gas > T_dust), so the detected lines carry information about gas heating mechanisms (e.g., UV/X‑ray irradiation) and the efficiency of neutral‑neutral water formation pathways. Because the emission originates well above the mid‑plane, Herschel cannot directly locate the snow line; the snow line truncates the radial extent of Region 1, but the observable lines are insensitive to that truncation.

The paper also discusses observational strategies. Low‑excitation lines (with upper‑level energies < 100 K) have contributions from the cold belt and can, in principle, hint at the presence of icy material, but their fluxes are weak and blended with background. High‑excitation lines are bright and spatially compact, making them ideal tracers of the warm surface layer. The authors suggest that combining Herschel spectra with high‑resolution interferometric data from ALMA or future JWST observations could break degeneracies and map the vertical distribution of water more precisely.

In summary, the study demonstrates that rotational water emission in Herbig Ae disks is dominated by a thin, warm, UV‑irradiated surface layer rather than the massive icy mid‑plane reservoir. Accurate modeling requires a full 2‑D thermo‑chemical structure and non‑LTE radiative transfer. Herschel will be sensitive to the conditions and radial extent of this warm layer, providing valuable constraints on gas heating, water chemistry, and the decoupling of gas and dust temperatures in the planet‑forming region of intermediate‑mass stars.