El arco iris

Vista de los arcos primario y secundario, en Durango (Vizcaya). Al fondo el monte Oiz.

Roger Bacon en 1266 midió el ángulo que forman los rayos del arco iris con la luz solar incidente. Descartes demostró que el arco iris primario está formado por los rayos que penetran en una gota refractándose, se reflejan una vez en su superficie interna y salen de la gota refractándose de nuevo. El arco iris secundario está formado por los rayos que penetran en la gota y se reflejan dos veces en su superficie interna. La figura de abajo, es la explicación de Descartes a la formación del arco primario y secundario.

En la figura de abajo, se muestra el arco primario y secundario producidos en una simulación con ordenador, no se ha tenido en cuenta que el arco iris secundario es menos intenso que el primario.

Observamos también que el arco secundario tiene una anchura angular mayor que el primario. Los ángulos de desviación de la luz monocromática se calculan en los siguientes apartados

Indice de refracción

La refracción, como observó Newton con el prisma, depende del color. La luz de longitudes de onda diferentes se desvían según ángulos ligeramente distintos, esta dependencia entre el índice de refracción y la frecuencia de la luz se denomina dispersión.

Longitud de onda, nm Indice de refracción del agua a 20º C
670.8 1.3308
643.8 1.3314
589.3 1.3330
486.4 1.3371
404.7 1.3428

Fuente: Koshkin N. I., Shirkévich M. G.. Manual de Física Elemental. Editorial Mir 1975. pág. 210

Tomando como modelo la fórmula de Cauchy para los índices de refracción del agua en la región visible del espectro, ajustamos los datos de la tabla a la función no lineal

y=a+ b x 2

x=[670.8,643.8,589.3,486.4,404.7];  %longitud de onda en nm
y=[1.3308,1.3314,1.3330,1.3371,1.3428]; %índice de refracción del agua
hold on
plot(x,y,'bo','markersize',4,'markerfacecolor','b')
f_ajuste =@(a,x) a(1)+a(2)./x.^2;          
a0=[1.3 3000];  %valor inicial de los parámetros
af=nlinfit(x,y,f_ajuste,a0);
%representa la función 
g=@(x) f_ajuste(af,x);
fplot(g,[400,700],'r')
hold off
xlabel('\lambda (nm)')
ylabel('n')
grid on
title('Indice de refracción del agua')
>> af =   1.0e+03 *    0.0013    3.0843

n=1.324+ 3084.3 λ 2 λ en nm

Trayectoria de un rayo de luz monocromático

Cuando un rayo de luz incide en la superficie de una gota de forma esférica, una parte del rayo se refleja y la otra parte se transmite (1). La luz transmitida pasa al interior de la gota donde cambia de dirección debido a la refracción. Cuando llega a la superficie interna de la gota (2), una parte se refleja en su superficie y la otra parte se transmite al exterior cambiando de dirección debido a la refracción y así, sucesivamente.

El rayo incidente forma un ángulo i con la dirección radial, se refracta acercándose a la normal, formando un ángulo r, con la dirección radial tal que

sini=n·sinr

El rayo cambia de dirección un ángulo (i-r)

Se refleja k veces en la superficie de la gota, cambiando de dirección (π-2r) en cada reflexión

Se transmite al exterior de la gota. El ángulo incidente es ahora r y el ángulo refractado forma un ángulo i con la dirección radial. El rayo cambia de dirección un ángulo (i-r)

La desviación total del rayo, es

D=(i-r)+k(π-2r)+ (i-r)=kπ+2i-2(k+1)r

Donde k es el número de veces que el rayo se refleja en el interior de la gota

Creamos un script de MATLAB, para tazar la trayectoria de un rayo de luz de una determinada longitud de onda. Se introduce el ángulo incidente en grados, la longitud de onda en nanómetros y el número de reflexiones internas (una o dos). El programa nos proporciona el ángulo de desviación en grados en la esquina superior derecha

angulo=60;  %angulo incidente en grados
lonOnda=671; % en nm=10^-9 m
k=1; % número de reflexiones

%cálculo
n=1.324+3084.3/lonOnda^2; %índice de refracción
%dibuja la gota
R=1; % radio de la gota
ang=0:360;
x=R*sind(ang);
y=R*cosd(ang);
hold on
%marco
xx=[-3*R/2,3*R/2,3*R/2,-3*R/2];
yy=[3*R/2,3*R/2,-3*R/2,-3*R/2];
fill(xx,yy,'k')
fill(x,y,[0.75 0.75 0.75]) %gota de color gris claro, silver

%trazado de rayos
angIncidente=angulo*pi/180;
y1=R*sin(angIncidente);
x1=-R*cos(angIncidente);
line([-3*R/2,x1],[y1,y1],'color',[1 0 0]);
angRefractado=asin(sin(angIncidente)/n); %ángulo refractado
line([0 x1],[0,y1],'color','w')
for i=0:k
    m=-tan(angIncidente-angRefractado+i*(pi-2*angRefractado));
    b=y1-m*x1;
    x2=-x1-2*m*b/(m^2+1);
    y2=m*x2+b;
    line([0 x2],[0,y2],'color','w');
    line([x1,x2],[y1,y2],'color',[1 0 0]);
    x1=x2; y1=y2;
end
delta=-k*pi-2*angIncidente+2*(k+1)*angRefractado;
m=tan(delta);
b=y2-m*x2;
x2=x1+R*cos(delta)/2;
y2=y1+R*sin(delta)/2;
line([x1,x2],[y1,y2],'color',[1 0 0]);	
texto=sprintf('%3.2f\n',-delta*180/pi);
text(R,R,texto,'color','w')

axis equal
xlim([-3*R/2,3*R/2])
ylim([-3*R/2,3*R/2])
hold off

Modificamos el código del script para trazar un haz de rayos paralelos al eje horizontal de la gota. Cambiando el valor del número de reflexiones internas k=2, trazamos la trayectoria que sigue la luz de una determinada longitud de onda para formar el arco secundario

lonOnda=671; % en nm=10^-9 m
k=1; % número de reflexiones

%cálculo
n=1.324+3084.3/lonOnda^2; %índice de refracción
%dibuja la gota
R=1; % radio de la gota
ang=0:360;
x=R*sind(ang);
y=R*cosd(ang);
hold on
%marco
xx=[-3*R/2,3*R/2,3*R/2,-3*R/2];
yy=[3*R/2,3*R/2,-3*R/2,-3*R/2];
fill(xx,yy,'k')
fill(x,y,[0.75 0.75 0.75]) %gota de color gris claro, silver

%trazado de rayos paralelos al eje horizontal
for y1=0:0.05:0.95
    angIncidente=asin(y1/R);
    x1=-R*cos(angIncidente);
    line([-3*R/2,x1],[y1,y1],'color',[1 0 0]);
    angRefractado=asin(sin(angIncidente)/n); %ángulo refractado
         for i=0:k
            m=-tan(angIncidente-angRefractado+i*(pi-2*angRefractado));
            b=y1-m*x1;
            x2=-x1-2*m*b/(m^2+1);
            y2=m*x2+b;
            line([x1,x2],[y1,y2],'color',[1 0 0]);
            x1=x2; y1=y2;
        end
    delta=-k*pi-2*angIncidente+2*(k+1)*angRefractado;
    m=tan(delta);
    b=y2-m*x2;
    x2=x1+R*cos(delta)/2;
    y2=y1+R*sin(delta)/2;
    line([x1,x2],[y1,y2],'color',[1 0 0]);	
end
axis equal
axis off   %elimina los ejes
xlim([-3*R/2,3*R/2])
ylim([-3*R/2,3*R/2])
hold off

Rayo cartesiano

Fijando el número k de veces que el rayo se refleja internamente. Para distintos valores del ángulo de incidencia i obtenemos distintos valores del ángulo de desviación D.

Los arcos se forman para un ángulo de incidencia θc tal que el ángulo de desviación D alcanza el valor mínimo. Entonces la densidad de los rayos que dejan la gota es máxima.

d D d i = 2 2 ( k + 1 ) d r d i cos i = n cos r d r d i cos r = 1 sin 2 r = 1 1 n 2 sin 2 i = 1 1 n 2 + 1 n 2 cos 2 i d D d i = 2 2 ( k + 1 ) cos i n 2 1 + cos 2 i

Poniendo dD/di=0, se despeja el ángulo de incidencia i.

cos 2 i = n 2 1 k ( k + 2 )

Este es el ángulo de incidencia del denominado rayo cartesiano

Para determinar si las direcciones de los arcos corresponden a un máximo o a un mínimo, volvemos a derivar el ángulo de desviación D, respecto del ángulo de incidencia i.

d 2 D d i 2 = 2 ( k + 1 ) ( n 2 1 ) sin i ( n 2 1 + cos 2 i ) 3 / 2

La derivada segunda es positiva ya que n>1.

Arco primario

El arco primario se forma cuando k=1, hay una sola reflexión interna.

En color azul claro, el rayo cuyo ángulo de incidencia da lugar a una desviación mínima.

Creamos un script para representar el ángulo de desviación D en función del ángulo de incidencia i. Fijamos el número de reflexiones internas k=1 (arco primario) y la longitud de onda λ en nanómetros

k=1;  %número de reflexiones internas
lonOnda=671;  %longitud de onda en nm
n=1.324+3084.3/lonOnda^2; % índice de refracción
x=55:0.1:65;
%x es el ángulo incidente, 
%el ángulo refractado, sin(x)=n·sin(r), r=asin(sin(x)/n)
D=(k*pi+2*(x*pi/180)-2*(k+1)*asin(sin(x*pi/180)/n))*180/pi;
plot(x,D)
xlabel('i')
ylabel('D')
title('Angulo de desviación')
grid on
[xmin,nmin]=min(D);
fprintf('El mínimo de desviación se produce para el 
ángulo %2.2f y vale %3.2f\n',x(nmin), xmin)
El mínimo de desviación se produce para el 
        ángulo 59.50 y vale 137.61

En la figura, se muestra la desviación D en función de ángulo de incidencia, para la luz de color rojo, λ=671 nm, el valor mínimo se obtiene para i=59.5º, el ángulo de desviación vale Dm=137.6º

Arco secundario

El arco primario se forma cuando k=2, hay dos reflexiones internas.

En color azul claro, el rayo cuyo ángulo de incidencia da lugar a una desviación mínima.

Modificamos el script anterior, para mostrar el ángulo de desviación D en función del ángulo incidente i cuando el número de reflexiones internas k=2

k=2;  %número de reflexiones internas
lonOnda=671;  %longitud de onda en nm
n=1.324+3084.3/lonOnda^2; % índice de refracción
x=65:0.1:75;
%x es el ángulo incidente, 
%el ángulo refractado, sin(x)=n·sin(r), r=asin(sin(x)/n)
D=(k*pi+2*(x*pi/180)-2*(k+1)*asin(sin(x*pi/180)/n))*180/pi;
plot(x,D)
xlabel('i')
ylabel('D')
title('Angulo de desviación')
grid on
[xmin,nmin]=min(D);
fprintf('El mínimo de desviación se produce para el 
ángulo %2.2f y vale %3.2f\n',x(nmin), xmin)
El mínimo de desviación se produce para el 
        ángulo 71.90 y vale 230.33

En la figura, se muestra la desviación D en función de ángulo de incidencia, para la luz de color rojo, λ=671 nm, el valor mínimo se obtiene para i=71.9º, el ángulo de desviación vale Dm=230.3º

Angulo de desviación mínima

Las fórmulas que precisamos para calcular el ángulo de desviación mínimo son las siguientes:

Arco, k Rojo, λ=671 nm, Violeta, λ=404 nm, ΔD
i Dm i Dm
1 59.5º 137.6º 58.8º 139.3º 1.7º
2 71.9º 230.3º 71.5º 233.4º 3.1º

En la última columna, tenemos la anchura angular ΔD de los arcos primario y secundario. El arco secundario es casi el doble de ancho que el primario.
El ángulo de la zona oscura de Alejandro Δφ es la diferencia entre el ángulo del rojo en el arco primario visto por el observador terrestre φ1=180-137.6=42.4º y el rojo del arco secundario visto por el mismo observador φ2=230.3-180=50.3º. Δφ=φ2-φ1=7.9º. Véase la tercera figura de esta página.

Intensidad y polarización de la luz

Las fórmulas de Fresnel nos dan los coeficientes de reflexión y de transmisión de las ondas electromagnéticas en la separación entre dos medios no conductores.

R N = ( n 1 cos θ i n 2 cos θ t n 1 cos θ i + n 2 cos θ t ) 2 T N = 4 n 1 n 2 cos θ i cos θ t ( n 1 cos θ i + n 2 cos θ t ) 2 R P = ( n 2 cos θ i + n 1 cos θ t n 2 cos θ i + n 1 cos θ t ) 2 T P = 4 n 1 n 2 cos θ i cos θ t ( n 2 cos θ i + n 1 cos θ t ) 2

Fuente: Lorrain P. Corson D. L. Campos y Ondas Electromagnéticas. Selecciones Científicas (1972) pág. 554

Vamos a calcular las intensidades relativas de los rayos que dan lugar al arco primario y secundario.

Un rayo de luz no polarizada de longitud de onda  λ, incide sobre una gota de agua de índice de refracción n2=n, que a su vez depende de la longitud de onda. El índice de refracción del aire es constante e igual a n1=1. Parte de la luz incidente se refleja y otra parte se refracta. Si el ángulo de incidencia es θi=i y el ángulo que forma la luz refractada con la normal es θt=r. Los coeficientes se expresan

R N = ( cosincosr cosi+ncosr ) 2 T N =1 R N R P = ( ncosi+cosr ncosi+cosr ) 2 T P =1 R P

Vamos a calcular los coeficientes para la luz de color rojo, λ=671 nm, n=1.331

Arco primario (k=1)

El ángulo de incidencia que corresponde a la desviación mínima es i=59.54º

El ángulo r se calcula aplicando la ley de Snell, sini=n·sinr

En (1) tenemos que una parte de la luz se transmite y otra parte se refleja.

La luz atraviesa la gota y en (2) tenemos que una parte se refleja en su superficie interior y otra se transmite.

Los datos ahora son: n2=1, n1=n. El ángulo de incidencia es θi=r y el ángulo que forma la luz refractada con la normal es θt=i.

R ' N = ( ncosrcosi ncosr+cosi ) 2 = R N R ' P = ( cosr+ncosi cosr+ncosi ) 2 = R P

que son las mismas fórmulas que hemos obtenido anteriormente

La luz atraviesa la gota y en (3) tenemos que una parte se refleja en su superficie interior y otra se transmite.

Los datos ahora son: n2=1, n1=n. El ángulo de incidencia es θi=r y el ángulo que forma la luz refractada con la normal es θt=i.

La relación T3N/T3P=24

Total, I=T3N+T3P=915·10-4

Arco secundario (k=2)

El ángulo de incidencia que corresponde a la desviación mínima es i=71.92º

El ángulo r se calcula aplicando la ley de Snell, sini=n·sinr,

En (1) tenemos que una parte de la luz se transmite y otra parte se refleja.

La luz atraviesa la gota y en (2) tenemos que una parte se refleja en su superficie interior y otra se transmite.

La luz atraviesa la gota y en (3) tenemos que una parte se refleja en su superficie interior y otra se transmite.

La luz atraviesa la gota en (4) tenemos que una parte se refleja en su superficie interior y otra se transmite.

La relación T4N/T4P=9

Total, I=T4N+T4P=390.4·10-4

La luz proveniente del arco iris primario está mucho más polarizada que la proveniente del arco secundario.

En general, después de k reflexiones internas la intensidad del rayo trasmitido es

T N = ( 1 R N ) 2 R N k = ( 1 ( cosincosr cosi+ncosr ) 2 ) 2 ( cosincosr cosi+ncosr ) 2k T P = ( 1 R P ) 2 R P k = ( 1 ( ncosi+cosr ncosi+cosr ) 2 ) 2 ( ncosi+cosr ncosi+cosr ) 2k T= T N + T P

Los valores de TN y TP para los arcos supernumerarios (k>2) son muy pequeños.

Actividades

Se introduce

Se pulsa el botón titulado Nuevo

El programa interactivo calcula el ángulo de desviaciación D, y los coeficientes IP e IN que hay que multiplicar por 10-4. La magnitud relativa de estos coeficientes nos da una idea acerca de la polarización de la luz y su suma de la reducción de la intensidad de la luz.

Se sugiere al lector que para el color violeta λ=404 nm, calcule



En el siguiente programa, se muestra las trayectorias de todos los rayos monocromáticos, para sus ángulos de desviación mínima, respectivos. Se muestran los valores numéricos de las desviaciones mínimas de los rayos cuya longitud de onda corresponde a los extremos del espectro visible: rojo 780 nm, violeta 380 nm.

Referencias.

Whitaker R. J. Physics of the rainbow. The Physics Teacher. Vol 12, May 1974, pp. 283-286.

Hendry A. W. A tripe rainbow?. The Physics Teacher Vol 41, November 2003, pp. 460-463.

Nussenzveig H. M. Teoría del arco iris. Investigación y Ciencia, nº 9, Junio 1977, págs. 82-94.

Walker J. D., Multiple rainbows from single drops of water and other liquids. Am. J. Phys. Vol. 44 nº 5, May 1976, pp. 421-433