Estructura Atómica de La Materia
Estructura Atómica de La Materia
Estructura Atómica de La Materia
LA ESTRUCTURA ATÓMICA
La estructura atómica
1. Espectros atómicos. Series del Hidrógeno.
1.1. Espectros de emisión y absorción
2. Revisión del modelo de Rutherford.
2.1. Evolución de los modelos atómicos
2.2. Limitaciones del modelo de Rutherford
3. Teoría cuántica de Planck
4. Efecto fotoeléctrico
5. Modelo de Bohr.
5.1. Determinación de la energía y del radio
5.2.Determinación de la constante de Rydberg
6. Inconvenientes del modelo de Bohr.
6.1.Modelo de Bohr-Sommerfeld
6.2. Origen de los números cuánticos l,m y s
6.3.Principio de Exclusión de Pauli
7. Modelo mecanocuántico del átomo.
7.1. Teoría Corpúsculo ondulatoria de De Broglie
7.2. Principio de Incertidumbre de Heisenberg
7.3. Orbitales atómicos.
7.4.Los números cuánticos en el modelo mecanocuántico
10
11
1. Espectros atómicos. Series del Hidrógeno.
Espectros de emisión
Por espectro de emisión se conoce a aquellos espectros obtenidos a partir
de la descomposición de la radiación emitida por una sustancia. Estos
espectros dependen de las condiciones a las que se haya sometido la
sustancia.
Espectros de absorción
Si sobre una sustancia se hace incidir radiación
compuesta por un gran número de longitudes de
onda y posteriormente se registra el espectro de
dicha radiación una vez atravesada la sustancia
observaremos que aparecen unas rayas o bandas
negras correspondientes a las longitudes de onda
absorbidas por la sustancia.
Este espectro recibe el nombre de espectro de
absorción y se puede comprobar que para una
misma sustancia las rayas o bandas negras que
aparecen en este espectro corresponden con las
rayas o bandas que aparecen en su espectro de
emisión (fig 1.2).
Fig.1.2. Espectro s
11
12
1.2. Espectro atómico del hidrógeno
Ya Ångström(1814-1874) había descubierto en 1853 en el espectro
solar, una raya característica del hidrógeno que llamó " (alfa). Después en
el estudio realizado sobre la radiación emitida por una lámpara de descarga
del hidrógeno se observó la presencia de tres rayas más con unas
determinadas longitudes de onda 8 (lambda) que se podían obtener con el
esquema del montaje (fig.1.3):
1 = ⎛⎜ 4 ⎞⎟ ⎛⎜ 1 − 1 ⎞⎟ = Ry ⎛⎜ 1 − 1 ⎞⎟
λ ⎝ b ⎠ ⎝ 22 n 2 ⎠ ⎝2
2
n2 ⎠
12
13
Ejemplo 1.
1.Calcular la frecuencia emitida cuando el electrón excitado del hidrógeno cae del nivel 4 al nivel 2. ¿A qué serie
correspondería en el espectro del hidrógeno? ¿Qué raya produciría?.
7 -1 8 -1
DATOS: Ry=1,09.10 m c=3.10 ms
Se emplea la fórmula de Balmer-Rydberg, 1/8 = Ry(1/nf2 – 1/ni2), dado que en el enunciado se precisa el nivel de salida y el
de llegada, que indica que se trata de la serie de Balmer (saltos al nivel 2) y a la segunda raya (de 4 a 2).
1/8 = 1,09.107m-1 (1/22 – 1/42) = 2,04.106 m-1
Determinación de frecuencia: < = c/8 < = c(1/8) = (3.108ms-1)(2,04.106 m-1) = 6,13.1014 s-1.
13
14
.
- Casi toda la masa del átomo y toda su carga positiva está concentrada en
un determinado espacio, el núcleo responsable del desvío de las partículas
α
QQ '
Fig.1.6.Detalle de la dispersión Fe = k
r2
Como las cargas son numéricamente iguales aunque de diferente signo e
iguales a la del electrón e, tendremos que:
e2 mv 2
Fe = k = ma N =
r2 r
14
15
3. Teoría cuántica de Planck
Se va a hacer un resumen y un recuerdo de lo explicado en este punto
en el Heurema I
Hechos previos:
En 1859 Kirchhoff, define la superficie absorbente ideal de
radiación, que a su vez se convertiría en emisor ideal. Lo llama cuerpo
negro. Podría realizarse en una esfera hueca recubierta por dentro de negro
de humo y con un orificio por que penetra o sale radiación que será la
radiación del cuerpo negro(fig.1.7.).
La teoría:
En diciembre de 1900, Planck (1858-1947), que había sido
alumno de Kirchhoff, presenta en la Academia de Ciencias de Berlín, la 14
memoria sobre la radiación del cuerpo negro. Justificaba las anomalías con
un pequeño postulado que llamó “hipótesis de los cuantos”. Fig 1.8. Catástrofe ultravioleta
Por lo tanto la energía total será un múltiplo entero del cuanto energético
(o cantidad energética más pequeña). El valor que calculó Planck para h
fue de 6,55,10-34 Js. Actualmente se toma como 6,626.10-34Js
15
16
4. Efecto fotoeléctrico
Hechos previos:
En 1872, el profesor Stoliétov de la universidad de Moscú descubre que
entre dos láminas metálicas separadas por vacío y conectadas a una pila
se producía una corriente cuando eran iluminadas por una lámpara de
vapor de mercurio, corriente eléctrica que se interrumpía al apagarse la luz.
13 años después, Hertz y Hallwachs llaman a este hecho: efecto
fotoeléctrico, esto es electricidad producida por la luz.
Contradicciones:
El efecto fotoeléctrico no dependía de la intensidad de la iluminación como
parecería lógico sino del “color” de la luz que iluminaba la placa metálica, o
sea de la frecuencia luminosa. Separando la placa hasta una distancia
prudencial la producción de corriente se mantenía, y ocurría
instantáneamente. Sólo para frecuencias mínimas específicas, dependía de
la intensidad.
Una vez que se descubre el electrón, y que se explica dicho efecto como
debido a la acción de las ondas luminosas de una determinada frecuencia
sobre los electrones del metal, ya a principios de siglo XX, Lenard (1862-
1947), calcula que la energía que necesitaban las ondas electromagnéticas
para ionizar los electrones de una película de sodio, implicaba una
En todo el trabajo
iluminación continuada de una bujía durante 11 días, sin embargo el efecto
Einstein refiriéndose a
los corpúsculos lumi- fotoeléctrico se producía instantáneamente.
nosos se expresó en
términos de “cuantos Explicación de Einstein:
de luz”, que 20 años La radiación monocromática se comporta como un medio discreto formado
después serían bauti- por cuantos de energía, hν (toma la teoría de Planck, en la que nadie
zados como fotones creía). Así, cuando se ilumina una lámina metálica los electrones salen con
por el químico ameri- una determinada velocidad, y por lo tanto tienen una energía cinética (½
cano Lewis.
mv2) que debe ser igual a la energía del cuanto energético, hν menos la
energía de extracción del electrón. La ecuación se formula así:
1
2 mv 2 = hν − Φ 0
La energía o trabajo de extracción que efectúa la radiación, también llamada
“función trabajo” será en virtud de la misma teoría:
Φ 0 = hν 0
siendo ν0, lo que se conoce como frecuencia umbral
Como experimentalmente se determinó la energía de extracción y la energía
cinética del electrón, se pudo comprobar así la teoría de Planck.
La explicación de Einstein, justifica plenamente los hechos experimentales,
la luz roja por muy cerca que se pusiera de la placa metálica, no conseguía
arrancarle electrones, ello era debido a que su frecuencia era muy pequeña,
y por lo tanto la energía cuántica era inferior al trabajo de extracción.
Ejemplo de Aplicación 2:
1.La energía mínima requerida para arrancar el electrón externo del sodio es de 2,3eV.Cuando se ilumina la superficie del
-34 8 -1
metal con luz amarilla de 5890Å ¿se producirá emisión fotoelectrónica?. h=6,63.10 Js. c=3.10 ms .
La energía del cuanto de luz amarilla será E=hc/λ =6,63.10-34 Js 3.108 ms-1/5,89.10-7m=3,40.10-19J.
-19 -19
Como 1eV=1,6.10 J, el trabajo de extracción será 3,68.10 J, inferior a la energía del cuanto por lo que no se producirá
emisión fotoelectrónica.
16
17
5.Modelo de Bohr.
En 1918, justificará Bohr la cuantización
La teoría de Bohr, para el átomo de hidrógeno se suele explicar partiendo del momento angular del electrón del
de 3 principios fundamentales denominados POSTULADOS DE BOHR, que hidrógeno partiendo del principio de
son los siguientes: correspondencia, nombre que se origina a
través de la correspondencia de una
1º. El átomo de hidrógeno está constituido por un núcleo, de carga positiva y órbita clásica no estacionaria a una
masa 1836 veces la del electrón que gira alrededor de aquél en una distancia grande del núcleo, hasta otra
órbita circular de radio r no radiante, denominada estacionaria. Este órbita no radiante en la que el electrón
modelo se puede extender a los átomos monoelectrónicos, sustituyendo estaría en un estado estacionario.
la carga del núcleo por Ze, siendo Z su número de protones o número
El trabajo mínimo para pasar al
atómico. electrón de una órbita a otra, debe
corresponder aproximadamente a la
2º. La condición de cuantización o fijación de la energía en cada órbita para energía cuantizada por Planck= hν.
que sea un estado estacionario es que el momento angular del electrón Teniéndolo en cuenta, calculará
sea un número entero de veces el cuanto de acción de Planck ( h/2π): las condiciones del camino mínimo, y
⎛ h ⎞ considerará que el radio de una órbita
mv.r = n⎜ ⎟ deberá ser un múltiplo entero de veces n,
⎝ 2π ⎠ la distancia mínima ∆r. Así:
siendo n = 1, 2, 3 … según el estado estacionario a que se haga
referencia. Para el H en su estado de mínima energía o fundamental
sería n=1. h es la constante de Planck = 6,63.10-34 J.s
e2
Dado que k = mv 2 ,
r Demostrando la expresión que se atribuye
Si en el segundo postulado se despeja v, se se eleva al cuadrado y se a su segundo postulado
sustituye se obtendrá el valor de r dependiendo de la única variable n que
se denominará número cuántico principal o de Bohr.
2 2
⎛ h ⎞⎛ 1 ⎞ ⎛ h ⎞ ⎛ 1 ⎞
v = n⎜ ⎟ ⎜ ⎟ v2 = n2 ⎜ ⎟ ⎜ ⎟ La clave de la teoría de
⎝ 2π ⎠ ⎝ mr ⎠ ⎝ 2π ⎠ ⎝ mr ⎠ Bohr, fue el considerar la
frecuencia de radiación
similar a la de oscilación
de los electrones alrede-
ke 2 ⎡ 2 ⎛ h ⎞2 ⎛ 1 ⎞2 ⎤ ⎡⎛ h ⎞ 2 ⎛ 1 ⎞⎤ dor del núcleo
= m ⎢n ⎜ ⎟ ⎜ ⎟ ⎥ r = n ⎢⎜ ⎟ ⎜
2
2 ⎟⎥
r ⎣⎢ ⎝ 2π ⎠ ⎝ mr ⎠ ⎦⎥ ⎣⎢⎝ 2π ⎠ ⎝ kme ⎠⎥⎦
17
18
Para dicho cálculo se emplearán las fórmulas y sistemas actuales y no los
-11
empleados por Bohr. Para n=1, r=52,9.10 m = 0,529D que es lo que se
toma como radio de Bohr y constituye un parámetro de referencia en la
teoría atómica, con el nombre bohr que constituye una longitud de 52,9pm.
1 = ⎛⎜ 2,18.10 ⎞⎛ 1 1 ⎞
−18
⎟⎟ ⎜⎜ 2 − 2 ⎟⎟ (1.11)
λ ⎜⎝ hc ⎠ ⎝ nf ni ⎠
Al sustituir los parámetros constantes quedaría:
1 = 1,09557.10 7 ⎛⎜ 1 − 1 ⎞⎟ m −1
λ ⎜ nf 2 ni 2 ⎟
⎝ ⎠
que coincide con la fórmula de Balmer dada por Rydberg, de forma que su
constante sería Ry=109557cm-1, con una pequeña desviación respecto a la
dada por Rydberg, Ry=109677cm-1 en 1900, de forma meramente empírica.
Fig.1.9.Saltos entre niveles La energía en los diferentes saltos cuantizados, variaría según la
gráfica adjunta, observándose el por qué los saltos al estado estacionario
con n=1, nivel 1, denominado K, liberan un cuanto con un contenido
energético mayor, por eso se emite en el ultravioleta (UV), sin embargo,
en los saltos electrónicos a partir del nivel 3, o M, la emisión energética es
menor, produciéndose en el infrarrojo(IR), y la única visible (por eso se
18
19
descubrió antes), fue la serie de Balmer, con saltos al nivel 2 o L.
Véase las fig 1.9 y 1.10
2. Determinar la velocidad que le correspondería al electrón en su órbita normal, según la teoría de Bohr. ¿Y si estuviera
excitado, en el nivel n=3?
DATOS: me =9,1.10-31 kg. Energía =- 2,18.10-18 /n2 J
2 2 2
La energía total correspondiente al electrón es –e /2r según se ha visto en el desarrollo teórico. Como se conoce que e /r
2
= mv /r,
½mv2 = e2/2r = 2,18.10-18 /n2 J , de lo que:
v= √ (2. 2,18.10-18/. 9,1.10-31 kg) / n = 2,19. 106 /n ms-1
Desarrollo: aplicando n=1 (estado fundamental) o n=3 (excitado)
v1 = 2,19. 106 /1 ms-1 = 2,19. 106 ms-1 v3 = 2,19. 106 /3 ms-1= 0,73. 106 ms-1
19
20
6.1.Modelo de Bohr-Sommerfeld
20
21
k l nombre Nº de e
1 0 s (sharp =neta, definida) 2
2 1 p (principal = principal) 6 Tabla 2
3 2 d (diffuse = difusa) 10
4 3 f (fundamental = fundamental) 14
Número cuántico s
En 1926, Uhlenbeck y Goudsmit, encuentran el cuarto número cuántico
sugerido por Pauli, debido al giro del electrón sobre sí mismo y
denominado por eso de espín, con valores de +½, - ½ del cuanto de
acción. Esta cuantización del momento angular intrínseco justificaba la
aparición en los espectros de grupos de dos rayas muy juntas (dobletes)
aún en ausencia de campo magnético externo.
Si n=3, l=0,1 y 2. Para l=0, m=0 y s= ½ y –½ .O sea que en esta órbita podría haber dos electrones con sus cuatro números
cuánticos distintos. Por lo mismo para l=1, m=+1,0 y –1 y como s= ½ y –½, implica 6 electrones y para l=2, m=+2,+1,0,-1,-2 y
como s= ½ y –½, implica 10 electrones. En total 18.
21
22
22
23
que había descubierto el electrón. De estos experimentos se dedujo que
el electrón que era un corpúsculo, se comportaba también como una onda,
puesto que daba lugar a imágenes de difracción.
Ejemplo 5
introducirá un error en la medida del fenómeno, y por lo tanto una sin que se modifique el orden de
incertidumbre. magnitud, que es lo importante.
Si se quisiera observar la posición de un electrón, necesitaríamos Estas relaciones se aplican
un dispositivo óptico que lo iluminara con una longitud de onda del tamaño actualmente en la resolución de
del observable, como la longitud de onda es inversamente proporcional a la problemas.
frecuencia (fórmula A1.3), y ésta directamente proporcional a la energía (ley
de Planck, fórmula A1.2). Este electrón recibiría tal cantidad de energía que
produciría una imprecisión en la determinación de su velocidad.
23
EL MODELO ATÓMICO DE BOHR EN SU PUNTO JUSTO
La mayoría de los textos de Química enfocan el modelo del átomo de Bohr, sin que sus autores vayan a la fuente original, esto es
los trabajos que Bohr publicó en junio, septiembre y noviembre de 1913, sobre la constitución de átomos y moléculas, en el Philo-
sophical Magazine. Si los hubieran leído, estoy plenamente seguro que no pondrían en boca de Bohr, cosas que él jamás estable-
ció. Lo que ocurre que es cada autor toma como base lo escrito por otro y así sucesivamente sin ir a las fuentes.
Vamos a intentar reducir y comentar el artículo que Niels Bohr envió a la revista citada el 5 de abril de 1913, base de toda su teo-
ría.Consta de cinco partes. Una introducción. Consideraciones generales: Unión de electrones por núcleos positivos. Emisión de
espectros de rayas. Continuación de consideraciones generales y Estado permanente de un sistema atómico.
En el enfoque que hace parte del modelo de Rutherford, su maestro en Mánchester, al que anticipó el borrador del trabajo el año
anterior.
Sitúa al electrón de carga –e, y masa m, describiendo órbitas elípticas cuyo eje mayor es 2a con una frecuencia de revolución T,
alrededor de un núcleo con carga E, y con masa mucho mayor que la del electrón, cuya velocidad a su vez es muy pequeña compa-
rada con la de la luz. Supone que en principio no hay radiación de energía y que el electrón describirá órbitas elípticas estaciona-
rias. La frecuencia de revolución T, y el eje mayor de la órbita 2a dependerán de la cantidad de energía W que debe suministrarse
al sistema para alejar al electrón una distancia infinitamente grande del núcleo.
Al aplicar el equilibrio de fuerzas llega a la conclusión (tomada de trabajos de Rutherford y colaboradores) que:
3
⎛ 2⎞ W 2 eE
ω = ⎜⎜ ⎟⎟ y 2a =
⎝ π ⎠ eE m W
4W 2 2W 3
=
τ 2 h 2 π 2 m( eE ) 2
⎛ 2π 2 me 2 E 2 ⎞⎛ 1 ⎞
simplificando y despejando W = ⎜⎜ ⎟⎟⎜ 2 ⎟
⎝ h2 ⎠⎝ τ ⎠
eE ⎛ h2 ⎞ 2
y como 2 a = , llevando a esta expresión el valor obtenido de W, quedaría: 2a = ⎜⎜ 2 ⎟⎟τ
W ⎝ 2π mEe ⎠
haciendo en las expresiones anteriores e=E = 4,7.10-10 ; e/m = 5,31.1017 y h = 6,5.10-27, se obtiene:
2 a = 1,1. 10-8 cm; T= 6,2.1015 ( 1/s ) y W/e =13 volt2
Valores que son del orden de magnitud de las dimensiones lineales de los átomos, las frecuencias ópticas y los potenciales de
ionización.
La parte segunda del artículo hace referencia al espectro del hidrógeno, y la explicación de las series de Bal-
mer y Paschen.
⎛ 2π 2 me 4 ⎞⎛ 1 ⎞
Wτ = ⎜⎜ ⎟⎟⎜ 2 ⎟ y la cantidad de energía emitida por el paso del sistema de un estado correspondiente a J =J1 , a otro
⎠⎝ τ ⎠
2
⎝ h
⎛ 2π 2 me 4 ⎞⎛ 1 1 ⎞
correspondiente a J2 , será : Wτ 2 − Wτ 1 = ⎜⎜ ⎟⎟⎜⎜ 2 − 2 ⎟⎟
⎠⎝ τ 2 τ 1 ⎠
2
⎝ h
Si suponemos que la radiación es homogénea y que la cantidad de energía emitida es h<, siendo < la fre-
cuencia de la radiación, obtendremos :
Wτ 2 − Wτ 1 = hν
⎛ 2π 2 me 4 ⎞⎛ 1 1 ⎞
de lo que ν = ⎜⎜ ⎟⎟⎜⎜ 2 − 2 ⎟⎟ expresión que constituye la ley que relaciona las rayas del espectro del hidrógeno.
⎠⎝ τ 2 τ 1 ⎠
3
⎝ h
Si hacemos J2 = 2, obtenemos la serie de Balmer. Si J2 = 3, obtendremos la serie observada por Paschen en la región infrarroja.
Escribiendo J2 = 1 y J2 = 4, 5,. obtendremos series situadas en el extremo ultravioleta e infrarrojo , no observadas pero cuya
existencia debe esperarse.
La concordancia no es solo cualitativa sino también cuantitativa.
Haciendo e=4,7.10-10 ; e/m=5,31.1017 y h=6,5.10-27, obtendremos ⎛⎜ 2π me ⎞⎟ = 3,1.1015
2 4
⎜ h3 ⎟
⎝ ⎠
El valor experimental observado para este factor es 3,290.1015”
La tercera parte del artículo, hace referencia a las consideraciones generales y enuncia de otra forma sus
ideas
“....Representando por M el momento angular del electrón en torno al núcleo, y para una órbita circular,
T
π M = siendo T la frecuencia de revolución y T la energía cinética del electrón. Además T = W, como
ω
ω
ω , por lo que se obtiene que τh h h
W =τ h π M = 2 =τ . Si hacemos M 0 = , se tendrá que M = τ M 0
2 ω 2 2π
Por consiguiente si postulamos que la órbita del electrón es circular en los estados estacionarios, se puede expresar que el mo-
mento angular del electrón en torno al núcleo en un estado estacionario del sistema, es un múltiplo entero de un valor universal
independiente de la carga del núcleo. La posible importancia del momento angular en la discusión de los sistemas atómicos en
relación con la teoría de Planck, es destacada por Nicholson”.
1
El primero al que se le ocurrió aplicar la teoría de Planck al átomo, fue a Arthur Eric Haas, estudiante de doctorado de la universidad de Viena,
en 1910, refiriéndolo al modelo de Thomson y generando la expresión hL= e2/a. Lo llevó a la conferencia Solvay de 1911, donde lo leyó Nichol-
son, astrónomo de Oxford que andaba buscando una forma de explicar los espectros de la luz procedentes de la corona solar, donde había identi-
ficado a dos nuevos elementos el protofluor y el nebulio. Nicholson cuantificó en 1912, el momento angular del electrón L=n(h/2B). Bohr se
inspiró en los trabajos de Nicholson, como así lo reconoce en este artículo.
2
Thomson en 1911, había calculado experimentalmente el potencial de ionización del hidrógeno, dando un valor de 11 voltios.
3
El trabajo de Bohr consta de 36 páginas, y en el resumen que se reseña se ha destacado lo más importante, referido a lo que se enseña sobre la
teoría de Bohr en la Química de Enseñanza Secundaria y Bachillerato, en ninguno de los cuales se hace referencia a la hipótesis fundamental o
sea que la frecuencia de radiación es la mitad de la frecuencia de revolución.
2
EL ÁTOMO MECANO CUÁNTICO
Obtención de la ecuación de Schrödinger Ecuación de Schrödinger y Orbitales atómicos
Se parte de la ecuación de una onda
ψ= Asen( ωt –kx) La incertidumbre introducida por Heisenberg, hace que se
en la cual Ψ es la elongación, A la amplitud, ω busque un modelo matemático que determine la probabilidad de
frecuencia angular, k el número de onda, y x la encontrar al electrón en una zona determinada del átomo.
posición respecto al centro de la perturbación.
Derivándola dos veces respeto a la posición: El tratamiento matemático considerando su condición ondulatoria
lo hace en 1926, Schrödinger (1887-1961) que convierte la ecuación de
d2ψ/dx2 = - k2Asen(ωt-kx)= - k2ψ una onda, en la ecuación general de ondas para el electrón:
3
4
orbital n l m
En la figura 2, se observa la génesis de lo que se dibuja como orbital atómico 3s 3 0 0
2px, obtenida a partir de las soluciones de la ecuación de Schrödinger para los 3px 3 1 1
valores indicados en coordenadas rectangulares (visión tridimensional), y elevándolas 3py 3 1 -1
al cuadrado(densidad de probabilidad) y proyectando las curvas de igual densidad de 3pz 3 1 0
probabilidad sobre en plano XY en unidades de radio de Bohr o bohrs y por fin lo que 3dx2-y2 3 2 2
se dibuja como OA 2px, es la línea de contorno que abarca el 90% de la probabilidad 3dxy 3 2 -2
acumulativa. 3dxz 3 2 1
3dyz 3 2 -1
3dz2 3 2 0
Para l=2,
Existirán dos planos nodales perpendiculares entre sí, y el orbital forma 4 lóbulos (orbital tipo d), parecido a un trébol
de 4 hojas y así sucesivamente, formándose los diferentes orbitales atómicos que se dan para n=3 ( tabla 1).
De esta forma la distribución espacial de los diferentes OA, para n=3, en sus diferentes valores de , l y m, dan
lugar a los OA cuyos contornos deberás dibujar
En el caso de orbitales d (l=2) habría 5 valores de m, y a cada uno correspondería una distribución
espacial de probabilidad. Obsérvense los mapas tridimensionales que hacen referencia al 3dz2 y al 3dxy y las
correspondientes superficies de contorno de densidades de probabilidad y la envolvente del 90%, que es lo
que se dibuja como OA (fig.3 y 4).
4
5
PROYECCIÓN DE LA DENSIDAD
DE PROBABILIDAD PROYECCIÓN SOBRE EL PLANO XY
SOBRE EL PLANO XZ DE LA DENSIDAD DE PROBABILIDAD
CURVAS DE NIVEL DE PROBABILIDAD CONSTANTE
PLANOS
EJE Z NODALES
EJE X
EJE X EJE Y
Los OA d dados en las figuras son los dos tipos perfectamente diferenciados ya que el 3dxz, es igual que el
3dxy aunque orientado sobre otros ejes, y los 3dyz y 3dx2-y2, también son similares aunque están orientados sobre
los ejes de coordenadas xy y no entre ellos.
FUNCIÓN DE ONDA
Representando funciones ψ
2px Si se toman las funciones R y no sus cuadrados,
10.00 aparecerán zonas con valores negativos (por debajo
EJE Y
del plano), y la forma no es tan estilizada como las
5.00 distribuciones de probabilidad.
0.00 EJE X
En los p aparecerá un lóbulo con el signo positivo y
-5.00 otro con el negativo como se puede apreciar en la
ZONA NEGATIVA ZONA POSTIVA
GRIS MAS OSCURO GRIS MAS CLARO figura 5
VALORES MAS NEGATIVOS VALORES LAS POSTIVOS
-10.00
-20.00 -15.00 -10.00 -5.00 0.00 5.00 10.00 15.00 20.00
5
Fig.5.. ψ 2px
6
En los d, surgirán dos lóbulos con signo -, como se puede 20.00
-5.00
-10.00
-15.00
7
Orbitales s con ayuda de la hoja de cálculo EXCEL
Las ecuaciones que se suelen encontrar sobre los orbitales están dadas en coordenadas
polares
3 Zr Zr
1 ⎛Z ⎞ 2 − ao ⎛
1 Zr ⎞ − 2a o
Ψ1s = ⎜ ⎟⎟ e Ψ2s = ⎜ 2 − ⎟⎟e
π ⎜⎝ a o ⎠ 32π ⎜⎝ ao ⎠
3 − Zr
⎛Z
1 ⎞2 ⎛ 18Zr 2Z 2 r 2 ⎞ 3a o
Ψ3s = ⎜ ⎟⎟ ⎜⎜ 27 − + ⎟⎟e
81 3π ⎜⎝ a o ⎠ ⎝ ao a o2 ⎠
Los orbitales s no tienen componentes angulares por lo que es posible realizar los
cálculos con una hoja adecuada de cálculo como por ejemplo Excel.
Para el átomo de hidrógeno Z=1 y utilizando para r unidades del radio de bohr resulta
r
−
Ψ1s = 0,564 * e −r ; Ψ = 0,1(2 − r ) e 2
2s
−r
Ψ3s = 4.10 −3
(27 − 18r + 2r ) 2
e3
Orbital Ψ1s
orbital 1s
0,6
0,5
0,4
0,3
Ψ
0,2
0,1
0
0 1 2 3 4 5 6
r en radios de bohr
Fig.1
La función de onda disminuye exponencialmente y a cinco radios de bohr es
prácticamente nula.
P = ∫ Ψ 2 dV
se puede hacer otra representación gráfica. Con centro en el núcleo consideramos una
capa esférica de radio r y espesor dr, el volumen de esa capa esférica es
dV = 4π r 2 dr
orbital 1s
0,6
0,5
0,4
2
4*π *r *Ψ
2
0,3
0,2
0,1
0
0 1 2 3 4 5 6
r en radios de bohr
Fig.2
b) el área bajo la curva mide la probabilidad de encontrar al electrón. Dado que la curva
se extiende hasta el infinito esa área vale 1 o sea probabilidad 100%.
En general la representación de la probabilidad se limita a 0,9 (90%). Si medimos el
área comprendida entre r = 0 y r ≈ 3 radios de bohr, el área vale aproximadamente 0,9.
En consecuencia el orbital es una esfera centrada en el núcleo y cuyo radio es
aproximadamente igual a 3 radios de bohr indica que en ese volumen la probabilidad de
encontrar al electrón es el 90%.
Si en vez de dibujar el orbital tratamos de representar la probabilidad de encontrar al
electrón, mediante un diagrama de puntos, dibujamos un círculo cuyo radio es
aproximadamente 3 radios de bohr (con mayor precisión 2,6 radios de bohr) en una
imagen de puntos como la de la figura 3.
Fig.3
Nota
Este dibujo se hizo a escala.
Originalmente se tomó 8mm=1 Bohr, haciendo circunferencias concéntricas, cada
milímetro, después de fueron insertando punto a punto proporcionalmente al indicativo
de densidad de probabilidad radial. El perímetro del OA, está hecho a partir del 90%,
o sea 2,66 bohr, o sea 21 mm aproximadamente.
Si imprime este dibujo esa escala varía, pero siempre encontrará que la región de
máxima probabilidad se encuentra a 1 bohr de distancia del núcleo.
Orbital Ψ2s
r
−
Dando valores a r en la función Ψ = 0,1(2 − r )e 2 se obtiene la curva de la figura 4
2s
orbital Ψ 2s
0,25
0,2
0,15
Ψ
0,1
0,05
-0,05
0 2 4 6 8 10 12
r, en radios de bohr
Fig.4
La figura 4 nos indica que la función de onda entre r =0 y r =2 es positiva y a partir de
r =2 negativa, siendo nula en r =2. Prácticamente a unos 12 radios de bohr la función
de onda es nula.
La representación espacial de la función de onda es una esfera pero debido a que se
anula su valor en r =2 y además toma valores positivos y negativos no es posible hacer
una representación tridimensional clara.
Si dibujamos la función 4π r 2 Ψ frente a r obtenemos la figura 5
orbital 2s
0,2
0,19
0,18
0,17
0,16
0,15
0,14
0,13
2
0,12
4π *r Ψ
0,11
2
0,1
0,09
0,08
0,07
0,06
0,05
0,04
0,03
0,02
0,01
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
r en radios de bohr
Fig.5
Si se comparan las figuras 2 y 5, se observa que el primer máximo del orbital 2s está
más cerca del núcleo que el máximo del orbital 1s.
Fig.6
Cuando r =2 radios de bohr aparece una zona en blanco (en realidad es una corona
circular de cierto espesor, que depende de la escala que se ha utilizado para dibujar)
que corresponde a la región donde la probabilidad de encontrar al electrón es nula. Si se
pusiesen juntas las representaciones bidimensionales de los orbitales 1s y 2s y se
respetasen sus tamaños, el orbital 2s tendría unas tres veces mayor tamaño que el 1s,
abarcando en ambos casos la misma probabilidad del 90%.
Orbital Ψ3s
−r
Dando valores a r en la función Ψ3s = 4.10 −3
(27 − 18r + 2r )
2
e3 se obtiene la curva de
la figura 7
0,12
0,1
0,08
0,06
Ψ
0,04
0,02
0
-0,02 0 5 10 15 20 25
-0,04
r en radios de bohr
Fig.7
0,11
0,1
0,09
0,08
0,07
2
4π r Ψ
0,06
2
0,05
0,04
0,03
0,02
0,01
0
0 1 2 3 4 5 6 7 8 9 10 11 1213 14 15 16 17 18 19 2021 22 23 24
r en radios de bohr
Fig.8
El área bajo la curva mide la probabilidad de encontrar al electrón. Si nos limitamos a
una probabilidad del 90%, el área bajo la curva comprendida entre r =0 y
aproximadamente r = 20 radios de bohr.
Si se comparasen los tamaños de los orbitales para el 90 % de probabilidad el orbital 2s
sería aproximadamente 3 veces mayor que el 1s y el 3s aproximadamente 6 veces
mayor
probabilidad orbital 1s
1,1
1
P= 1-(exp(-2r))*(1+2r+2r )
2
0,9
0,8
0,7
0,6
0,5
0,4
0,3
0,2
0,1
0
0 0,2 0,4 0,6 0,8 1 1,2 1,4 1,6 1,8 2 2,2 2,4 2,6 2,8 3 3,2 3,4 3,6 3,8 4
r en radios de bohr
Fig.9
Una probabilidad del 90 % de encontrar al electrón es tal como indica la gráfica r = 2,6
radios de bohr
Para el orbital 2s
( ) ( )
= ∫ 0,125 * 4r 2 + r 4 − 4r 3 e − r dr = ∫ 0,5r 2 + 0,125r 4 − 0,5r 3 e − r dr =
= 1 − e −r * (1 + r + 0,5r 2
+ 0,125r 4 )
Con la hoja de cálculo Excel representamos P frente a r (fig. 10)
probabilidad orbital 2s
1,1
1
0,9
0,8
0,7
0,6
P
0,5
0,4
0,3
0,2
0,1
0
0 1 2 3 4 5 6 7 8 9 10 11 12
r en radios de bohr
Fig.10
Una probabilidad del 90 % de encontrar al electrón es tal como indica la gráfica de la
figura 10, r ≈ 9 radios de bohr
Para el orbital 3s
⎛ −r ⎞ 2 −2 r
P=∫
⎜
( ) ⎟
(
⎜ 4.10 −3 27 − 18r + 2r 2 e 3 ⎟ * 4πr 2 dr = 2.10 − 4 27 − 18r + 2r 2
∫ )
2
*r 2
* e 3 dr =
⎝ ⎠
−2r
= 1− e 3 ⎛1 + 2 2 4
r + r2 + r4 +
8 5
r +
8 6⎞
⎜ r ⎟
⎝ 3 9 81 729 8561 ⎠
1,1
1
0,9
0,8
Probabilidad
0,7
0,6
0,5
0,4
0,3
0,2
0,1
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22
r en radios de bohr
Fig.11
Una probabilidad
del 90 % de
encontrar al
electrón es tal
como indica la
gráfica de la
figura 11, r
≈ 19,5 radios de
bohr.
Si se dibujan en
dos dimensiones
los tres orbitales
comparando sus
tamaños se obtie-
ne la fig 12.
Comparación
entre orbitales s y
con el átomo de
Bohr
Fig.12
ORBITALES ATÓMICOS HIDROGENOIDES p
La enseñanza del tema de orbitales atómicos en la Química de Enseñanzas Medias, siempre es muy
problemática, por la interpretación matemática de la función de onda, de su cuadrado y la densidad de
probabilidad radial. Para intentar aclararlo, vamos a aplicar dichos conceptos a los orbitales atómicos 2p y
3p. Se han elegido precisamente estos orbitales porque siempre dan lugar a preguntas curiosas por parte del
alumno, que todavía no sabe disociar el concepto de órbita del de orbital y el de densidad de probabilidad.
Muchas veces hemos escuchado comentarios tan curiosos referidos a los orbitales p como que “los electrones
hacen ochos junto al núcleo”, que ¿cómo era posible que un electrón pasara por el núcleo desde una parte del
orbital a otra? Para intentar que comprendan esto, diremos que los orbitales atómicos son consecuencia de un
tratamiento matemático, partiendo de la ecuación de una onda, al cuadrado de su característica fundamental o
función de onda se le atribuyó un significado físico asociándolo a la densidad de probabilidad de encontrar
al electrón, a una determinada distancia del núcleo atómico. Así como se ha visto en el tema de OAs, la
máxima densidad de probabilidad radial para la función R1s, coincidía con el radio de Bohr para el
hidrógeno en su estado fundamental.
Los factores numéricos k que aparecen multiplicando ambas funciones se determinan con las condiciones de
normalización, que implican que la probabilidad total de encontrar al electrón en el espacio es 1.
Así, aplicando ∫ r 2 (k i ψ i ) dr = 1 , se obtiene que k1=0,204 y k2=0,0202.
2
0,1
Fig. 2
0,075
0,05
OA 3p
0,025
Si se repite el método y aplica a la
0 ecuación de R3p se obtienen las
0 1 2 3 4 5 6 7 8 9 10 11 12 siguientes gráficas:
radio/ bohr
0,085
0,08
0,075
0,07
0,065
0,06
0,055
0,05
0,045
0,04
0,035
0,03
0,025
0,02
0,015
Fig. 3
0,01
0,005
0
-0,005
-0,01 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22
-0,015
-0,02
-0,025
-0,03
-0,035
-0,04
-0,045
radio/ bohr
2 2 2 2
Densidad de probabilidad radial r Ψ 3p (r Ψ /r)
0,11
0,1
0,09
0,08
0,07
0,06
Fig. 4
0,05
0,04
0,03
0,02
0,01
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22
radio /bohr
En la figura 3 el máximo de R3p en función del radio está a 1,8 bohr (hay un efecto de penetración o
acercamiento al núcleo al aumentar n), pero en la figura 4, se observa que el primer máximo de la densidad
de probabilidad radial tiene lugar a unos 3 bohr de distancia al núcleo, mientras que el segundo lo hace a
unos 12 bohr. La función se anula para valores de r = 0 y 6 bohr. Este último valor indica que a esa
distancia existe una superficie nodal, esto es, una superficie en la cual la densidad radial de probabilidad es
cero.
La representación de la función que
Probabilidad acumulativa 2p y 3p da la probabilidad acumulativa (PA)
con la distancia al núcleo para los
1 OAp realizada mediante cálculo
0,95
0,9 numérico a partir de la suma de los
0,85
0,8
0,75
PA 2p
PA 3p valores r2R2 multiplicada por dr, se
0,7 muestra en la fig 5.
0,65
0,6
0,55 Si se considera el OA como la
PA
0,5
0,45 superficie que encierra el 90% de
0,4
0,35
0,3
densidad de probabilidad radial, se
0,25 tendrá que para el OA 2p, la
0,2
0,15
0,1 distancia de extensión es de 8,1 bohr
0,05
Todo 0lo anterior hace referencia a las funciones radiales. (si fuera el 99% sería 11 bohr),
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 mientras que para el OA 3p, es de
OA p sobre los ejes de coordenadas cartesianas.
r(bohr) 18,5 bohr (para el 99% sería de 24,5
bohr).
Fig. 5
Sin embargo para orbitales determinados por condiciones cuánticas en las que l>0, la componente radial se
ve influenciada por una componente angular que condicionará la forma del orbital y modificará las zonas de
máxima probabilidad. Ahora existirán tres variables: r, y los ángulos N y 2, que forma r con los ejes .
De esta forma la función R debe reunir ambas componentes, tomando similares valores para cada par (n,l).
Según el valor empleado de m (1, 0, y -1) se orientará sobre diferente eje de coordenadas.
Para representar el OA de forma completa, se usará un programa de dibujo tridimensional, en ejes
cartesianos, lo que implicará pasar de coordenadas polares a rectangulares con las conversiones necesarias.
En este caso x= r sen Ncos 2, y = r sen N sen 2, y z= r cos 2 y r2 = x2+y2+z2
Las funciones características de los orbitales atómicos 2p toman las expresiones siguientes para el caso del
hidrógeno Z=1,y en función del radio de Bohr:
x2 + y2 + z2
−
Ψ(2,1,0) = 0,1ze 2
que corresponde al OA 2pz
x2 + y2 + z2
−
Ψ(2,1,1) = 0,1xe 2
que corresponde al OA 2px
x2 + y2 + z2
−
Ψ(2,1,−1) = 0,1 ye 2
que corresponde al OA 2py
Se va a representar como ejemplo la función R2px. El programa matemático empleado, genera una malla
superficial, en la que cada punto (x,y) corresponde valor de ψ. Como se trabaja sobre el plano XY, se hace
z=0.
Dado que los valores de ψ pueden ser
positivos o negativos, la malla
superficial contendrá “montañas ψ >0
y “abismos”, ψ <0. Los abismos
corresponden a valores negativos y en
los dibujos se representan en la escala
de grises hacia el negro; cuanto más
oscuro, más negativo, mientras que los
valores positivos se toman en la escala
de grises hacia el blanco, cuanto más
blanco más positivo (fig.6). La
proyección sobre el plano XY,
conforma unas superficies de nivel,
que determinarán la forma del orbital
atómico 2p, con dos lóbulos bien
diferenciados, debido a la aparición de
un plano nodal, puesto que l=1, que
será la característica de todos los OAp
Fig. 6
Fig. 7
Si se eleva la función R al cuadrado se obtiene la densidad de probabilidad que se representa en las figuras 8
y 9. Esta última se puede comparar con la gráfica de la función de onda 2px, que se da en la figura 7 y
observar cómo se modifica el tamaño de la línea de contorno y las isolíneas correspondientes. En este caso
cada punto de la malla superficial indicará una densidad de probabilidad, que al proyectarse sobre el plano
base (el XY para el OA 2px), producirá las curvas de nivel de probabilidad en las que R2 es constante.
La línea de contorno en el
plano a que se haga referencia
o superficie límite (si se
tratara de volúmenes),
indicará la zona del espacio
que encierra una densidad de
probabilidad acumulativa del
90% que es lo que se dibuja
habitualmente como orbital
atómico.
En la figura 8, se observa la
génesis del orbital atómico 2px,
obtenida a partir de las
soluciones de la ecuación de
Schrödinger para los valores in-
dicados en coordenadas rectan-
gulares (visión tridimensional),
la visión de los diagramas de
contorno con la distribución
probabilística acumulativa sobre
en plano XY, y por fin lo que se
dibuja generalmente como OA
2px, esto es la línea de contorno
que encierra el 90%. El máximo
tal como que se observaba en la
gráfica de R/r, tiene lugar a los
2 bohr.
Fig. 8
Fig. 9
La figura 10, corresponde a la
gráfica tridimensional de la
densidad de probabilidad para
el OA 3px, y a través de ella se
puede comparar con la obtenida
para el OA 2px.
Fig. 10
Fig. 11
Si se hiciera girar la línea de
contorno externa dada sobre el
eje X se produciría la
representación espacial del
orbital atómico 3px.
Fig. 12
Hasta ahora se han visto las funciones de onda y las densidades de probabilidad sobre el eje X (m=1). Se
trabajaría de la misma forma sobre otro eje. Si se representa la función R 2py (m=-1), las gráficas serían las
dadas en las figuras 13 y 14. Como se observa, la forma es la misma y lo único que varía es la orientación
sobre el mencionado eje.
Fig. 13 Fig. 14
ORBITALES ATÓMICOS HIDROGENOIDES tipo d, (primera parte)
Si los orbitales atómicos p tenían bastante complejidad ¿qué decir de los d? El alumno de enseñanzas
medias, dibuja un trébol de cuatro hojas, que toma como modelo de este orbital y se queda tan
tranquilo, sin enterarse de lo que está haciendo.
La condición cuántica que lo determina es l=2 y las funciones radiales que se deben representar
corresponden a valores de R que para n=3, y para n=4, tienen las siguientes ecuaciones:
r
−
r
⎛ r ⎞ −
Ψ3d = 0,00902r e 2 3
y Ψ 4d = 0,007⎜1 − ⎟ r 2 e 4
⎝ 12 ⎠
Si se representan en una hoja de cálculo se obtienen las gráficas, correspondientes a las figuras 1 y 2
0 ,0 5
0 ,0 4
0 ,0 3
0 ,0 2
0 ,0 1
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22
radio /bohr
Fig. 1
0,029
0,027
0,025
0,023
0,021
0,019
0,017
0,015
0,013
0,011
0,009
0,007
0,005
0,003
0,001
-0,001
-0,003
-0,005 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42
-0,007
-0,009
-0,011
-0,013
-0,015
radio/ bohr
Fig. 2
En la función R 4d, se observa que para r=0 y para r=12, la función se anula. Este último valor va a
generar una superficie nodal (superficie interior cuya densidad de probabilidad es nula)
Si las funciones dadas se elevan al cuadrado y se multiplican por r2, se obtiene la densidad de
probabilidad radial (fig. 3 y 4).
Densidad de probabilidad radial r2 Ψ 2 3d
0,125
0,12
0,115
0,11
0,105
0,1
0,095
0,09
0,085
0,08
0,075
0,07
0,065
0,06
0,055
0,05
0,045
0,04
0,035
0,03
0,025
0,02
0,015
0,01
0,005
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22
radio/ bohr
Fig. 3
0,07
0,065
0,06
0,055 hasta aquí el 90%
0,05
0,045
0,04
0,035
0,03
0,025 hasta aquí el 99%
0,02
0,015
0,01
0,005
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42
radio /bohr
Las propiedades acumulativas revelan que hasta el 90%, las distancias máximas se separación son de
16 y 30 bohr respectivamente (fig. 5)
Fig. 4
Probabilidad acumulativa 3d y 4d
1
0,95
0,9 PA 3d
0,85 PA 4d
0,8
0,75
0,7
0,65
0,6
0,55
PA
0,5
0,45
0,4
0,35
0,3
0,25 Fig. 5
0,2
0,15
0,1
0,05
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45
r(bohr)
La contribución angular a la función de onda total, es en este caso muy grande, y depende de los
valores que tome el número cuántico m. Para m=0, se obtiene la función R dz2. Para valores de m=1,-
1 sería R dxz y R dyz y si m = 2, y -2 ; R xy y R d (x2-y2).
Otra de las preguntas características de los alumnos de enseñanzas medias respecto a los nombres de
los orbitales d, es ¿De dónde sale el orbital 3dxy? ¿En qué se parece y diferencia del 3dxz? ¿Y del
4dxy?
Para explicarlo no queda mas remedio que ir a las ecuaciones de definición considerando la
contribución radial y angular y que al ser representadas darán los diferentes OA d.
Serían:
x2 + y2 x2 + z2 z2 + y2
− − −
Ψ 3dxy = 0,00985 xye 3
; Ψ3dxz = 0,00985xze 3
; Ψ3dyz = 0,00985 yze 3
Como se aprecia, las tres funciones son similares, lo único diferenciante es el plano en que se
representan en función de las variables que intervienen. Así obtendríamos los tres orbitales d, de la
misma forma, el 3dxy, 3d xz y 3dyz, con planos nodales coincidentes con dichos ejes. Obsérvese en
las figuras 6 y 7, la forma de la función superficie para el orbital 3dxy y el detalle de las curvas de
nivel, siempre en unidades bohr de distancia atómica, en los ejes correspondientes.
Fig. 6
Las zonas mas oscuras corresponden a
valores de R3d negativos (valles), mientras
que las más claras corresponden a valores
positivos (montañas), en la gráfica de
superficie.
Fig. 7
Lo cual indica un punto nodal, cuando el factor binómico se anule. Este punto nodal engendrará una
superficie nodal, o sea una superficie interior en la que la probabilidad de encontrar al electrón es nula,
que será el factor diferenciante en la forma del orbital frente al 3dxy.
La representación de la función de onda, y el detalle de las curvas de nivel, se dan en las figuras 10 y
11. A través de ellas se puede diferenciar perfectamente ambos orbitales atómicos.
También los dibujos permiten diferenciar, los planos nodales de las superficies nodales. Los primeros
como dependen única y exclusivamente del valor de l, que es 2 para los orbitales d, ya sean 3d o 4d,
es un elemento de analogía, mientras que las superficies nodales que dependen de n-(l+1), es
diferenciante. Para el 3d, las superficies nodales serían 3-(2+1)= 0, mientras que para el 4d; 4-
(2+1)=1.
En las gráficas dadas, siguiendo el criterio general presentado en estos trabajos, a los valles, valores
negativos, les corresponden las zonas mas oscuras mientras que a las cumbres, valores positivos, les
corresponden las zonas mas claras.
Tampoco la forma de las curvas de nivel es la misma. El acotamiento al 90%, en zonas de muy poca
pendiente, hace que aparezcan eliminadas zonas superficiales que modifican el trazado habitual.
Fig. 11
Fig. 10
Fig. 13
Fig. 12
En estas gráficas, los valores de la función al estar
elevados al cuadrado siempre van a ser positivos, y por
lo tanto no existe el criterio de tonalidades de grises
anterior.
En la distribución espacial de lo que se considera como
orbital atómico, haciendo girar la figura plana por el eje
X, las superficies nodales interiores, quedan tapadas, y
no se apreciarían.
ORBITALES ATÓMICOS HIDROGENOIDES tipo d, 2 (continuación)
x2 + z2
⎛ z 2 + x 2 ⎞⎟ 2
x +z
2 2
− −
Ψ 3dz = 0,002843( 2 z − x )e
2 2 2
; Ψ 4dz = 0,0022 1 −
3 2 ⎜ ( 2 z − x )e 4
2
⎜ 12 ⎟
⎝ ⎠
Como podemos apreciar existe en ambas, la participación de la componente z2, que da nombre al
orbital, procedente de la contribución angular del orbital atómico. Sus representaciones y detalles, se
dan en las figuras 1 y 2 (OA 3d) y 3 y 4(OA4d).
Fig. 1 Fig. 2
Fig. 3
Fig. 4
La observación de los dibujos anteriores, evita cualquier tipo de aclaración. Los criterios de diseño,
son los empleados en la primera parte del tema (OAd1), o sea los valores positivos (montañas), más
claros, y los negativos (valles), más oscuros. El reforzamiento sobre el eje z, de las zonas de
probabilidad es lo más destacado.
Las funciones elevadas al cuadrado, nos dan, las densidades de probabilidad que se van a comparar
siguiendo la misma distribución anterior. Véanse las figuras 5,6, 7 y 8.
Fig. 6
Fig. 5
Fig. 7
Fig. 8
Pasemos a los OA 3d(x2-y2) y 4d(x2-y2). Las ecuaciones de ambos son:
x2 + y2
−
Ψ 3d ( x − y ) = 0,004925( x − y )e
2 2 2 2 3
y
⎛ y2 + x2 ⎞ 2 −
x2 + y2
Ψ 4d ( x − y ) = 0,003824⎜1 −
2 2 ⎟( x − y 2 ) e 4
⎜ 12 ⎟
⎝ ⎠
Como se aprecia, el factor que da nombre al orbital, surge en ambas ecuaciones, procedente de la
contribución angular. Sus representaciones gráficas y las curvas de nivel correspondientes se dan en
las figuras 9,10, 11 y 12. Las comparaciones se pueden hacer a partir de estos dibujos, que emplean los
mismos criterios de dibujos anteriores.
Fig. 10
Fig. 9
Fig. 12
Fig. 11
Por fin las densidades de probabilidad respectivas
serían las indicadas en las figuras 13,14 ,15 y 16:
Fig. 14
Fig. 13
Fig. 16
Fig. 15
ORBITALES ATÓMICOS HIDROGENOIDES tipo f ( primera parte)
Por lo general en el bachillerato, los profesores dejan al margen de sus explicaciones los orbitales
f. Describen los d, estudian los problemas que encierran, y se marginan los f, calificándolos como muy
complejos. Se argumenta que tampoco se van a explicar los elementos cuyos electrones completan esos
orbitales, aunque se estudien otros que sí los han rellenado. Al ser unos orbitales internos, no superponen
con otros orbitales y por lo tanto no dan lugar a enlaces. Por eso no son importantes.
Es evidente que se supone una pérdida de tiempo explicar los orbitales f en bachillerato, porque
realmente su descripción es compleja. Se definen a partir de unas condiciones cuánticas n=4, l=3, que
implican la existencia de 3 planos nodales en cada orbital. Esto provoca dificultades al dibujarlos.
r
−
La función radial 4f, tiene por ecuación Ψ4 f = 0,000651 r e 3 4
0,06
0,055
0,05
0,045
0,04
0,035
0,03
0,025
0,02
0,015
0,01
0,005
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35
radio/ bohr
Fig. 1
0,07
0,065
0,06 68%
0,055
90%
0,05
0,045
0,04
0,035
0,03
99%
0,025
0,02
0,015
0,01
0,005
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35
radio/ bohr
Fig. 2
A 16 bohr, existe un máximo de densidad de probabilidad, que como se ve en las gráficas de
probabilidades acumulativas abarca el 68%.
Probabilidad acumulativa 4f
1
0,95
0,9
0,85
0,8
0,75
0,7
0,65
0,6
0,55
PA
0,5
0,45
0,4
0,35
0,3
0,25
0,2
0,15
0,1
0,05
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36
r(bohr)
Fig. 3
Si comparamos con la gráfica 4d (orbitales atómicos d, primera parte), observamos que en la 4d, el
máximo tenía lugar a 22 bohr, mientras que en el 4f ocurre a 16 bohr ; es más interno. En el 4d, existe un
punto nodal a 12 bohr, mientras que en el 4f no existe, dado que n-(l+1) = 0.
El 90% de probabilidad se encuentra encerrado en un volumen de radio 22 bohr, mientras que en el 4d, el
radio que encerraba dicho porcentaje era de 30 bohr.
Donde se dan las mayores diferencias es en la función angular, y por lo tanto en la forma del propio
orbital. Ya se ha explicado que depende de los valores de m, que en este caso pueden ser siete: ±3, ±2, ±1
y 0. Por lo tanto existen 7 orbitales 4f.
Estudiaremos primero aquellos orbitales 4f condicionados por m = 0, m = ±1.
En el primer caso el orbital atómico recibe el nombre de 4fz3. Suele recibir este nombre debido al término
con mayor exponente que aparece en la función o con el producto de las variables que lo definen, como
va a ocurrir en la mayoría de los OA 4f.
La función angular es, en coordenadas rectangulares:
5 z 2 − 3r 2
F ( x, y, z ) = 0,167 z 3
, que modificada teniendo en cuenta que r 2 = x 2 + y 2 + z 2 , y
r
después de incorporar la parte radial común a todas, produce una función R4f(4,3,0) que al representarse
mediante un programa matemático, genera el siguiente gráfico superficial (para reducir una dimensión y
representar las características más sobresalientes del orbital, proyectamos sobre los ejes “y” y “x” el corte
de la función de ondas cuando y=z).
En las figuras 4 y 5, se dibujan la función superficial (fig 4) y su proyección o mapa de contorno (fig.5).
15.00
10.00
5.00
Eje X
0.00
-5.00
CURVAS DE NIVEL
Proyección sobre el plano XZ -10.00
-15.00
-20.00
PLANOS
NODALES -25.00
-30.00
-25.00 -20.00 -15.00 -10.00 -5.00 0.00 5.00 10.00 15.00 20.00 25.00
Fig. 5
5.00
Eje X
0.00
-5.00
-10.00
-15.00
CURVAS DE NIVEL
-20.00
-25.00
-30.00
-25.00 -20.00 -15.00 -10.00 -5.00 0.00 5.00 10.00 15.00 20.00 25.00
Eje Z
Fig. 6 Fig. 7
Línea de contorno
90% de probabilidad La línea de contorno marca el
PLANOS
dibujo del OA 4f. Como se
NODALES observa, no sería muy difícil de
Fig. 5 dibujar para un alumno de
Enseñanzas Medias.
Con las condiciones cuánticas n = 4, l = 3, m = ±1, se tienen los denominados orbitales 4fx3 y el 4fy3, que
son exactamente iguales aunque orientados en los ejes dominantes (los que le dan el nombre).
La funciones angulares respectivas en coordenadas rectangulares son:
5z 2 − r 2 5z 2 − r 2
F ( x, y, z ) = 0,0342 x y F ( x, y, z ) = 0,0342 y
r3 r3
Operando como en el caso anterior, produce los siguientes gráficos para R4f(4,3,±1) (fig 8 y 9).
CONDICIONES
CUÁNTICAS
n=4 ; l=3 ; m=+1,-1
CURVAS DE NIVEL
Proyección sobre el plano XZ
Eje Z
PLANOS
NODALES
Eje X
Fig. 8
Eje Z
20.00
15.00
10.00
5.00
Este orbital atómico se conoce con el
0.00
nombre de 4fx3.
-5.00 La representación del OA 4fy3, sólo
Eje X
cambiaría el eje x por el y, al tener en
-10.00
cuenta la simplificación de ejes.
-15.00
-20.00
-25.00
-25.00 -20.00 -15.00 -10.00 -5.00 0.00 5.00 10.00 15.00 20.00 25.00
Fig. 9
Al elevar al cuadrado las funciones anteriores, y representarlas, se obtienen las figuras 10 y 11. La 11
corresponde a las curvas de nivel de densidad de probabilidad, cuya línea de contorno definiría el dibujo
del orbital atómico.
DENSIDAD DE PROBABILIDAD 4fx3
Curvas de nivel
25.00
DENSIDAD DE PROBABILIDAD 4fx3 Eje Z
CONDICIONES
20.00
CUÁNTICAS
n=4; l=3; m=+1,-1
15.00
10.00
5.00
Eje X
0.00
-10.00
-15.00
-20.00
-25.00
-25.00 -20.00 -15.00 -10.00 -5.00 0.00 5.00 10.00 15.00 20.00 25.00
Línea de contorno
90% de probabilidad
Fig. 11
PLANOS
NODALES
Fig. 10
Fig. 5
Fig. 5
ORBITALES ATÓMICOS HIDROGENOIDES tipo f (continuación)
Siguiendo la descripción de los OA 4f, se van a representar los condicionados por los números cuánticos
n=4, l=3, m= ±2. Estos OA 4f se denominan en función de las funciones angulares que los representan:
4fz(x2-y2) y 4fxyz.
Las funciones angulares respectivas en coordenadas rectangulares son:
x2 − y2
(
Fz x − y2 2
) = 0,04324 z
r3
y Fxyz = 0,04324
2 xyz
r3
En las figuras se dibujan la función superficial (fig 1) y su proyección o mapa de contorno (fig.2) que
corresponden al primer orbital el 4fz(x2-y2)
15.00
10.00
5.00
0.00
-5.00
CURVAS DE NIVEL
-10.00
PLANOS -15.00
NODALES
-20.00
-25.00
-30.00
-20.00 -15.00 -10.00 -5.00 0.00 5.00 10.00 15.00 20.00
Fig. 2
Fig. 1
Eje Z
Eje Z
Curvas de nivel
Curvas de nivel
Línea de contorno
Probabilidad 95% Línea de contorno
99% de Probabilidad
PLANOS
PLANOS NODALES
NODALES
Fig. 3 Fig. 4
20.00
20.00
20.00
15.00
15.00
10.00
10.00 10.00
5.00
5.00
0.00 0.00
0.00
-5.00
-5.00
-10.00
-10.00
-10.00
-15.00 -20.00
-15.00
-20.00
-20.00
-30.00
-25.00
-25.00
-20.00 -15.00 -10.00 -5.00 0.00 5.00 10.00 15.00 20.00 -20.00 -10.00 0.00 10.00 20.00
-20.00 -15.00 -10.00 -5.00 0.00 5.00 10.00 15.00 20.00
La línea de contorno marca el dibujo del OA 4f. Como se observa, solo aparecen los 3 planos nodales
para probabilidades acumulativas a partir del 95%.
Siguiendo el mismo procedimiento a partir del orbital R4fxyz se obtiene la función:
( x2 + y2 )
−
Ψ ( x, y, z ) = 0,00002814 ( 2 xy )e 2 4
,
Para visualizar los lóbulos del orbital, se nuevo se eligirán las condiciones de contorno para el 90 y 99%
de probabilidad (fig 8 y 9)
ORBITAL ATÓMICO 4f xyz ORBITAL ATÓMICO 4f xyz
Curvas de nivel Curvas de nivel
90% de Probabilidad 99% de Probabilidad
40.00 40.00
Eje Y Eje Y
30.00 30.00
20.00 20.00
10.00 10.00
0.00 0.00
Eje X Eje X
-10.00
-10.00
-20.00
-20.00
-30.00
-30.00
-40.00
-40.00 -40.00 -30.00 -20.00 -10.00 0.00 10.00 20.00 30.00 40.00
-40.00 -30.00 -20.00 -10.00 0.00 10.00 20.00 30.00 40.00
Fig. 8 Fig. 9
Como la gráfica del cuadrado de la función no presenta ninguna variación frente a la del orbital, se obvia
su representación.
Con las condiciones cuánticas n = 4, l = 3, m = ±3, se tienen los orbitales 4fx(x2-3y2) y el 4fy(3x2-y2), que
tienen la misma forma aunque orientados hacia sus respectivos ejes dominantes (los que le dan el
nombre).
x2 − 3y 2 3x 2 − y 2
( )
Fx x 2 − 3 y 2 = 0,0265 x
r3
y ( )
Fy 3x 2 − y 2 = 0,0265 y
r3
Operando como en los casos anteriores, obtenemos los siguientes gráficos (fig 10 y 11) para el orbital
R4f (x(x2-3y2))
ORBITAL ATÓMICO 4f x(x2-3y2)
ORBITAL ATÓMICO 4f x(x2-3y2) Curvas de nivel
CONDICIONES 25.00
CUÁNTICAS Eje Y
n=4; l=3; m=+3,-3
20.00
15.00
10.00
5.00
CURVAS DE NIVEL
Proyección sobre el plano XY 0.00
Eje X
-5.00
-10.00
-15.00
PLANOS
NODALES
-20.00
-25.00
-30.00 -25.00 -20.00 -15.00 -10.00 -5.00 0.00 5.00 10.00 15.00 20.00 25.00 30.00
Fig. 11
Fig. 10
Al elevar al cuadrado las funciones anteriores, y representarlas, se obtienen las figuras 12 y 13. La 13
corresponde a las curvas de nivel de densidad de probabilidad, cuya línea de contorno define el dibujo del
orbital atómico.
10.00
5.00
0.00
Curvas de nivel Eje X
-5.00
-10.00
-15.00
-20.00
Línea de contorno
90% de probabilidad
-25.00
-30.00 -25.00 -20.00 -15.00 -10.00 -5.00 0.00 5.00 10.00 15.00 20.00 25.00 30.00
PLANOS
NODALES
Fig. 13
Fig. 12
El orbital que nos queda es simétrico del anterior, con la única permutación de los ejes X e Y . Véanse las
figuras 14,15 16 y 17
DENSIDAD DE PROBABILIDAD 4fy(3x2-y2)
ORBITAL ATÓMICO 4fy(3x2-y2) Curvas de nivel
30.00
CONDICIONES Eje Y
CUÁNTICAS 25.00
n=4; l=3; m=+3,-3
20.00
15.00
10.00
5.00
Eje X
CURVAS DE NIVEL 0.00
Proyección sobre el plano XY
-5.00
PLANOS -10.00
NODALES
-15.00
-20.00
-25.00
-30.00
-25.00 -20.00 -15.00 -10.00 -5.00 0.00 5.00 10.00 15.00 20.00 25.00
Fig. 14
Fig. 15
20.00
15.00
10.00
5.00
Eje X
Curvas de nivel 0.00
-5.00
-10.00
-15.00
-20.00
Línea de contorno
90% de probabilidad
-25.00
PLANOS
NODALES
-30.00
-25.00 -20.00 -15.00 -10.00 -5.00 0.00 5.00 10.00 15.00 20.00 25.00
Fig. 17
Fig. 16
ORBITALES HIBRIDOS sp
La enseñanza del tema de orbitales híbridos (OH) en la Química de Enseñanzas Medias está llena de
tópicos que trataremos de resolver y aclarar.
En primer lugar, la forma. Aparecen con un lóbulo muy grande y otro muy pequeño, debido a la
combinación lineal de los orbitales s y p, pero ¿tienen siempre esa forma?. ¿Es igual que se mezclen un
orbital s y uno p que un s y dos o tres orbitales atómicos p?. ¿Tiene la misma forma el orbital híbrido
producido por la mezcla de un 2s y un 2p, que un 3s y un 3p?. Siempre se dibujan iguales, pero no lo son.
Nunca se habla de la densidad de probabilidad en las zonas que delimitan los lóbulos ¿dónde están las curvas
de los niveles de probabilidad?
Para aclararlo usaremos el mismo sistema empleado en el tratamiento de los orbitales atómicos p, d y f,
trabajando con funciones de onda en coordenadas cartesianas.
Las funciones de onda que se combinarán para formar dos orbitales híbridos sp son:
ψ1 =
1
(s + p x ) y ψ2 =
1
(s − px ) (1)
2 2
Ahora bien, la función de onda 2s como se ha visto tiene por ecuación general:
3
Z 'r
⎛ 1 ⎞ ⎛ Z' ⎞2 ⎛ Z ' r ⎞ − 2 a0
ψ 2s =⎜ ⎟ * ⎜⎜ ⎟⎟ * ⎜⎜ 2 − ⎟*e (2)
⎝ 32π ⎠ ⎝ a 0 ⎠ ⎝ a 0 ⎟⎠
mientras que la 2p
5 5
Z 'r Z 'r
⎛ 1 ⎞ ⎛ Z' ⎞2 − ⎛ Z' ⎞2 −
ψ 2 px =⎜ ⎟ * ⎜⎜ ⎟⎟ * xe 0 = 0.1 * ⎜⎜ ⎟⎟ * xe 2 a0 (3)
2a
⎝ 32π ⎠ ⎝ a0 ⎠ ⎝ a0 ⎠
Hasta ahora se habían usado para el hidrógeno, esto es en orbitales atómicos hidrogenoides en los que Z’=1.
Sin embargo los orbitales atómicos que normalmente se hibridan corresponden a átomos característicos,
como son C, N, etc, y por lo tanto habrá que introducir en las ecuaciones 2 y 3, el valor de Z’ que es la carga
nuclear efectiva de cada átomo, esto es la carga nuclear modificada por los parámetros de Slater. Este hecho
es el factor que modificará esencialmente el tamaño del orbital híbrido.
En nuestro trabajo tomaremos para el C, Z’=3,15. Para el N, Z’=3,8. Para el O, Z’=4,45 y para el F, Z’=5,1
( ))
⎛
⎜
Ψsp = 0,385 2 − 3,15 x + y e
⎜
2 2
(−
2 ⎟
⎟
+ 1,232 xe
−
2
(4)
⎝ ⎠
Esta expresión correspondería a una de las combinaciones lineales; bastaría con cambiar el signo para tener
la otra.
Para representar el OH, se usará el programa de dibujo tridimensional que se ha empleado en trabajos
anteriores obteniéndose para el carbono, el gráfico superficial de la fig.1 , con el corte de líneas de nivel, que
se puede apreciar mejor en la fig.2.
El programa genera una malla superficial, en la que cada punto (x, y) corresponde valor de ψ. Como se
trabaja sobre el plano XY, se hace z=0.
Fig. 2
Fig.1
Dado que los valores de ψsp pueden ser positivos o negativos, la malla superficial contendrá “montañas ψ
>0 y “abismos”, ψ <0. Los abismos corresponden a valores negativos y en los dibujos se representan en la
escala de grises hacia el negro; cuanto más oscuro, más negativo, mientras que los valores positivos se toman
en la escala de grises hacia el blanco, cuanto más blanco más positivo. La proyección sobre el plano XY,
conforma unas superficies de nivel, que determinarán la forma del orbital híbrido sp para el carbono, con dos
lóbulos bien diferenciados, separados por un plano nodal(fig.2). Como se puede apreciar, el mayor desnivel
se encuentra en el lóbulo pequeño, en contra de lo que se cree.
Si ahora se eleva al cuadrado la funciones anteriores, para buscar las densidades de probabilidad, se
obtendrán, gráficos bien diferentes. Los cortes son parecidos, con una mayor diferencia de densidades en
función de la distancia. El lóbulo pequeño en contra de lo que se cree es el que tiene mayor densidad de
probabilidad, pero debido a su proximidad al núcleo del átomo presentaría mayores inconvenientes en la
superposición.
Fig. 4
Fig. 3
Se van a comparar ahora estos orbitales y densidades de probabilidad con los obtenidos en la hibridación de
un 3s y un 3p, para el átomo “hermano” del carbono; el silicio, a fin de observar si su forma es similar o no.
La combinación lineal es la misma, pero no así las funciones que en este caso serían:
3 − Zr
⎛Z
1 ⎞2 ⎛ 18Zr 2Z 2 r 2 ⎞ 3a o
Ψ3s = ⎜ ⎟⎟ ⎜⎜ 27 − + ⎟⎟e (5)
81 3π ⎜⎝ a o
y
⎠ ⎝ ao a o2 ⎠
5 5
⎛ 2 ⎞ ⎛ Z ⎞2 ⎛ Zr ⎞ ⎛ Z ⎞2 ⎛ Zr ⎞
Zr Zr
− −
ψ 3 px = ⎜⎜ ⎟⎟ * ⎜⎜ ⎟⎟ * ⎜⎜ 6 − ⎟⎟ * xe 3a0 = 0.01 * ⎜⎜ ⎟⎟ * ⎜⎜ 6 − ⎟⎟ * xe 3a0 (6)
⎝ 81 π ⎠ ⎝ a0 ⎠ ⎝ a0 ⎠ ⎝ a0 ⎠ ⎝ a0 ⎠
La carga nuclear efectiva Z para el Si es 4,28. Sustituyendo estos valores en 4 y 5, y llevándo las ecuaciones
a la 1, dan lugar a la función R3sp, que se representa como en los casos anteriores obteniéndose:
Fig. 7
Fig. 6
Si se eleva la función R3sp al cuadrado se obtiene la densidad de probabilidad que se representa en las
figuras 8 y 9.
Fig. 9
Fig. 8
Las comparaciones demuestran las diferencias entre estos orbitales híbridos de elementos de la misma
familia, no solamente en el tamaño, como es lógico, sino también en la forma, como corresponde a la
combinación lineal de funciones diferentes. Como podemos apreciar, la función R3sp (fig.7), lleva incluida
la forma del sp anterior, surgiendo una superficie nodal interior, como ocurría en los OA sin hibridar. Los
dos planos nodales se aprecian mucho mejor en las figuras 8 y 9, de las densidades de probabilidad, al
delimitar esta al 90%.
Se tratará ahora de ir hacia la derecha en el sistema periódico para observar la variación en la forma y el
tamaño de los mismos orbitales híbridos, para el N y F. Como se ha dicho se modificará la carga nuclear
efectiva, en la fòrmula de la función de onda para el carbono y se sustituirá por el valor correspondiente a
cada átomo. Las representaciones gráficas corresponderán a las figuras
Fig. 11
Fig. 10
Fig. 12 Fig. 13
Se observa muy poca diferencia de tamaño entre los orbitales híbridos sp del carbono y del nitrógeno, pero al
saltar dos posiciones en el sistema periódico y pasar al flúor, la diferencia de tamaños se hace mas
considerable. Mientras en el flúor el lóbulo grande avanza hasta los 2 bohr, en el carbono hasta los 3,25 y en
el nitrógeno hasta los 3 bohr.
ORBITALES HIBRIDOS sp2
Las funciones de onda que se combinarán para formar dos orbitales híbridos sp2 son la 2s, 2px y 2py en las
siguientes combinaciones lineales:
1 2 1 1 1 1 1 1
Ψ1 = s+ px Ψ2 = s− px + py Ψ3 = s− px − py (1)
3 3 3 6 2 3 6 2
Ahora bien, la función de onda 2s como se ha visto tiene por ecuación general:
3
Z 'r
⎛ 1 ⎞ ⎛ Z' ⎞2 ⎛ Z ' r ⎞ − 2 a0
ψ 2s =⎜ ⎟ * ⎜⎜ ⎟⎟ * ⎜⎜ 2 − ⎟*e (2)
⎝ 32π ⎠ ⎝ a 0 ⎠ ⎝ a 0 ⎟⎠
⎝ 32π ⎠ ⎝ a 0 ⎠ ⎝ a0 ⎠
y la 2py
5 5
⎛ 1 ⎞ ⎛ Z ' ⎞2 ⎛ Z ' ⎞2
Z 'r Z 'r
− −
ψ 2 py =⎜ ⎟ * ⎜⎜ ⎟⎟ * ye 0
2a
= 0.1 * ⎜⎜ ⎟⎟ * ye 2 a0 (4)
⎝ 32π ⎠ ⎝ 0⎠a ⎝ a0 ⎠
siendo como en desarrollos anteriores, Z’ la carga nuclear efectiva, y a0 = 1, al referirse a unidades bohr.
Sustituyendo 2,3 y 4 en 1, se obtendrán las funciones de onda de los orbitales híbridos a representar
Como se trata de comparar orbitales híbridos , siempre se representará para cada átomo, la función R1
En el C ,
3,15 (x 2 + y 2 ) ⎞ 3,15 (x 2 + y 2 )
( ))
⎛
⎜
Ψ1sp 2 = 0,335 2 − 3,15 x + y e
⎜
2 2
−
( 2 ⎟
⎟
+ 1,4335 xe
−
2
⎝ ⎠
Fig. 4
Fig. 3
Fig. 6
Fig. 5
Fig. 7
Si ahora se trabaja con las funciones de onda 3s, 3px y 3py, a fin de formar los orbitales híbridos de
elementos del mismo grupo pero del periodo siguiente, se deberá considerar la variación de la carga nuclear
efectiva, que para el Si es 4,28, para el P, 4,886, para el S, 5,48 y para el Cl, 6,116
Fig. 10
Fig. 11
La observación detallada de los dibujos permite establecer las diferencias entre estos orbitales y los de
los elementos que encabezan sus respectivos grupos, tanto en el lóbulo enlazante como en el no apto para la
superposición. Por lo tanto he aquí aclarado uno de los tópicos mas generalizados en los libros de textos; no
se pueden dibujar los orbitales híbridos sp2 de diferentes átomos como si fueran iguales, ya que no lo son; ni
en la forma ni en el tamaño.
ORBITALES HIBRIDOS sp3
Para terminar el tema de orbitales híbridos s y p, se tratará de los sp3. Las funciones de onda que se
combinarán en principio para formar cuatro orbitales híbridos sp3 son la 2s, 2px , 2py y 2pz en las
siguientes combinaciones lineales:
1 3 1 2 1 1 1 1 1
Ψ1 = s+ pz ; Ψ2 = s+ px − pz ; Ψ3 = s− px + py − pz y
2 2 2 3 12 2 6 2 12
1 1 1 1
Ψ4 = s− px − py + pz (1)
2 6 2 12
Ahora bien, la función de onda 2s como se ha visto tiene por ecuación general:
3
Z 'r
⎛ 1 ⎞ ⎛ Z' ⎞2 ⎛ Z ' r ⎞ − 2 a0
ψ 2s =⎜ ⎟ * ⎜⎜ ⎟⎟ * ⎜⎜ 2 − ⎟⎟ * e (2)
⎝ 32π ⎠ ⎝ a 0 ⎠ ⎝ a 0 ⎠
x2 + y2 + z2 x2 + y2 + z2
− −
mientras que la 2px : Ψ(2,1,1) = 0,1xe 2
; la 2py : Ψ( 2,1,−1) = 0,1 ye 2
x2 + y 2 + z2
−
y la 2pz : Ψ(2,1,0) = 0,1ze 2
(3)
Como en el tema anterior, se tomará para el C, Z’=3,15. Para el N, Z’=3,8. Para el O, Z’=4,45 y para el F,
Z’=5,1
Como se trata de comparar orbitales híbridos siempre se representará para cada átomo, la función R1 de sus
orbitales híbridos sp3
La representación aplicando el programa matemático de gráficas superficiales, tomando el eje z como x tal
como en temas anteriores, viene dada por las figuras 1 a 8
Fig. 2
Fig. 1
Fig. 3 Fig. 4
Fig. 6
Fig. 5
Fig. 7 Fig. 8
Como se observa especiamente en los dibujos de curvas de nivel, el avance del extremo del lóbulo enlazante
varía desde el carbono que llega hasta los 3,8 bohr, hasta el flúor que alcanza los 2,5 bohr. El lóbulo no
enlazante varía también desde -3 bohr hasta -1,7 bohr.
Las densidades de probabilidad que se obtienen elevando al cuadrado las funciones de onda, tienen un
aspecto algo diferente, aunque varían desde el carbono al flúor de la misma manera como se puede observar
en las gráficas de las figuras 9 y 10.
Fig. 10
Fig. 9
Si ahora se trabaja con las funciones de onda 3s, 3px y 3py, a fin de formar los orbitales híbridos de
elementos del mismo grupo pero del periodo siguiente, se deberá considerar la variación de la carga nuclear
efectiva, que para el Si es 4,28, para el P, 4,886, para el S, 5,48 y para el Cl, 6,116.
Se presentan las curvas de nivel del primero y último para observar su variación; corresponden a las figuras
11 y 12. En estos orbitales los dos lóbulos están mucho menos diferenciados que en los anteriores.
Fig. 11 Fig. 12
Ahora se va a establecer la comparación total entre los diferentes orbitales híbridos spX, y su variación,
haciendo un resumen de los tres trabajos presentados, tomando como referencia los híbridos en el C
Fig. 14 Fig. 15
Fig. 13
Se aprecia perfectamente que tanto el lóbulo no enlazante como el enlazante van aumentanto en superficie,
especialmente el negativo, por ello la forma externa se va modificando. Los dos lóbulos están mucho mas
diferenciados en el sp, que en el sp3.
Orbitales moleculares I
Introducción a laTeoría de orbitales moleculares. Aplicación a moléculas homonucleares sencillas.
La teoría de orbitales moleculares (OM) para describir un enlace covalente nace unida a la de
orbitales atómicos (OA), y con su misma filosofía. Es elaborada por Hund, Lennart-Jones y
Mulliken, entre 1928 y 1932.
La función de onda que representa al OM, a partir de la cual se determina su energía, se puede
calcular de muchas maneras pero la más general se obtiene por combinación lineal de sus orbitales
atómicos (CLOA), o sea por suma o resta de las funciones atómicas (idea de Mulliken), dado que al
aproximarse los OA de un átomo X a los de otro igual y superponerse, los electrones están controlados
por la interacción nuclear de uno u otro núcleo la mayoría del tiempo. La suma produce los OM
enlazantes (OME), mientras que la resta da lugar a los OM antienlazantes (OMAE). Ahora bien para
realizarla deberán cumplirse los requisitos inherentes a toda la teoría cuántica del enlace covalente:
superposición de los OA de energía y simetría semejante, así como los característicos de la teoría
atómica: Principio de Exclusión de Pauli y leyes de Hund.
Los OA, venían caracterizados por dos números cuánticos n y l, que determinaba su forma ( l=0,
OA s, l=1 OA p, l=2 OA d). En los OM, se conservan los números cuánticos n y l que corresponden a las
funciones atómicas que se combinan, incrementado en el número cuántico λ característico del OM, que
cuantifica el momento angular molecular respecto al eje de enlace de la molécula, y que también
determina la forma, debiendo tomar valores equivalentes al módulo de los que tomaba m como se aprecia
en la tabla siguiente. La comparación entre los números cuánticos atómicos de los electrones de X y
moleculares de X2, implica:
X X2
n l m OA Nºe λ OM NºeE NºeAE
1 0 0 1s 2 0 σ 2 2
2 0 0 2s 2 0 σ 2 2
1 1,0,-1 2p 6 0 σ 2 2
1 π 4 4
3 0 0 3s 2 0 σ 2 2
1 1,0,-1 3p 6 0 σ 2 2
1 π 4 4
2 2,1,0,-1,-2 3d 10 0 σ 2 2
1 π 4 4
2 δ 4 4
Tabla 2.2.3.
En el caso de la molécula de hidrógeno (H2), el
CLOA de las ψ individuales produciría la función de onda
molecular φ = aψA ± b ψB , considerando los hidrógenos
A y B, y a y b los coeficientes de participación de los
respectivos ψ atómicos en la función de onda molecular,
que en este caso serán iguales al serlo los OA. La suma
produce el OM enlazante que se denominaría σ1sE y la
resta, el antienlazante σ1sAE de mayor energía. En este caso
los dos electrones 1s que la forman cubrirían el OM de
menor energía, o sea el σ1sE que genéricamente se
denomina con las siglas inglesas de HOMO (orbital
molecular ocupado de mayor energía), mientras que el
σ1sAE quedaría desocupado de electrones, y sería el
LUMO (orbital molecular vacío de menor energía) (fig.1)
Fig.1
Los números cuánticos (n,l, λ, s) de los electrones en dicho OM, serían (1,0,0,½ ) y (1,0,0,-½ ), o sea
siguiendo el Principio de Exclusión de Pauli, tendría capacidad para dos electrones (general para cada
OM).
El nivel energético que proporcionaría el σ1sE sería de 337 kJ/mol, mejor que la obtenida a través de la
TEV (Teoría del enlace de valencia) de 303kJ/mol, pero todavía lejos de su valor experimental,
458kJ/mol.
3
⎛⎜ − (( x − x1)2 *+ y* y ) ⎞⎟
Zr
⎛ 1 ⎞ ⎛ Z ⎞ 2 − a0
Conocemos que para el átomo de hidrógeno ψ 1s = ⎜ ⎟ ⎜ ⎟
* * e = 0,56 * e ⎝ ⎠
⎝ π ⎠ ⎝ a0 ⎠
Teniendo en cuenta que Z ( carga nuclear efectiva ) =1
Y que el radio covalente experiental del H2 es 37pm, y que por lo tanto la distancia entre las
posciones de los núcleos de H deberá ser 74pm, en función del radio de Bohr (53 pm) será 1,49 bohr. Por
lo tanto las funciones de los OM a representar serán las combinaciones lineales :
⎛ −⎛⎜ (( x + 0 , 7 )2 + y 2 ) ⎞⎟ ⎞
0,56⎜⎜ e −( (( x −0,7 ) + y ) ) ± e ⎝
2
⎟⎟
2
⎠
⎝ ⎠
Las gráficas obtenidas de los OM y las curvas de nivel, serán fig 2,3,4 y 5, referidas como
siempre a distancias bohr
Fig.2 Fig.3
Fig.4 Fig.5
Si elevamos las funciones al cuadrado para determinar las densidades de probabilidad se obtienen
las gráficas respectivas que copmo se aprecia difieren poco de la descripción de los OM ( fig 6, 7, 8 y 9).
Fig.6
Fig.7
Fig.8 Fig.9
Tópicos de la descripción de OM
Es muy curiosa la descripción gráfica de los orbitales moleculares (OM) en los textos en
Enseñanzas Medias. Al aplicar el método CLOA (combinación lineal de orbitales atómicos), siempre se
toman los orbitales atómicos hidrogenoides, aunque la molécula de hidrógeno no excitada no forme
enlaces con combinaciones lineales superiores al 1s. Sin embargo en los dibujos apreciamos la forma de
los orbitales moleculares F2s, F2p e incluso los B2p propocionados por las combinaciones lineales de los
OA del hidrógeno que son sólo teóricas e inexistentes en la práctica.
Este es el primer tópico que surge en los textos, el segundo es creer que los OM moleculares del
mismo tipo tienen la misma forma o son aparentemente iguales.
En este trabajo vamos a comparar los OM, formados por combinaciones lineales de átomos no
hidrogenoides, comparados con los hidrogenoides ( que aparecen en los textos).
Si aplicamos todo lo anterior a la posible molécula de Li2, la función de onda átomica generalizada
3
Z 'r
⎛ 1 ⎞ ⎛ Z ' ⎞ 2 − a0
ψ 1s = ⎜⎜ ⎟⎟ * ⎜ ⎟ * e =
⎝ π ⎠ ⎝ a0 ⎠
vendría modificada no sólo por los valores de r ( la distancia de separación aumenta) sino sobre todo por
los de de Z’ (carga nuclear efectiva), con lo cual el segundo factor ya no es 1 , como en el caso del H
Con lo que la combinación lineal para el Li2, de orbitales 1s, proporciona una función
⎛ −1, 3*⎛⎜ (( x + 2 , 53 )2 + y 2 ) ⎞⎟ ⎞
0,56 * 1,48⎜⎜ e −1,3*( (x − 2,53 + y ) ) ± e ⎝ ⎠⎟
2 2
⎟
⎝ ⎠
Si la elevamos al cuadrado para calcular la densidad de probabilidad y compararla con la del H, para estos
mismos OA, se obtendrían las gráficas siguientes
Fig.11
Fig.10
Fig.12 Fig.13
¿Qué significa esto?. Pues sencillamente que al fijar los núcleos de los átomos de Li, a la distancia
internuclear experimental, los OA 1s, no superponen, y por lo tanto no se puede hacer la combinación
lineal. Por eso al no existir esta, no hay enlace entre ellos; son OM no enlazantes. Si creemos que a lo
mejor, al aumentar la carga nuclear efectiva y disminuir la distancia de enlace, avanzando hacia la
derecha en el sistema periódico a lo mejor conseguimos la combinación lineal, estamos equivocados,
como puede verse al aplicar la teoría a la combinación lineal de las funciones 1s para la molécula de
oxígeno ( fig.14)
Fig.14
Esto ocurriría con todos los elementos del segundo período y por lo tanto con los demás.
Orbitales moleculares II
Tópicos en la descripción de OM
En este trabajo vamos a comparar los OM F2s, formados por combinaciones lineales de átomos no
hidrogenoides comparados con los hidrogenoides (que aparecen en los textos).
3 ⎛
⎜− (( x − x1) + ( y* y )) ⎞⎟
2
((x − x1) )
Z 'r
⎛ ⎞ ⎛ Z' ⎞ ⎛ Z'r ⎞ − ⎜⎜ ⎟⎟
= 0,1 * ⎛⎜ 2 − + ( y * y ) ⎞⎟ * e
1 2 2
ψ 2 s = ⎜⎜ ⎟⎟ * ⎜⎜ ⎟⎟ * ⎜⎜ 2 − ⎟*e ⎝ ⎠
2 a0 2
⎝ 32π ⎠ ⎝ a 0 ⎠ ⎝ a 0 ⎟⎠ ⎝ ⎠
Por lo que aplicándolo a dos átomos de hidrógeno cuyas funciones se sumarían o restarían y
tomando una distancia de 1,5 radios de bohr, sería:
⎛ (( x −0,75 ) + ( y* y )) ⎞⎟
2
La combinación lineal produciría los OM cuyas gráficas y curvas de nivel serían los dados en las
fig 1 y 2, en combinaciones enlazantes y 3 y 4, en antienlazantes.
Fig.1 Fig.2
Fig.3 Fig.4
Sin embargo esta misma combinación lieal, si cambiamos las condiciones de contorno,
extendiéndolas hasta los 10 bohr, como se hará con la mayoría de los elementos del segundo período,
producirían unos gráficos muy diferentes como se vé en las fig 5 y 6.
Fig.5
Fig.6
Fig.9
Como se observa los límites de la densidad de probabilidad son inferiores a los del OM
Al trasladarnos hacia la derecha a través del sistema periódico a otros elementos del segundo período,
desde el litio hasta los correspondiente del oxígeno ( fig 10 y 11) y del flúor ( fig 12 y 13) percibiremos
en los gráficos delas curvas de nivel, la disminución del tamaño del OM F2s,. Esta disminución, va de
unos límites de "10 bohr en el F2s del Li2 hasta los " 3 bohr para el correspondiente del F2.
Fig.10 Fig.11
Fig.13
Fig.12
Como es lógico esta disminución del tamaño todavía es mas acusada en las curvas de densidad de
probabilidad, para OM F2s del mismo tipo, como las de la hipotética molécula de carbono (C2), y la de O2
( fig 14 y 15)
Fig.14 Fig.15
Orbitales moleculares III
Tópicos de la descripción de OM F 2p
1.El modelo hidrogenoide
Si se aplica el sistema empleado en OM1 y OM2, veremos primero los OM hidrogenoides, que se
representarán a partir de la combinación lineal de OA 2px. Para ello se tomará la función atómica
5 5
Z 'r Z 'r
⎛ 1 ⎞ ⎛ Z' ⎞2 − ⎛ Z' ⎞2 −
ψ 2 px = ⎜ ⎟ * ⎜⎜ ⎟⎟ * xe 0 = 0.1 * ⎜⎜ ⎟⎟ * xe 0 (1)
2a 2a
⎝ 32π ⎠ ⎝ a 0 ⎠ ⎝ a0 ⎠
Aplicándola al H; Z’=1, a0=1, y la posición de los núcleos sobre el eje X, en unidades bohr en x+0,75 y
x-0,75.
Obteniéndose las dos funciones enlazante y antienlazante:
( ( x − 0 , 75 )*( x − 0 , 75 ) − y* y ) ( ( x + 0 , 75 )*( x + 0 , 75 ) − y* y )
− −
Φ σ 2 px = 0.1 * ( x − 0,75) * e 2
± 0.1 * ( x + 0,75) * e 2
OM SIGMA 2px
OM SIGMA 2px
ENLAZANTE ANTIENLAZANTE
PROYECCIÓN SOBRE
EL PLANO XY
PLANOS
NODALES
Fig.1 Fig. 2
En las proyecciones sobre el plano XY, se aprecian perfectamente en el fig 1, los dos planos nodales
procedente de los OA 2px, que hereda el OM. En cambio en el caso del OM antienlazante (fig.2). se
produce el plano nodal molecular y los dos atómicos, o sea que existen tres planos nodales.
2. Moléculas homonucleares de átomos no hidrogenoides
Si en la ecuación 1, se introducen para cada elemento los parámetros de Z’, y las distancias
correspondientes al radio covalente de cada elemento, para una hipotética molécula Li2, tendríamos unas
funciones moleculares:
⎛ 1, 32* (( x − 2 , 53 )2 + y 2 ) ⎞ ⎛ 1, 3* (( x + 2 , 53 )2 + y 2 ) ⎞
⎡ − ⎜⎜ ⎟ −⎜⎜ ⎟⎤
⎢ ⎟ ⎟⎟ ⎥
φ ' = 0,19 * ⎢( x − 2,53).e ± ( x + 2,53).e
⎜ 2 ⎟ ⎜ 2
⎝ ⎠ ⎝ ⎠
⎥
⎢⎣ ⎥⎦
Fig.3 Fig.4
Fig.6
Fig.5
Fig.7 Fig.8
Fig. 9 Fig.10
Se puede comparar el tamaño del OM desde el Li2 al N2, relacionando las figuras 6 y 10, y 4 y 8.
Casi en el extremo del período , el O2, presenta unas funciones de onda molecular F2px
Fig.13 Fig.14
Orbitales moleculares IV
5 5
Z 'r Z 'r
⎛ 1 ⎞ ⎛ Z' ⎞2 − ⎛ Z' ⎞2 −
ψ 2 py =⎜ ⎟ * ⎜⎜ ⎟⎟ * ye 0 = 0.1 * ⎜⎜ ⎟⎟ * ye 2 a0
2a
⎝ 32π ⎠ ⎝ a 0 ⎠ ⎝ a0 ⎠
Aplicada al hidrógeno Z’ carga nuclear efectiva =1, a0=1, como en los casos anteriores, tomará la
expresión anterior.
La CLOA proporciona estas dos funciones
⎛ (( x − 0 , 75. )2 + y 2 ) ⎞ ⎛ (( x + 0 , 75 )2 + y 2 ) ⎞
⎡ −⎜⎜ ⎟ −⎜⎜ ⎟⎤
⎢ ⎜ 2 ⎟
⎟ ⎜ 2 ⎟⎟ ⎥
φ = 0,1⎢ y * e ⎝ ⎠
± y * .e ⎝ ⎠
⎥
⎢ ⎥
⎣ ⎦
Representadas con un programa de gráficos superficiales nos da la combinación enlazante (fig.1) y la
antienlazante (fig.2). Se aprecian perfectamente en los dibujos los planos nodales procedentes de los OA
p y los planos nodales generados por la combinación lineal antienlazante
OM 2PI py OM 2PI py
ENLAZANTE ANTIENLAZANTE
PROYECCIÓN SOBRE
EL PLANO XY PROYECCIÓN SOBRE
EL PLANO XY
Fig.1 Fig.2
Si ahora aplicamos a los OM de otros átomos del segundo período, por ejemplo, para el carbono, la
función OM tomará, al realizar la combinación lineal, teniendo en cuenta que Z’=3,25, y la distancia
media interatómica es 2,92 bohr, la siguiente expresión:
⎛ 3, 25* (( x −1, 46 )2 + y 2 ) ⎞ ⎛ 3, 25* (( x +1, 46 )2 + y 2 ) ⎞
⎡ − ⎜⎜ ⎟ − ⎜⎜ ⎟⎤
⎢ ⎜ 2 ⎟
⎟ ⎜ 2 ⎟⎟ ⎥
φ = 1,9 ⎢ y * e ⎝ ⎠
± y * .e ⎝ ⎠
⎥
⎢ ⎥
⎣ ⎦
Fig.3 Fig.4
Fig.5 Fig.6
Si queremos comparar estos OM, con los del mismo tipo para el Si2, especie molecular del
“elemento hermano” en el tercer periodo, se deberá emplear la función:
5 5
⎛ 2 ⎞ ⎛ Z ⎞2 ⎛ Zr ⎞ ⎛ Z ⎞2 ⎛ Zr ⎞
Zr Zr
− −
ψ 3 py = ⎜⎜ ⎟ * ⎜ ⎟ * ⎜ 6 − ⎟ * y * e 3 a0
⎟ ⎜a ⎟ ⎜ = 0.01* ⎜⎜ ⎟⎟ * ⎜⎜ 6 − ⎟⎟ * y * e 3a0
⎝ 81 π ⎠ ⎝ 0⎠ ⎝ a0 ⎟⎠ ⎝ a0 ⎠ ⎝ a0 ⎠
Fig.7
Fig.10
Fig.9
Comparando la fig 4 y la 8, por un lado, y la 5 con la 10, se pueden establecer las diferencias entre los
OM de moléculas diatómicas de elementos del mismo grupo pero de distinto período. Se observa que
tanto en la combinación enlazante como en la antienlazante, aparece una superficie nodal interior,
heredada de la que surge en los OA respectivos.
Si pasamos a otros elementos como el nitrógeno y el oxígeno, avanzando hacia la derecha del sistema
periódico (deberá disminuir el radio y volumen atómico), veremos como influye este hecho en el tamaño
del OM. Así las curvas de nivel de los OM antienlazantes y enlazantes, se representan en las figuras
11, 12 , 13 y 14
Fig.11 Fig.12
Fig.13 Fig.14
El efecto de la modificación del radio atómico tiene menos incidencia al avanzar hacia la derecha en el
período 3, tal como se puede apreciar al considerar el OM enlazante del azufre diatómico, con su curva de
nivel (fig. 15 y 16) y compararlo con el del silicio diatómico (fig 9 y 10).
Fig.16
Fig.15
48
AMPLIACIÓN TEMÁTICA
NÚCLEO ATÓMICO
Introducción.
A partir del descubrimiento de los rayos X por Röntgen, en 1895 la
búsqueda de rayos con similares características, propició el de la radiac-
tividad natural, por Becquerel y el matrimonio Curie, en 1889.
De los átomos salían tres tipos de radiaciones, con diferentes
características:
La α con carga positiva, masa similar al He y velocidad media de 2.106ms-1
La β con carga negativa, masa similar al electrón y velocidad media de
2,9.108ms-1, muy próxima a la velocidad de la luz.
La γ sin carga ni masa, con las propiedades de una radiación
electromagnética de frecuencia elevada y por lo tanto de gran energía.
Fig 1.23. Reglas radiactivas Cuando 7 años más tarde Rutherford, descubre el núcleo de los
átomos, donde se concentraba la masa y la carga positiva, ya se supo el
lugar de donde procedían las radiaciones; el núcleo atómico, e incluso se
En 1913 Soddy llamará isótopos a llegó a la conclusión de la existencia de la partícula neutra, formada por un
los núcleos de átomos del mismo protón y un electrón fuertemente unidos, que justificaría la emisión β.
elemento con diferente número
másico, como ocurre en la tercera El descubrimiento del protón en 1914 y del neutrón 18 años
regla.
después permitieron averiguar la composición nuclear e interpretar el origen
de las radiactividad.
60
48
49
La emisión de radiación γ que acompaña a todos los procesos anteriores, se Isótopo Aplicaciones
debe a la reorganización de protones y neutrones dentro del núcleo , que en 226
Ra Pararrayos
un modelo nuclear en capas (similar al electrónico) promoverían saltos entre 192
Ir Industria, gammagrafia
niveles con emisión de energía en forma de rayos γ. 133
Ba Medida de espesores
85
Kr Ionizadores
137
Radiactividad artificial. Cs Medida de densidades
60
El matrimonio Federico Joliot, e Irene Curie, hija de la descubridora de la Co Medida de niveles
82
radiactividad natural, descubrió en 1934, la radiactividad artificial que será Br Transporte de fluidos
60
fundamental para la creación de núcleos radiactivos “a la carta” o Co Obstrucciones en tuberías
137
Cs Irradiación de alimentos
radioisótopos, con innumerables aplicaciones médicas e tecnológicas. Ver 14
C Datación orgánica
tabla 18
O Fugas en embalses
131
I Medicina, tiroides
Parece evidente que un núcleo estable si se le bombardea en un reacción 198
Medicina, hígado
Au
nuclear con protones a gran velocidad, núcleos de deuterio 2H1, partículas 197
Medicina, riñón
Hg
alfa o incluso neutrones, se le puede sacar de la franja de estabilidad y por 51
Cr Medicina, intestino
lo tanto, artificialmente convertirlo en un emisor de radiación. Eso fue lo que 75
Se Medicina, páncreas
hicieron los Jolito-Curie con el aluminio como blanco de partículas α, hasta 32
P Medicina, sangre, betaterapia
convertirlo en un silicio inestable, que a su vez se convirtió en fósforo. 64
Cu Medicina, cerebro
47
Ca Medicina, huesos
99
Las radiaciones emitidas debido a la gran sensibilidad de los detectores de Tc Medicina, corazón
60
radiactividad (contadores de Geiger, cámaras de niebla etc), se pueden Co Bomba de cobalto, cáncer
localizar fácilmente lo que aumenta su aplicabilidad.
2
denominó L, en función de su apellido. Estas durezas dependían del peso atómico de
cada elemento irradiado, y por lo tanto podían servir de criterio sistematizador. Su
trabajo recibirá como premio, el Nobel de Física de 1917. Sin embargo continuaba
sin conocerse su naturaleza que unos postulaban corpuscular como los Bragg, en
Adelaida (Australia), y otros ondulatoria.
Hasta que en 1912, Von Laue del Instituto de Física teórica de Munich, con
En 1899, Haga y Win ya
la colaboración de Friedrich y Knipping, consiguen difractarlos; eran de naturaleza habían observado el
ondulatoria. La muestra empleada era de vitriolo de cobre (sulfato cúprico), y la ensanchamiento de un
imagen obtenida después de 12 horas de irradiación, fue bastante difusa. Sin haz de rayos X al pasar a
embargo usando como blanco la blenda de cinc (sulfuro de cinc), los puntos través de rendijas de
producidos eran mucho más simétricos y perfectamente ordenados: era el primer milésimas de milímetro
“lauegrama”. El 8 de junio de 1912, el profesor Sommerfeld, director del Instituto (había una difracción)
comunica en la Academia de Ciencias de Baviera, el descubrimiento de Von Laue y
colaboradores que reciben el Nobel de Física de 1914.
¿Cómo se explican los lauegramas?.
La explicación y justificación la van a realizar a partir de 1912, William Bragg y su
hijo Lawrence. Estos rayos tenían una longitud de onda muy pequeña, y por lo tanto
eran muy energéticos, y se difractaban a través de la materia, dado que su longitud
de onda era del mismo orden que el tamaño de los cristales, produciéndose una
T1.2. Difracción de rayos X dispersión al atravesar el cristal (iones o átomos como puntos negros en fig T1.2)).
3
Ley de Moseley
En 1913, unos meses más tarde que Bohr presentara sus primeras teorías
Moseley murió a los 27 años, atómicas, Moseley, que también había sido discípulo de Rutherford, imbuido por las
poco después de enunciar su ideas de Barkla, que relacionó la energía y frecuencia de la radiación X con el peso
ley, en 1915, durante la atómico del elemento que actuaba como blanco y que la emitía, atribuye el hecho a
primera guerra mundial , en que al aumentar el peso atómico también lo hacía el número de cargas positivas del
el sitio de Gallípoli, luchando núcleo, que denomina número atómico Z, y que relaciona en la radiación K
contra los turcos específica de él, con la frecuencia, obteniendo una gráfica que proporciona la
Fig.T1.3. Ley de Mosele
fórmula denominada actualmente ley de Moseley:
14
10
0
40 45 50 55 60 65 70 75 80
Z
ν
R 3
La pendiente de la gráfica = = 0,866 (2)
Z 2
Siendo R la constante de Rydberg, </c el número de onda de la radiación y
Z el número atómico. Esta ley era experimental y sólo fue explicada mas tarde, y se
hizo con los elementos descubiertos hasta aquella época. Incluso faltaba escandio,
pues no logró separarlo para emplearlo como anticátodo.
4
La interpretación de la ley de Moseley, fue realizada por J.J.Thomson, el
descubridor del electrón, quien sugirió que la excitación de electrones interiores, por
colisión con los del haz externo (rayos catódicos),producían vacantes que eran
cubiertas por los de niveles superiores, que al caer a los inferiores (K, L, etc),
producían la emisión (fig.T1.4). Por lo tanto, el espectro de rayos X se convertía en
la radiografía de los niveles energéticos del átomo que actuaba como blanco. Y lo
mas importante, la frecuencia emitida dependía del cuadrado del número atómico Z
del elemento. Por eso conociendo experimentalmente la frecuencia, se podía calcular
Z, y ver si correspondía a un elemento conocido o a uno nuevo ( así se descubrieron
algunos elementos químicos y se hicieron algunos falsos descubrimientos)
despejando
ν
c
=ν =
hc
( )
2,18.10 −18 Z 2 µ 3
4
(4)
1 =ν ⎛ 2π 2 k 2 me 4 ⎞⎛ 1 1 ⎞ ⎛ 1 1 ⎞
=⎜ ⎟⎟ ⎜⎜ 2 − 2 ⎟⎟ = R⎜⎜ 2 − 2 ⎟⎟ (6)
λ c ⎜ h 3c
⎝ ⎠ ⎝ nf ni ⎠ ⎝ nf ni ⎠
pues bien esta fórmula que es la de Balmer también se convierte en la de Moseley
con las consideraciones hechas antes: Espectro de RX producido por un tránsito del
nivel L(n=2), al K(n=1). En el H estos tránsitos correspondía a la serie de Lymen
con frecuencias en el UL, en elementos de Z mayor, la energía desprendida era
mayor y la frecuencia correspondía a los rayos X. con lo cual esta ley sirvió para
comprobar la teoría de Bohr.
Actualmente la ley de Moseley se expresa: 1/8=νK /c = AR(Z-F)2 (7)
siendo Z-F, la carga nuclear efectiva, F la constante de pantalla, y A el factor
correspondiente al salto electrónico
5
Espectroscopía: Determinación experimental de constante de
Rydberg.
Se necesita para esta serie de experimentos que sólo se pueden hacer en el
laboratorio, tubos de descarga de mercurio y de hidrógeno, una fuente de 5000V o
un carrete de Rumhkkorff, una red de difracción de 600 rayas/mm, cinta métrica y
soportes.
El primer paso será conocer la constante de la red de difracción que
tenemos, conociendo la longitud de onda 8 de una de las rayas que apreciamos al
observar el espectro del mercurio con la red de difracción, en este caso la verde
manzana de 546,074 nm. Al aplicar la ley de Bragg para la difracción n8 =2dsen2,
para n=1, encontramos que 2d=8/sen2 . (EX1)
Para ello se monta tal como indica la figura, la lámpara de Hg, uniendo sus
terminales a la salida del carrete de Rumhkkorff (aprox. 5000V de ddp). Se apaga la
luz del laboratorio, y entre las múltiples rayas que se observan a la derecha del tubo,
se selecciona la verde, de longitud de onda conocida. Se superpone a ella un poste
para medir la distancia L y se enciende la luz.
Se miden las distancias D y L, y se calcula el sen2= L/(√(D2+L2), al
sustituir en la fórmula (EX1), se determina 2d, que en este caso da 1662nm.
D=1,61m L=0,56m 2d=1662nm
El valor dado por las tablas para este salto es de 656nm, por lo tanto el error que se
comete es despreciable, pese a que las medidas están hechas con cinta métrica, y la
observación puede llevar consigo un error accidental según la posición del ojo sobre
la red.
47