REVISTA CUBANA DE FÍSICA Vol. 22, No.

1, 2005

NUEVA FORMULACIÓN DEL PRINCIPIO DE HUYGENS

PARA ONDAS PLANAS
Julio A. Romero*, Centro de Estudios para la Educación Superior (CEPES), Universidad de La Habana, Cuba
Luis Hernández**, Facultad de Física, Universidad de La Habana, Cuba

RESUMEN
En el presente artículo se presenta una nueva formulación matemática del Principio de Huygens para
ondas planas, el cual es descrito a través de una ecuación vectorial integral. Se demuestra que las
soluciones de los problemas que pueden plantearse con la nueva formulación satisfacen la ecuación
homogénea de Helmholtz. Una expresión para la densidad de energía de las ondas emitidas por las
fuentes secundarias es obtenida. Para ilustrar la nueva formulación se realiza el cálculo exacto de la
intensidad de la luz difractada por una abertura circular a lo largo del eje perpendicular a su centro,
encontrándose una ecuación analítica exacta que describe la intensidad de la onda difractada. A conti-
nuación se hace un análisis comparativo entre los resultados aquí obtenidos y los correspondientes a la
teoría de Fresnel. Se realiza un análisis similar para una placa zonal.

Palabras clave: Principio de Huygens, ondas planas, Teoría de Fresnel.

ABSTRACT
The article presents results corresponding to a new mathematical formulation of Huygens’ Principle for
plane waves which is described by a vectorial integral equation. We showed that the solutions of the
integral equation fulfill the Helmoltz’s equation. Subsequently the density of energy of the secondary
waves was determined. The task is illustrated by the diffraction due to a circular aperture. We found an
exact analytical equation that describes relative intensity of the diffracted light as a function of the
distance from the center of the circular aperture. For the first maxima and minima, the intensity behavior
is similar to the one predicted by Fresnel’s theory. However, for distances near to the aperture, the
difference is remarkable. An extension of the theory is made to study the diffraction through the zonal
plates. It must also be pointed out that contrary to previous methods; the present formulation is dealt
with a vectorial analysis without recurring to any approximations.

Key words: Huygens’ Principle, plane waves, Fresnel theory.

INTRODUCCIÓN La difracción que tiene lugar a través de una

abertura circular iluminada por una onda con-
Christiaan Huygens dio a conocer el principio que vergente ha recibido especial atención durante
lleva su nombre en 1678, y puede decirse, que este años [1]. Por supuesto, un problema particular
principio le debe mucho a la confianza que depositó puede ser resuelto mediante métodos numéricos
su autor en la teoría ondulatoria de la luz. Pero no pero las teorías simplificadas impiden entender las
sería hasta el año 1818, que Augustin Jean Fresnel tendencias y predecir comportamientos. Partiendo
planteara una fórmula matemática para hacer un uso de la teoría simple de Debye [2] varias mejoras han
cuantitativo de este principio. sido introducidas al estudio de la difracción que
incluye la extensión a sistemas ópticos de gran
La aproximación escalar, para explicar la luz apertura [3], considerar los efectos de la polariza-
difractada a través de las integrales de Fresnel- ción [4] y sistemas ópticos con un número finito de
Kirchhoff, es relativamente simple y una herramienta zonas de Fresnel [5]. Por otra parte, recientemente
muy útil para describir la propagación de la luz se han publicado varios trabajos que estudian la
coherente. Con esta aproximación es posible resolver aproximación de Fresnel para el cálculo de la
un buen número de problemas actuales de la óptica intensidad difractada fuera del eje de simetría [6.7].
con un significado físico real. Sin embargo, la aproxi- Sin embargo, en todos estos trabajos se utiliza la
mación de Fresnel está restringida a la difracción bajo aproximación escalar.
pequeños ángulos donde es excluida la naturaleza
vectorial del campo eléctrico, la cual hay que tener en En un trabajo preliminar anterior [8] presentamos
cuenta cuando la difracción ocurre bajo grandes una nueva formulación matemática del Principio de
ángulos. Huygens para ondas planas, el cual es descrito

*
E-mail: julioarr@cepes.uh.cu
**
luisman@ff.oc.uh.cu;

67
 r r
mediante una ecuación vectorial integral. Hemos E( z, t ) = E 0 ei( wt + k ( z−R ) + α ) (4)
demostrado que las soluciones de los problemas
que pueden plantearse con la nueva formulación
satisfacen la ecuación homogénea de Helmholtz y El vector del campo eléctrico en el punto p, se
las condiciones de Sommerfeld. En el presente obtiene evaluando en z = 0. Teniendo en cuenta el
trabajo se reporta una expresión para la densidad de carácter transversal de las ondas electromagnéticas,
energía de las ondas emitidas por las fuentes asumiendo que Ep(t) representa cualquiera de las
secundarias. Para ilustrar la nueva formulación se componentes transversales y que Eo es su amplitud
realiza el cálculo exacto de la intensidad de la luz correspondiente, se tiene:
difractada por una abertura circular a lo largo del eje
perpendicular a su centro, encontrándose una ecua- Ep(t) = E0ei(wt - kR + α) (5)
ción analítica exacta que describe la intensidad de la
onda difractada. Para los primeros máximos y Sustituyendo E en la ecuación 2 por la relación 5,
mínimos de intensidad de la luz difractada, su el valor de cada componente transversal del vector
comportamiento es idéntico al predicho por la teoría del campo eléctrico se obtiene integrando 2 por todo
de Fresnel, sin embargo para distancias cercanas a la el plano z = R. Para determinar la componente z del
abertura su diferencia es sustancial. La teoría también vector del campo eléctrico de las ondas secundarias,
es utilizada para el estudio de las placas zonales se asume que éstas son también transversales.
encontrándose que la posición de los máximos no Considerando lo expresado anteriormente y el hecho
se altera, pero si el ancho medio de los picos. de que el origen de coordenadas del sistema de
Igualmente se observa el efecto de la difracción referencia no tiene que coincidir con el punto donde
disminuyendo la intensidad de los máximos al se está evaluando el vector del campo eléctrico, se
acercarnos al plano de la placa zonal. llega a una ecuación vectorial relevante, que es la
base de la nueva formulación del Principio de
2. ECUACIÓN INTEGRAL Huygens para las ondas planas:
r r
r r ⎛ e m ik r ′− r ⎞
Aplicando la segunda fórmula de Green a un
dominio en forma de pirámide recta de base circular
r r r
( r r
) r′ − r
dE( r ′, r , t ) = ds′ × Ep ( t ) × r r r r
∂ ⎜ r r
r ′ − r ∂ r ′ − r ⎜ 2π r ′ − r
⎟ (6)
⎟
de altura R y pasando al límite cuando el radio de la ⎝ ⎠
base del cono tiende a infinito, se obtiene una
función escalar acotada, continuamente derivable y donde
que es además solución de la ecuación homogénea
r r
de Helmholtz. Esta función toma la forma de: E p ( t ) = E 0 e ±i( wt + α ) (7)

E Re ±ik R ∂ ⎛⎜ e mik r′ ⎞⎟ ds′

F(A, E) = −
2π ∫
A
∂ r ′ ⎜⎝ r ′ ⎟⎠ r ′
(2) El vector del campo eléctrico dado por la ecuación 7
corresponde al campo de una onda plana en el
punto de observación, y en general, tomará una u
otra forma en función de las condiciones de contorno
f(plano, E) = E (3)
del problema de difracción específico que se
considere. Es fácil comprobar que la integral de la
En la relación, f(A,E) es el valor de la función
ecuación 6 para cualquier área A del plano de
escalar en el punto p que se encuentra situado a difracción satisface la ecuación homogénea de
una distancia R > 0 del plano de integración, A es el Helmholtz.
área de integración sobre el plano y r’ es la distancia
entre el punto p y cualquier punto p’ del área de 4. ONDAS SECUNDARIAS
integración. Integrando sobre todo el plano se
obtiene la relación 3. A partir de la ecuación 6 y considerando esféricos
los frentes de ondas secundarios, la densidad de
3. PRINCIPIO DE HUYGENS energía de las ondas esféricas secundarias viene
dada por:
Se considera un régimen estacionario de ondas r 2
ε E0 (t) ⎛
electromagnéticas planas monocromáticas y un
medio homogéneo, isótropo y lineal; además el frente W=
8 π 2r 2 Γ 2 ⎝
1⎞
r ⎠
(
⎜ k 2 + 2 ⎟ 1 + Cos 2 (φ) (8) )
de onda plano con fase α en el instante inicial se
encuentra a la distancia R de un punto p y a partir de
dicho instante se mueve acercándose al punto en En esta última ecuación Γ es el número de fuentes
cuestión. Colocando el origen de coordenadas del secundarias por unidad de área, ε es la constante
sistema de referencias en dicho punto, entonces, el dieléctrica del medio y φ es el ángulo que forma la
vector del campo eléctrico de la onda plana dirección de propagación de la onda esférica
monocromática, viene dado por: secundaria con la dirección de propagación de la
onda plana, como puede apreciarse en la Figura 1.

68
 En la ecuación puede observarse que existe El vector del campo eléctrico de la onda difractada
simetría axial respecto a la dirección de propagación se determina integrando la ecuación 6 por el área de
de la onda plana, además, la densidad de energía de la abertura circular:
las ondas secundarias disminuye con φ y con 1/r2.
r v
E( r , t ) =
−1
2π
[ (r r
)
r
ê z ρ 0E p ( t ) I1 + zE p ( t )I2 ] (9)

r r r
r r′
donde E p ( t ) ⊥ ê z ρ 0 = ρ 0 k = k ê z r2 = ρ2 + z2
r′

π
e m ik ρ 02 + r 2 − 2ρ 0 ρ Cos(ψ )
I1 = 2 ∫
0 ρ 02 + r 2 − 2ρ 0 ρ Cos(ψ )
Cos(ψ ) dψ (10)

2 2
π ρ0 + r − 2ρ 0 ρ Cos( ψ )
∂ ⎛⎜ e m ik R ⎞
I2 = 2 ∫0
∫
r
∂ R ⎜⎝ R
⎟
⎟
⎠

⎛ ⎞
⎜ ρ Cos(ψ ) ⎟
⎜⎜1 + ⎟ dR dψ (11)
⎝ R − r + ρ Cos (ψ ) ⎟⎠
2 2 2 2

Figura 1. Esquema referido a la ecuación 8, para describir Las ecuaciones 9-11 permiten expresar un resul-
la propagación de las ondas secundarias. tado analítico para la intensidad relativa I/I0 de la
r r onda difractada a lo largo del eje de simetría axial de
El ángulo que forman los vectores k y r es φ.
la abertura:

5. DIFRACCIÓN POR UNA ABERTURA CIRCULAR

I z2 2z ⎛ ⎞
= 1+ 2 − Cos ⎜ k ⎛⎜ z 2 + ρ 02 − z ⎞⎟ ⎟
I0 z + ρ 02 z 2
+ ρ 02 ⎝ ⎝ ⎠ ⎠
Supongamos que tenemos una onda plana mono-
cromática incidiendo de forma normal sobre un (12)
plano opaco dotado de una abertura circular de radio
ρo. Coloquemos el origen de coordenadas O sobre el En la Figura 3 puede verse representado los
centro de la abertura y tomemos el eje z en la gráficos de intensidad relativa dados por la ecuación
dirección y sentido de propagación de la onda, como 12 en función de la distancia para una longitud de
muestra la Figura 2. onda λ = 633 nm y ρ0 = 1 mm, valores que se toman
como ejemplo para ilustrar los resultados. Obsérvese
que para valores de z >> ρ0 la ecuación 12 se
reduce a la expresión clásica de acuerdo a la teoría
de Fresnel:

⎛ ⎡ kρ 2 ⎤ ⎞
I = 2 I0 ⎜1 − Cos ⎢ 0 ⎥ ⎟ (13)
⎜ ⎟
⎝ ⎣⎢ 2z ⎦⎥ ⎠

De tal manera que la posición de los máximos y

los mínimos a lo largo del eje coincide con los
predichos por la teoría de Fresnel para los primeros
órdenes. Para distancias cercanas o menores al
valor de z = ρ0, la diferencia entre los resultados que
se obtienen a partir de la ecuación 12 y la teoría de
Fresnel (ecuación 13) se hacen evidentes, como se
muestra en la Gráfica 3. Este hecho puede también
Figura 2. observarse en la forma analítica de las envolventes
Pantalla opaca con una abertura circular de la ecuación 12:
de radio ρ0.

69
 2 2
r r 2n +1 m ik ρl + z
( −1)l e
E( z, t ) = z Ep ( t ) ∑
l= 0 ρl2 + z 2
(17)

donde ρl son los radios de los anillos. La intensidad

relativa para un número de anillos n+1 viene ahora
dada por:

I
2n +1
((
( −1)l+ j z 2Cos k w l − w j ))
I0
= ∑l=0
wl − w j
(18)
j=0

donde wl = ρl2 + z 2 (19)

La relación que guardan los radios de los anillos ρl

Figura 3. Intensidad relativa de la luz difractada con la posición del máximo principal z0 se obtiene de
en función de la distancia al centro la ecuación 18:
de la abertura. La figura principal muestra
el máximo de primer orden 2
⎛ 1λ ⎞
en tanto las gráficas insertadas revelan ρl2 = ⎜ ⎟ + ρ 02 + l λ z 02 + ρ 02 0 ≤ l ≤ 2n + 1 (20)
el comportamiento para otros valores de z/ρ0. ⎝ 2 ⎠

⎛ ⎞
2
⎛ ⎞
2 Para el caso de ρ0 = 0, tomando x = z/z0, y = λ/z0,
⎜ z ⎟ I ⎜ z ⎟ sustituyendo en 19 y redefiniendo wl se obtiene:
⎜⎜1 − ≤ ≤ 1+ ⎟⎟ (14)
z + ρ 02
2 ⎟⎟ I0 ⎜⎜ z 2 + ρ 02
⎝ ⎠ ⎝ ⎠ 2
⎛ 1λ ⎞
wl = ⎜ ⎟ + ly + x 2 0 ≤ l ≤ 2n +1 (21)
⎝ 2 ⎠
Estas últimas se obtienen directamente de la
ecuación 12 bajo las condiciones: Teniendo en cuenta estos últimos cambios, la
ecuación 18 se transforma en la siguiente relación:
k ⎛⎜ z 2 + ρ 02 − z ⎞⎟ = 2 n π
⎝ ⎠
I
2n +1
(
( −1)l+ j z 2Cos k( w l − w j ) )
I0
= ∑ wlw j
0 ≤ l ≤ 2n +1 (22)
k ⎛⎜ z 2 + ρ 02 − z ⎞⎟ = (2 n + 1) π
l= 0
j =0
⎝ ⎠

ρ0
n≤ (15)
λ

Para ρo << z, estas se convierten en las condi-

ciones de mínimo y máximo respectivamente, de la
teoría de Fresnel:

ρ02 ρ 02
z≈ z≈ (16)
2n λ λ(2 n + 1)

6. DIFRACCIÓN POR UNA PLACA ZONAL

Considerando una placa zonal, un juego de n+1

anillos concéntricos que se alternan entre opacos
y transparentes (ver Figura 4) e integrando la
ecuación 6 por el área de los anillos transparentes
se obtiene, a lo largo del eje de simetría axial, una
expresión para el vector campo eléctrico de la onda
difractada por la placa zonal, a saber:
Figura 4. Croquis de la placa zonal.

70
 Evaluando la ecuación 21 se obtiene la intensidad Las condiciones iniciales y de contorno, para los
de la onda difractada por la placa zonal a lo largo del problemas que pueden plantearse con la ayuda de
eje axial en función del número de anillos. La Figura 5 esta nueva formulación del Principio de Huygens, se
muestra esta dependencia para λ = 633 nm y una establecen por medio de la ecuación 7 y por la
posición del máximo principal fijada en 100 mm. integral de esta última por el área correspondiente al
Obsérvese que los máximos se van estrechando en problema en cuestión. También se satisfacen las
la medida que el número de anillos se incrementa condiciones de Sommerfeld, garantizándose de esta
pero que no cambian de posición, mostrando como forma la unicidad del problema tratado.
era de esperarse un comportamiento análogo a una
red de difracción. Nótese también en la Figura 5-d Por otro lado, partiendo de la nueva formulación y
que los máximos de intensidad disminuyen su considerando esféricos los frentes de ondas de las
valor en la medida que su posición se acerca a la fuentes secundarias, se obtuvo una expresión para
placa zonal. la densidad de energía espacial de las ondas
secundarias (ecuación 8), realizándose un análisis
7. CONCLUSIONES del comportamiento de esta función.

Encontrada una función escalar (ecuación 2) y Al considerar los casos de la placa zonal y de la
sobre la base de las características propias de las abertura circular, se demostró que con la ayuda de
ondas electromagnéticas en medios homogéneos, esta nueva formulación puede llegarse a resultados
isótropos y lineales, se obtuvo por primera vez la analíticos exactos, demostrándose de esta forma su
ecuación vectorial (6), consecuente con el Principio sencillez. En los casos considerados como ejemplos
de Huygens. Esta ecuación vectorial describe se obtuvieron de forma exacta las intensidades de los
correctamente la propagación de los frentes de campos eléctricos. Se hizo un análisis de la forma
ondas planos en el medio en cuestión, en plena funcional de estos campos y se compararon los resul-
concordancia con dicho principio. Se demostró que tados con los de la teoría de Fresnel. Los resultados
las soluciones de los problemas que pueden válidos de la vieja formulación de Fresnel fueron corro-
plantearse con su ayuda, satisfacen la ecuación borados, en tanto que se obtuvieron nuevos resulta-
homogénea de Helmholtz, como corresponde a dos más acorde a la realidad en los casos donde la
regímenes estacionarios. teoría de Fresnel falla en sus predicciones.

Figura 5. Las gráficas muestran la dependencia de la intensidad en función de la posición (x = z/z0,

siendo z0 la posición del máximo principal). (a) Para una abertura circular, (b) para cuatro anillos,
(c) para nueve anillos. En la figura (d) se observa el patrón para distancias cercanas a la placa zonal
para nueve anillos, en donde los máximos disminuyen su valor.

71
 REFERENCIAS

[1] STAMNES, J.J. (1986): “Waves in focal regions”. Hilger, Bristol.

[2] BORN, M. and E. WOLF (1983): “Principles in Optics” 6th ed. Pergamon, Oxford.

[3] HOPKINS, H.H. (1943): “The Airy disc formula for systems of high relative aperture”. Proc. Phys.
Soc. London 55, 116.

[4] RICHARDS, B. and E. WOLF (1959): “Electromagnetic diffraction in optical systems”. Proc. R.
Soc. London, Ser. A 253, 358.

[5} LI, Y. and E. WOLF (1984): “Three-dimensional intensity distribution near the focus in systems
of different Fresnel numbers”, J. Opt. Soc. Am. A1, 801.

[6] MIELENZ, KLAUS D. (1998): “Algorithms for a Fresnel Diffraction at Rectangular and Circular
Apertures”, J. Res. Natl. Inst. Stand. Technol. 103, 497.

[7] SHEPPARD, C.l.R. and I.I. COOPER (2004): “Fresnel diffraction by a circular aperture with off-
axis illumination and its use in deconvolution of microscope images”, J. Opt. Soc. Am. A21,
540-545.

[8] ROMERO, J.; O. CALZADILLA y L. HERNÁNDEZ (2000): “Principio de Huygens para las ondas
planas monocromáticas y difracción por una abertura circular”, Revista Cubana de Física
17, 10.

72