Atkins hizo un buen trabajo explicando lo que es la dependencia paramétrica de la función de las ondas electrónicas (y la energía electrónica) de las coordenadas nucleares, pero no mencionó que en la derivación de la aproximación de Born-Oppenheimer esta dependencia se supone que es continua y diferenciable y que tanto el primer como el segundo derivado de estas cantidades con respecto a las coordenadas nucleares son en general no cero. En particular, para el sistema descrito en el texto, tenemos ∂ψ∂Zj≠0,∂2ψ∂Z2j≠0.
Ahora, cuando la solución de la forma (8.3) se sustituye en la ecuación de Schrödinger (8.2), el término que implica Te es trivial: ya que la función de la onda nuclear χ no es una función de las coordenadas electrónicas, es sólo una constante al diferenciarse con respecto a ellas, por lo que obtenemos Te(ψχ)=χTeψ. Pero el término que implica TN no se transforma trivialmente de manera similar, TN(ψχ)≠ψTNχ ya que ambos ψ y χ dependen de las coordenadas nucleares. Más bien, aplicando el regla del producto dos veces, tenemos \begin {alinear} \frac { \partial ^2}{ \partial Z_j^2} ( \psi \chi ) &= \frac { \partial }{ \partial Z_j} \left ( \frac { \partial }{ \partial Z_j} ( \psi \chi ) \right ) \\ &= \frac { \partial }{ \partial Z_j} \left ( \psi \frac { \partial \chi }{ \partial Z_j} + \chi \frac { \partial \psi }{ \partial Z_j} \right ) \\ &= \psi \frac { \partial ^2 \chi }{ \partial Z_j^2} + 2 \frac { \partial \psi }{ \partial Z_j} \frac { \partial \chi }{ \partial Z_j} + \chi \frac { \partial ^2 \psi }{ \partial Z_j^2} \ ~ - , \end {alinear} para que TN(ψχ)=−ψ∑j=1,2ℏ22mj∂2χ∂Z2j−∑j=1,2ℏ22mj(2∂ψ∂Zj∂χ∂Zj+χ∂2ψ∂Z2j), donde el primer término no es nada más que ψTNχ y los restos se designan como W .