The Cartan decomposition of a semisimple Lie algebra is pivotal in both theoretical exploration and practical application. It generalizes the classical spectral, polar, and singular value decompositions found in linear algebra. This paper explores an advanced refinement of the Cartan decomposition of \(\mathfrak {su}(2^{n})\) , arising from the parameterization of a Hamiltonian for quantum simulation. Specifically, given an AI-type Cartan decomposition \({\mathfrak {g}} = {\mathfrak {k}} \oplus {\mathfrak {p}}\) of any Lie subalgebra \({\mathfrak {g}}\) in \(\mathfrak {su}(2^{n})\) , where \({\mathfrak {p}} = \widetilde{{\mathfrak {p}}} \oplus {\mathfrak {h}}\) and \({\mathfrak {h}}\) is the maximal abelian subalgebra in \({\mathfrak {p}}\) , it is revealed that both the subalgebra \({\mathfrak {k}}\) and the corresponding \(\widetilde{{\mathfrak {p}}}\) (which is generally not a subalgebra) can further be partitioned, respectively, as the direct sums of equal-dimensional commutative Lie subalgebras. With respect to the involution \(\theta (g)=-g^{\top }\) , any Lie subalgebra \({\mathfrak {g}}\) in \(\mathfrak {su}(2^{n})\) with nondegenerate Cartan decomposition can thus be fully decomposed as the orthogonal direct sum of \({\mathfrak {h}}\) and pairs of commutative Lie subalgebras over \({\mathfrak {k}}\) and \(\widetilde{{\mathfrak {p}}}\) , all of which, except \({\mathfrak {h}}\) , share the same dimension.