Gravitational trapping and ram pressure trapping of ultracompact and hypercompact H II regions

Observationally, early H II regions are classified by size into ultracompact and hypercompact configurations. It remains unclear whether these phases are long-lived or transient. Understanding the physical processes that stall H II region growth may help to solve the so-called lifetime problem: the observation of more compact H II regions than expected from theory. Utilizing two-dimensional, axially symmetric radiation hydrodynamic simulations of young expanding H II regions, including the phase of early star and disk formation, we seek to better understand the trapping of H II regions. Trapping forces include gravity and ram pressure, which oppose forces such as thermal pressure expansion, radiation pressure, and centrifugal force. Without radiation pressure, the H II region remains gravitationally trapped in the ultracompact phase indefinitely. With radiation pressure, the H II region escapes gravitational trapping but experiences ram pressure trapping on larger scales. For initial mass reservoirs with high central density, no trapping occurs, while a less steep density gradient yields clear trapped phases. Hypercompact trapped phases exhibit a so-called flickering variation in H II region radius, in agreement with observations of stalling and even contraction over small time scales. With radiation pressure, low-density reservoirs experience both gravitational and ram pressure trapping, while high-mass reservoirs undergo only the latter.

💡 Research Summary

The paper tackles the long‑standing “lifetime problem” of compact H II regions—why ultracompact (UC) and hypercompact (HC) H II regions are observed far more frequently than simple expansion models predict. Using two‑dimensional, axisymmetric radiation‑hydrodynamic simulations that follow a massive molecular core from gravitational collapse through protostar and disk formation to the early expansion of an ionized region, the authors include all the major forces: thermal pressure of the ionized gas, stellar gravity, self‑gravity of the gas, centrifugal forces, radiation pressure, and the ram pressure of infalling material. The numerical framework couples the PLUTO MHD code with the Makemake radiation‑transfer module, the Sedna photo‑ionization solver, and the Haumea self‑gravity solver.

A suite of models explores different core masses (≈30–120 M⊙) and initial density gradients (ρ ∝ r⁻ᵖ with p = 1.5 or 2.0). In runs without radiation pressure, the ionized bubble remains trapped by stellar gravity at radii ≲0.01 pc (the HC regime), confirming the analytical binding radius concept of Keto (2003). When radiation pressure is turned on, the gravitational trap is broken and the H II region begins to expand, but the inflow of gas from larger scales exerts a ram pressure that can halt further growth. This “ram‑pressure trapping” is strongest for shallow density profiles (p ≈ 1.5) and leads to a prolonged UC phase (0.05–0.1 pc) lasting tens of thousands of years. Steeper density profiles (p ≈ 2.0) produce rapid infall that overwhelms the ram pressure, allowing the region to expand more freely.

A particularly striking result is the reproduction of “flickering” in the HC phase: the radius oscillates on timescales of years to decades, driven by episodic accretion bursts and the formation of dense clumps in the surrounding flow. This behavior matches multi‑epoch radio observations that have recorded both brightening and dimming of UC/H II sources. The simulations also show that high‑mass reservoirs experience only ram‑pressure trapping, whereas low‑mass, low‑density clouds can exhibit both gravitational and ram‑pressure trapping.

Overall, the study demonstrates that (1) gravitational trapping dominates only when radiation pressure is neglected; (2) realistic radiation pressure eliminates the gravitational trap but introduces a ram‑pressure barrier that can sustain UC‑size H II regions for long periods; and (3) the combination of these mechanisms naturally explains the observed abundance and variability of compact H II regions, offering a quantitative resolution to the lifetime problem. Future work extending the models to three dimensions and incorporating stellar winds or proper motion will further refine the connection between theory and observation.