Line cooling (Λ_line; dashed lines) and photoelectric heating (Γ_peh; solid lines) rates as a function of radius at the time when the central density is n(H₂) = 10⁵ cm⁻³ in simulations M_TD, M_10, and M_0 (top panel). Radial distributions of selected cooling molecules at the same time (bottom panel). In this panel, solid, dashed, and dotted lines represent simulations M_TD, M_10, and M_0, respectively, but they overlap almost perfectly.