In B-O approximation, one of the basic assumptions is that the total many-body wavefunction can be expanded as the following: $$\Psi(\bf{r},\bf{R})=\sum_n\phi_n(\bf{R})\psi_n(\bf{r;\bf{R}})$$ where $\phi_n(\bf{R})$ is the wavefunction for the ions alone and the index $n$ can be thought of as the index for energy levels for the ions. $\psi(\bf{r};\bf{R})$ is the electron wavefunction for fixed ion position $\bf{R}$. The electron wavefunction $\psi$ depends parametrically on $\bf{R}$ and we can assume that electrons adjust themselves very quickly to the ion movement due to the huge difference in their velocities.
Now my question is why we can assign the same index $n$ to the electron part of the wavefunction? My reasoning is the following: For each fixed $\bf{R}$, we have a fixed background potential for the electrons, then we can in principle solve for electronic energy eigenstates, indexed by $m$. This $m$, under the aforementioned conditions, should be independent from the level index $n$ of the ions. This seems to invalidate the above expansion.