# Potencial de Morse

El potencial de Morse, nombrado así en honor al físico Philip M. Morse, es un modelo conveniente para la energía potencial de una molécula diatómica. Es una mejor aproximación para la estructura vibracional de la molécula que el oscilador armónico cuántico porque incluye explícitamente los efectos de la ruptura del enlace, tales como la existencia de estados no enlazados. También da cuenta de la anarmonicidad de los enlaces reales y de la probabilidad de la transición no-cero para los sobretonos y las bandas de combinación. El potencial de Morse también puede usarse para modelar otras interacciones como la que se da entre un átomo y una superficie. Debido a su simplicidad (sólo tres parámetros ajustables), no es utilizada en la espectroscopía moderna. Sin embargo, su forma matemática inspiró el potencial MLR (Morse de largo alcance, Morse/Long-range), el cual es la función de energía potencial más popular para explicar los datos espectroscópicos.

## Función de energía potencial

La función de la energía potencial de Morse tiene la siguiente forma:
$V'(r)=D_{e}(1-e^{-\alpha (r-r_{e})})^{2}$
Donde $r$ es la distancia entre los átomos, $r_{e}$ es la distancia de enlace al equilibrio, $D_{e}$ es la profundidad del pozo (definida en relación con los átomos disociados) y $\alpha$ es una variable que controla el ancho del potencial (mientras más pequeño sea $\alpha$, más grande es el pozo). La energía de disociación del enlace puede ser calculada al restar la energía de punto cero $E_{0}$ de la profundidad del pozo. La constante de fuerza (rigidez) del enlace puede encontrarse mediante una expansión de Taylor de $V'(r)$ cuando $r=r_{e}$ para la segunda derivada de la función de energía potencial, de al cual puede demostrarse que el parámetro $\alpha$ es:
$\alpha ={\sqrt {\frac {k_{e}}{2D_{e}}}}$
Donde $k_{e}$ es la constante de fuerza en el mínimo del pozo.
Debido a que el cero de la energía potencial es arbitrario, la ecuación del potencial de Morse puede ser reescrita en cualquier número de formara agregando o restando un valor constante. Cuando se utiliza para modelar la interacción átomo-superficie, la energía cero puede ser redefinida para que el potencial tenga la forma:
$V'(r)=V'(r)-D_{e}=D_{e}(1-e^{-\alpha (r-r_{e})})^{2}-D_{e}$
Lo cual, normalmente se escribe como:
$V(r)=D_{e}(e^{-2\alpha (r-r_{e})}-2e^{-\alpha (r-r_{e})})$
Donde $r$ es la coordenada perpendicular a la superficie. Esta forma aproxima el cero a un valor de $r$ infinito e iguala $-D_{e}$ en su mínimo. Esto muestra que el potencial de Morse es la combinación de un término repulsivo de corto alcance (el primero) y un término atractivo de largo alcance (el segundo), análogo al potencial de Lennard-Jones.

## Estados vibracionales y energías

De la misma forma que oscilador armónico cuántico, las energías y eigenestados del potencial de Morse pueden encontrarse utilizando operadores. Una aproximación involucra aplicar el método de factorización al Hamiltoniano.

### Estados estacionarios

Para escribir los estados estacionarios en el potencial de Morse, por ejemplo las soluciones $\Psi _{n}(r)$ y $E_{n}$ de la siguiente ecuación de Schrödinger:
$\left(-{\frac {h^{2}}{2m}}{\frac {\partial ^{2}}{\partial r^{2}}})+V(r)\right)\Psi _{n}(r)=E_{n}\Psi _{n}(r)$
es conveniente introducir las variables:
$x=ar\|x_{e}=ar_{e}\|\lambda ={\frac {\sqrt {2mD_{e}}}{a\hbar }}\|\epsilon _{n}={\frac {2m}{a^{2}\hbar ^{2}}}E_{n}$
Por lo tanto, la ecuación de Scrhödinger, toma la forma:
$\left(-{\frac {\partial ^{2}}{\partial x^{2}}}+V(x)\right)\Psi _{n}(x)=\varepsilon _{n}\Psi _{n}(x)$
y
$V(x)=\lambda ^{2}\left(e^{-2(x-x_{e})}-2e^{-(x-x_{e})}\right)$

### Eigenvalores y eigenestados

Y sus eigenvalores y eigenestados pueden ser escritos:
$\varepsilon _{n}=-\left(\lambda -n-{\frac {1}{2}}\right)^{2}$
donde
$n=0,1,...,\left[\lambda -{\frac {1}{2}}\right]$
con $x$ representando la mayor integral menor a $x_{e}$
$\Psi _{n}(z)=N_{n}z^{\lambda -n-{\frac {1}{2}}}e^{-{\frac {1}{2}}z}L_{n}^{(2\lambda -2n-1)}(z)$
donde:
$z=2\lambda e^{-(x-x_{e})}$
$N_{n}=\left[{\frac {n!(2\lambda -2n-1)}{\Gamma (2\lambda -n))}}\right]^{\frac {1}{2}}$
y $L_{n}^{(a)}(z)$ es un polinomio de Laguerre generalizado de la forma:
$L_{n}^{(a)}(z)={\frac {z^{-a}e^{z}}{n!}}{\frac {d^{n}}{dz^{n}}}(z^{n+a}e^{-z})={\frac {\Gamma (a+n+1)/\Gamma (a+1)}{n!}}_{1}F_{1}(-1n,\alpha +),z)$

### Expresión analítica

También existe la siguiente expresión analítica para elementos de la matriz del operador coordinado:
$\left\langle \Psi _{m}|x|\Psi _{n}\right\rangle ={\frac {2(-1)^{m-n+1}}{(m-n)(2N-n-m)}}{\sqrt {\frac {(N-n)(N-m)\Gamma (2N-m+1)m!}{\Gamma (2N-n+1)n!}}}.$
que es válido siempre que $m>n$ y $N=\lambda -1/2$.

### Energía

Las eigenenergías de las variables iniciales tienen la forma:
$E_{n}=h\nu _{0}(n+1/2)-{\frac {[h\nu _{0}(n+1/2)]^{2}}{4D_{e}}}$
donde n es el número cuántico vibracional y $\nu _{0}$ tiene unidades de frecuencia. Esta última está relacionada matemáticamente a la masa de la partícula ($m$) y a las constantes de Morse mediante la expresión:
$\nu _{0}={\frac {a}{2\pi }}{\sqrt {\frac {2D_{e}}{m}}}$

### Espaciamiento de niveles energéticos

Mientras que el espaciamiento de energía entre niveles vibracionales en el oscilador armónico cuántico es constante (con un valor de $h\nu _{0}$), la energía entre niveles adyacentes en el oscilador de Morse disminuye conforme aumenta $\nu$. Por lo tanto, este espaciamiento entre niveles de Morse puede describirse con la expresión:
$E_{n+1}-E_{n}={\frac {h\nu _{0}-(n+1)(h\nu _{0})^{2}}{2D_{e}}}$
Esta tendencia coincide con la anarmonicidad de moléculas reales. Sin embargo, esta ecuación falla por valores arriba de $n_{m}$ donde $E(n_{m}+1)-E(n_{m})$ posee valores negativos o cero. Específicamente:
$n_{m}={\frac {2D_{e}-h\nu _{0}}{h\nu _{0}}}$
Esta falla se debe al número finito de niveles enlazados en el potencial de Morse, donde existe un nivel máximo $n_{m}$ que aún permanece enlazado y para energías superiores, el enlace se rompe y todos los niveles de energía posibles son permitidos pues la ecuación para $E_{n}$ cuantizada ya no es válida.
Por debajo de ese nivel máximo $n_{m}$, $E_{n}$ es una buena aproximación para el valor vibracional verdadero en moléculas diatómicas que no rotan. De hecho, el espectro molecular real generalmente se ajusta a la forma:
${\frac {E_{n}}{hc}}=\omega _{e}(n+1/2)-\omega _{e}\chi _{e}(n+1/2)^{2}$
en el que las constantes $\omega _{e}$ y $\omega _{e}\chi _{e}$ pueden relacionarse directamente con los parámetros del potencial.
Como puede observarse si se realiza el análisis dimensional, la última ecuación utiliza notación espectroscópica (por razones históricas) en las que $\omega _{e}$ representa el número de onda según la ecuación $E=hc\omega$, en lugar de utilizar frecuencia angular, según la ecuación $E=\hbar \omega$