TÓPICOS ESPECIALES EN FÍSICA BFI04 B – 2020 1

Las funciones de onda del átomo de hidrógeno

Arnold S. Ore a, Jeanpiere Retamozo a, Steef M. Diaz a, Rafael B. Saldaña a, Gustavo J.
Muñoz a
a
Facultad de Ingeniería Química y Textil, Universidad Nacional de Ingeniería, Lima, Perú

INFORMACIÓN DEL
ARTÍCULO
RESUMEN
Historia del artículo El siguiente trabajo tiene como fin dar a conocer la
aplicación de la Ecuación de Schrödinger al átomo de
Hidrógeno. Explicaremos a detalle el desarrollo de la
Recibido 12 de Junio de 2020
ecuación de Schrödinger en tres dimensiones para luego por
sustitución de variables a coordenadas esféricas, debido a
Palabras Clave: que le átomo de hidrógeno se considera esférico, distribuir
Ecuación de Schrödinger las tres funciones: polar, radial y azimutal. Posteriormente,
tridimensional explicaremos el método para poder hallar la Probabilidad
Separación de variables Radial y Angular descritas por la combinación de los
Ecuación azimutal números cuánticos n, m, ml.
Ecuación polar
Ecuación radial
Probabilidad radial
ABSTRACT
Proabibilidad angular
The following work aims to publicize the application of the
Keywords: Schrödinger equation to the Hydrogen atom. We will
Three-dimensional Schrödinger explain in detail the development of the Schrödinger
equation equation in three dimensions and then by substituting
Variable separation variables to spherical coordinates, because the hydrogen
Azimuth equation atom is considered spherical, distributing the three
Polar equation
Radial equation functions: polar, radial and azimuth. Later, we will explain
Radial probability the method to be able to find the Radial and Angular
Angular chance Probability described by the combination of the quantum
numbers n, m, ml.

llamo principio de incertidumbre lo cual fue

1. Antecedentes fundamental para el desarrollo de la mecánica
cuántica.
Se sabe que la mecánica clásica no puede
explicar el comportamiento de los electrones del La mecánica cuántica hace el uso de números
átomo de Bohr, Heinsenberg se percató que no se complejos. Por ejemplo, la función de onda(Ψ)
podía determinar la posición y velocidad del que al elevarla al cuadrado nos da la probabilidad
electrón con su mecánica matricial lo cual decía de encontrar a la partícula en un cierto intervalo,
que es imposible conocer la posición y velocidad al ser una probabilidad nos da un número
del electrón al mismo tiempo, a este postulado lo positivo.
 TÓPICOS ESPECIALES EN FÍSICA BFI04 B – 2020 1

Por otro lado, tenemos al Hamiltoniano (Ĥ) que ℏ2 2

𝛨=− (𝑃 + 𝑃𝑦2 + 𝑃𝑧2 ) + 𝑈(𝑥,𝑦,𝑧)
es un operador de la mecánica cuántica, 2𝑚 𝑥
comparando con la mecánica clásica, tiene
operadores asociados con la energía cinética y
potencial. En seguida se sustituye las cantidades dinámicas
𝑃𝑥 , 𝑃𝑦 , 𝑃𝑧 y Η por sus operadores diferenciales.
Con el Hamiltoniano (Ĥ) y la función de
onda(Ψ), tenemos la ecuación de Schrödinger
independiente del tiempo
ℏ2 𝑑2 𝑑2 𝑑2
ĤΨ=EΨ − ( 2+ + ) + 𝑈(𝑥,𝑦,𝑧)
2𝑚 𝑑𝑥 𝑑𝑦 2 𝑑𝑧 2
𝑑
De aquí partiremos para hablar sobre la ecuación = 𝑖ℏ
𝑑𝑡
de Schrödinger en 3 dimensiones.
2. Ecuación de Schrödinger en tres Dimensiones
De acuerdo con lo explicado sobre la ecuación Debido que en primer lugar la función de onda
de Schrödinger hasta el momento; solo son para (ᴪ) es una función del tiempo y de las 3
problemas unidimensionales, haciendo la coordenadas espaciales (x, y, z) se interpreta de
comparación entre una partícula newtoniana que la siguiente manera:
se mueve a lo largo de la recta. El modelo de
línea recta es adecuado para algunas 𝜑 = 𝜑(𝑥,𝑦,𝑧)
aplicaciones; sin embargo, para entender la
estructura atómica, necesitamos una Operando cada termino sobre la función de
generalización tridimensional. onda, se obtiene la ecuación de Schrödinger para
el sistema:
Primeramente, la función de energía potencial
depende también de las 3 coordenadas, y se ℏ2 𝑑2 𝜑(𝑥,𝑦,𝑧,𝑡) 𝑑2 𝜑(𝑥,𝑦,𝑧,𝑡)
− ( +
puede escribir como: 2𝑚 𝑑𝑥 2 𝑑𝑦 2
2
𝑑 𝜑(𝑥,𝑦,𝑧,𝑡)
𝐾𝑒 2 + )
𝑈(𝑥,𝑦,𝑧) = − 𝑑𝑧 2
√𝑥 2 + 𝑦 2 + 𝑧 2 + 𝑈(𝑥,𝑦,𝑧) 𝜑(𝑥,𝑦,𝑧,𝑡)
𝑑𝜑(𝑥,𝑦,𝑧,𝑡)
= 𝑖ℏ
𝑑𝑡
Si la partícula se puede mover en tres
dimensiones, entonces el momento lineal de (Ecuación General Tridimensional de
dicha partícula tiene tres componentes Schrödinger)
(𝑃𝑥 , 𝑃𝑦 , 𝑃𝑧 ), y su energía cinética es:

𝑃𝑥2 𝑃𝑦2 𝑃𝑧2 2.1 Ecuación de Schrödinger para el átomo de

𝐾= + + Hidrógeno
2𝑚 2𝑚 2𝑚
Una vez descrito la versión tridimensional de
En el desarrollo de la ecuación de Schrödinger dicha ecuación y que la función energía
en el sistema de tres dimensiones se procederá potencial es esférica simétrica, tal que solo
con la expresión del hamiltoniano mecano depende de la distancia “r”, al origen de
clásico: coordenadas.
 TÓPICOS ESPECIALES EN FÍSICA BFI04 B – 2020 1

1 𝑒2 Entonces las coordenadas rectangulares están

𝑈(𝑟) = − relacionadas con las esféricas por:
4𝜋𝜀0 𝑟
𝑥 = 𝑟 sin 𝜃 cos ∅
Luego, en la sección anterior notamos que la
ecuación de onda independiente del tiempo de
𝑦 = 𝑟 sin 𝜃 sin ∅
una partícula sometida a un potencial 𝑈(𝑟) en
tres dimensiones.
𝑧 = 𝑟 cos 𝜃
2
ℏ
− ∇2 𝜑(𝑟) + 𝑈(𝑟) 𝜑(𝑟) = 𝐸𝜑(𝑟) Donde el operador laplaciano en coordenadas
2𝑚 esféricas es:
1 𝑑 2 𝑑 1 𝑑 𝑑
∇2 = 2
(𝑟 )+ 2 (sin 𝜃 )
Debido a que el potencial 𝑈(𝑟) tiene simetría 𝑟 𝑑𝑟 𝑑𝑟 𝑟 sin 𝜃 𝑑𝜃 𝑑𝜃
esférica, vamos a utilizar coordenadas polares 1 𝑑2
+ 2
esféricas. 𝑟 sin 𝜃 𝑑𝜑2

Luego de transformar al operador hamiltoniano

Donde:
de coordenadas rectangulares a esféricas, la
r, el radio vector ecuación de Schrödinger de estado estacionario
𝜃, el ángulo polar resulta ser:
∅, el ángulo azimutal

ℏ2 1 𝑑 2 𝑑𝜑 1 𝑑 𝑑𝜑
− [ 2 (𝑟 )+ 2 (sin 𝜃 )
2𝑚 𝑟 𝑑𝑟 𝑑𝑟 𝑟 sin 𝜃 𝑑𝜃 𝑑𝜃
1 𝑑2 𝜑
+ 2 ] + 𝑉(𝑟) 𝜑
𝑟 sin 𝜃 𝑑𝜑2
= 𝐸𝜑
(Ecuación de Schrödinger en 3 dimensiones del
átomo de Hidrógeno)

Figura 2.1.1 Coordenadas polares

esféricas r,𝜃, ∅ del electrón y sus
coordenadas cartesianas.
 TÓPICOS ESPECIALES EN FÍSICA BFI04 B – 2020 1

Separación de variables
Ahora si sustituimos la expresión anterior de la Ahora bien, tomando la ecuación de la parte
función de onda en la ecuación de Schrödinger, derecha e igualamos a la constate 𝑙(𝑙 + 1) y la
tenemos: otra constante 𝑚𝑙 aparece al aplicar el método de
separación a la ecuación de armónico esféricos.
𝑑 𝑑𝜃 1 𝑑 2∅
ℏ2 𝜃∅ 𝑑 2 𝑑𝑅 𝑑𝜃
(sin 𝜃 𝑑𝜃) + 𝑠𝑖𝑛 2 𝜃∅ 𝑑∅
+ 𝑙(𝑙 + 1) = 0 … (2)
− [ (𝑟 )
2𝑚 𝑟 2 𝑑𝑟 𝑑𝑟
𝑅∅ 𝑑 𝑑𝜃
+ 2 (sin 𝜃 ) 𝑑2 ∅
𝑟 sin 𝜃 𝑑𝜃 𝑑𝜃
𝑑∅2
= −𝑚𝑙 2 , ∅ = 𝑎𝑒 𝑖𝑚𝑙∅ + 𝑏𝑒 −𝑖𝑚𝑙 ∅
𝑅𝜃 𝑑2 ∅
+ 2 2 ] + 𝑈(𝑟) 𝑅𝜃∅ , 𝐵(∅ + 2𝜋) = 𝐵(∅)
𝑟 𝑠𝑖𝑛 𝜃 𝑑∅2
= 𝐸𝑅𝜃∅ 𝑒 𝑖𝑚𝑙 ∅ = 𝑒 𝑖𝑚𝑙 ∅ . 𝑒 𝑖2𝜋𝑚𝑙 , 𝑚 𝜖 𝑧
Finalmente, la ecuación resulta:

En esta ecuación la derivada parcial de R con 1 𝑑 𝑑𝜃 𝑚𝑙 2

(sin 𝜃 ) + (𝑙(𝑙 + 1) − )𝜃 = 0
respecto a r se escribirá como la derivada total. sin 𝜃 𝑑𝜃 𝑑𝜃 𝑠𝑖𝑛2 𝜃
Luego multiplicando esta ecuación por 𝜑 =
2𝜇𝑟2
De la ecuación (2) ahora, podemos separarlos de
∶ manera similar de dos partes, una dependiente
ℏ2 𝑅∅𝜃
sólo del ángulo polar 𝜃 y la otra sólo del ángulo
azimutal ∅.Multiplicando por 𝑠𝑖𝑛2 𝜃 y
1 𝑑 2 𝑑𝑅 2𝜇
(𝑟 ) − 2 (𝑈(𝑟) − 𝐸) agrupando en términos:
𝑅 𝑑𝑟 𝑑𝑟 ℏ
1 ∅ 𝑑 𝑑𝜃
=− [ (sin 𝜃 )]
𝜃∅ sin 𝜃 𝑑𝜃 𝑑𝜃
𝜃 𝑑2 ∅ 1 𝑑2 ∅ sin 𝜃 𝑑 𝑑𝜃
+ 2
= 𝑙 (𝑙 + 1)𝑠𝑖𝑛2 𝜃 − (sin 𝜃 )
𝑠𝑖𝑛2 𝜃 𝑑∅2 ∅ 𝑑∅ 𝜃 𝑑𝜃 𝑑𝜃

Se observa que se crea una forma de ecuación en Debido a que cualquier solución de estas dos
la cual aquellos términos dependen de la ecuaciones llevará un valor correspondiente de 𝑙
coordenada r no depende ni de 𝜃 ni de ∅. como un parámetro, para fines convenientes
También vemos que la parte izquierda de esta elegimos una nueva constante 𝑚𝑙 2, por lo tanto,
ecuación no depende de las variables de la parte una de las ecuaciones resultantes será:
derecha, esto quiere decir que la única forma de 𝑑2 ∅
satisfacer la igualdad es que ambas partes sean 𝑑∅2
+ 𝑚𝑙 2 ∅ = 0 … (3)
igual a una constante, por lo tanto, dicha
La ecuación de Schrödinger para el átomo de
constante de separación debe ser 𝑙(𝑙 + 1) de
hidrogeno idealizado ha sido separado en tres
modo que se obtiene 2 ecuaciones, estas son:
ecuaciones, cada una dependiente de una sola
coordenada.
ℏ2 𝑑 𝑑𝑅 ℏ𝑙(𝑙+1)
− 2𝜇𝑟2 𝑑𝑟 (𝑟 2 𝑑𝑟 ) + ( 2𝜇𝑟2
+𝑈(𝑟) ) 𝑅 = 𝐸𝑅 … (1)
 TÓPICOS ESPECIALES EN FÍSICA BFI04 B – 2020 1

3. Ecuación Azimutal De esta nomenclatura si los valores de ml en la

En la separación de la ecuación de Schrödinger ecuación 2 se elevan al cuadrado, son los valores
para el átomo de hidrogeno la ecuación azimutal característicos propios de las varias funciones
es: características de la ecuación 1, como también
son las soluciones permitidas de la ecuación
𝑑2 ∅ 𝑑2 ∅
azimutal (𝑑∅2 + 𝑚𝑙 2 ∅ = 0 ). La constante ml es
+ 𝑚𝑙 2 ∅ = 0
𝑑∅2
el primer ejemplo de un numero cuántico en el
Describe el comportamiento de la función onda a tema del átomo de hidrógeno. Posteriormente se
medida que el sistema gira alrededor del eje z. su le llamara número cuántico magnético.
forma es similar al de la ecuación del oscilador
4. Ecuación polar
armónico, ya que tiene entre sus soluciones dos
funciones periódicas reales y una compleja. La ecuación que describe el comportamiento
angular de la función de onda alrededor del
∅ = 𝐴𝑠𝑒𝑛(𝑚𝑙∅) electrón es la ecuación polar:
∅ = 𝐴𝑐𝑜𝑠(𝑚𝑙∅) 𝑚𝑙 2 1 𝑑 𝑑∅
− (sen 𝜃 ) = 𝑙(𝑙 + 1)
sen2 𝜃 𝜃 sen 𝜃 𝑑𝜃 𝑑𝜃
∅ = 𝐴𝑒 𝑖𝑚𝑙∅
Es bien conocida en el estudio de los armónicos,
Si el electrón es girado una vuelta completa con una clase de soluciones calificadas como
alrededor del eje z, esta alcanza el mismo valor polinomios asociados de Legendre representado
para la solución ø cuando se alcance de nuevo el por 𝑃𝑙𝑚𝑙 (cos 𝜃) .Estas son funciones del cos 𝜃 y
mismo ángulo ø.
dependen de dos constantes paramétricas 𝑚𝑙 y 𝑙.
F(Φ)=F(Φ+2nπ) Ya que 𝑚𝑙 solo puede tomar valores enteros
positivos, negativos y cero, de las propiedades
conocidas de los polinomios de Legendre las
soluciones dadas existen solamente cuando 𝑙 es
entero y cuando es mayor o igual al valor
absoluto de 𝑚𝑙 .
5. Ecuación Radial
Esta ecuación especifica el comportamiento de la
función de onda en su dependencia de la distancia
al protón, la ecuación es:
𝑑 𝑑𝑅 2𝑚𝑟2 𝑒2
(𝑟 2 )+ (𝐸 + ) = 𝑙 (𝑙 + 1)𝑅
𝑑𝑟 𝑑𝑟 ℎ2 4𝜋𝜖𝑟

El conjunto resultante de soluciones a la ecuación

Figura 3.1 Giro del electrón alrededor del eje z de onda radial, las funciones características
radiales, están dadas por: 𝑅𝑛,𝑙 = 𝑒 −𝑛𝑟 𝑟 𝑙 𝐿𝑛,𝑙 (𝑟).
Se cumple en las funciones de la ecuación
siempre y cuando ml ø sea un múltiplo entero de Donde n es cualquier numero entero positivo
2𝜋. Incluye también al cero como rotaciones en distinto de cero y l es el numero cuántico orbital
direcciones opuestas, los valores permitidos de de las otras ecuaciones, las soluciones existen
ml, son: solamente cuando n 𝑛 ≥ (𝑙 + 1)
ml= 0, ±1, ±2, ±3…
 TÓPICOS ESPECIALES EN FÍSICA BFI04 B – 2020 1

6. La función de onda completa Cada combinación permitida de estos tres

Del método de la separación de variables, números define un estado del átomo y cada uno
logramos ver que la función de onda en de estos tiene una función característica
coordenadas esféricas está representada por: correspondiente.

Ψ=Ψ(r, 𝜃, Φ) Definiendo la probabilidad por unidad de

Es decir, la función de onda es dependiente volumen o densidad de probabilidad como:
ahora de 𝜃, Φ y r, si suponemos que esta 𝑑𝑝 ∗
= 𝜓𝑛,𝑙,𝑚𝑙 𝜓𝑛,𝑙,𝑚𝑙
ecuación tiene soluciones expresables como 𝑑𝑣
Si tomamos un diferencial de volumen tal que:
producto de funciones de las tres coordenadas, o
𝑑𝑣 = 𝑟 2 𝑠𝑒𝑛𝛳 𝑑𝑟𝑑𝛳𝑑∅
sea:
pasando a multiplicar al otro término
Ψ(r, 𝜃, Φ) = R(r)𝜃(𝜃)Φ(∅) … (1)
obtenemos:
Normalizando cada una de las tres ecuaciones
∗
𝑑𝑝 = 𝜓𝑛,𝑙,𝑚𝑙 𝜓𝑛,𝑙,𝑚𝑙 𝑑𝑣
componentes, en el intervalo correspondientes a
cada una de ellas, y multiplicándolas como en la
ecuación (1) obtenemos la ecuación de onda Pero de la función de onda completa sabemos

completa, se obtiene: que:

Ψ(r, 𝜃, Φ) = 𝑅𝑛,𝑙 𝑌𝑚,𝑙 Ψ(n,l,ml) = 𝑅𝑛,𝑙 𝑌𝑚,𝑙 … función radial

Donde: 𝑅𝑛,𝑙 : función radial Ψ*(n,l,ml) = 𝑅 ∗𝑛,𝑙 𝑌 ∗𝑚,𝑙 … función angular

(armónicos esféricos).
𝑌𝑚,𝑙 : función angular
Por lo tanto, operando y agrupando:
𝑑𝑝 = 𝑅 ∗𝑛,𝑙 𝑅𝑛,𝑙 𝑌 ∗𝑚,𝑙 𝑌𝑚,𝑙 𝑑𝑣
7. La función de onda del átomo de hidrógeno
La ecuación de Schrödinger para el átomo de 𝑑𝑝 = 𝑅 ∗𝑛,𝑙 𝑅𝑛,𝑙 𝑟 2 𝑑𝑟𝑌 ∗𝑚,𝑙 𝑌𝑚,𝑙 𝑠𝑒𝑛𝛳 𝑑𝛳𝑑∅

hidrógeno tiene como solución general: Donde notamos que:

Ψ(n, 𝑙, ml) = 𝑅𝑛,𝑙 𝜃𝑙,𝑚𝑙 Φ𝑚𝑙 dP (r): 𝑅 ∗𝑛,𝑙 𝑅𝑛,𝑙 𝑟 2 𝑑𝑟

Donde: 𝑅𝑛,𝑙 : 𝑝𝑎𝑟𝑡𝑒 𝑟𝑎𝑑𝑖𝑎𝑙 dP (ϴ, Φ): 𝑌 ∗𝑚,𝑙 𝑌𝑚,𝑙 𝑠𝑒𝑛𝛳

𝜃𝑙,𝑚𝑙 : 𝑝𝑎𝑟𝑡𝑒 𝑝𝑜𝑙𝑎𝑟 De la condición de normalización:

Φ𝑚𝑙 : 𝑝𝑎𝑟𝑡𝑒 𝑎𝑧𝑖𝑚𝑢𝑡𝑎𝑙 ∫ 𝜓 ∗ 𝜓 𝑑𝑣 = 1

Donde n, m y ml son respectivamente los
En el átomo de Hidrógeno:
números cuántico principal, orbital y magnético.
∫ 𝑅 ∗𝑛,𝑙 𝑅𝑛,𝑙 𝑟 2 𝑑𝑟𝑌 ∗𝑚,𝑙 𝑌𝑚,𝑙 𝑠𝑒𝑛𝛳𝑑𝛳𝑑∅ = 1
Estos toman valores como:
n = 1, 2, 3, … Desarrollando la integral:
l = 0, 1, 2, …, (n-1); n valores ∞ 𝛱 2𝛱
∫ 𝑟 2 𝑅 ∗𝑛,𝑙 𝑅𝑛,𝑙 𝑑𝑟 ∫ ∫ 𝑌 ∗𝑚,𝑙 𝑌𝑚,𝑙 𝑠𝑒𝑛𝛳𝑑𝛳𝑑∅ = 1
ml = 0, ±1, ±2, …, ±l; (2l+1) valores 0 0 0
 TÓPICOS ESPECIALES EN FÍSICA BFI04 B – 2020 1

Donde:
La probabilidad radial viene dada por:
∞
𝑃𝑛,𝑙 (𝑟) = ∫0 𝑟 2 𝑅 ∗𝑛,𝑙 𝑅𝑛,𝑙 𝑑𝑟
La probabilidad angular viene dada por: El cual da como resultado los siguientes
𝛱 2𝛱 polinomios:
𝑃𝑙,𝑚𝑙 (𝛳, ∅) = ∫ ∫ 𝑌 ∗𝑚,𝑙 𝑌𝑚,𝑙 𝑠𝑒𝑛𝛳𝑑𝛳𝑑∅
0 0

Tabla 7.1 Valores para las funciones de onda del átomo de

hidrógeno

Utilizando el polinomio para distintos valores de

n y l, siguiendo la normalización, se obtiene:

Tabla 8.1 Valores de funciones radiales normalizadas

para el átomo de hidrógeno

8. Probabilidad Radial
De la parte radial de la solución general se
obtiene:
𝑅𝑛,𝑙 = 𝑒 −𝑛𝑟 𝑟 𝑙 𝐿𝑛,𝑙 (𝑟)
Donde Ln,l son llamados los polinomios
asociados a Laguerre, el siguiente representado
por: Donde ao: radio de Bohr

Es la probabilidad total de encontrar el electrón

en un elemento diferencial de capa esférica, es
decir, a una distancia del núcleo entre r y r+dr, es
decir, dentro de un diferencial de volumen
esférico.
 TÓPICOS ESPECIALES EN FÍSICA BFI04 B – 2020 1

Figura 8.1 Caparazón esférico del átomo de

hidrógeno de radio r y diferencial dr.

El volumen será: dv = 4Πr2dr, por lo que la

probabilidad será igual a:

dP(r) = |𝑅𝑛,𝑙 |2 4Πr2dr

Se puede concluir lo siguiente:
• La probabilidad de encontrar al electrón
muy cerca al núcleo es pequeña, porque
r2 es pequeño.
• La distancia del electrón al núcleo en el 9. Probabilidad Angular
máximo de las curvas de distribución Es la probabilidad de encontrar al electrón en un
radial en las que solo hay un máximo (1s, parámetro delimitado por los parámetros
2p, 3d, …) es exactamente igual al radio angulares 𝜃, ϕ. Se denota de la siguiente manera:
de la correspondiente órbita calculada
por la teoría de Bohr (n2ao/Z).
 TÓPICOS ESPECIALES EN FÍSICA BFI04 B – 2020 1

10. Referencias
• YOUNG, HUGH D. y FREEDMAN, ROGER
A. (2013). Física Universitaria, 41, 1372-
1377.
• VIRGILIO A., CLYDE L. COWAN, B.J.
GRAHAM. (1999). Curso de Física
Moderna.México: Oxford University
Press.
• LEVINE, IRA N. (2001). QUÍMICA
CÚANTICA. 5𝑎 edición,
Madrid, PEARSON EDUCACIÓN, S. A.
Figura 9.1 Diagramas polares de la dependencia • RAYMOND A., SERWAY Y JOHN W.
direccional de las densidades de probabilidad JEWETT, Jr. (2006). Física para ciencias e
angular para el átomo con un electrón para l=3; ingeniería con Física Moderna, 42, 1224-
1230.
ml= 0, ±1, ±2, ±3.
• Alonso, M., FINN, Física, Addison-Wesley,
Tabla 9.1 Funciones de onda para el átomo de 1967, Vol. III
hidrógeno • http://hyperphysics.phy-
astr.gsu.edu/hbasees/quantum/sch3D.html
• Acosta, Virgilio, Cowan,B. J GRAHAM- Curso
de Física Moderna , pag. 213,214
• Acosta Virgilio, Curso de Física moderna,
pág. 239 – 243.