Beyond the Hohenberg−Kohn density functional theory for the ground state, it has been established that the Hamiltonian matrix for a finite number (N) of lowest eigenstates is a matrix density functional. Its fundamental variable�the matrix density D(r)�can be represented by, or mapped to, a set of auxiliary, multiconfigurational wave functions expressed as a linear combination of no more than N 2 determinant configurations. The latter defines a minimal active space (MAS), which naturally leads to the introduction of the correlation matrix functional, responsible for the electronic correlation effects outside the MAS. In this study, we report a set of rigorous conditions in the Hamiltonian matrix functional, derived by enforcing the symmetry of a Hilbert subspace, namely the subspace invariance property. We further establish a fundamental theorem on the correlation matrix functional. That is, given the correlation functional for a single state in the N-dimensional subspace, all elements of the correlation matrix functional for the entire subspace are uniquely determined. These findings reveal the intricate structure of electronic correlation within the Hilbert subspace of lowest eigenstates and suggest a promising direction for efficient simulation of excited states.