From forced collapse to H ii region expansion in Mon R2: Envelope density structure and age determination withHerschel

Abstract

Context. The surroundings of H ii regions can have a profound influence on their development, morphology, and evolution. This paper explores the effect of the environment on H ii regions in the MonR2 molecular cloud.

Aims. We aim to investigate the density structure of envelopes surrounding H ii regions and to determine their collapse and ionisation expansion ages. The Mon R2 molecular cloud is an ideal target since it hosts an H ii region association, which has been imaged by the Herschel PACS and SPIRE cameras as part of the HOBYS key programme.

Methods. Column density and temperature images derived from Herschel data were used together to model the structure of H ii bubbles and their surrounding envelopes. The resulting observational constraints were used to follow the development of the Mon R2 ionised regions with analytical calculations and numerical simulations.

Results. The four hot bubbles associated with H ii regions are surrounded by dense, cold, and neutral gas envelopes, which are partly embedded in filaments. The envelope’s radial density profiles are reminiscent of those of low-mass protostellar envelopes. The inner parts of envelopes of all four H ii regions could be free-falling because they display shallow density profiles: ρ(r) ∝ r− q with . As for their outer parts, the two compact H ii regions show a ρ(r) ∝ r-2 profile, which is typical of the equilibrium structure of a singular isothermal sphere. In contrast, the central UCH ii region shows a steeper outer profile, ρ(r) ∝ r-2.5, that could be interpreted as material being forced to collapse, where an external agent overwhelms the internal pressure support.

Conclusions. The size of the heated bubbles, the spectral type of the irradiating stars, and the mean initial neutral gas density are used to estimate the ionisation expansion time, texp ~ 0.1 Myr, for the dense UCH ii and compact H ii regions and ~ 0.35 Myr for the extended H ii region. Numerical simulations with and without gravity show that the so-called lifetime problem of H ii regions is an artefact of theories that do not take their surrounding neutral envelopes with slowly decreasing density profiles into account. The envelope transition radii between the shallow and steeper density profiles are used to estimate the time elapsed since the formation of the first protostellar embryo, tinf ~ 1 Myr, for the ultra-compact, 1.5−3 Myr for the compact, and greater than ~6 Myr for the extended H ii regions. These results suggest that the time needed to form a OB-star embryo and to start ionising the cloud, plus the quenching time due to the large gravitational potential amplified by further in-falling material, dominates the ionisation expansion time by a large factor. Accurate determination of the quenching time of H ii regions would require additional small-scale observationnal constraints and numerical simulations including 3D geometry effects.