We study properties of asymptotically free vectorial gauge theories with gauge groups G = SO(Nc) and G = Sp(Nc) and N f fermions in a representation R of G, at an infrared (IR) zero of the beta function, αIR, in the non-Abelian Coulomb phase. The fundamental, adjoint, and rank-2 symmetric and antisymmetric tensor fermion representations are considered. We present scheme-independent calculations of the anomalous dimensions of (gauge-invariant) fermion bilinear operators γψ ψ,IR to O(∆ 4 f ) and of the derivative of the beta function at αIR, denoted β ′ IR , to O(∆ 5 f ), where ∆ f is an N f -dependent expansion variable. It is shown that all coefficients in the expansion of γψ ψ,IR that we calculate are positive for all representations considered, so that to O(∆ 4 f ), γψ ψ,IR increases monotonically with decreasing N f in the non-Abelian Coulomb phase. Using this property, we give a new estimate of the lower end of this phase for some specific realizations of these theories.