SIMULACIÓN DE LA RADIACIÓN DE CUERPO NEGRO - LEY DE WIEN

Willy W. K. Menacho N. ∗

∗ Carrera de Física, Universidad Mayor de San A n d r é s , La Paz, Bolivia

(e-mail: [email protected])

Resumen

Se hizo uso de un simulador de espectros de cuerpos negros (PHET) logrando ver las gráficas de cuerpo negro
y notando como va cambiando el espectro en relación a la temperatura. Descubrimos la relación que existe
entre la temperatura de los cuerpos negros y la longitud de onda máxima a su respectiva irradiación (ley de
Wien) determinando mediante una regresión lineal la constante de Wien la cual nos dio como resultado
b=0,00289824 ±4,11552E-07 [mK] con una discrepancia del 0.016% , realizamos un programa capaz de
integrar la función de Planck, verificando nuestras respuestas con las mediciones del simulador logrando una
discrepancia menor al 5%, también se logró la observación de las estrellas con 9602[K] para la estrella vega
y para la estrella Aldebarán 3910[K] de temperatura.

1. Introducción. Según Sears, ZemansKy (2005), un cuerpo negro es

una superficie ideal que absorbe todas las longitudes
El nombre de cuerpo negro obedece al hecho de que de onda de la radiación electromagnética que le
es un cuerpo que absorbe perfectamente toda la llegan, también es el mejor emisor posible de la
radiación que incide sobre ´el. Sin embargo, fue radiación electromagnética de cualquier longitud de
Kirchhoff quien sostuvo que un cuerpo que absorbe onda. La distribución espectral de la radiación
perfectamente energía electromagnética es un cuerpo procedente de un cuerpo negro se obtuvo
que también puede emitir energía electromagnética. experimentalmente y responde a una gráfica en forma
Con ´esta definición un cuerpo negro es aquel que de campana como la mostrada en la figura 1.
absorbe toda la radiación que le llega a todas las
longitudes de onda y la radiación que ´el emite es solo
función de la temperatura y de la frecuencia de la
onda. Obviamente no existe ningún objeto con tales
características, es decir, es una idealización como lo
es un gas ideal. Sin embargo sí existen cuerpos que se
aproximan a la definición del cuerpo negro.

El término de radiación se refiere a la emisión

continua de energía desde la superficie de cualquier
cuerpo, esta energía se denomina radiante y es
transportada por las ondas electromagnéticas que
viajan en el vacío a la velocidad de 3x108 [m/s].
Las ondas de radio, las radiaciones infrarrojas, la luz
visible, la luz ultravioleta, los rayos X y los rayos Figura 1. Distribución espectral de la radiación para
gamma, constituyen las distintas regiones del diferentes temperaturas. Fuente: Halliday D, Resnick
espectro electromagnético. R. Física Vol. 2. Tercera Edición. Compañía Editorial
Continental, S.A. México. (1996). Pág. 479.
En 1900 existían dos fórmulas, ninguna de las cuales
se podían ajustar a las curvas para la gama entera de Los osciladores, de origen electromagnético, pueden
longitudes de onda. La primera, debida originalmente tener solamente ciertos valores de energía discretos,
a Lord Rayleigh, deducida luego independientemente determinados por entero. En = nhv, donde n es un
por Einstein y modificada por James Jeans, se número.
desarrolló rigurosamente de su base clásica. Los osciladores pueden absorber o emitir energía en
Desafortunadamente falla por completo en el ajuste cantidades que son múltiplos discretos del cuanto
de las curvas, no pasando ni siquiera por un máximo, fundamental de energía el cual está dado por: ΔE = hv
ajustándose a las curvas bastante bien en la región de
las longitudes de onda muy largas, mientras que no se La segunda hipótesis de Planck establece que la
ajustaba a la realidad para la región de ondas cortas energía de los osciladores está cuantizada. La energía
(Eisberg R. 1974). de un oscilador de frecuencia sólo puede tener ciertos
valores que son 0, hv,2hv, 3hv, …,nhv. (Alonso M.,
Wilhelm Wien también dedujo una expresión teórica Finn E. 1971).
para la radiancia espectral. La Ley de Wien, obtenida En 1900 Max Planck, pudo deducir una función, que
a partir de los principios de la Física clásica, coincide hoy se llama ley de radiación de Planck, que
con los resultados experimentales sólo en la zona de concordaba muy bien con las curvas experimentales
las ondas cortas, pasa por un máximo, mientras que de distribución de intensidad de radiación para un
para la zona de las ondas largas estaba en cuerpo negro, la cual permite obtener la emitancia
contradicción con los resultados. espectral, que corresponde a la potencia por unidad de
área por unidad de longitud de onda, formulada por la
El problema antes mencionado acerca de las ecuación:
dificultades de la Ley de Wien y la fórmula de 2hc 2 f 3 1
Rayleigh y Jeans se denominó “Catástrofe I(f, T) = (1)
λ5 hc
Ultravioleta” y represento una etapa de crisis para la e λ k B T −1
Física de entonces. Max Planck buscando reconciliar
estas dos leyes de la radiación, realizó una Dónde: I(f, T) es la radiancia espectral de un cuerpo,
interpolación entre ellas que dio por resultado el i.e., la potencia por unidad de ángulo sólido y por
ajuste de los datos para todas las longitudes de onda. unidad de área normal a la propagación, de frecuencia
f de radiación por unidad de frecuencia en equilibrio
Planck era un experto en termodinámica y mecánica térmico a la temperatura T, h es la constante de
estadística, reconoció en 1900 que los Planck, k B es la constante de Boltzmann, c es la
descubrimientos empíricos sobre la distribución de velocidad de la luz en el vacío, f es la frecuencia de la
energía de la radiación negra, sólo se pueden entender radiación electromagnética y T es la temperatura
cuantitativamente cuando se introduce una de las absoluta del cuerpo. La relación entre la longitud de
hipótesis ajenas fundamentalmente a la física clásica: onda, λ la frecuencia, f y la rapidez es c= λf.
La hipótesis de los cuantos. Un sistema radiante no
puede intercambiar energía en porciones arbitrarias La longitud de onda λm en la que se produce el
con el campo de radiación, sino solamente en máximo de emisión viene dada por la ecuación (2),
múltiplos enteros del cuanto de energía hv, donde v toque describe la ley de Wien
es la frecuencia de la radiación y h una nueva
λmT = 2.898 ∗ 10−3 [mK] (2)
constante de la naturaleza h=6.6261*10-34[Js]. La
hipótesis cuántica se justifica por sus éxitos: ella
Y la potencia a la cual una superficie emite energía
proporciona exactamente la distribución de energía
radiante en todas las longitudes de onda, por unidad
medida de la radiación negra.
de área está dada por la ley de Stefan-Boltzmann:
Históricamente, el nacimiento de la Mecánica
Cuántica, se sitúa en el momento en el que Max W
Planck explica el mecanismo que hace que los átomos P(T) = 5.67 ∗ 10−8 T 4 [ ] (3)
m2 K4
radiantes produzcan la distribución de energía
observada. Max Planck sugirió en 1900 que:
En esta simulación se pretende analizar la relación 3. Resultados y discusiones.
que existe entre la temperatura emitida por los
cuerpos negros y la longitud máxima de su respectiva Se tomaron los datos del simulador PHET, los
radiancia espectral, corroborando la ley de Wien. datos obtenidos se encuentran en la tabla 1.

2. Materiales y método. Tabla 1. Datos obtenidos simulador PHET.

Fuente: elaboración propia.
Para este experimento necesitaremos: Simulador n λ [um ] T[K]
PHEThttps://phet.colorado.edu/sims/html/blackbody 1 0,5000 5800
-spectrum/latest/blackbody-spectrum_en.html. Una 2 0,4830 6000
computadora con el paquete office Excel, paquete
científico anaconda versión 2.7, la aplicación 3 0,4460 6500
Android Sky Map y un celular. 4 0,4140 7000
5 0,3860 7500
Nos familiarizamos con el uso del simulador 6 0,3690 7850
PHET sobre el Espectro de un Cuerpo Negro en el 7 0,3410 8500
cual tiene la Densidad de Potencia Espectral o
8 0,3220 9000
Radiancia Espectral, [MW/m2/m] en función de la
9 0,3050 9500
longitud de onda,  [m].
10 0,2910 9950
Una vez hecho esto pasaremos a tomar por lo 11 0,2900 10000
menos 20 datos del simulador para determinar la 12 0,2760 10500
relación que existe entre la longitud de onda máxima 13 0,2630 11000
de cada distribución y la temperatura del cuerpo 14 0,4170 6950
negro, con los cuales se realizaran el ajuste
correspondiente indicando las respectivas constantes 15 0,5570 5200
y sus errores, mediante una regresión lineal a partir de 16 0,5850 4950
la ecuación (2), de la forma: 17 0,9660 3000
18 1,2330 2350
λm = A + B/T (4) 19 1,4140 2050
20 1,4860 1950
Donde A y B tienen que cumplir las siguientes
condiciones, A=0 y B=hc/as Kb, es decir B será el
valor de la constante de Wien y A estará ligado al Realizaremos el análisis estadístico paso a paso
error y al desfase de los datos que se esperan seas casi con el uso de las siguientes fórmulas para nuestra
nulos. correspondiente regresión lineal:

Observaremos dos estrellas, una roja y otra azul, y=A+Bx (5)

indicando cuál de ellas es más caliente, con ayuda de
la aplicación Sky Map localizaremos dichas estrellas A=1/D [− (Σx) (Σxy) + (Σx2) (Σy)] (6)
y compararemos la información teórica que se tiene
de estas para corroborar nuestra observación. B=1/D [− (Σx) (Σy)+ n(Σxy)] (7)

Finalmente realizaremos un programa en el lenguaje D= − (Σx)2+n(Σx2) (8)

Python, con ayuda del interfaz del paquete científico
Anaconda, el cual tendrá como fin integrar R= [ n(Σxy)−(Σx)(Σy)] .
numéricamente la función de Planck, con los cuales [{n(Σx2)−(Σx)2}{n(Σy2)−(Σy)2}]1/2 (9)
realizaremos una comparación con los datos
obtenidos del simulador y realizaremos su respectivo Para el cálculo de las respectivas constantes y para el
análisis. análisis de error y usaremos las siguientes formulas:
ΔB= [n (Ơ2/D)]1/2 (10)
 brindó sobre su temperatura fueron: 9602[K] para
ΔA= ΔB [(Σx2/n]1/2 (11) la estrella vega y para la estrella Aldebarán
3910[K] corroborando que sus temperaturas están
Ơ2= [1/(n−2)] Σ (y−A−Bx)2 (12) ligadas a la longitud de onda en la que emiten,
nótese también que a simple vista puede aparecer
Para la determinación de la constante de Wien afectos de aberración del color debido a la gran
realizamos la regresión lineal de acuerdo a la distancia y al giro de las estrellas conocido como
ecuación (4), la cual nos dio la siguiente grafica efecto Doppler.
(véase figura 2).
Realizamos el programa y comparamos los
1,60E-06 resultados con los datos obtenidos por el
1,40E-06 simulador PHET con los de nuestro programa los
1,20E-06 cuales se muestran en la siguiente tabla 3:
λ [µm]

1,00E-06
8,00E-07 Tabla 3. Comparación de datos obtenidos por el
6,00E-07 y = 0,0029x
simulador PHET y nuestro programa. Fuente:
4,00E-07 R² = 1 elaboración propia.
2,00E-07
0,00005 0,00025 0,00045
1/T [1/K]
N T [K]   m DISC
mP PROG. [%]
Figura 2. Grafica longitud de onda vs inverso de HET
la temperatura. Fuente: elaboración propia. 1 5800 0,5000 0,4772 4,78
2 6000 0,4830 0,4620 4,35
Ahora realizamos las sumatorias necesarias para 3 6500 0,4460 0,4279 4,07
cumplir con las ecuaciones 5-12 se obtuvieron los
siguientes datos mostrados en la tabla 2 para las 4 7000 0,4140 0,3985 3,76
constantes y sus respectivos errores: 5 7500 0,3860 0,3728 3,41
6 7850 0,3690 0,3568 3,32
Tabla 2. Obtención de las constantes y sus 7 8500 0,3410 0,3303 3,13
correspondientes errores. Fuente: elaboración
propia. 8 9000 0,3220 0,3125 2,95
A B D 9 9500 0,3050 0,2965 2,78
-1.124E-10 0.002898236411 6,71198E-06 10 9950 0,2910 0,2835 2,59
ΔA ΔB Des Ơ 11 10000 0,2900 0,2821 2,73
5.81582E-38 4,16289E-16 5,6842E-20 12 10500 0,2760 0,2690 2,54
13 11000 0,2630 0,2571 2,25
Donde podemos notar como previamente
14 6950 0,4170 0,4012 3,79
habíamos afirmado que la constante A es tan
pequeña que se puede decir que es casi nula, así 15 5200 0,5570 0,5295 4,94
como el resultado de la constante B nos da el valor 16 4950 0,5850 0,5548 5,16
de la constante de Wien. 17 3000 0,9660 0,8854 8,35
Una vez terminado el análisis de los datos 18 2350 1,2330 1,1048 10,40
obtenidos por el simulador PHET, se observó en 19 2050 1,4140 1,2470 11,81
fecha 12/06/2020 a horas 11:20 la estrella vega 20 1950 1,4860 1,3036 12,28
color azulado y como estrella rojiza se observó la
estrella Aldebarán con ayuda de la aplicación Sky
Notemos que a medida que la temperatura
Map para conocer sus nombres Se comparó con
desciende el error de nuestro programa aumenta
la biblioteca de Stellarium y los datos que nos
pero la discrepancia con el valor dado por el Comparando con el valor teórico para la constante
simulador PHET son más que aceptables al ser de Wien de 0.002897771955 [mK], la
estas casi todas menores al 5%, pero a medida que discrepancia que obtenemos es del 0.016%, estas
la longitud de onda crece el error va aumentando. discrepancias tan bajas se deben a que los
Podemos notar esto al analizar las gráficas cálculos del simulador PHET son exactos y
brindadas por el simulador y nuestro programa precisos. Se logró elaborar un programa capaz de
como se ven en la figura 4 y 5 con las realizar gráficas de cuerpo negro y va cambiando el
temperaturas de 10000K Y 8000K en naranja y espectro en relación a la temperatura con una
azul respectivamente. discrepancia de hasta 5% respecto al simulador
PHET, pero a medida que nuestro programa va
aumentando la longitud de onda también aumenta el
error, finalmente se logró corroborar los datos de las
estrellas observadas, con la diferencia que al
observarlas existe corrimiento de la luz por efecto
Doppler, llegando así a que esto influya en nuestra
decisión de elección de dichas estrellas, pero se logró
obtener los datos de sus temperaturas y su respectiva
emisión siendo estas : 9602[K] para la estrella vega
y para la estrella Aldebarán 3910[K]
corroborando una vez más el efecto de cuerpo
negro ligado la longitud de onda emitida por las
Figura 4. Grafica brindada por el programa estrellas con su temperatura.
PHET.Fuente:https://phet.colorado.edu/sims/htm
l/blackbody-spectrum/latest/blackbody- 5. Bibliografía.
spectrum_en.html
R. Bustos, Guía experimental I. templete para
trabajos cusco 2015.pdf. 2020
1

1160

Serway R.A., Moses C.J., Moyer C.A., Modern

Physics, primera edición, Brooks Cole, 2004.

Halliday D, Resnick R. Física Vol. 2. Tercera

Edición. Compañía Editorial Continental, S.A.
México. (1996). Pág. 479.

https://phet.colorado.edu/sims/html/blackbody-
Figura 5. Grafica brindada por el programa spectrum/latest/blackbody-spectrum_en.html
realizado. Fuente: elaboración propia.
https://micro.magnet.fsu.edu/primer/lightandcolor/sp
Notemos al comparar la figura 4 y 5 que eedoflight.html
efectivamente la discrepancia es baja y los
resultados ofrecidos con nuestro programa son https://www.fisica.unam.mx/personales/romero/TER
aceptables aunque con un error un poco más MO2017/Radiacion-Cuerpo-Negro-LIBRO.pdf
amplio, se adjunta en la parte de bibliografía el
link respectivo del programa elaborado. https://www.uv.es/zuniga/3.2_Propagacion_de_err
ores.pdf.
4. Conclusiones y recomendaciones.
https://es.wikipedia.org/wiki/Ley_de_desplazamient
Se logró realizar una regresión lineal de nuestros o_de_Wien
datos obteniendo como constante de Wien con un
valor de (0.002898236411± 4,1628E-16) [mK].
https://es.wikipedia.org/wiki/Historia_Cuerpo_Negr
o

https://www.anaconda.com/products/individual

Link programa realizado para descarga:

https://drive.google.com/file/d/1ODfOuePvlXo_ZnR
UeeQ3N_ahDsFDr8dq/view?usp=shari