VSI-stability maps showing the ratio of the relaxation time to the critical cooling time defined in the local stability criterion (Eq. (7)). From left to right, t_rel/t_crit values are computed assuming only radiative cooling (t_rel = t_cool), including dust collisions (), and including as well molecular emission with (Eq. (13)). The top and bottom row correspond to the nominal and dust-depleted cases, respectively. Also shown in the leftmost panels is the location of the domain used in our Rad-HD simulations (white dashed box). The gas opacities reported in Freedman et al. (2008), Malygin et al. (2014, 2017) for solar abundances correspond to .