Chemistry · Libro 4 · Bachelor Year 3

Química universitaria — tercer año

Química universitaria — tercer año · Bachelor Year 3

1Mecánica cuántica para químicos: sistemas modelo

Sobre una mesa de laboratorio hay tres pequeños matraces, cada uno con una disolución de un colorante de cianina en etanol: el primero es magenta, el segundo azul, el tercero de un azul verdoso intenso que absorbe en el rojo lejano. Las tres moléculas están construidas igual —dos anillos nitrogenados idénticos unidos por una cadena de átomos de carbono— y solo difieren en la longitud de esa cadena, de dos en dos átomos de carbono. Al alargar la cadena, el color se desplaza hacia el rojo. Un químico de 1949 explicó la serie con uno de los problemas más sencillos de la mecánica cuántica: un electrón libre de moverse a lo largo de una línea de longitud fija, una “partícula en una caja”. Este capítulo monta la maquinaria que hay detrás de esa explicación —operadores, valores propios, la ecuación de Schrödinger— y resuelve los sistemas modelo que la química usa a diario: la caja, el oscilador armónico de un enlace que vibra, el rotor rígido de una molécula que gira, el átomo de hidrógeno y la barrera que una partícula ligera puede atravesar sin tener la energía necesaria para subirla.

Lo que ya sabes

El volumen del primer año mostró que la energía de un átomo está cuantizada (niveles de energía, estados fundamental y excitados, las rayas del hidrógeno) y que un fotón transporta E=hν=hc/λE = h\nu = hc/\lambda. El volumen del segundo año describió un electrón mediante una función de onda ψ\psi cuyo cuadrado es una densidad de probabilidad, separó los orbitales del hidrógeno en partes radial y angular y contó sus nodos; tomó como dadas las funciones de onda del hidrógeno. Aquí se deducen, en la sección 5. De la física: la cantidad de movimiento p=mvp = mv y la longitud de onda de de Broglie λ=h/p\lambda = h/p.

Disoluciones de tres colorantes de cianina de una misma familia. Cada par de átomos de carbono añadido a la cadena desplaza la absorción hacia el rojo.
Disoluciones de tres colorantes de cianina de una misma familia. Cada par de átomos de carbono añadido a la cadena desplaza la absorción hacia el rojo.

1.1 Operadores, valores propios y la ecuación de Schrödinger

La mecánica clásica asigna a una partícula una posición y una cantidad de movimiento en cada instante. La mecánica cuántica le asigna un estado, la función de onda ψ\psi, y sustituye cada magnitud medible por una operación efectuada sobre ψ\psi.

Definición 1.1 (Operador, función propia, valor propio)

Un operador A^\hat A es una regla que transforma una función en otra función: A^f=g\hat A f = g. Es lineal si A^(λf+μg)=λA^f+μA^g\hat A(\lambda f + \mu g) = \lambda\hat Af + \mu\hat Ag. Una función no nula ff tal que A^f=af\hat Af = a f para un número aa es una función propia de A^\hat A, y aa es el valor propio correspondiente.

Ejemplo 1.2 (Dos operadores: posición y cantidad de movimiento)

En una dimensión, el operador posición multiplica por xx, x^f=xf\hat x f = xf, y el operador cantidad de movimiento deriva, p^f=−iℏ  ⁣df/ ⁣dx\hat p f = -\iu\hbar\, \dd f/\dd x, con ℏ=h/2π\hbar = h/2\pi. La función eikx\eu^{\iu kx} es una función propia de p^\hat p con valor propio ℏk\hbar k: una onda de longitud de onda λ=2π/k\lambda = 2\pi/k tiene cantidad de movimiento h/λh/\lambda, la relación de de Broglie. La función sin⁡kx\sin kx no es función propia de p^\hat p (su derivada es un coseno), pero sí lo es de p^2=−ℏ2  ⁣d2/ ⁣dx2\hat p^2 = -\hbar^2\,\dd^2/\dd x^2, con valor propio ℏ2k2\hbar^2k^2.

Definición 1.3 (Operador hermítico, valor esperado)

Se escribe ⟨f∣g⟩=∫f∗g  ⁣dτ\langle f|g\rangle = \int f^*g\,\dd\tau para dos funciones que se anulan en los límites de la región. Un operador es hermítico si ⟨f∣A^g⟩=⟨A^f∣g⟩\langle f|\hat Ag\rangle = \langle \hat Af|g\rangle para todas esas ff y gg. Para un estado normalizado ψ\psi (⟨ψ∣ψ⟩=1\langle\psi|\psi\rangle = 1), el valor esperado de A^\hat A es ⟨A⟩=⟨ψ∣A^ψ⟩\langle A\rangle = \langle\psi|\hat A\psi\rangle: la media de muchas medidas efectuadas sobre sistemas preparados de forma idéntica.

Las reglas de trabajo de la química cuántica —sus postulados— pueden enunciarse ahora de un tirón: toda magnitud medible está representada por un operador hermítico; una medida da uno de sus valores propios; y si ψ\psi es función propia de A^\hat A, toda medida sobre ella da el mismo valor, su valor propio.

Teorema 1.4 (Operadores hermíticos)

Los valores propios de un operador hermítico son reales. Dos funciones propias con valores propios distintos son ortogonales: ⟨f1∣f2⟩=0\langle f_1|f_2\rangle = 0.

Demostración. Sea A^f=af\hat Af = af. Entonces ⟨f∣A^f⟩=a⟨f∣f⟩\langle f|\hat Af\rangle = a\langle f|f\rangle y ⟨A^f∣f⟩=a∗⟨f∣f⟩\langle\hat Af|f\rangle = a^*\langle f|f\rangle. El operador es hermítico, luego ambas cantidades son iguales, y ⟨f∣f⟩>0\langle f|f\rangle > 0: a=a∗a = a^*. Sean ahora A^f1=a1f1\hat Af_1 = a_1f_1 y A^f2=a2f2\hat Af_2 = a_2f_2 con a1≠a2a_1 \ne a_2, ambos reales. Entonces ⟨f1∣A^f2⟩=a2⟨f1∣f2⟩\langle f_1|\hat Af_2\rangle = a_2\langle f_1|f_2\rangle y ⟨A^f1∣f2⟩=a1⟨f1∣f2⟩\langle\hat Af_1| f_2\rangle = a_1\langle f_1|f_2\rangle; su igualdad da (a2−a1)⟨f1∣f2⟩=0(a_2 - a_1) \langle f_1|f_2\rangle = 0, de donde ⟨f1∣f2⟩=0\langle f_1|f_2\rangle = 0. ∎

Un valor medido es un número real, y por eso los observables deben ser hermíticos; la ortogonalidad es lo que hace de los orbitales de un átomo, o de los niveles de una caja, un conjunto limpio de estados independientes.

Definición 1.5 (Conmutador)

El conmutador de dos operadores es [A^,B^]=A^B^−B^A^[\hat A,\hat B] = \hat A\hat B - \hat B\hat A. Los operadores conmutan si su conmutador es nulo.

Ejemplo 1.6 (El conmutador canónico)

Para cualquier función ff, x^p^f=−iℏxf′\hat x\hat pf = -\iu\hbar xf' y p^x^f=−iℏ(f+xf′)\hat p\hat xf = -\iu\hbar(f + xf'). La diferencia es [x^,p^]f=iℏf[\hat x,\hat p]f = \iu\hbar f, luego [x^,p^]=iℏ[\hat x,\hat p] = \iu\hbar: la posición y la cantidad de movimiento no conmutan. Esta única relación sustenta el principio de incertidumbre y, más abajo, todo el espectro del oscilador armónico.

Proposición 1.7 (Operadores que conmutan)

Si A^\hat A y B^\hat B conmutan y ff es una función propia de A^\hat A con un valor propio no degenerado aa, entonces ff es también función propia de B^\hat B.

Demostración. A^(B^f)=B^A^f=a(B^f)\hat A(\hat Bf) = \hat B\hat Af = a(\hat Bf): la función B^f\hat Bf es una función propia de A^\hat A con el mismo valor propio aa, o es nula. Como el valor propio no es degenerado, B^f\hat Bf es un múltiplo de ff, digamos bfbf. ∎

El operador de la energía es el hamiltoniano. Para una partícula de masa mm en un potencial VV, la energía clásica p2/2m+Vp^2/2m + V se convierte en un operador.

Definición 1.8 (Operador hamiltoniano)

El operador hamiltoniano de una partícula de masa mm que se mueve en un potencial V(r)V(\mathbf r) es

H^=−ℏ22m∇2+V(r),∇2=∂2∂x2+∂2∂y2+∂2∂z2,\hat H = -\frac{\hbar^2}{2m}\nabla^2 + V(\mathbf r), \qquad \nabla^2 = \frac{\partial^2}{\partial x^2} + \frac{\partial^2}{\partial y^2} + \frac{\partial^2}{\partial z^2},

y, para varias partículas, la suma de sus términos cinéticos y de todas sus energías potenciales.

Teorema 1.9 (La ecuación de Schrödinger)

Los estados de energía definida de un sistema —sus estados estacionarios— son las funciones propias de su hamiltoniano, y sus energías permitidas son los valores propios:

H^ψ=Eψ,\hat H\psi = E\psi,

la ecuación de Schrödinger independiente del tiempo. Una ψ\psi aceptable es univaluada, continua y de cuadrado integrable (puede normalizarse).

Estatus. Es un postulado, justificado por sus consecuencias: los niveles de todos los sistemas resueltos a continuación concuerdan con la espectroscopia. La ecuación dependiente del tiempo, de la que esta es el caso estacionario, pertenece a la física. ∎

Son las condiciones de contorno las que cuantizan: una ecuación diferencial tiene soluciones para cualquier EE, pero solo un conjunto discreto de ellas se mantiene finito y continuo. La química se reduce entonces a elegir un potencial, resolver y leer los valores propios.

Método 1.10 (Resolver un modelo unidimensional)

  1. Escribir el potencial V(x)V(x) y el hamiltoniano; dividir el espacio en regiones donde VV sea sencillo.
  2. Resolver −ℏ22mψ′′+Vψ=Eψ-\frac{\hbar^2}{2m}\psi'' + V\psi = E\psi en cada región: senos y cosenos donde E>VE > V, exponenciales reales donde E<VE < V.
  3. Imponer las condiciones: ψ=0\psi = 0 donde VV es infinito; ψ\psi y ψ′\psi' continuas en un escalón finito; ψ→0\psi \to 0 en el infinito.
  4. Las condiciones solo se cumplen para ciertos valores de EE: esos son los niveles.
  5. Normalizar cada función propia; contar sus nodos como comprobación (el nn-ésimo nivel de un problema unidimensional tiene n−1n - 1 nodos).

1.2 La partícula en una caja

Definición 1.11 (Partícula en una caja, modelo del electrón libre)

Una partícula en una caja se mueve libremente (V=0V = 0) entre x=0x = 0 y x=Lx = L y no puede salir (V=∞V = \infty fuera). El modelo del electrón libre de una molécula conjugada trata sus electrones π\pi como partículas independientes en una caja cuya longitud es la de la cadena conjugada.

Teorema 1.12 (Niveles de la caja)

Los niveles y las funciones propias normalizadas de una partícula de masa mm en una caja de longitud LL son

En=n2h28mL2,ψn(x)=2L sin⁡nπxL,n=1,2,3,…E_n = \frac{n^2h^2}{8mL^2}, \qquad \psi_n(x) = \sqrt{\frac2L}\,\sin\frac{n\pi x}{L}, \qquad n = 1, 2, 3, \dots

Demostración. Dentro, ψ′′=−k2ψ\psi'' = -k^2\psi con k2=2mE/ℏ2k^2 = 2mE/\hbar^2, luego ψ=Asin⁡kx+Bcos⁡kx\psi = A\sin kx + B\cos kx. Fuera ψ=0\psi = 0 y la continuidad da ψ(0)=0\psi(0) = 0, de donde B=0B = 0, y ψ(L)=0\psi(L) = 0, de donde sin⁡kL=0\sin kL = 0 con A≠0A \ne 0: kL=nπkL = n\pi, con nn un entero positivo (n=0n = 0 da ψ=0\psi = 0; los nn negativos repiten las mismas funciones). Entonces E=ℏ2k2/2m=n2π2ℏ2/2mL2=n2h2/8mL2E = \hbar^2k^2/2m = n^2\pi^2\hbar^2/2mL^2 = n^2h^2/8mL^2. Por último, ∫0LA2sin⁡2(nπx/L)  ⁣dx=A2L/2=1\int_0^L A^2\sin^2(n\pi x/L)\,\dd x = A^2L/2 = 1 da A=2/LA = \sqrt{2/L}. ∎

Tres rasgos se trasladan a todo sistema ligado: la energía más baja no es nula (la partícula nunca puede estar en reposo), los niveles se separan al crecer nn y se apiñan al agrandarse la caja (E∝1/L2E \propto 1/L^2): una caja grande es casi clásica.

Los cuatro primeros niveles de una partícula en una caja, con cada función dibujada sobre su propio nivel (en trazo discontinuo). Izquierda: _n, con n - 1 nodos dentro de la caja. Derecha: la densidad de probabilidad | _n|2. Los cuatro primeros niveles de una partícula en una caja, con cada función dibujada sobre su propio nivel (en trazo discontinuo). Izquierda: _n, con n - 1 nodos dentro de la caja. Derecha: la densidad de probabilidad | _n|2.
Los cuatro primeros niveles de una partícula en una caja, con cada función dibujada sobre su propio nivel (en trazo discontinuo). Izquierda: ψn\psi_n, con n−1n - 1 nodos dentro de la caja. Derecha: la densidad de probabilidad ∣ψn∣2|\psi_n|^2.

Proposición 1.13 (La caja cúbica)

En una caja cúbica de arista LL, ψ=ψnx(x)ψny(y)ψnz(z)\psi = \psi_{n_x}(x)\psi_{n_y}(y)\psi_{n_z}(z) y E=(nx2+ny2+nz2) h2/8mL2E = (n_x^2 + n_y^2 + n_z^2)\,h^2/8mL^2. Ternas distintas con la misma suma de cuadrados dan la misma energía.

Demostración. Con V=0V = 0 dentro, H^\hat H es una suma de tres operadores unidimensionales, uno por coordenada. El producto de tres funciones propias unidimensionales es una función propia de la suma, con la suma de los tres valores propios; se anula en cada cara, como se exige. ∎

Definición 1.14 (Niveles degenerados)

Las funciones propias linealmente independientes con el mismo valor propio son niveles degenerados; su número es la degeneración gg de esa energía.

Ejemplo 1.15 (Degeneración y simetría)

El nivel nx2+ny2+nz2=6n_x^2 + n_y^2 + n_z^2 = 6 del cubo procede de (2,1,1)(2,1,1), (1,2,1)(1,2,1) y (1,1,2)(1,1,2): g=3g = 3. Si se estira ligeramente la caja según zz, la tercera función se aparta de las otras dos: la degeneración procedía de la simetría del cubo. El mismo vínculo entre simetría y degeneración organiza los Capítulos 4 y 5.

Colorantes conjugados

En un colorante de cianina simétrico, una cadena de jj átomos une los dos nitrógenos (j=5j = 5, 7, 9 para los tres colorantes de la escena inicial), y sus electrones π\pi están deslocalizados a lo largo de ella. Una cadena de jj átomos lleva N=j+1N = j + 1 electrones π\pi: uno por carbono, más el par no enlazante de un nitrógeno, compartido por todo el catión.

Proposición 1.16 (Longitud de onda de absorción del electrón libre)

Si NN electrones π\pi (un número par) llenan los niveles de una caja de longitud LL, dos por nivel, la absorción de mayor longitud de onda, del nivel lleno más alto n=N/2n = N/2 al siguiente, está en

λ=8mcL2h(N+1).\lambda = \frac{8mcL^2}{h(N+1)}.

Demostración. ΔE=EN/2+1−EN/2=[(N/2+1)2−(N/2)2] h2/8mL2=(N+1) h2/8mL2\Delta E = E_{N/2+1} - E_{N/2} = [(N/2+1)^2 - (N/2)^2]\,h^2/8mL^2 = (N+1)\,h^2/8mL^2, y λ=hc/ΔE\lambda = hc/\Delta E. ∎

Ejemplo 1.17 (La caja desnuda falla, una caja alargada funciona)

Se toma cada enlace de la cadena igual al enlace C–C del benceno, l=139.7 pml = 139.7\,\mathrm{pm}, de modo que la cadena desnuda mide (j−1)l(j - 1)l. Para los tres colorantes (N=6N = 6, 8, 10) la fórmula da 147, 257 y 374 nm, muy por debajo de los máximos medidos, 524, 603 y 711 nm en etanol. Los electrones no se detienen en los nitrógenos: se extienden por los anillos. Si se alarga la caja una longitud δ\delta en cada extremo y se ajusta δ\delta con el colorante central, δ=222 pm\delta = 222\,\mathrm{pm}, alrededor de enlace y medio. El mismo δ\delta predice entonces 474 y 732 nm para los otros dos, con un error del 10 % y del 3 %. Una sola longitud ajustable explica una familia de colores (figura siguiente; los valores se calculan en el script de la figura).

Máximos de absorción de los tres colorantes de cianina frente al modelo del electrón libre: la cadena desnuda y la caja alargada en  en cada extremo, con  ajustado sobre el colorante central.
Máximos de absorción de los tres colorantes de cianina frente al modelo del electrón libre: la cadena desnuda y la caja alargada en δ\delta en cada extremo, con δ\delta ajustado sobre el colorante central.

En el laboratorio — Medir una serie de colorantes

Se pesan los tres colorantes (unos pocos miligramos), se disuelven en etanol y se diluyen hasta que la absorbancia en el máximo quede entre 0.3 y 1. Se registra el espectro de cada uno entre 400 y 800 nm en una cubeta de 1 cm frente a un blanco de etanol puro, y se lee λmax⁡\lambda_{\max} en la cima de la banda principal. Los colorantes de cianina son sólidos que manchan e irritan: guantes, y pesada en vitrina de gases.

1.3 El oscilador armónico

Cerca del fondo de cualquier pozo de potencial suave, V≈12k(x−xe)2V \approx \frac12 k(x - x_e)^2: las pequeñas vibraciones de un enlace, de un átomo en un cristal o de una molécula en una jaula son todas aproximadamente armónicas.

Definición 1.18 (Oscilador armónico, constante de fuerza, masa reducida, energía de punto cero)

Un oscilador armónico es una partícula de masa mm en el potencial V=12kx2V = \frac12 kx^2; kk es la constante de fuerza, y ω=k/m\omega = \sqrt{k/m} su frecuencia angular clásica. Una molécula diatómica de masas atómicas m1m_1, m2m_2 vibra como una sola partícula de masa reducida μ=m1m2/(m1+m2)\mu = m_1m_2/(m_1 + m_2) en el potencial de su enlace. La energía del nivel más bajo de un oscilador es su energía de punto cero.

Definición 1.19 (Operadores escalera)

Con x0=ℏ/mωx_0 = \sqrt{\hbar/m\omega}, los operadores escalera son

a^=12(x^x0+ix0p^ℏ),a^†=12(x^x0−ix0p^ℏ).\hat a = \frac{1}{\sqrt2}\Big(\frac{\hat x}{x_0} + \frac{\iu x_0\hat p}{\hbar}\Big), \qquad \hat a^\dagger = \frac{1}{\sqrt2}\Big(\frac{\hat x}{x_0} - \frac{\iu x_0\hat p}{\hbar}\Big).

Teorema 1.20 (Niveles del oscilador armónico)

Los niveles de un oscilador armónico son

Ev=(v+12)ℏω=(v+12)hν,v=0,1,2,…,E_v = \big(v + \tfrac12\big)\hbar\omega = \big(v + \tfrac12\big)h\nu, \qquad v = 0, 1, 2, \dots,

equidistantes en hνh\nu, con la energía de punto cero 12hν\frac12h\nu.

Demostración. A partir de [x^,p^]=iℏ[\hat x,\hat p] = \iu\hbar se calcula [a^,a^†]=1[\hat a,\hat a^\dagger] = 1 y H^=ℏω(a^†a^+12)\hat H = \hbar\omega(\hat a^\dagger\hat a + \frac12). Sea N^=a^†a^\hat N = \hat a^\dagger\hat a. Entonces [N^,a^]=−a^[\hat N,\hat a] = -\hat a y [N^,a^†]=a^†[\hat N,\hat a^\dagger] = \hat a^\dagger: si N^ψ=νψ\hat N\psi = \nu\psi, entonces N^(a^ψ)=(ν−1)a^ψ\hat N(\hat a\psi) = (\nu - 1)\hat a\psi y N^(a^†ψ)=(ν+1)a^†ψ\hat N(\hat a^\dagger\psi) = (\nu + 1)\hat a^\dagger\psi; a^\hat a baja el valor propio en una unidad, a^†\hat a^\dagger lo sube. Pero ν=⟨ψ∣a^†a^ψ⟩/⟨ψ∣ψ⟩=∥a^ψ∥2/∥ψ∥2≥0\nu = \langle\psi|\hat a^\dagger\hat a\psi\rangle/\langle\psi|\psi \rangle = \|\hat a\psi\|^2/\|\psi\|^2 \ge 0: el descenso debe detenerse, lo que solo ocurre en una función con a^ψ0=0\hat a\psi_0 = 0, de valor propio ν=0\nu = 0. Los valores propios de N^\hat N son, por tanto, v=0,1,2,…v = 0, 1, 2, \dots, y los de H^\hat H son (v+12)ℏω(v + \frac12)\hbar\omega. ∎

Proposición 1.21 (Las primeras funciones de onda)

En la coordenada reducida y=x/x0y = x/x_0,

ψ0∝e−y2/2,ψ1∝2y e−y2/2,ψ2∝(4y2−2)e−y2/2,ψ3∝(8y3−12y)e−y2/2,\psi_0 \propto \eu^{-y^2/2}, \quad \psi_1 \propto 2y\,\eu^{-y^2/2}, \quad \psi_2 \propto (4y^2 - 2)\eu^{-y^2/2}, \quad \psi_3 \propto (8y^3 - 12y)\eu^{-y^2/2},

los polinomios de Hermite multiplicados por una gaussiana. ψv\psi_v es par para vv par e impar para vv impar.

Demostración. a^ψ0=0\hat a\psi_0 = 0 se escribe yψ0+ ⁣dψ0/ ⁣dy=0y\psi_0 + \dd\psi_0/\dd y = 0, cuya solución es la gaussiana. Cada función siguiente es a^†\hat a^\dagger aplicado a la anterior, a^†=12(y− ⁣d/ ⁣dy)\hat a^\dagger = \frac{1}{\sqrt2}(y - \dd/\dd y), que multiplica el polinomio por 2y2y y le resta su derivada: 1→2y→4y2−2→8y3−12y1 \to 2y \to 4y^2 - 2 \to 8y^3 - 12y (salvo factores). Cambiar yy por −y-y cambia el signo de yy y de  ⁣d/ ⁣dy\dd/\dd y, de modo que a^†\hat a^\dagger invierte la paridad en cada paso. ∎

El potencial armónico, sus cuatro primeros niveles y sus funciones de onda, cada una dibujada sobre su nivel. Las funciones de onda llegan más allá de los puntos de retroceso clásicos (y = ±1 para v = 0, puntos rojos), donde una partícula clásica se detendría.
El potencial armónico, sus cuatro primeros niveles y sus funciones de onda, cada una dibujada sobre su nivel. Las funciones de onda llegan más allá de los puntos de retroceso clásicos (y=±1y = \pm1 para v=0v = 0, puntos rojos), donde una partícula clásica se detendría.

Ejemplo 1.22 (Energías de punto cero de HX2\ce{H2} y DX2\ce{D2})

Los números de onda armónicos son ω~e=4401.2 cm−1\tilde\omega_e = 4401.2\,\mathrm{cm}^{-1} para HX2\ce{H2} y 3115.5 cm−13115.5\,\mathrm{cm}^{-1} para DX2\ce{D2}, en la razón 1.4131.413, próxima a 2\sqrt2: mismo enlace y misma constante de fuerza, masa reducida doble. Las energías de punto cero armónicas son la mitad, 2200.6 y 1557.8 cm−11557.8\,\mathrm{cm}^{-1}, 26.3 kJ/mol26.3\,\mathrm{kJ}/\mathrm{mol} y 18.6 kJ/mol18.6\,\mathrm{kJ}/\mathrm{mol}. Una molécula nunca puede perder esta energía; la diferencia entre isótopos es el origen de los efectos isotópicos de los Capítulos 11 y 12.

1.4 El rotor rígido

Una molécula diatómica también gira en torno a su centro de masas. Tratada como dos masas a una distancia fija rr, es un rotor rígido de momento de inercia I=μr2I = \mu r^2, la masa reducida a la distancia rr del eje.

Definición 1.23 (Rotor rígido, armónicos esféricos)

Un rotor rígido es un cuerpo de forma fija que gira libremente; para una molécula lineal, H^=L^2/2I\hat H = \hat L^2/2I, donde L^2\hat L^2 es el operador del cuadrado del momento angular. Sus funciones propias, que solo dependen de los ángulos θ\theta y ϕ\phi, son los armónicos esféricos Yl,m(θ,ϕ)Y_{l,m}(\theta,\phi), indexados por l=0,1,2,…l = 0, 1, 2, \dots y m=−l,…,lm = -l, \dots, l; para una molécula, ll se escribe JJ.

El rotor en un anillo (rotación en un plano, en torno a un eje fijo) se resuelve primero; contiene lo esencial.

Teorema 1.24 (El anillo)

Una partícula de momento de inercia II que gira en un plano tiene los niveles Em=m2ℏ2/2IE_m = m^2\hbar^2/2I con m=0,±1,±2,…m = 0, \pm1, \pm2, \dots; todos los niveles salvo m=0m = 0 son doblemente degenerados.

Demostración. La única coordenada es el ángulo ϕ\phi, y H^=−(ℏ2/2I)  ⁣d2/ ⁣dϕ2\hat H = -(\hbar^2/2I)\, \dd^2/\dd\phi^2. Las soluciones de ψ′′=−(2IE/ℏ2)ψ\psi'' = -(2IE/\hbar^2)\psi son eimϕ\eu^{\iu m\phi} con m2=2IE/ℏ2m^2 = 2IE/\hbar^2. Una función univaluada debe tomar el mismo valor en ϕ\phi y en ϕ+2π\phi + 2\pi: e2πim=1\eu^{2\pi\iu m} = 1, luego mm es entero. Las funciones con mm y −m-m tienen la misma energía y son independientes. ∎

Teorema 1.25 (Niveles del rotor rígido)

L^2YJ,m=J(J+1)ℏ2 YJ,m\hat L^2Y_{J,m} = J(J+1)\hbar^2\,Y_{J,m} y L^zYJ,m=mℏ YJ,m\hat L_zY_{J,m} = m\hbar\, Y_{J,m}, de modo que los niveles de un rotor rígido lineal son

EJ=ℏ22I J(J+1),J=0,1,2,…,E_J = \frac{\hbar^2}{2I}\,J(J+1), \qquad J = 0, 1, 2, \dots,

cada uno con degeneración gJ=2J+1g_J = 2J + 1.

Demostración. Admitido a este nivel. ∎

La demostración general, mediante operadores escalera del momento angular, corresponde a cursos más avanzados; he aquí una comprobación. En coordenadas esféricas,

L^2=−ℏ2[1sin⁡θ∂∂θ(sin⁡θ∂∂θ)+1sin⁡2θ∂2∂ϕ2].\hat L^2 = -\hbar^2\Big[\frac{1}{\sin\theta}\frac{\partial}{\partial\theta} \Big(\sin\theta\frac{\partial}{\partial\theta}\Big) + \frac{1}{\sin^2\theta}\frac{\partial^2}{\partial\phi^2}\Big].

Para Y1,0∝cos⁡θY_{1,0} \propto \cos\theta, el corchete vale 1sin⁡θ(−2sin⁡θcos⁡θ)=−2cos⁡θ\frac{1}{\sin\theta}(-2\sin\theta \cos\theta) = -2\cos\theta, luego L^2Y1,0=2ℏ2Y1,0=1(1+1)ℏ2Y1,0\hat L^2Y_{1,0} = 2\hbar^2Y_{1,0} = 1(1+1)\hbar^2 Y_{1,0}. Para Y1,±1∝sin⁡θ e±iϕY_{1,\pm1} \propto \sin\theta\,\eu^{\pm\iu\phi} el mismo cálculo vuelve a dar 2ℏ22\hbar^2, y L^z=−iℏ ∂/∂ϕ\hat L_z = -\iu\hbar\,\partial/\partial \phi da ±ℏ\pm\hbar.

Proposición 1.26 (Armónicos esféricos reales)

Para l=1l = 1, las combinaciones Y1,0Y_{1,0}, (Y1,−1−Y1,1)/2(Y_{1,-1} - Y_{1,1})/\sqrt2 e i(Y1,−1+Y1,1)/2\iu(Y_{1,-1} + Y_{1,1})/\sqrt2 son reales y proporcionales a z/rz/r, x/rx/r e y/ry/r; para l=2l = 2 la misma construcción da funciones proporcionales a (3z2−r2)(3z^2 - r^2), xzxz, yzyz, xyxy y x2−y2x^2 - y^2, divididas por r2r^2.

Demostración. Y1,±1∝∓sin⁡θ e±iϕY_{1,\pm1} \propto \mp\sin\theta\,\eu^{\pm\iu\phi} con el convenio de signos habitual, de modo que su diferencia dividida por 2\sqrt2 es ∝sin⁡θcos⁡ϕ=x/r\propto\sin\theta\cos\phi = x/r y la suma ponderada con i\iu es ∝sin⁡θsin⁡ϕ=y/r\propto\sin\theta\sin\phi = y/r; cos⁡θ=z/r\cos\theta = z/r. Toda combinación de funciones propias degeneradas sigue siendo función propia de L^2\hat L^2 (con el mismo ll), aunque ya no de L^z\hat L_z. El caso l=2l = 2 es el mismo álgebra con sin⁡2θ e±2iϕ\sin^2\theta\,\eu^{\pm2\iu\phi} y sin⁡θcos⁡θ e±iϕ\sin\theta\cos\theta\,\eu^{\pm\iu\phi}. ∎

Son las partes angulares de los orbitales pp y dd dibujados en el volumen del segundo año: sus nombres pxp_x, dxyd_{xy}, dz2d_{z^2} son las formas cartesianas que acaban de encontrarse.

Los primeros niveles de un rotor rígido lineal, E_J = J(J+1)\, 2/2I, con sus degeneraciones 2J + 1. Las separaciones crecen como 2(J+1): es la razón de las rayas equidistantes del .
Los primeros niveles de un rotor rígido lineal, EJ=J(J+1) ℏ2/2IE_J = J(J+1)\,\hbar^2/2I, con sus degeneraciones 2J+12J + 1. Las separaciones crecen como 2(J+1)2(J+1): es la razón de las rayas equidistantes del Capítulo 6.

Ejemplo 1.27 (Niveles de rotación de HCl\ce{HCl})

Los espectroscopistas escriben EJ=hcB J(J+1)E_J = hcB\,J(J+1) con la constante de rotación B=ℏ/(4πcI)B = \hbar/(4\pi cI) en números de onda. Para HX35X2235Cl\ce{H^{35}Cl}, B=10.593 cm−1B = 10.593\,\mathrm{cm}^{-1}: los primeros niveles están en 0, 21.19, 63.56 y 127.12 cm−1127.12\,\mathrm{cm}^{-1}, todos muy por debajo de la energía térmica a temperatura ambiente, kT/hc=207 cm−1kT/hc = 207\,\mathrm{cm}^{-1}. Hay muchos niveles de rotación poblados; de vibración, solo uno (ω~e=2991 cm−1\tilde\omega_e = 2991\,\mathrm{cm}^{-1}).

1.5 El átomo de hidrógeno y el efecto túnel

El átomo de hidrógeno, resuelto

Un electrón de carga −e-e en torno a un núcleo de carga +e+e tiene V=−e2/4πε0rV = -e^2/ 4\pi\varepsilon_0r. Con la masa reducida μ\mu del electrón y el protón,

H^=−ℏ22μ∇2−e24πε0r,∇2=1r2∂∂r(r2∂∂r)−L^2ℏ2r2.\hat H = -\frac{\hbar^2}{2\mu}\nabla^2 - \frac{e^2}{4\pi\varepsilon_0r},\qquad \nabla^2 = \frac{1}{r^2}\frac{\partial}{\partial r}\Big(r^2\frac{\partial}{ \partial r}\Big) - \frac{\hat L^2}{\hbar^2r^2}.

La parte angular de ∇2\nabla^2 es el operador del rotor: las soluciones se separan en la forma ψ=R(r) Yl,m(θ,ϕ)\psi = R(r)\,Y_{l,m}(\theta,\phi), y RR obedece a la ecuación radial

−ℏ22μ1r2 ⁣d ⁣dr(r2 ⁣dR ⁣dr)+[l(l+1)ℏ22μr2−e24πε0r]R=ER.-\frac{\hbar^2}{2\mu}\frac{1}{r^2}\frac{\dd}{\dd r}\Big(r^2\frac{\dd R}{\dd r} \Big) + \Big[\frac{l(l+1)\hbar^2}{2\mu r^2} - \frac{e^2}{4\pi\varepsilon_0r} \Big]R = ER.

La rotación añade un término centrífugo l(l+1)ℏ2/2μr2l(l+1)\hbar^2/2\mu r^2 que mantiene a los electrones con l>0l > 0 alejados del núcleo.

Teorema 1.28 (Niveles del átomo de hidrógeno)

Los niveles ligados de un átomo hidrogenoide de carga nuclear ZeZe son

En=−μme Z2hcR∞n2,n=1,2,…,hcR∞=mee48ε02h2,E_n = -\frac{\mu}{m_e}\,\frac{Z^2hcR_\infty}{n^2}, \qquad n = 1, 2, \dots, \quad hcR_\infty = \frac{m_ee^4}{8\varepsilon_0^2h^2},

con l=0,…,n−1l = 0, \dots, n - 1; las funciones radiales 1s1s y 2p2p son R10∝e−Zr/aR_{10} \propto \eu^{-Zr/a} y R21∝r e−Zr/2aR_{21} \propto r\,\eu^{-Zr/2a}, con a=(me/μ)a0a = (m_e/\mu)a_0.

Demostración parcial. Se toma Z=1Z = 1 y se ensaya R=e−r/aR = \eu^{-r/a} con l=0l = 0. Entonces R′=−R/aR' = -R/a y 1r2(r2R′)′=R/a2−2R/(ar)\frac{1}{r^2}(r^2R')' = R/a^2 - 2R/(ar). La ecuación radial se convierte en −ℏ22μ(1a2−2ar)−e24πε0r=E-\frac{\hbar^2}{2\mu}\big(\frac{1}{a^2} - \frac{2}{ar}\big) - \frac{e^2}{4\pi \varepsilon_0r} = E para todo rr: los términos en 1/r1/r se cancelan si a=4πε0ℏ2/μe2a = 4\pi\varepsilon_0\hbar^2/\mu e^2, y entonces E=−ℏ2/2μa2=−μe4/8ε02h2E = -\hbar^2/2\mu a^2 = -\mu e^4/8\varepsilon_0^2h^2, el nivel n=1n = 1. Con R=re−r/bR = r\eu^{-r/b} y l=1l = 1, la misma sustitución deja términos en 1/r21/r^2 (que se cancelan con el término centrífugo), en 1/r1/r (que se cancelan si b=2ab = 2a) y una constante, E=−ℏ2/2μb2E = -\hbar^2/2\mu b^2, la cuarta parte de la anterior: el nivel n=2n = 2. El caso general, un polinomio multiplicado por e−r/na\eu^{-r/na} cuya serie debe terminar para seguir siendo normalizable, corresponde a cursos más avanzados. ∎

Funciones radiales del átomo de hidrógeno, obtenidas de las soluciones anteriores (cada una escalada a un máximo próximo a 1). R_20 tiene un nodo radial en r = 2a; R_21 se anula en el núcleo, expulsada por el término centrífugo.
Funciones radiales del átomo de hidrógeno, obtenidas de las soluciones anteriores (cada una escalada a un máximo próximo a 1). R20R_{20} tiene un nodo radial en r=2ar = 2a; R21R_{21} se anula en el núcleo, expulsada por el término centrífugo.

Ejemplo 1.29 (La energía de ionización del hidrógeno)

hcR∞=13.6057 eVhcR_\infty = 13.6057\,\mathrm{eV} y μ/me=1/(1+me/mp)=0.999456\mu/m_e = 1/(1 + m_e/m_p) = 0.999456, luego −E1=13.598 eV-E_1 = 13.598\,\mathrm{eV}: la energía de ionización medida del hidrógeno, hasta la última cifra dada en el volumen del primer año. La masa reducida cambia la cuarta cifra significativa.

Efecto túnel

Definición 1.30 (Efecto túnel, probabilidad de transmisión)

Se llama efecto túnel al paso de una partícula a través de una región en la que su energía es inferior a la energía potencial, prohibido en mecánica clásica. La probabilidad de transmisión TT de una barrera es la fracción de partículas incidentes que se encuentran más allá de ella.

Proposición 1.31 (Barreras delgadas y gruesas)

Dentro de una barrera de altura V0>EV_0 > E, ψ\psi es una combinación de e±κx\eu^{\pm\kappa x} con κ=2m(V0−E)/ℏ\kappa = \sqrt{2m(V_0 - E)}/\hbar. Para una barrera rectangular de anchura aa,

T=[1+sinh⁡2(κa)4ε(1−ε)]−1,ε=EV0,T = \Big[1 + \frac{\sinh^2(\kappa a)}{4\varepsilon(1-\varepsilon)}\Big]^{-1}, \qquad \varepsilon = \frac{E}{V_0},

y para una barrera gruesa (κa≫1\kappa a \gg 1), T≈16ε(1−ε) e−2κaT \approx 16\varepsilon(1-\varepsilon) \,\eu^{-2\kappa a}.

Demostración parcial. En la barrera, la ecuación de Schrödinger se escribe ψ′′=κ2ψ\psi'' = \kappa^2\psi, cuyas soluciones son las dos exponenciales. Empalmar ψ\psi y ψ′\psi' en las dos paredes con las ondas e±ikx\eu^{\pm\iu kx} del exterior da cuatro ecuaciones lineales; su solución, una página de álgebra, es la fórmula de TT, que se admite. Para κa\kappa a grande, sinh⁡κa≈12eκa\sinh\kappa a \approx \frac12\eu^{\kappa a} domina y da la forma de barrera gruesa. ∎

El factor decisivo es e−2κa\eu^{-2\kappa a}, y κ∝m\kappa \propto \sqrt m: el efecto túnel es cosa de las partículas más ligeras, los electrones primero, luego los protones, y mucho menos los deuterones.

Probabilidad de transmisión a través de una barrera modelo de 0.40 eV de altura y 50 pm de anchura. A la mitad de la altura de la barrera, el protón pasa unas 58 veces más a menudo que el deuterón (escala logarítmica).
Probabilidad de transmisión a través de una barrera modelo de 0.40 eV de altura y 50 pm de anchura. A la mitad de la altura de la barrera, el protón pasa unas 58 veces más a menudo que el deuterón (escala logarítmica).

Observación 1.32 (Dónde se manifiesta el efecto túnel)

La molécula de amoniaco se invierte como un paraguas, con su nitrógeno atravesando el plano de los tres hidrógenos, por efecto túnel a través de una barrera baja; las transferencias de protón y de átomo de hidrógeno en las enzimas muestran efectos isotópicos cinéticos mucho mayores de lo que predice el Capítulo 12 sin efecto túnel; y el microscopio de efecto túnel obtiene imágenes de átomos individuales sobre una superficie gracias a la corriente de electrones que atraviesan por efecto túnel un hueco de vacío de unas décimas de nanómetro, una corriente que se divide por diez por cada 0.1 nm adicional.

Historia — Schrödinger, 1926

Schrödinger en 1933.

En cuatro artículos escritos en la primera mitad de 1926, Erwin Schrödinger sustituyó las reglas cuánticas del átomo de Bohr por una ecuación de ondas, la resolvió para el átomo de hidrógeno, el oscilador y el rotor, y recuperó los niveles que la espectroscopia había medido. Los enteros que Bohr había impuesto a mano aparecían ahora por sí solos, con la misma naturalidad que el número de nodos de una cuerda que vibra. Compartió el premio Nobel de Física de 1933 con Paul Dirac.

1.6 Ejercicios

Ejercicio 1.1 ★

Un electrón está confinado en una caja de longitud (a) 1.00 nm1.00\,\mathrm{nm}, (b) 0.50 nm0.50\,\mathrm{nm}. Calcúlese E1E_1 en julios y en electronvoltios, y la longitud de onda de la transición n=1→2n = 1 \to 2.

Solución

Solución de Ejercicio 1.1.

E1=h2/8meL2E_1 = h^2/8m_eL^2. (a) L=1.00 nmL = 1.00\,\mathrm{nm}: E1=6.02×10−20 J=0.376 eVE_1 = 6.02 \times 10^{-20}\,\mathrm{J} = 0.376\,\mathrm{eV}; ΔE12=3E1\Delta E_{12} = 3E_1 y λ=hc/3E1=1099 nm\lambda = hc/3E_1 = 1099\,\mathrm{nm}, en el infrarrojo cercano. (b) L=0.50 nmL = 0.50\,\mathrm{nm}: E1E_1 es cuatro veces mayor, 2.41×10−19 J2.41 \times 10^{-19}\,\mathrm{J} =1.50 eV= 1.50\,\mathrm{eV}, y λ=275 nm\lambda = 275\,\mathrm{nm}, en el ultravioleta.

Ejercicio 1.2 ★

Enumérense los cinco niveles más bajos de una partícula en una caja cúbica, en unidades de h2/8mL2h^2/8mL^2, con sus degeneraciones.

Solución

Solución de Ejercicio 1.2.

nx2+ny2+nz2=3n_x^2 + n_y^2 + n_z^2 = 3 (1,1,1): g=1g = 1; 66 (2,1,1 y permutaciones): g=3g = 3; 99 (2,2,1): g=3g = 3; 1111 (3,1,1): g=3g = 3; 1212 (2,2,2): g=1g = 1.

Ejercicio 1.3 ★

Calcúlense los conmutadores [x^2,p^][\hat x^2,\hat p] y [p^2,x^][\hat p^2,\hat x] haciéndolos actuar sobre una función f(x)f(x).

Solución

Solución de Ejercicio 1.3.

[x^2,p^]f=−iℏx2f′+iℏ(x2f)′=2iℏxf[\hat x^2,\hat p]f = -\iu\hbar x^2f' + \iu\hbar(x^2f)' = 2\iu\hbar xf, luego [x^2,p^]=2iℏx^[\hat x^2,\hat p] = 2\iu\hbar\hat x. [p^2,x^]f=−ℏ2[(xf)′′−xf′′]=−2ℏ2f′[\hat p^2,\hat x]f = -\hbar^2[(xf)'' - xf''] = -2\hbar^2f', luego [p^2,x^]=−2iℏp^[\hat p^2,\hat x] = -2\iu\hbar\hat p.

Ejercicio 1.4 ★

Con B=10.593 cm−1B = 10.593\,\mathrm{cm}^{-1} para HX35X2235Cl\ce{H^{35}Cl}, dense las energías (en cm−1\mathrm{cm}^{-1}) y las degeneraciones de los niveles J=0J = 0 a 3, y el momento de inercia de la molécula.

Solución

Solución de Ejercicio 1.4.

EJ/hc=BJ(J+1)E_J/hc = BJ(J+1): 0, 21.19, 63.56 y 127.12 cm−1127.12\,\mathrm{cm}^{-1}, con g=1g = 1, 3, 5, 7. I=h/(8π2cB)=2.643×10−47 kg m2I = h/(8\pi^2cB) = 2.643 \times 10^{-47}\,\mathrm{kg}\,\mathrm{m}^{2} (con cc en cm/s\mathrm{cm}/\mathrm{s}).

Ejercicio 1.5 ★★

Para el estado ψn\psi_n de una partícula en una caja, calcúlense ⟨x⟩\langle x\rangle y ⟨x2⟩\langle x^2\rangle. Evalúese ⟨x2⟩\langle x^2\rangle para n=1n = 1 y compárese con el valor clásico L2/3L^2/3 de una partícula con igual probabilidad de estar en cualquier punto.

Solución

Solución de Ejercicio 1.5.

∣ψn∣2|\psi_n|^2 es simétrica respecto de L/2L/2, luego ⟨x⟩=L/2\langle x\rangle = L/2. Con sin⁡2u=12(1−cos⁡2u)\sin^2u = \frac12(1 - \cos2u) y dos integraciones por partes, ⟨x2⟩=L2(13−12n2π2)\langle x^2\rangle = L^2\big(\frac13 - \frac{1}{2n^2\pi^2}\big). Para n=1n = 1, ⟨x2⟩=0.283L2\langle x^2\rangle = 0.283L^2, menos que el valor clásico L2/3L^2/3: el estado fundamental se acumula en el centro. Al crecer nn se recupera el valor clásico.

Ejercicio 1.6 ★★

Calcúlense las energías de punto cero armónicas de HX2\ce{H2} y DX2\ce{D2} en kJ/mol\mathrm{kJ}/\mathrm{mol} a partir de ω~e=4401.2 cm−1\tilde\omega_e = 4401.2\,\mathrm{cm}^{-1} y 3115.5 cm−13115.5\,\mathrm{cm}^{-1}, y su diferencia. ¿Qué razón de números de onda cabe esperar de las masas reducidas, y cuánto se acerca a ella la medida?

Solución

Solución de Ejercicio 1.6.

1 cm−1=0.011 963 kJ/mol1\,\mathrm{cm}^{-1} = 0.011\,963\,\mathrm{kJ}/\mathrm{mol}. EPC(HX2\ce{H2}) =2200.6×0.011963=26.33 kJ/mol= 2200.6 \times 0.011963 = 26.33\,\mathrm{kJ}/\mathrm{mol}, EPC(DX2\ce{D2}) =18.63 kJ/mol= 18.63\,\mathrm{kJ}/\mathrm{mol}; diferencia: 7.69 kJ/mol7.69\,\mathrm{kJ}/\mathrm{mol}. Misma constante de fuerza, luego ω~∝μ−1/2\tilde\omega \propto \mu^{-1/2} y la razón esperada es μD/μH=1.4137\sqrt{\mu_D/\mu_H} = 1.4137 (masas atómicas 2.01410 y 1.00783); la razón medida es 1.4127, a un 0.07 % (la pequeña diferencia procede de los electrones, que no siguen a los núcleos perfectamente).

Ejercicio 1.7 ★★

La función ϕ=Nx(L−x)\phi = Nx(L - x) es una aproximación burda del estado fundamental de la caja. Normalícese y calcúlese después su solapamiento ⟨ψ1∣ϕ⟩\langle\psi_1|\phi\rangle con el verdadero estado fundamental.

Solución

Solución de Ejercicio 1.7.

∫0Lx2(L−x)2  ⁣dx=L5/30\int_0^Lx^2(L - x)^2\,\dd x = L^5/30, luego N=30/L5N = \sqrt{30/L^5}. Entonces

⟨ψ1∣ϕ⟩=2L30L5∫0Lx(L−x)sin⁡πxL  ⁣dx=60L3⋅4L3π3=460π3=0.99928.\langle\psi_1|\phi\rangle = \sqrt{\tfrac2L}\sqrt{\tfrac{30}{L^5}}\int_0^Lx(L - x)\sin\tfrac{\pi x}{L}\,\dd x = \frac{\sqrt{60}}{L^3}\cdot\frac{4L^3}{\pi^3} = \frac{4\sqrt{60}}{\pi^3} = 0.99928 .

La parábola es en un 99.86 % estado fundamental (el cuadrado del solapamiento).

Ejercicio 1.8 ★★

Demuéstrese que p^=−iℏ  ⁣d/ ⁣dx\hat p = -\iu\hbar\,\dd/\dd x es hermítico para funciones que se anulan en los extremos de un intervalo, y que x^p^\hat x\hat p no lo es.

Solución

Solución de Ejercicio 1.8.

⟨f∣p^g⟩=−iℏ∫f∗g′  ⁣dx=−iℏ[f∗g]+iℏ∫f∗′g  ⁣dx=∫(−iℏf′)∗g  ⁣dx=⟨p^f∣g⟩\langle f|\hat pg\rangle = -\iu\hbar\int f^*g'\,\dd x = -\iu\hbar[f^*g] + \iu\hbar\int f^{*\prime}g\,\dd x = \int(-\iu\hbar f')^*g\,\dd x = \langle\hat pf|g\rangle, pues el corchete se anula en los extremos. Para x^p^\hat x\hat p, usando que x^\hat x y p^\hat p son ambos hermíticos: ⟨f∣x^p^g⟩=⟨x^f∣p^g⟩=⟨p^x^f∣g⟩\langle f|\hat x\hat pg\rangle = \langle\hat xf|\hat pg\rangle = \langle\hat p\hat xf|g\rangle, y p^x^=x^p^−iℏ≠x^p^\hat p\hat x = \hat x\hat p - \iu\hbar \ne \hat x\hat p: no es hermítico (su forma simetrizada 12(x^p^+p^x^)\frac12(\hat x\hat p + \hat p\hat x) sí lo es).

Ejercicio 1.9 ★★

Modélese el hexa-1,3,5-trieno como seis electrones π\pi en una caja de longitud 6l6l, l=139.7 pml = 139.7\,\mathrm{pm} (cinco enlaces y medio enlace en cada extremo). Predígase la longitud de onda de su primera absorción. ¿Es coloreada la molécula?

Solución

Solución de Ejercicio 1.9.

L=6×139.7=838 pmL = 6 \times 139.7 = 838\,\mathrm{pm}, N=6N = 6: λ=8mecL2/(7h)=331 nm\lambda = 8m_ecL^2/(7h) = 331\,\mathrm{nm}. La absorción está en el ultravioleta: el hexatrieno es incoloro.

Ejercicio 1.10 ★★★

Para la barrera modelo de la figura (V0=0.40 eVV_0 = 0.40\,\mathrm{eV}, a=50 pma = 50\,\mathrm{pm}) en E=V0/2E = V_0/2, calcúlese κ\kappa para un protón y para un deuterón y, con la fórmula de barrera gruesa, la razón TH/TDT_H/T_D. Compárese con la razón exacta, 58.

Solución

Solución de Ejercicio 1.10.

κ=2m(V0−E)/ℏ\kappa = \sqrt{2m(V_0 - E)}/\hbar con V0−E=0.20 eVV_0 - E = 0.20\,\mathrm{eV}: κH=9.82×1010 m−1\kappa_H = 9.82 \times 10^{10}\,\mathrm{m}^{-1}, κD=1.388×1011 m−1\kappa_D = 1.388 \times 10^{11}\,\mathrm{m}^{-1}. Los prefactores se cancelan en la razón: TH/TD=e2a(κD−κH)=e4.06=58T_H/T_D = \eu^{2a(\kappa_D - \kappa_H)} = \eu^{4.06} = 58, el valor exacto: la barrera es gruesa para ambos (κa=4.9\kappa a = 4.9 y 6.9).

Ejercicio 1.11 ★★★

Trátense los seis electrones π\pi del benceno como partículas en un anillo de radio 139.7 pm139.7\,\mathrm{pm} (la distancia C–C es igual al radio de un hexágono regular). Llénense los niveles y calcúlese después la longitud de onda de la transición del nivel lleno más alto al nivel vacío más bajo.

Solución

Solución de Ejercicio 1.11.

Niveles Em=m2ℏ2/2meR2E_m = m^2\hbar^2/2m_eR^2: m=0m = 0 aloja dos electrones, m=±1m = \pm1 cuatro. El nivel lleno más alto es ∣m∣=1|m| = 1, el vacío más bajo ∣m∣=2|m| = 2: ΔE=3ℏ2/2meR2=9.38×10−19 J\Delta E = 3\hbar^2/2m_eR^2 = 9.38 \times 10^{-19}\,\mathrm{J}, λ=212 nm\lambda = 212\,\mathrm{nm}, en el ultravioleta, como la fuerte absorción del benceno.

Ejercicio 1.12 ★★★

Para el estado 1s1s del hidrógeno, ψ=(πa3)−1/2e−r/a\psi = (\pi a^3)^{-1/2}\eu^{-r/a}. Calcúlense los valores esperados ⟨r⟩\langle r\rangle y ⟨V⟩\langle V\rangle, y compruébese que ⟨V⟩=2E1\langle V\rangle = 2E_1. (Úsese ∫0∞rne−βr ⁣dr=n!/βn+1\int_0^\infty r^n\eu^{-\beta r}\dd r = n!/\beta^{n+1}.)

Solución

Solución de Ejercicio 1.12.

⟨r⟩=4ππa3∫0∞r3e−2r/a ⁣dr=4a3⋅3!(2/a)4=32a\langle r\rangle = \frac{4\pi}{\pi a^3}\int_0^\infty r^3\eu^{-2r/a}\dd r = \frac{4}{a^3}\cdot\frac{3!}{(2/a)^4} = \frac32a, 79.4 pm79.4\,\mathrm{pm} para a=a0a = a_0. ⟨1/r⟩=4a3⋅1!(2/a)2=1a\langle 1/r\rangle = \frac{4}{a^3}\cdot\frac{1!}{(2/a)^2} = \frac1a, luego ⟨V⟩=−e2/(4πε0a)\langle V\rangle = -e^2/(4\pi\varepsilon_0a). Como E1=−e2/(8πε0a)E_1 = -e^2/(8\pi\varepsilon_0a) (a partir de a=4πε0ℏ2/μe2a = 4\pi\varepsilon_0\hbar^2/\mu e^2), ⟨V⟩=2E1\langle V\rangle = 2E_1, y la energía cinética media es −E1-E_1.

1.7 Problema: ¿cuánto mide un colorante?

Problema 1.1

Problema de fin de semana — el modelo del electrón libre de tres colorantes de cianina: por qué falla una cadena desnuda, qué transiciones están permitidas, las escalas de energía de una molécula y la longitud de caja que se ajusta

Tres colorantes de cianina simétricos absorben a 524, 603 y 711 nm en etanol; sus cadenas conjugadas, entre los dos nitrógenos y contándolos, tienen j=5j = 5, 7 y 9 átomos. Se toma cada enlace igual a l=139.7 pml = 139.7\,\mathrm{pm}, me=9.109×10−31 kgm_e = 9.109 \times 10^{-31}\,\mathrm{kg}, h=6.626×10−34 J sh = 6.626 \times 10^{-34}\,\mathrm{J}\,\mathrm{s}, c=2.998×108 m/sc = 2.998 \times 10^{8}\,\mathrm{m}/\mathrm{s}, k=1.381×10−23 J/Kk = 1.381 \times 10^{-23}\,\mathrm{J}/\mathrm{K}, 1 eV=1.602×10−19 J1\,\mathrm{eV} = 1.602 \times 10^{-19}\,\mathrm{J}.

Parte I — Una cadena desnuda.

  1. Explíquese por qué una cadena de jj átomos lleva N=j+1N = j + 1 electrones π\pi, y dese NN para cada colorante.
  2. ¿Qué niveles de la caja son el lleno más alto y el vacío más bajo para el primer colorante?
  3. Demuéstrese que la primera absorción está en λ=8mcL2/h(N+1)\lambda = 8mcL^2/h(N+1).
  4. Dese la longitud desnuda L0=(j−1)lL_0 = (j - 1)l de cada cadena.
  5. Calcúlese λ\lambda para cada colorante con L=L0L = L_0.
  6. Compárese con los valores medidos. ¿En qué sentido se equivoca el modelo y qué indica eso sobre LL?

Parte II — Qué transiciones están permitidas. La intensidad de una transición n→n′n \to n' es proporcional a ∣⟨ψn∣x−L/2∣ψn′⟩∣2|\langle \psi_n|x - L/2|\psi_{n'}\rangle|^2.

  1. Demuéstrese que ψn\psi_n es simétrica respecto del centro de la caja para nn impar y antisimétrica para nn par.
  2. Dedúzcase que la integral se anula cuando n+n′n + n' es par.
  3. ¿Está permitida la transición HOMO →\to LUMO de cada colorante?
  4. Para el segundo colorante, ¿está permitida la transición del HOMO al nivel situado encima del LUMO? ¿Y la del nivel situado debajo del HOMO al LUMO?
  5. En el modelo, ¿a qué longitud de onda (como fracción de la primera) estaría la siguiente transición permitida desde el HOMO?
  6. ¿Por qué una regla así, deducida solo de la simetría, sobrevive a lo burdo del modelo?

Parte III — Las escalas de energía de una molécula.

  1. Conviértase la absorción del segundo colorante en una energía en eV.
  2. La vibración de HCl\ce{HCl} tiene ω~e=2991 cm−1\tilde\omega_e = 2991\,\mathrm{cm}^{-1}: conviértase su cuanto en eV.
  3. La rotación de HCl\ce{HCl} tiene B=10.59 cm−1B = 10.59\,\mathrm{cm}^{-1}: conviértase la separación entre J=0J = 0 y J=1J = 1 en meV.
  4. Calcúlese kTkT en meV a 298 K298\,\mathrm{K}.
  5. ¿Qué tipos de excitación están poblados a temperatura ambiente?
  6. ¿En qué regiones del espectro se observan las transiciones electrónicas, vibracionales y rotacionales?

Parte IV — La caja que se ajusta.

  1. A partir de los 603 nm medidos, calcúlese la longitud LL que debe tener la caja del segundo colorante.
  2. Escribiendo L=L0+2δL = L_0 + 2\delta, dedúzcase δ\delta.
  3. Con este δ\delta, predígase λ\lambda para el primer colorante; dese el error relativo.
  4. La misma pregunta para el tercer colorante.
  5. Interprétese δ\delta: ¿adónde van los electrones π\pi más allá de los nitrógenos?
  6. Exprésese δ\delta en longitudes de enlace y enúnciese el resultado del problema: la prolongación efectiva de la caja en cada extremo.
Solución

Solución de Problema 1.1.

1. Cada uno de los jj átomos aporta un orbital p; los j−1j - 1 carbonos y un nitrógeno aportan un electrón cada uno y el otro nitrógeno su par no enlazante, ya que la carga del catión está repartida por la cadena: N=j+1N = j + 1, es decir, 6, 8 y 10. 2. Seis electrones llenan n=1n = 1, 2, 3: HOMO n=3n = 3, LUMO n=4n = 4. 3. ΔE=EN/2+1−EN/2=(N+1)h2/8mL2\Delta E = E_{N/2+1} - E_{N/2} = (N+1)h^2/8mL^2 y λ=hc/ΔE=8mcL2/h(N+1)\lambda = hc/\Delta E = 8mcL^2/h(N+1). 4. L0=4lL_0 = 4l, 6l6l, 8l8l: 559, 838 y 1118 pm1118\,\mathrm{pm}. 5. 147, 257 y 374 nm374\,\mathrm{nm}. 6. Todas demasiado cortas (por factores de 3.6, 2.3 y 1.9): la caja del modelo es demasiado pequeña, ya que λ∝L2\lambda \propto L^2; los electrones se mueven sobre una distancia mayor que la cadena comprendida entre los nitrógenos. 7. Con el origen en el centro, ψn∝cos⁡(nπu/L)\psi_n \propto \cos(n\pi u/L) para nn impar y sin⁡(nπu/L)\sin(n\pi u/L) para nn par (u=x−L/2u = x - L/2): funciones par e impar de uu. 8. Si n+n′n + n' es par, ψn\psi_n y ψn′\psi_{n'} tienen la misma paridad; su producto multiplicado por uu es impar y su integral sobre un intervalo simétrico es nula. 9. El HOMO N/2N/2 y el LUMO N/2+1N/2 + 1 tienen n+n′=N+1n + n' = N + 1, impar: siempre permitida. 10. Segundo colorante: HOMO 4, LUMO 5. 4→64 \to 6: suma par, prohibida. 3→53 \to 5: suma par, prohibida. 11. La siguiente permitida desde n=4n = 4 es 4→74 \to 7, con un ΔE\Delta E mayor en un factor (49−16)/(25−16)=3.67(49 - 16)/(25 - 16) = 3.67: a λ/3.67\lambda/3.67, 164 nm164\,\mathrm{nm} para el segundo colorante, en pleno ultravioleta. 12. Solo depende de la simetría del potencial respecto del centro de la cadena, que comparten los colorantes simétricos reales. 13. E=hc/λ=3.294×10−19 J=2.06 eVE = hc/\lambda = 3.294 \times 10^{-19}\,\mathrm{J} = 2.06\,\mathrm{eV}. 14. hcω~e=0.371 eVhc\tilde\omega_e = 0.371\,\mathrm{eV}. 15. 2hcB=4.21×10−22 J=2.63 meV2hcB = 4.21 \times 10^{-22}\,\mathrm{J} = 2.63\,\mathrm{meV}. 16. kT=25.7 meVkT = 25.7\,\mathrm{meV}. 17. Solo los niveles de rotación (2hcB≪kT2hcB \ll kT); las excitaciones vibracionales (15kT15kT) y electrónicas (80kT80kT) prácticamente no están pobladas. 18. Electrónicas: visible y ultravioleta; vibracionales: infrarrojo; rotacionales: microondas (e infrarrojo lejano). 19. L=λh(N+1)/8mcL = \sqrt{\lambda h(N+1)/8mc} con N=8N = 8: L=1283 pmL = 1283\,\mathrm{pm}. 20. δ=(1283−838)/2=222 pm\delta = (1283 - 838)/2 = 222\,\mathrm{pm}. 21. L=559+445=1003 pmL = 559 + 445 = 1003\,\mathrm{pm}, λ=474 nm\lambda = 474\,\mathrm{nm}, un 9.5 % por debajo de 524 nm. 22. L=1118+445=1562 pmL = 1118 + 445 = 1562\,\mathrm{pm}, λ=732 nm\lambda = 732\,\mathrm{nm}, un 2.9 % por encima de 711 nm. 23. Hacia los anillos aromáticos que llevan los nitrógenos: el sistema π\pi no termina en los átomos de nitrógeno, y la caja efectiva incluye parte de cada anillo. 24. δ/l=222/139.7=1.6\delta/l = 222/139.7 = 1.6: la caja se prolonga unos 222 pm222\,\mathrm{pm}, enlace y medio, más allá de cada nitrógeno, y con esa sola longitud el modelo reproduce los tres colores con un error inferior al 10 %.

Términos definidos en este capítulo

Ver los 852 términos del glosario