Asking for theorems might be overkill, in the absence of a firm grasp of the construction, and they would hinder such a grasp, rather than enhance it. Nevertheless, 30 years ago, technicolorists beat such theorems to death. In any case, chiral symmetry breaking dynamics particularities are not well-understood to the point of excluding potential condensation outcomes.
What I'll do here is illustrate the point and answer your questions for the simplest model that I assume you learned on week one of your course, based on just SU(2) and its generators, the half of the prosaic three Pauli matrices.
So consider your chiral model bifundamental,
$$
\Phi =\phi_0 {\mathbb 1}+ \vec{\phi} \cdot \vec{\sigma},
$$
whose columns transform as $SU(2)_L$ doublets on the left, and its rows like $SU(2)_R$ doublets on the right.
A v.e.v. (condensate) matrix $\langle \Phi\rangle=M$ may or may not be left invariant under parts of the two SU(2) s. So your problem amounts to which v.e.v. s are left intact by which subgroups of the full six generator $SU(2)_L\times SU(2)_R$ group,
$$
e^{i\vec{\theta}_L \cdot \vec{\sigma}/2} M e^{-i\vec{\theta}_R \cdot \vec{\sigma}/2} =M ~, ~~\Longrightarrow \\
\vec{\theta}_L \cdot \vec{\sigma} ~M - M ~\vec{\theta}_R \cdot \vec{\sigma}=0 ~.
$$
You already appreciated that the 3 axials are gone and broken, and at most the three diagonal generators ($SU(2)_V$ so, $\theta_L=\theta_R$) may survive by virtue of the condensate
$$
M= v {\mathbb 1},
$$
since the identity commutes with all three diagonal (vector) generators,
$$
[\vec{\theta}\cdot \vec{\sigma} ,{\mathbb 1}]=0 ~.
$$
What happens if we complicate M and tack onto the identity a further piece
$$
M= \vec{v}\cdot \vec{\sigma} ~? ~~~~\Longrightarrow\\
[\vec{\theta}\cdot \vec{\sigma} , \vec{v}\cdot \vec{\sigma}]=2i(\vec{\theta}\times \vec{v} )\cdot \vec{\sigma}~.
$$
For the two θ directions (generators) perpendicular to $\vec v$, the cross product survives, so the r.h.s. cannot vanish, so these two generators are gone and broken. But for the direction parallel to $\vec v$, the corresponding generator survives, as the cross product collapses and the r.h.side vanishes to ensure invariance. So, indeed, the residual diagonal group is a vector U(1).
You cannot kill this U(1) with a condensate, as you know multiplication of the v.e.v. by a complex number and its conjugate will always preserve it. (So, forget about the singlet outcome.)
You should be able to trivially generalize this to SU(3) which has several nifty SU(2) subgroups, as well as U(1) s, etc. It's really straightforward.
(Hint: what does a term proportional to $\lambda_8$ do? This actually obtains in our real world: condensation, trailing vector strangeness explicit breaking; cf. §4.1.2 of Scherer 2002. Of course, in real life, isospin $\lambda_3$ is broken as well, and all three condensates differ.)
NB in response to question's edit . No, the full matrix you wrote down does not break you down to a singlet... it can be written as an identity plus a linear combination of Gell-Mann matrices, so, then a generator in the Lie algebra of SU(3) : it then defines an unbroken direction commuting with it. Construction of potentials dialing-in specific triplets (~ diagonal 3×3 M s) is standard and routine.