Lowering energy in molecular orbital mixing
The orbital $\psi$ is obtained as a mixture of the orbitals {$\phi_i$}

$$\psi= c_1 \phi_1 + c_2 \phi_2$$

choosing the {$c_i$} through energy minimization.

Once obtained $\psi$, if you try to mix it with the initial {$\phi_i$} to obtain a new function, lets say $\mu$, you'll have:

$$\mu = c_3 \psi + c_4 \phi_1 + c_5 \phi_2$$ 

(notice that the coefficients don't need to be the same)

or equally

$$\mu = c_3 (c_1 \phi_1 + c_2 \phi_2) + c_4 \phi_1 + c_5 \phi_2$$ 
$$\mu = c_3 c_1 \phi_1 + c_3c_2 \phi_2 + c_4 \phi_1 + c_5 \phi_2$$ 

$$\mu = (c_3 c_1 + c_4) \phi_1 + (c_3c_2+c_5) \phi_2$$ 

that has the same functional form that 
$$\psi= c_1 \phi_1 + c_2 \phi_2$$

so as the criteria is the energy minimization it should be no surprising that we can not get better energy that making
$(c_3 c_1 + c_4) = c_1$ and $(c_3c_2+c_5)=c_2$.

The idea behind this is that you want to obtain the "better" function $\psi$, but you are restricted to a set of functions of the form: $c_1 \phi_1 + c_2 \phi_2$.  You impose this restriction because it only requires to obtain two coefficients, that is far simpler than obtaining an arbitrary function (defined with infinite images). 
One restriction is that the equation $c_i \phi_i = c_k \phi_k$  can only be true if $c_i = c_k = 0$.

You can expand your search to a larger set of functions, for example: $c_1 \phi_1 + c_2 \phi_2 + c_3 \phi_3$, but it requires that the set still be linearly independent, that is, you can not get any of the functions with a linear combination of the others, (i.e. it is not possible to get $c_1 \phi_1  =  c_2 \phi_2 + c_3 \phi_3$ if at least one $c_i$ is different from 0). If not, as in the discussion above, the space of function will contain the same functions so the inclusion of the new function is in some sense redundant or useful.

You can add infinite linearly independent functions, but, it does not mean that the energy will be always get lower.