We consider the critical behavior at the threshold of black hole formation for the five dimensional vacuum Einstein equations satisfying the cohomogeneity-two triaxial Bianchi IX ansatz. Exploiting a discrete symmetry present in this model we predict the existence of a codimension-two attractor. This prediction is confirmed numerically and the codimension-two attractor is identified as a discretely self-similar solution with two unstable modes.Introduction. Since the pioneering work of Choptuik on the collapse of self-gravitating scalar field [1], the nature of the boundary between dispersion and black hole formation in gravitational collapse has been a very active research area (see [2] for a review). One of the most intriguing aspects of these studies is the occurrence of discretely self-similar critical solutions. Discrete self-similarity (DSS) means invariance under rescalings of space and time variables by a constant factor e ∆ where ∆ is a number, usually called the echoing period. Critical solutions possessing this curious symmetry (with different echoing periods) have been found for several selfgravitating matter models (massless scalar field [1], YangMills field [3], σ-model [4] and few others) and recently also in the vacuum gravitational collapse in higher dimensions [5,6].Our current understanding of DSS solutions is very limited in comparison with continuously self-similar (CSS) solutions. In the case of spherical symmetry the CSS ansatz leads to an ODE system which can be handled analytically and sometimes even rigorous proofs of existence are feasible. For example, the existence of a countable family of CSS solutions was proved in the Einstein-sigma model [7] and the ground state of this family was identified as the critical solution (which was known previously from numerical studies of critical collapse). In contrast, the DSS ansatz leads to a 1 + 1 PDE eigenvalue problem which seems intractable analytically. Although this eigenvalue problem can be solved numerically, as was done by Gundlach for two models (scalar field [8] and ), the numerical iterative procedure requires a good initial seed in order to converge. Thus, Gundlach's method is efficient in validating and refining DSS solutions which are already known from direct numerical simulations but it is not useful in searching for new solutions.In this paper we consider the critical collapse for the five dimensional vacuum Einstein equations satisfying the cohomogeneity-two triaxial Bianchi IX ansatz and provide heuristic arguments and numerical evidence for the the existence of a DSS solution with two unstable modes. This is a continuation of our studies in [5] where we have shown the existence of a critical DSS solution with one unstable mode and the associated type II critical behavior in this model. On the basis of our result it is tempting to conjecture that the critical DSS solution is a ground state of a countable family of DSS solutions with increasing number of unstable modes.Background. Our starting point is the cohomogeneitytwo sym...