$\textbf{General Question}:$ We may write the nonrelativistic Hamiltonian for a molecule as a sum of five terms:
$$H = T_N(\textbf{R}) + T_e(\textbf{r}) + V_{eN}(\textbf{r},\textbf{R}) + V_{NN}(\textbf{R})+V_{ee}(\textbf{r}),$$
where $\textbf{R}$ is the set of nuclear coordinates and $\textbf{r}$ is the set of electronic coordinates.(Here we have ignored the spin-orbit effects.)
$\textbf{Problem}:$ Unfortunately, the $V_{eN} (\textbf{r};\textbf{R})$ term prevents us from separating H into nuclear and electronic parts, which would allow us to write the molecular wavefunction as a product of nuclear and electronic
terms, $\Psi(\textbf{r};\textbf{R} ) = \Psi(\textbf{r})\Psi(\textbf{R})$.
$\textbf{Solution}:$ We thus introduce the Born-Oppenheimer approximation, by which we conclude that this nuclear and electronic separation is approximately correct.The term $V_{eN}(\textbf{r};\textbf{R})$ is large and cannot be neglected; however, we can make the $\textbf{R}$ dependence parametric, so that the total wavefunction is given as $\Psi(\textbf{r};\textbf{R})\Psi(\textbf{R})$. The Born-Oppenheimer approximation
rests on the fact that the nuclei are much more massive than the electrons, which allows us to say
that the nuclei are nearly fixed with respect to electron motion. We can fix $\textbf{R}$, the nuclear configuration, at some value $\textbf{R}_a$, and solve for the electronic wavefunction $\Psi(\textbf{r};\textbf{R}_a)$, which depends
only parametrically on $\textbf{R}$. If we do this for a range of R, we obtain the potential energy curve along which the nuclei move.
$\textbf{Mathematical Details}:$ You could refer this link,which will support a nice derivation for the mathematical statement of this approximation.
Hope it helps.