Partícula en un anillo



La partícula en un anillo es ejemplo sencillo de sistema cuántico con propiedades interesantes. El modelo reproduce las características hipotéticas de una partícula libre que pudiera moverse de manera solamente a lo largo de un anillo (espacio topológico homeomorfo a S1) y de manera uniforme. Además en modelo aquí presentado ha encontrado aplicación en explicar la regla de Hückel sobre la estabilidad de los hidrocarburos aromáticos.

Tabla de contenidos

Descripción cuántica del sistema

Suponemos una partícula que se mueve libremente a lo largo de un anillo, las coordenadas de posición angular sobre el anillo y coordenadas cartesianas son:

x = R\cos(\varphi) \quad,\quad y = R\sin(\varphi)

donde R^2 = x^2 + y^2\,. Los operadores de momento lineal vienen dados por:

P_x=-i\hbar \frac{\partial}{\partial x}
P_y=-i\hbar \frac{\partial}{\partial y}

Buscamos la forma funciona de la energía clásica:

E_{clasica}=T+V=\frac{1}{2}mv^2+0=\frac{1}{2}m(V_x^2+V_y^2)=\frac{1}{2m}(P_x^2+P_y^2)


Operador hamiltoniano

Si llamamos \tilde{\psi}(x,y) a la función de onda de un estado estacionario del sistema, el operador hamiltoniano es aquel que cumple:

\hat{H}\tilde{\psi}=E\tilde{\psi}

Donde E es el valor de la energía del estado, que por ser estacionario estará perfectamente bien definida. Ahora vamos a procurar transformar nuestra expresión de coordenadas cartesianas: \tilde{\psi}(x,y) \, a coordenadas polares:\psi(r,\varphi):

\hat{H}_\varphi = -\frac{\hbar^2}{2m} \left(\frac{d^2}{dx^2}+\frac{d^2}{dy^2} \right) = -\frac{\hbar^2}{2m}\left(\frac{1}{R^2}\frac{\partial^2}{\partial \varphi^2} \right) =  -\frac{\hbar^2}{2I} \frac{d^2}{d\varphi^2}

Buscaremos las funciones propias del Hamiltoniano que son las que satisfacen la ecuación de Schrödinger independiente del tiempo:

-\frac{\hbar^2}{2I} \frac{d^2\psi}{d\varphi^2}=E\psi\Rightarrow \frac{d^2\psi}{d\varphi^2}=\frac{-2IE}{\hbar^2}\psi \qquad (1)

Soluciones de la ecuación de Schrödinger

Los posibles estados estacionarios del sistema son las soluciones de la ecuación (1) anterior. Por otro lado, cualquier estado no estacionario será combinación de estados estacionarios de diferente energía. Como candidatos canónicos para representar los estados estacionarios debo escoger funciones propias del hamiltoniano, que por tanto deben ser solución de ecuación (1). En un sistema cuántico, puede pasar que existan diversos estados estacionarios con un mismo valor de la energía, esto es lo que va a suceder en el caso de la partícula en un anillo. Cuando eso sucde se dice que un estado de energía concreto presenta degeneración (un término poco explicativo que se introdujo por motivos históricos relacionados con el átomo de hidrógeno, pero que ha sido mantenido a pesar de ser poco explicativo).

Puede verificarse fácilmente que las funciones trigonométricas seno y el coseno cumplen esos requisitos para ser mis funciones candidatas, aunque simplificar calculos, escogeremos exponenciales complejas que dependan de mi coordenada angular. Así buscamos candidatos a solución de la forma:

\psi_1(\varphi) = e^{+k\varphi} \qquad \psi_2(\varphi) = e^{-k\varphi}


Substituyendo esas funciones candidatas en la ecuación (1), se obtiene el valor necesario de k para que cualquiera de las dos sea solución:

(1) \ \frac{d^2\psi_1}{d\varphi^2} = k^2 e^{k\varphi}=\frac{-2IE}{\hbar^2} e^{k\varphi}\longrightarrow k^2=\frac{-2IE}{\hbar^2}\Rightarrow k=\sqrt{\frac{2IE}{\hbar^2}} i


(2) \ \frac{d^2\psi_2}{d\varphi^2} = k^2 e^{-k\varphi}=\frac{-2IE}{\hbar^2} e^{-k\varphi}\longrightarrow k^2=\frac{-2IE}{\hbar^2}\Rightarrow k=\sqrt{\frac{2IE}{\hbar^2}} i


Por lo tanto obtenemos nuestras funciones:

\psi_1(\varphi)=e^{i \sqrt{\frac{2IE}{\hbar^2}} \varphi} \qquad \psi_2(\varphi)=e^{-i \sqrt{\frac{2IE}{\hbar^2}} \varphi}


Mi función (o vector) de estado, por tanto, debe de cumplir que sea una combinación lineal de ambas funciones:

\psi(\varphi)=Ae^{i \sqrt{\frac{2IE}{\hbar^2}} \varphi}+Be^{-i \sqrt{\frac{2IE}{\hbar^2}} \varphi}

Número cuántico principal

Para simplificar, definimos una constante matemática n que vamos a llamar simplemente número cuántico principal como:

n=\sqrt{\frac{2IE_n}{\hbar^2}}


De nuevo, tendremos que imponer una condición para que mi función se comporte bien (well-behaviour function), en este caso, mi función tiene que ser continua en todos sus puntos y al dar una vuelta completa de mi partícula en el anillo, me tiene que dar el mismo valor, lo que nos lleva a imponer la siguiente condición de periodicidad:

\psi(\varphi)= \psi(\varphi+2\pi)


De hecho esa condición de periodicidad, no se da para cualquier valor del número cuántico n. Si estamos interesados sólo en los estados estacionarios que sí cumplen la condición de periodicidad y por tanto sí representan adecuadamente las restricciones físicas del problema debemos examinar que valores de n satisfacen la condición de periodicidad:

Ae^{i n \varphi}+Be^{-i n\varphi} = Ae^{i n (\varphi+2\pi)}+Be^{-i n(\varphi+2\pi)} \,
Ae^{i n \varphi}+Be^{-i n\varphi} = Ae^{i n \varphi}e^{i n 2\pi}+Be^{-i n\varphi}e^{-i n 2\pi} \,
Ae^{i n \varphi}(1-e^{i n 2\pi})+Be^{-i n\varphi}(1-e^{-i n 2\pi}) =0 \,
\ \ \ \{(1-e^{-i n 2\pi})= e^{-in2\pi}(e^{i n 2\pi}-1)= -e^{-in2\pi}(1-e^{i n 2\pi})\}
Ae^{i n \varphi}(1-e^{i n 2\pi})+Be^{-i n\varphi}(-e^{-i n 2\pi})(1-e^{i n 2\pi})= 0
(1-e^{i n 2\pi})(Ae^{in\varphi}-Be^{-i n\varphi}e^{-i n 2\pi}) =0
(1-e^{i n 2\pi})= 0\Rightarrow e^{in 2\pi} = 1
cos(2n\pi)+ i sen(2n\pi)= 1 \ \ \ cos(2n\pi)= 1\ \ \ sen(2n\pi)= 0


Y de la última ecuación concluimos que sólo los valores enteros de n satisfacen la ecuación, los posibles valores del número cuántico principal son n = 0, 1, 2, 3... y n = -1, -2, -3... Es interesante notar que como consecuencia de exigir la condición periodicidad la energía del sistema está cuantizada:

\sqrt{\frac{2IE_n}{\hbar^2}}= n\Rightarrow E_n=\frac{\hbar^2}{2I}\ n^2


Degeneración

El número cuántico n puede tomar, en este caso, el valor 0, porque no anula a mi función en el espacio. Y por otra parte tenemos que para dos números cuánticos que sean iguales y opuestos, se ve que obtenemos la misma energía (al estar n al cuadrado): a esto se le llama que los estados de energía definida son degenerados ya que existen varios estados de la misma energía. Sin embargo, esos estados de misma energía no son del todo idénticos como veremos, puesto que sobre ellos podemos medir otras magnitudes físicas diferentes de la energía y podemos ver que arrojan valores diferentes, lo cual significa que existe un procedimiento físico para distinguir entre estados "degenerados" de la misma energía. Esto se puede comprobar introduciendo el momento angular.


Momento angular

A continuación veremos el momento angular de la partícula, es decir, vamos a ver si el operador \hat{L_z} se comporta correctamente, partiendo de las expresiones clásicas:

\vec{L}=\vec{r}\times \vec{p} \qquad L_z = xP_y-yP_x

Construimos el correspondiente operador cuántico:

\hat{L_z}=\frac{\hbar}{i} \ (x\frac{\partial}{\partial y}-y\frac{\partial}{\partial x} )= \frac{\hbar}{i}\frac{\partial}{\partial \varphi}
\hat{L_z}\psi=l_z \psi
\hat{L_z}\psi=\frac{\hbar}{i}\frac{\partial}{\partial \varphi} \ (Ae^{im\varphi}+Be^{-in\varphi} )=\frac{\hbar}{i}(inAe^{in\varphi}-inBe^{-in\varphi} )=\hbar n(Ae^{in\varphi}-Be^{-in\varphi})\not = L_z \psi

No es una función propia. Pero yo tengo infinitas funciones que determinan mi estado, como tal podré escoger un subconjunto propio de ellas en este caso:

Si \ \ B=0: \ \ \ \ L_z\psi=\frac{\hbar}{i}\frac{\partial}{\partial \varphi} A e^{in\varphi}=\hbar n A e^{in\varphi}=\hbar n\psi

Y compruebo que efectivamente existe un subconjunto que es propio, en este caso sí es función propia de Lz:

\psi(\varphi)=Ae^{in\varphi} \ \ / \ l_z=n\hbar

Con esta magnitud acabo de romper la degeneración de la energía, pudiendo distinguir dos estados que tienen distinto valor de momento angular.

Normalización de los estados propios de momento angular

Por último, para poder determinar el valor de A, y que nuestra función en el espacio valga 1, normalizamos a todo el espacio:

\int_0^{2\pi}|\psi|^2 d\varphi=1
|\psi|^2=\psi^*\psi=(A e^{-in\varphi})^*Ae^{in\varphi}=A^* e^{-in\varphi}Ae^{in\varphi}
|A|^2\int_0^{2\pi}d\varphi=1\longrightarrow A=\frac{1}{\sqrt{2\pi}}

Con la cual la función normaliada salvo fase es:

\psi(\varphi)=\frac{1}{\sqrt{2\pi}}e^{in\varphi}


Estados estacionarios generales del sistema

Puede comprobarse que cualquier otro estado estacionario del sistema tiene la forma:

\psi(\varphi)=\frac{1}{\sqrt{2\pi}} \left( \cos\alpha e^{i(n\varphi+\beta)} + \sin\alpha e^{-i(n\varphi+\beta)} \right)


Este estado tiene la propiedad interesante de que a pesar de que tiene una energía bien definida, su momento angular Lz no está bien definido, sino que una medida de esa magnitud con una probabilidad p1 da el valor +nh/2π y con una probabilidad p2 da el valor -nh/2π, cumpliéndose además:

p_1 = |\cos \alpha|^2 \qquad p_2 = |\sin \alpha|^2

Aplicación a los hidrocarburos aromáticos

En química orgánica, los hidrocarburos aromáticos como el el benceno y otros, contienen estructuras en forma de anillo formado por cinco o seis átomos de carbono. Los experimentos muestran que estos compuestos químicos son extraordinariamente estables, debido a que de acuerdo con la discusión anterior los electrones se comportan como si estuvieran girando en ambas direcciones y están altamente deslocalizados.

De acuerdo con el cálculo cuántico presentado anteriormente, rellenar todos los niveles de energía hasta el nivel n-ésimo requiere 2·(2n+1) electrones (donde el factor 2 inicial procede del hecho de que los electrones tienen dos posibles valores de espín). Esa es precisamente la Regla de Hückel que afirma que un exceso de 4n+2 electrones en un anillo de Kekulé produce un compuesto aromático excepcionalmente estable.

 
Este articulo se basa en el articulo Partícula_en_un_anillo publicado en la enciclopedia libre de Wikipedia. El contenido está disponible bajo los términos de la Licencia de GNU Free Documentation License. Véase también en Wikipedia para obtener una lista de autores.
Su navegador no está actualizado. Microsoft Internet Explorer 6.0 no es compatible con algunas de las funciones de Chemie.DE.