CAPÍTULO 7

RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

En el capítulo anterior se examinaron las soluciones de las ecuaciones de Maxwell

en ausencia de fuentes y se estudió el carácter de la propagación de las ondas
electromagnéticas en diversas situaciones. En el presente capítulo se examina la generación
o radiación de ondas electromagnéticas. Para ello se requiere estudiar inicialmente el
problema formal de resolver las ecuaciones para los potenciales electromagnéticos en
presencia de fuentes. Después de mostrar cómo resolver en forma general este problema,
nuestra atención se concentrará en la descripción de los sistemas radiantes más simples
(dipolo, cuadrupolo, antenas sencillas, carga puntual que se mueve arbitrariamente).

§ 7.1. POTENCIALES DE UN SISTEMA DE CARGAS.

Se examinan las soluciones de las ecuaciones para los potenciales escalar y vectorial en
presencia de fuentes, las cuales, para simplificar la discusión se supondrán en el vacío,
aunque los resultados que se obtendrán se pueden generalizar fácilmente a un medio lineal,
isotrópico y homogéneo caracterizado por una permitividad dieléctrica  y una
permeabilidad magnética . En el vacío, los potenciales escalar y vectorial satisfacen las
ecuaciones

 2 1 2 
r ,t     r ,t  ,
1
   2 2 
(7.1.1)
 c t  0
 2 1 2 
   2  Ar , t     0 j r , t  . (7.1.2)
 c  t 2 

Se supone que entre estos potenciales tiene lugar la condición de Lorentz

1 
divA  0. (7.1.3)
c2  t

y se considera que a partir de cierto instante t=0 las cargas se mueven de una manera
arbitraria, y que para t<0 las cargas eran estacionarias y por lo tanto se puede considerar
que para t<0 (r,t)=0, j(r,t)=0. Esto implica que en el instante t=0, que es cuando se
“conecta” la perturbación del sistema, los campos E y B son nulos. Dado que
A
E  grad   y B=rotA, en t=0 los potenciales deben satisfacer las condiciones
t
iniciales A(r,t)t=0=0, (r,t)t=0=0, A(r,t)/tt=0=0.
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

En el caso de distribuciones estacionarias de carga las ecuaciones (7.1.1) y (7.1.2) se

reducen a las ecuaciones de Poisson para campos electrostáticos y magnetostáticos y en el
caso de campos no estacionarios en ausencia de fuentes se reducen a la ecuación de onda.
Las relaciones (7.1.1) y (7.1.2) son ecuaciones no homogéneas lineales y sus
soluciones se pueden representar como la suma de dos términos, uno de los cuales
corresponde a la solución general de la respectiva ecuación homogénea (en ausencia de
fuentes). El otro término corresponde a una solución particular de la ecuación con fuente.
Dado que la estructura de (7.1.1) y (7.1.2) es semejante, se puede limitar nuestra
discusión a la ecuación para el potencial escalar. Para encontrar la respectiva solución
particular se procede de la siguiente manera: se divide todo el espacio en porciones
infinitesimales y se determina el potencial d creado en el punto P por la carga dq=idV´
situada dentro del elemento de volumen dV´. La geometría del problema se ilustra en la
Figura (7.1.1). Sea R=r-r´ la distancia entre el elemento de volumen considerado y el
punto P. Si R es diferente de cero, el potencial d satisface la ecuación

 2 1 2 
   2 2 d   0 .
 c t 

Figura (7.1.1)

De acuerdo con el principio de superposición, este potencial es creado por la carga

dq en el punto P como si no estuviesen presentes las cargas de los demás elementos de
volumen del sistema. Esto indica que el potencial d posee simetría central, es decir, d es
función exclusiva de R. En coordenadas esféricas se tiene entonces

1   2 (d)  1  2 (d)
R   0,
R 2 R  R  c 2 t 2

o bien

324
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

1 2 1  2 (d)
 
R(d)  2 0.
R R 2 c t 2

La solución de esta ecuación se representa en la forma

 ( R ,t )
d( R ,t )  .
R

Al reemplazar esta expresión en la ecuación anterior, se ve que la nueva función 

satisface la ecuación de D´Alembert

2 1 2
  0.
R 2 c 2 t 2

Esta ecuación se estudió con detalle en §7.1 y de acuerdo con los resultados allí
obtenidos, la función  tiene la forma (R,t)=F1(R-Vt)+ F2(R+Vt), la cual puede ser re-
escrita como

 R  R
 ( R ,t )  F1  t    F2  t   ,
 c  c

donde F1 y F2 son funciones que deben ser determinadas. Sin embargo, el sentido físico de
estas funciones se puede determinar fácilmente con base en el siguiente análisis: sea
R1=R+ct; en este caso, F1(t-R1/c)=F1(-R/c), es decir, en R1 la forma de la función F1
coincide con la forma que tiene en el instante inicial t=0. Se ve así que F1(t-R1/c) describe
una onda que se propaga en la dirección de crecimiento de R. Mediante un análisis similar
se concluye que F2(t+R1/c) describe una onda que se propaga en la dirección de
decrecimiento de R. Como se trata de buscar una solución particular de la ecuación (7.1.1),
se toma por ejemplo sólo la función F1 correspondiente a una onda que diverge del sistema
de cargas y se hace F2=0. Se obtiene entonces

 R
F1  t  
d  
c
.
R

Para determinar la función F1(R) obsérvese que R/c es el tiempo que gasta la
información en ir desde el elemento de volumen considerado al punto P (tiempo de
retardo). Este tiempo de retardo es despreciable cuando el punto P está cerca del elemento
dV´. Bajo estas condiciones el problema se reduce al caso del potencial coulombiano creado
por una carga dq(t)=(t)dV, es decir,

1  t dV
d  .
4 o R

325
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

 R
dV´ t 
 t dV´  R  c
Se ve que F1 t   y por lo tanto F1  t    . El potencial
4 o  c 4 o
escalar d toma así la forma

 R
  r´, t  dV ´
d 
1  c
. (7.1.4)
4 R

donde, al escribir los argumentos de los cuales depende la densidad de carga, se ha

retomado el hecho que el potencial d está determinado por las cargas situadas en la
vecindad del punto con radio vector r´.
De acuerdo con el principio de superposición, el potencial escalar total del sistema
se puede obtener integrando el anterior resultado con respecto al volumen del sistema. Se
obtiene:

 R
  r´,t  dV ´
1  c
(r , t ) 
4 0 
v
R
, R=r-r´. (7.1.5)

Esta expresión se puede escribir en forma abreviada como

1  t dV´
(r ,t ) 
4 0 
v
R
, (7.1.6)

donde los corchetes indican que el valor de la densidad de carga  debe tomarse en el
instante t-R/c.
La solución de la ecuación (7.1.2) para el potencial vectorial se obtiene fácilmente si
1
se hace el reemplazo    0 j . Se obtiene así
0

0 jt 
Ar ,t    dV´ . (7.1.7)
4 v
R

Las relaciones (7.1.6) y (7.1.7) determinan, respectivamente, los potenciales escalar

y vectorial creados por una distribución arbitraria de cargas y corrientes que varían en el
tiempo. Estos potenciales se denominan potenciales retardados.

326
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

Con base en estas expresiones se puede responder la siguiente pregunta: ¿Cómo

varían las leyes de Coulomb y de Biot-Savart cuando las distribuciones de carga y corriente
no son estacionarias?
Se examina inicialmente el caso del campo eléctrico, el cual puede calcularse de la
relación E=-grad-A/t. El cálculo de la derivada temporal del potencial vectorial es
A 1 1  j 
inmediato y da  
t 4 0 c 2 R  t 
dV´ . En el cálculo del gradiente del potencial

escalar es preciso tener en cuenta que, de acuerdo con (0.5.12),

    1 1
grad    grad     grad  ,
 R  R R
donde

grad(1/R)=eR/R2

y
    1  1   
grad     grad  t  R     e R .
 t   c  c  t 

Se tiene entonces la siguiente expresión para el campo eléctrico:

     j  
1     
E r , t    t  e   t  dV ´ .

4 0 v  R 2
 e R  R  (7.1.8)
cR c2 R 

 


El primer término de esta expresión se reduce a la conocida ley de Coulomb en

ausencia de efectos de retardo. Los demás términos corresponden al aporte asociado al
carácter no estacionario de la distribución de cargas y corrientes. Nótese que a grandes
distancias el campo eléctrico E difiere considerablemente de su carácter colombiano: en
1   
efecto, en la expresión subintegral de la relación (7.1.8), los términos eR y
cR  t 
1  j 
 2   , asociados a la evolución temporal de las distribuciones de carga y corriente,
c R  t 
son inversamente proporcionales a R, mientras que el término
  e
R , que prevalece en el
R2
caso estacionario, es inversamente proporcional a R2. Por lo tanto, a grandes distancias del
sistema de cargas y corrientes los términos dominantes son los asociados a las derivadas
temporales de las distribuciones de carga y corriente. A pequeñas distancias la situación es
al contrario: el término dominante es el coulombiano, lo cual es plausible intuitivamente,
ya que a pequeñas distancias del sistema los efectos de retardo son despreciables.

327
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

Examínese ahora el comportamiento de la inducción magnética B=rotA. En este

caso es preciso calcular

  j 1 1
rot   rot j    j   grad   ,
R R R
donde
 j   1  1  j 
rot j       grad  t  R      eR .
 t   c  c  t 

Así se obtiene la siguiente generalización de la ley de Biot-Savart:

  j  
  e 
0   j   e R  t  R
dV ´ .
Br , t  
4   R 2
 (7.1.9)
cR 
 
 

Se ve que a pequeñas distancias del sistema, el término dominante en la expresión

j e R
subintegral de (7.1.9) es y la inducción magnética toma la forma
R2
 j  eR
Br , t   0  dV ´ , la cual el corresponde a la ley de Biot-Savart (4.3.1).
4 R2
Sin embargo, a grandes distancias, donde los efectos de retardo son considerables, el
1  j 
término dominante en la expresión subintegral de (7.1.9) es  eR , el cual está
cR  t 
asociado al carácter no estacionario de la distribución de corrientes.
Así, a grandes distancias del sistema de cargas y corrientes los términos dominantes
deben su origen físico a los efectos de retardo en la propagación de los campos
electromagnéticos. En esta región, denominada zona de radiación, aparecen los efectos
físicos novedosos. Considérese con más detalle la estructura de los campos E y B en
1 1 1 1 R r  r´ r
dicha zona de radiación. Se puede escribir    , eR     er
R r  r´ r r R r  r´ r
y así,

     j        j  

1     t   t  
 
  t   t  

E r , t  
1
  2 R
4 0 v  R
e 
cR
e R  2
c R
dV ´
4 0 v  cR R c 2 R dV ´
e 
 

 
 
 

     j        j  

1   t   t  
 1 
  t   t  

 
4 0 v  c r  r´
e R 
c r  r´ 
2 dV ´
4 0 v  c r r c 2 r dV ´
e 
 

 
 
 


328
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

1 er 1    1 1 1  j   c       1  j 
  
4 0 r v c  t 
dV ´   2  
4 0 r v c  t 
dV ´ 0     dV ´er  0    dV ´ ,
4 r  v  t   4 r v  t 

  j     j  
  eR      eR 

Br , t   0  
 j  eR   t  dV ´ 0   t  dV ´
 
4  R 2 cR  4   cR 
   
   
  j     j  
    eR      er 
0   t  dV ´ 0   t  dV ´  0 1  j  dV ´  e .

4  c r  r´  
4  c r   4 cr   t 
  
 r

   
   

Es decir, en la zona de radiación los campos tienen la forma

0 c       1  j 
E r , t   
   dV ´er  0    dV ´ ,
4 r  v  t   4 r v  t 

  1  j  
Br , t    0  dV ´   er ,

 4 cr V  t 

Si se multiplica vectorialmente por er la expresión para el campo eléctrico se tiene

0 c        1 j    1 j 
E r , t   er   dV ´ er  er    0    dV ´  er    0    dV ´  er
4  
r  v  t    
 4 r  t    
 4 r  t  
  v   v 

Teniendo en cuenta que la última expresión de la derecha difiere en un factor c de la

expresión para la inducción magnética B en la zona de radiación, se obtiene entonces

er  E r , t   cBr , t  ,

Es decir, B y E son mutuamente perpendiculares. Además, de la expresión para la

inducción magnética B en la zona de radiación se nota que B es perpendicular a la
dirección er, a lo largo de la cual se transporta la energía electromagnética.
Es preciso reiterar que en la zona de radiación las magnitudes de los campos E y B son
inversamente proporcionales a la distancia r. Esto lleva a que en la zona de radiación se
defina un vector de Poynting S=EH cuya magnitud es inversamente proporcional al
cuadrado de la distancia, es decir, S1/r2. Recordando que la magnitud del vector de
Poynting es la energía electromagnética que fluye en unidad de tiempo a través de una
superficie transversal de área unitaria, por lo tanto la energía Sdf que fluye en unidad de

329
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

tiempo a través de un elemento de superficie df=r2der (siendo d el elemento de ángulo

sólido) será diferente de cero e independiente de la distancia al sistema:

1 

S  df   2 er   r 2 der  d . 
r 
Este resultado es de capital importancia, ya que indica que la potencia de la energía
electromagnética emitida por un sistema de cargas es diferente de cero en la zona de
radiación, independientemente de la distancia a la fuente.
En las siguientes secciones se estudiará con detalle la radiación de diversos sistemas
físicos. La metodología en todos los casos será la misma: se determina de manera
aproximada el potencial vectorial, se calcula la inducción magnética, el campo eléctrico y el
vector de Poynting. De allí se discute la distribución angular de la radiación emitida en
unidad de tiempo (potencia de radiación) y se determina la potencia total radiada por el
sistema.

§ 7.2. RADIACIÓN DIPOLAR.

Considérese un sistema no estacionario de cargas y corrientes distribuidos en una región

limitada del espacio. Sea d- el tamaño lineal característico del sistema y sea r la distancia
del punto P al origen de coordenadas. Se supone que el tamaño del sistema es tal, que tiene
lugar la desigualdad d<<r, lo cual nos permitirá despreciar el retardo de las ondas
electromagnéticas dentro del sistema. Se considerará además que la longitud  de la onda
radiada por el sistema es mucho mayor que d, es decir, >>d. Esta desigualdad implica que
los campos se pueden considerar homogéneos en la región interior del sistema. Si se tiene
en cuenta que la longitud de onda está relacionada con el período T de oscilación mediante
la expresión =cT, entonces la condición >>d es equivalente a v<<c, donde v=d/T es la
velocidad característica de una carga del sistema. En otras palabras, la condición >>d es
equivalente a considerar que el movimiento de las cargas es no relativista.
Se examinan los campos creados por este sistema en la zona de radiación << r. Esto
quiere decir que en el punto P prevalecen los campos electromagnéticos en forma de onda
plana y por lo tanto es suficiente determinar de manera adecuada el potencial vectorial A.
Se parte de la expresión exacta para el potencial vectorial retardado

 R
j  r´, t  
0  c
Ar , t    dV ´ , R=r-r´.
4 v
R

En el denominador de la expresión para el potencial retardado se puede tomar Rr.

En cuanto al argumento temporal de la densidad de corriente, éste se puede representar
aproximadamente como t-r/c. Se tiene así

330
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

0
Ar , t   j r´, t´dV ´ ,
4r v
(7.2.1)

donde se ha introducido el tiempo de retardo

r
t´ t  . (7.2.2)
c

Para entender el significado de la integral  jr´,t´dV ´ , se supone que el sistema

v

consta de un número finito de cargas puntuales. La densidad de corriente toma entonces la

forma

j(r´,t´)=(r´,t´)v(r´,t´)=iqi(r´-ri,t´)v(r´,t´)

y por lo tanto se obtiene

 jr´,t´dV ´   q  r´r vr´,t´dV ´   q   r´r vr´,t´dV ´ 

v i v
i i
i
i
v
i

 q   r´r vr´,t´dV ´   q vr , t´  q v t´  dt´  q r t´ .

d
i i i i i i i i
i v i i i

Por otro lado, para un sistema que consta de un número finito de cargas puntuales su
momento dipolar p es

pt '   r´ (r ' , t ' )dV ´    qi r´ r´ri dV ´   qi  r´ r´ri dV ´   qi ri t ' .
V i v i v i

Por lo tanto, en la aproximación considerada

dpt´
 jr´,t´dV ´
v
dt´
.

De esta forma, el potencial vectorial se representa a través de la primera derivada

temporal (calculada en el instante retardado t´) del momento dipolar del sistema:

 0 dpt´ 
A(r , t )  . (7.2.3)
4r dt´

Conociendo el potencial vectorial, se puede determinar la inducción magnética B

mediante la relación B=rotA. Se tiene:

331
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

 0  1 dpt´    0   1  dpt´  1 dpt´    0  dpt´  

B rot  = grad     rot  rot  .
4  r dt´  4   r  dt´ r dt´  4r  dt´ 

En esta cadena de relaciones se ha usado la identidad rot(a)=(rota)+(grad)a

1 1
(véase (0.5.13)) y se ha tenido en cuenta que    2 por lo que a grandes distancias
r r
del sistema se puede despreciar el primer término de la segunda igualdad en comparación
con el segundo (que es inversamente proporcional a r).
 dpt´  
Se debe ahora calcular rot  . Se tiene:
 dt´ 

 dpt´   dpt´ d t´ dpt´ d 2 pt´   r 

rot       t   ,
 dt´  r dt´ dt´ r dt´ dt´2 r  c 

es decir,

 dpt´  1 d pt´
2
rot   er,
 dt´  c dt´
2

y por lo tanto, la inducción magnética B toma la forma

 0 d 2 pt´ 
B  er . (7.2.4)
4rc dt´ 2

Conociendo la inducción magnética, se puede determinar el campo eléctrico

mediante la relación E=c[Ber]. Se obtiene entonces

  d 2 pt´ 
E   0  er   er . (7.2.5)
 4r dt´
2


Si se orienta el momento dipolar del sistema a lo largo del eje z (p=pez), las
anteriores expresiones se pueden escribir de la siguiente forma

 0 d 2 pt´
B sen e  ,
4rc dt´2
(7.2.6)
 0 d pt´
2
E sen e ,
4r dt´2

donde  es el ángulo formado por el radio-vector r y el momento dipolar p (Figura (7.2.1)).

332
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

Nótese que los campos de radiación están determinados por la segunda derivada del
momento dipolar.

Figura (7.2.1)

La parte del vector de Poynting que describe la radiación es:

 0  d 2 pt´   0  d 2 pt´ 
2 2

S=EH= 2 
 
 sen 2
  e  e   2 
  sen 2 e r . (7.2.7)
4r  c  dt´ 
2  
4r  c  dt´ 
2

Es de destacar que la magnitud de este vector de Poynting es inversamente

proporcional al cuadrado de la distancia entre el dipolo y el punto de observación P.
Considérese un elemento de área esférica con centro en el origen de coordenadas y radio r.
Este elemento se puede representar como df=r2der, donde d=sendd es el elemento
de ángulo sólido. Si se multiplica escalarmente el vector de Poynting por el elemento de
área df se obtiene entonces la potencia radiada (cantidad de energía en unidad de tiempo)
que pasa a través de la superficie transversal df: dP=Sdf=Sr2d. Por lo tanto la potencia
radiada en unidad de ángulo sólido por el sistema considerado en la aproximación dipolar
es

 0  d 2 pt´ 
2
dP
 2 
  sen 2 . (7.2.8)
d 4  c  dt´ 
2

En la Figura (7.2.2) se puede ver la representación polar de esta dependencia.

333
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

Figura (7.2.2)

Se ve que la potencia radiada es máxima en el plano ecuatorial (=/2) y se anula a lo largo

de la dirección del momento dipolar (=0)
Si se integra la expresión (7.2.8) sobre una superficie esférica de radio r, se obtiene
la potencia total radiada o intensidad total de la radiación:

0  d 2 pt´ 
2 2
P=Sdf=    sen  sen d  d ;
2

4 2 c  dt´2  0 0

la integral con respecto al ángulo polar  se calcula fácilmente introduciendo la variable


 
1
2
auxiliar u=cos. Se tiene entonces  sen 2 sen d   u 2 du  y por lo tanto la potencia
0 1
3
total radiada por el sistema de cargas cuyo momento dipolar evoluciona temporalmente es

  d 2 pt´ 
2

P  0   (7.2.9)
6c  dt´2 

En el instante t la potencia de la radiación depende del valor de d2p/dt2 en el instante

t-r/c y sólo a través de esta combinación dicha potencia depende de la distancia r al sistema.
Estos resultados se pueden ilustrar con algunos ejemplos concretos.
Radiación asociada al movimiento de una partícula no relativista. Considérese
una partícula cargada no relativista de masa m, la cual es acelerada bajo la acción de una
fuerza F. El momento dipolar de tal partícula es p=qr´ y por lo tanto d2p/dt´2=q(d2r´/dt´2)
=qa, donde a es la aceleración de la partícula. La potencia total radiada es igual a

  d 2 pt´ 
2
0 0 q 2 2
P  0    qa 2
 a t´ .
6c  dt´2  6c 6c

Se ve que P se determina por el cuadrado de la carga q (lo cual indica que la

potencia de la radiación, así como su intensidad, son independientes del signo de la carga) y
por el cuadrado de la aceleración. En particular, una carga que se mueva con velocidad

334
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

constante no puede emitir radiación. Este último resultado es consecuencia obvia del
principio de relatividad, ya que si una partícula se mueve con velocidad constante, entonces
se puede elegir un sistema de referencia con respecto al cual la partícula todo el tiempo está
en reposo. En tal sistema la partícula no puede emitir radiación ya que ello viola la ley de
conservación de la energía. Como tal resultado debe tener lugar en todos los sistemas
inerciales de referencia, se concluye que una partícula que se mueva con velocidad
constante no puede emitir radiación.
Radiación de una carga que se mueve en un campo magnético uniforme. Si una
partícula no relativista de masa m y carga q se mueve en un campo magnético uniforme,
ésta se encuentra sometida a la fuerza qvB y por lo tanto su aceleración será (q/m)vB. La
correspondiente potencia total radiada tendrá la forma

0 q 4
P v t´  B 2 .
6m c2

Dado que la potencia en este caso depende del cuadrado de la velocidad de la

partícula, esta radiación es demasiado pequeña en comparación con la del ejemplo anterior.
Radiación de una carga que oscila. Considérese una carga que efectúa
oscilaciones según la ley r=r0cos(t+). La aceleración de la carga será entonces
a=-2r0cos(t+) y por lo tanto la potencia total radiada será

0q 2 4 2 2
P  r0 cos t´   .
6c

Si se promedia con respecto al tiempo, se obtiene

0q 2 4 2
P   r0 .
12c

Se ve así que la potencia de la radiación emitida por un oscilador cargado es

proporcional a la cuarta potencia de la frecuencia de oscilación y al cuadrado de la amplitud
de las oscilaciones.
Radiación de un sistema compuesto por dos partículas. Considérese un sistema
compuesto por dos partículas de cargas q1, q2 y masas m1, m2. El momento dipolar de tal
sistema es p=q1r1+q2r2 y por lo tanto d2p/dt2=q1a1+q2a2. Si este es un sistema cerrado,
entonces la única fuerza F presente es la de interacción entre las partículas y las
aceleraciones toman la forma a1=F/m1, a2=-F/m2. Por lo tanto

d 2 p  q1 q 2 
  F
dt 2  m1 m2 

y la potencia total radiada toma la forma

335
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

2
 q q 
P  0  1  2  F 2 .
6c  m1 m2 

De aquí se obtiene un resultado fundamental: todo sistema que consta de dos

partículas iguales, o de dos partículas diferentes pero con igual relación carga-masa, no
puede radiar dentro de la aproximación dipolar. Para la radiación de tales sistemas es
preciso considerar aproximaciones de orden superior.
Radiación de una antena dipolar. Examínese el sistema representado en la Figura
(7.2.3). Este sistema se denomina antena dipolar. El dipolo oscilante se obtiene mediante
una corriente variable I=I0cos(t+) que fluye a lo largo del eje z. La respectiva carga
oscilante es q=Idt=(I0/)sen(t+) y por lo tanto el momento dipolar es
p=ez(I0l/)sen(t+) y su segunda derivada con respecto al tiempo toma la forma
d2p/dt2=-ez(I0l)sen(t+). Reemplazando en (7.2.9) se tiene

  d 2 pt´ 
2
  c  l 
2

P  0    0  lI 0sent´ 2  0   I 02sen 2 t´  .

6c  dt´ 
2
6c 6  c 

Al promediar con respecto al tiempo, se tiene sen 2 t´   1/ 2 , lo que lleva a

que la potencia media radiada por el sistema sea

0c  l  I 02
2

P    .
6  c  2

Figura (7.2.3)

336
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

Recuérdese que en un circuito de resistencia R por el que circula una corriente

variable
I=I0cos(t+), la potencia media disipada está dada por la fórmula de Joule-Lenz
I2
P  R 0 . Se ve así que la antena dipolar se caracteriza por una resistencia de radiación
2

 c  l  2 0  l 
2 2 2 2
2 l  l
Rrad  0       80    789   . (7.2.10)
6  c  3 0     

Dado que estos resultados son válidos en la aproximación dipolar, para la cual l<<,
la resistencia de radiación de la antena considerada es pequeña y por lo tanto la potencia
radiada también será pequeña. Por esta razón para aumentar la eficiencia de una antena se
busca que sus dimensiones sean del orden de la longitud de onda. Es claro que en este caso
pierde validez la aproximación dipolar. Una discusión más detallada se presenta en la
sección §7.5.

§ 7.3. RADIACIÓN CUADRUPOLAR Y MAGNETODIPOLAR.

Examínese de nuevo el sistema no estacionario de cargas y corrientes considerado en la

sección anterior. Se supone otra vez que el tamaño lineal característico del sistema d
cumple la condición d<<<<r donde r es la distancia del punto de observación P al origen
de coordenadas y  es la longitud de la onda radiada.
Se parte de la expresión exacta para el potencial vectorial retardado:

 R
j  r´, t  
0  c
Ar , t    dV ´ , R=r-r´.
4 v
R

y se toma en el denominador Rr. En el argumento de la densidad de corriente, a diferencia

de la aproximación dipolar, se determina con mayor exactitud el tiempo de retardo. Para
ello, téngase en cuenta que la condición d<<r permite expresar la magnitud R de la
siguiente forma:

R=r-r´=(r2-2rr´+r´2)1/2r-nr´, ner=r/r,

y por lo tanto
 R  n  r´  n  r´ j r´, t´
j  r´, t    j  r´, t´   j r´, t´  ,
 c  c  c t´
r
donde t´ t  .
c
El potencial vectorial es entonces

337
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

0  j r´, t´ 
A r , t     j r´, t´dV ´   n  r´
1
 dV ´ . (7.3.1)
4r  v cv t´ 

En la sección anterior se demostró que la integral  j r ,t´ dV´ no es otra cosa que la
v

derivada temporal del momento dipolar del sistema en el instante t´. Para ello se supuso que
el sistema consta de un número finito de cargas puntuales. Repitiendo este mismo
procedimiento con la segunda integral de la expresión (7.3.1), se obtiene

j r´, t´ 
  n  r´ dV ´  qa v a  n  ra  . (7.3.2)
v
t´ t´

La expresión va(nra) se transforma de la siguiente forma:

1 ra
v a  n  ra   v a  n  ra   v  n  ra    n  ra   v a  n  ra  ,
1 1 1
2 t´
a
2 2 2

o bien,

1  1  r 
va  n  ra   ra  n  ra   1 v a  n  ra   ra  n  a  =
2 t´ 2 2  t´ 
1 
 ra  n  ra   1 ra  v a  n .
2 t´ 2
Entonces (7.3.2) toma la forma

j r´, t´ 2 1 
  n  r´  q r  n  r   2 t´  q r  va  n .
1
dV ´ (7.3.3)
t´ 2 t´2
a a a a a
v

1 
q a ra  v a   n     r  v dV´   n  mn, donde m es el momento
1
Pero
2
 2 
magnético del sistema en el instante t´. De manera análoga, para la i-ésima componente de
la primera integral de la izquierda de (7.3.3) se obtiene

1 1 
2
 q a x ai n j x aj    x´ i x´ j dV´ n j .
2 

Por otro lado, de la definición del tensor cuadrupolar

Qik   3x' i x' k r'  ik ( r' )dV´ (véase (1.6.4)) se tiene que
2

V´

338
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

1 
 x' x' k dV´  Qik  ik  r'
2
i dV´
V´
3 3 V´

y por lo tanto

1 1 1 
2
 q a xai n j xaj  Qij n j    r ' 2 dV ni .
6 6 

Reemplazando estos resultados en la i-ésima componente de (7.3.3), se obtiene

ji r´, t´  mt´  1  Qij t´  1 2 

2

  n  r´ dV ´   n     r ' dV ni .

2
n
t´  t´  i 6 t´  6 t´
2 j 2
v 

Nótese que el último término de la derecha en esta expresión es proporcional al

vector unitario n por lo que dicho término no influye en las expresiones para los campos E
y B y por lo tanto se puede eliminar de la expresión para el potencial vectorial. Teniendo
esto en cuenta, la componente i-ésima de (7.3.1) queda

 0  pi t´  1  mt´   1  Qij t´  

2

Ai r , t      n  nj . (7.3.4)
4r  t´ c  t´  i 6c t´
2 

Con base en esta relación, la determinación de los campos E y B en la zona de radiación es
fácil.
Se ve que el potencial vectorial se compone de varios términos: El primero es
proporcional a la primera derivada temporal del momento dipolar y fue discutido con
detalle en la sección anterior. El segundo término depende de la primera derivada con
respecto al tiempo del momento dipolar magnético y recibe el nombre de término magneto-
dipolar. El tercer término está determinado por la segunda derivada del tensor cuadrupolar y
corresponde a la radiación cuadrupolar.
Para un sistema que consta de dos partículas cargadas la radiación magnetodipolar
es igual a cero: en efecto, en este caso se conserva el momento angular del sistema y por lo
tanto su derivada temporal se anula, lo que lleva a su vez a que se anule la derivada
temporal del momento magnético en virtud de la relación giromagnética. Sin embargo, la
radiación dipolar está presente para este sistema. Así mismo, si un sistema consta de
partículas con la misma relación carga masa, tanto la radiación dipolar como la
magnetodipolar son iguales a cero.
Se consideran ahora algunos sistemas concretos en los cuales se presentan las
radiaciones cuadrupolar y magnetodipolar.
Radiación de un cuadrupolo eléctrico lineal. Examínese la radiación producida
por un cuadrupolo eléctrico lineal formado por dos dipolos de polaridades opuestas, tal
como se muestra en la Figura (7.3.1). El momento cuadrupolar de este sistema es una
función periódica (sinusoidal) del tiempo.

339
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

Figura (7.3.1)

Se tienen así dos dipolos opuestos, uno centrado en z=-l/2 y otro centrado en z=l/2.
Por lo tanto el momento dipolar total de este sistema (al despreciar los efectos de retardo
dentro de éste) es nulo y en la zona de radiación (r>>) es preciso tomar en cuenta el aporte
asociado al tensor cuadrupolar. Para el sistema de cargas considerado la densidad
volumétrica de carga es  r´,t´  q z´l   q z´l   2q z´ x´  y´ ; dada la simetría,
el momento cuadrupolar es

Qij t´   2qt´  3l i l j  l 2 ij 
(compárese con el resultado (1.5.11) obtenido para un cuadrupolo lineal en la aproximación
estática). Por lo tanto el potencial vectorial toma la forma

 0  2 qt´ 
Ai r , t  
12cr t´ 2
3l i l j  l 2 ij n j ,

es decir,

 0  2 qt´ 
A r , t  
12cr t´ 2
3l  nl  l 2 n , (7.3.5)

y por lo tanto la inducción magnética B en el punto P es

0  1  2 qt´ 
B  rotA r , t   rot 
12c  r t´ 2
3l  er l  l 2er  .


 2 qt´
Usando la relación (0.5.13) rot(a)=(rota)+(grad)a con a 
t´2
 
3l  er l  l 2er y
1
  se obtiene
r

340
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

 1  2 qt´  1   2 qt´ 
rot  3l  e l  l 2

e r  rot  
3l  er l  l 2er   
 r t´  r  t´
2 r 2

 1   1  qt´ 
 
2
grad     3l  er l  l 2er 
 r   r t´
2

1 e
Dado que grad     r2 , el segundo término de la derecha en la última expresión se
r r
puede despreciar en la zona de radiación. Recordando, adicionalmente, que en la zona de
 t´  1 
radiación se tiene    er , entonces
r r t´ c t´

 1  2 qt´  1   2 qt´  1    2 qt´ 

rot  3l  e l  l 2

e r   rot   3l  e l  l 
2
e r    
3l  er l  l 2er 
 r t´  r  t´  r r  t´
2 r 2 r 2

1    qt´     qt´ 
3l  er l  l 2er    3l  er l  l 2er  .
2 2
11
 er  
r r  t´ 2
 rc t´  t´2


0  3qt´
Por lo tanto, B  l  er l  er  , o bien,
4c 2 r t´3

 0 l 2  3 qt´
B cossene . (7.3.6a)
4c 2 r t´3

El campo eléctrico de esta radiación se determina fácilmente de la relación

E=cBer. Se tiene así
 l 2  3 qt´
E 0 cossene . (7.3.6b)
4cr t´3

Se ve que los vectores E y B para la radiación del cuadrupolo eléctrico están

orientados de la misma manera que los correspondientes a la radiación del dipolo eléctrico.
Ello se debe a que los campos cuadrupolares se obtienen de la superposición de dipolos
eléctricos. Se recalca, sin embargo, dos diferencias fundamentales entre las radiaciones
dipolar y cuadrupolar: en la primera los campos de radiación están determinados por la
segunda derivada temporal de la carga y por lo tanto aumentan con el cuadrado de la
frecuencia de la radiación, mientras que en la radiación cuadrupolar los campos de
radiación están determinados por 3q/t´3 y por lo tanto aumentan como 3. Además, en la
radiación dipolar los campos de radiación se anulan a lo largo del eje polar (=0, ),
mientras que los campos de la radiación cuadrupolar se anulan tanto en los polos como en
el ecuador (=/2). Este último resultado es obvio, dado que en el plano ecuatorial los
dipolos crean campos de radiación de igual magnitud pero de direcciones opuestas.
Supóngase que la variación temporal del momento cuadrupolar del sistema
considerado se debe a que una corriente sinusoidal I(t´)=I0sen(t´+) actúa sobre el
sistema. Entonces 3q(t´)/t´3 =-2I(t´), de donde se obtiene para el vector de Poynting

341
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

1 1  0l 2  3qt´
2
1 
S EB  3 
cossen  e  e 
0 0 c  4r t´ 3

2
1 1  0 2l 2  0c  l  2
4

 I t´cossen  er  2 2   I t´cos sen er 

2 2
3 
0 c  4r  4  r  c 
0c  l  2 2
4

2 2   I 0 sen t´ cos sen er .

2 2

4  r  c 
Por lo tanto, el promedio temporal del vector de Poynting es

0c  l  I 02
4

S  2 2   sen 2 cos 2 er
4  r  c  2
y la potencia media radiada en unidad de ángulo sólido toma la forma

 c  l  I 2
4
dP
 S  r 2er  0 2   0 sen 2 cos 2 . (7.3.7)
d 4   c  2
El respectivo patrón de radiación se ilustra en la Figura (7.3.2). Nótese que el ángulo
de máxima emisión de radiación para el cuadrupolo lineal es 45º.

Figura (7.3.2)

La potencia total radiada por el sistema es entonces

0c  l  2 
4

P    I 0  sen  cos  d .
3 2

16  c  0

342
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

La integral es igual a 4/15. Recordando que 0c=(0/0)1/2, se puede escribir el

anterior resultado de la siguiente forma:

0  l  I 02
4
1
P    . (7.3.8)
30 0  c  2

Se ve que la resistencia de radiación del cuadrupolo considerado es

0  l   l 
4 4
1
Rrad     4.0  Ohm .
30 0  c   c 

Compárese la magnitud de la radiación cuadrupolar con la de la radiación dipolar. Se

recuerda que para esta última la potencia de radiación es

 0  l  I 02
2
1
P    .
6 0  c  2

Por lo tanto, la razón entre las potencias de la radiación cuadrupolar y dipolar es del orden
 l 
2

de   . Esto indica que en un sistema de cargas confinado en una región del espacio de
 c 
tamaño d<< la radiación cuadrupolar, en caso de estar presente, es despreciable si el
sistema tiene un momento dipolar diferente de cero.
Radiación de una espira circular. Considérese ahora una espira circular de radio a
que descansa en el plano xy, la cual está alimentada por una corriente sinusoidal I=I0
sen(t+). tal como se muestra en la Figura (7.3.3). Se propone discutir la radiación de este
sistema en un punto P situado a una distancia r del centro de la espira. El sistema descrito
se caracteriza por tener un momento dipolar magnético m=a2Iez.

343
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

Figura (7.3.3)

Para el potencial vectorial en el punto P se tiene:

0 mt´  a 2 dI t´  a 2 dI t´

A r , t    er  0 e z  er  0 e ,
4cr t´ 4cr dt´ 4cr dt´

donde t´=t-r/c, y r es la distancia entre el centro de la espira y el punto P.

De aquí se puede determinar fácilmente la inducción magnética

0 a 2  1 dI t´  0a  dI t´  0a d I t´

e  er  ,
2 2 2
B  rotA r , t   rot e   rot e   2
 r dt´  4cr  dt´  4c r dt´
2
4c
es decir,

 0 a 2 2 I 0
B sent´ sene . (7.3.9a)
4c 2 r

El correspondiente campo eléctrico E=cBer toma la forma

 0 a 2 2 I 0
E  sin t´ sine . (7.3.9b)
4cr

Nótese que la estructura de los campos de radiación dipolares magnéticos es similar

a la de los campos de radiación del dipolo eléctrico oscilante, pero con las direcciones de
los campos E y B intercambiadas.

344
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

Siguiendo el procedimiento estándar descrito con detalle para los casos de la

radiación dipolar y cuadrupolar, se puede obtener la potencia media radiada y la resistencia
de radiación de esta antena. Los respectivos cálculos se proponen como ejercicio. Los
resultados son

 c  a  2 2  a 
4 4

P  0  2  I 0  10.0   I 0 Watt
2
(7.3.10)
12  c   c 

para la potencia total P radiada por el sistema y

 a   a 
4 4

Rrad  20.0    197  

2

 c   c 

para la resistencia de radiación. Esta resistencia es proporcional a la cuarta potencia de la

frecuencia, resultado que contrasta con el dipolo oscilante, para el cual su resistencia es
proporcional a la segunda potencia de . Comparando con los resultados obtenidos para un
dipolo oscilante, se ve que la espira considerada es una antena menos eficiente.

§ 7.4. RADIACIÓN DE UNA ANTENA DE MEDIA ONDA.

En todos los sistemas descritos en las dos secciones anteriores se consideró que el tamaño
característico del sistema es mucho menor que la longitud de onda de la onda
electromagnética radiada. Se examina ahora un sistema en el cual tal restricción no tiene
lugar.
Considérese un conductor rectilíneo por el cual circula una corriente no homogénea
de la forma

I  I 0 cost cos(kz´) , (7.4.1)

donde k=/c es el número de onda. La respectiva distribución de carga en el sistema tiene

la forma q  I 0 /  sint cos(kz´ ) .
Esta distribución de corriente se logra en un dispositivo como el que se muestra en
la Figura (7.4.1): un conductor rectilíneo de longitud 2l=/2 se parte en dos segmentos
iguales y por él se hace circular una corriente variable I. Nótese que en los extremos del
sistema (z=l) la corriente I se anula y es máxima en el centro (z=0). Por este motivo, el
dispositivo descrito se denomina antena de media onda. En muchos casos se diseña un
mástil de longitud l situado en el suelo, que hace las veces de espejo; así, el mástil y su
imagen actúan como antena de media onda.

345
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

Figura (7.4.1)

Para aplicar los resultados obtenidos en la sección §7.2 para la antena dipolar, se
considera un elemento de antena de longitud z´ situado a una distancia z´ del origen. El
I 
momento dipolar asociado a este elemento es p= e z  0 sen t cos(kz´) z´. Por lo tanto, la
 
inducción magnética en el punto P es, de acuerdo con la fórmula (7.2.4),

0 d 2pt´  I
B   er   0 0 sent´coskz´z´ez  er .
4rc dt´ 2
4rc

donde t´=t-r´/c. Dado que r´ depende de la posición del elemento dipolar considerado, el
tiempo de retardo es función de z´. La inducción magnética total en P se obtiene al integrar
esta relación con respecto a z´.

 0I 0 l
  r´  
B = e z  e r  sen   t   cos(kz´)dz´ .
4rc l   c 

De acuerdo con la Figura (7.4.1), r´r-z´cos(). Por lo tanto, se puede escribir

 /2
0 I 0   r 
B = sine  sen   t    ucos cosu du ,
4r  / 2   c 

346
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

donde se ha introducido una nueva variable de integración adimensional u=kz´. Teniendo

en cuenta que

  r    r    r 
sen   t    ucos   sen   t  cosucos   senucos cos   t   ,
  c    c    c 

la expresión para la inducción magnética B toma la siguiente forma:

 /2
0 I 0   r 
B sine sen   t    cosucos cosu du  ,
4r   c    / 2
 /2
0 I 0   r 
 sine cos   t    senucos cosu du .
4r   c    / 2
 /2
Pero  senucos cosu du  0 , dado que la expresión subintegral
 / 2
senucos cosu

es una función impar de u y la integración se efectúa sobre límites simétricos. Por otro lado,
para la integral del término que no se anula en la expresión anterior de B se tiene:

 
 / 2l cos cos 
cos ucos cosu du  2  2 
2

 / 2
sen 

Así, la inducción magnética B creada por la antena considerada en el punto P es

  r    
cos   t   cos cos 
 I   c   2 
B = 0 0 e . (7.4.2a)
2r sen

El respectivo campo eléctrico E=cBer es

  r    
cos   t   cos cos 
 I   c   2 
E=  c 0 0 e . (7.4.2b)
2r sen

El vector de Poynting es

347
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

  r   
cos 2    t   cos 2  cos 
 0 cI
2
0   c  2 
S= E×H = eR . (7.4.3)
2r 
2
sen 2

Multiplicando escalarmente por el elemento de área df=err2sendd y promediando

con respecto al tiempo, se obtiene la potencia media radiada en unidad de ángulo sólido:

 
cos 2  cos  2
dP  c 2  I0 .
 0 2 (7.4.4)
d 2  sen 
2
2

El correspondiente diagrama de radiación se

presenta en la Figura (7.4.2) y se compara con el patrón de radiación de un dipolo elemental
(curva a trazos). El diagrama cualitativamente es semejante al de un dipolo, lo que se debe a
que este patrón se obtuvo como resultado de la superposición de dipolos elementales entre
los cuales la diferencia de fase es pequeña a lo largo de la antena. La diferencia básica entre
ambos sistemas radica en que la antena de media onda emite una mayor cantidad de
radiación en la región ecuatorial.

Figura (7.4.2)

Integrando (7.4.4) con respecto a las variables angulares, se obtiene la potencia

media radiada
2 
2  cos  cos 
 cI 2  d .
P  0 0 
2 2 0 sen

El valor de la integral es 1.21883 y para la resistencia de radiación de la antena de

media onda se obtiene

Rrad=73.13 . (7.4.5)

348
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

Arreglos de antenas. Los resultados obtenidos se pueden extender a sistemas más

complejos, que incluyen arreglos o agrupaciones de antenas. Para ilustrar la respectiva
metodología, se consideran dos antenas de media onda paralelas y separadas una distancia
d, tal como se muestra en la Figura (7.4.3). De manera intuitiva se puede anticipar los
siguientes resultados: si las antenas están excitadas en fase, el patrón de radiación será
equivalente al de un sistema con un momento dipolar neto diferente de cero. Si, al
contrario, las antenas están excitadas en oposición de fase, el momento dipolar total será
nulo y el respectivo patrón de radiación será equivalente al de un cuadrupolo oscilante.

Figura (7.4.3)

De acuerdo con el resultado (7.4.2a), las inducciones magnéticas creadas por cada
una de estas antenas en el punto P son

   
cos cos  cos cos 
I 2  cos   t  r1  e , B   0 I 0 2  cos   t  r2  e ,
B1   0 0        
2r sen 2r sen c  
2
  c   

Si  es el ángulo formado por el radio vector r y el eje x (Figura (7.4.3)), entonces

r1r-(d/2)cos , r2r+(d/2)cos.

Si se tiene en cuenta que er=sencosex+senseney+cosez, y que, por la definición

del producto escalar, cos= erex, se ve que cos=sencos. Por lo tanto

349
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

d d
r1  r  sencos , r2  r  sencos .
2 2

y por lo tanto

  r    r  d   kd 
cos   t  1, 2    cos   t    sencos    cos   sencos   ,
  c    c  2c   2 

 r
donde se ha introducido la fase     t   . Por lo tanto, las inducciones magnéticas de
 c
las antenas en el punto considerado P son:

 
cos cos 
0 I 0 2  cos   kd sencos  e ,
B1      
2r sen  2 
 
cos cos 
I 2  cos   kd sencos  e ,
B2   0 0    
2r sen  2 

La inducción magnética total B en el punto P es la superposición de los campos B1 y

B2, es decir, B=B1B2, donde el signo + se toma cuando las antenas están excitadas en fase
y el signo - cuando están excitadas en oposición de fase.

 
cos cos 
I 2  cos   kd sencos   cos   kd sencos  e ,
B 0 0      
2r sen   2   2 

Si se elige ahora la distancia d entre las antenas de tal forma, que kd= (es decir, d
es igual a media longitud de onda), entonces

    
cos cos  cos  cos sencos  (antenas en fase)
I 2 e  2 
B 0 0  .
r sen   
 sen  sen sencos  (antenas en oposición de fase)

 2 

Si las antenas están en fase, la dependencia angular de la magnitud de los campos de

 
radiación está determinada por el factor cos sencos  . A lo largo del eje x (=/2,
2 
=0) esta función se anula, lo que da una interferencia destructiva, mientras que a lo largo

350
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

del eje y (=/2, =/2) se obtiene interferencia constructiva que lleva a la formación de
un patrón de radiación tipo dipolo.
Un análisis semejante se puede realizar para el caso cuando las antenas están en
oposición de fase. Dicha actividad se propone en calidad de ejercicio.

§ 7.5. RADIACIÓN DIPOLAR A DISTANCIAS ARBITRARIAS DEL SISTEMA.

Al discutir la radiación de longitud de onda  correspondiente a un sistema de tamaño

característico d en la aproximación dipolar (d<<), se consideró que la distancia r era
mucho mayor que la longitud de onda (zona de radiación). Es preciso retomar la discusión
planteada en la sección §7.2 pero sin imponer la restricción <<r. Sin necesidad de realizar
cálculo alguno, se pueden anticipar algunas características esenciales de los campos creados
por un dipolo oscilante. Examínese, por ejemplo, el campo eléctrico. Si el punto de
observación se encuentra en la zona de radiación, el campo eléctrico deberá ser transversal
y tener la forma (7.2.6). Si, por el contrario, el punto de observación se encuentra cerca al
sistema (r<<), se pueden despreciar los efectos de retardo y el campo eléctrico respectivo
será equivalente al de un dipolo en el instante t.
Considérese ahora con más rigor esta situación. A pesar de que ya no se cumple la
condición <<r, la relación

 0 dpt´
A(r , t )  (7.5.1)
4r dt´

sigue siendo válida, ya que para su deducción sólo se usó el hecho de que el tamaño del
sistema es pequeño en comparación con la longitud de onda (véase (7.2.3)).
El potencial escalar  en la aproximación dipolar se puede hallar fácilmente de la
1 
condición de Lorentz divA  2  0 . Reemplazando (7.5.1), se tiene
c  t´

  dpt´  1     0  1 
div 0  2  div pt´  2 0
 4r dt´  c t´ t´   4r  c t´

Es decir, para el potencial vectorial (7.5.1) la condición de Lorentz toma la forma

   0  1 
div pt´  2    0 .
t´   4r  c 

La expresión entre corchetes es, por lo tanto, una magnitud que no evoluciona en el
tiempo; dado que se están considerando campos dinámicos, se puede considerar que
   1
div 0 pt´  2   0 . De aquí se obtiene, entonces, para el potencial escalar la
 4r  c
siguiente expresión:

351
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

1   1 
div pt´   grad  r   pt´  r divpt´ ,
1 1 1

4 0 r  4 0    

y a  pt´ . Si se
1
donde se ha usado la identidad div(a)=(grad)a+diva con  
r
obtiene en cuenta que grad(1/r)=-er/r2 y que

 t´ dpt´ dpt´   r  1 dpt´

divpt´   pt´     t      er ,
r r dt´ dt´ r  c  c dt´

se llega entonces a la siguiente expresión para el potencial escalar en la aproximación

dipolar:

1 1 1 dpt´
  pt´   er .
4 0  r 2
cr dt´ 

Si se dirige el eje z a lo largo de la dirección del momento dipolar, se obtiene

1 1 1 dpt´
  pt´  cos . (7.5.2)
4 0  r 2
cr dt´ 

El primer término en esta expresión tiene la estructura del término dipolar estático.
El segundo término (que es proporcional al cambio temporal del momento dipolar) varía
con el inverso de la primera potencia de la distancia r y es el que está asociado a la
radiación.
Para calcular los campos, se usan las relaciones

A
E  grad   , B = rot A.
t

Determínese inicialmente la inducción magnética B. Se tiene:

 0  1 dpt´  0   1  dpt´ 1  dpt´ 

B rot  grad  r   dt´  r rot dt´  ,
4  r dt´  4     

1 dpt´
donde se ha usado la identidad rot(a)=(rota)+(grad)a con   y a . El
r dt´
término inversamente proporcional a r es el que corresponde al campo de radiación y ya
fue determinado en la sección §7.2 (ec.(7.2.4)). El primer término, que contiene grad(1/r)
se calcula fácilmente. El resultado es entonces

352
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

0  1 dpt´ 1 d 2 pt´
B  2    er .
4  r dt´ rc dt´2 

Como el momento dipolar p está dirigido a lo largo del eje z, se obtiene

0  1 dpt´ 1 d 2 pt´
B  2  sene . (7.5.3)
4  r dt´ rc dt´2 

Obsérvese que la inducción magnética tiene componente azimutal a cualquier

distancia del sistema dipolar. El término inversamente proporcional al cuadrado de r es
dominante en la región cercana y corresponde a la ley de Biot-Savart. En el segundo
término, que es inversamente proporcional a r, se puede reconocer la inducción magnética
en la región de radiación (7.2.6).
Para determinar el campo eléctrico es necesario calcular A/t y grad. Para la
derivada temporal del potencial vectorial se tiene

A A t´  0 d 2 pt´ 1 1 d 2 pt´

  e  cose r - sene  . (7.5.4)
t t´ t 4r dt´2 4 0 c 2 r dt´2
z

En cuanto al gradiente del potencial escalar, es conveniente realizar el respectivo

cálculo en coordenadas esféricas. Se tiene

 1 
grad (r , )  er  e =
r r 
1    1 1 dpt´   1 1 dpt´   
=   2 pt´    cos e   p t´  sen e .
4 0  r  r 
r
cr dt´   r
3
cr dt´ 
2
 

En el cálculo de las derivadas es preciso tener en cuenta que

  t´    r  1 
  t     .
r t´ r t´ r  c  c t´

Se obtiene por lo tanto

1  2 2 dpt´ 1 d 2 pt´ 
grad     3 pt´  2  2 cose r 
4 0  r cr dt´ c r dt´2 
1 1 dpt´  
  3 pt´  2 sene  . (7.5.5)
r cr dt´  

Sumando las expresiones (7.5.4) y (7.5.5), se obtiene para el campo eléctrico el

siguiente resultado:

353
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

1  1 1 dpt´ 
E  3 pt´  2 2cose r 
4 0  r cr dt´ 
1 1 dpt´ 1 d 2 pt´  
  3 pt´  2
  2 sene  . (7.5.6)
c r dt´2 

r cr dt´ 

Nótese que en la componente radial del campo eléctrico desaparece el término que
es inversamente proporcional a r. Así, el término de radiación 1/r aparece sólo en la
componente transversal de E, lo que está en concordancia con el carácter transversal de los
campos de radiación. Los términos que son inversamente proporcionales a r3 corresponden
a los campos de un dipolo estático. Los términos que contienen la dependencia 1/cr2 son
totalmente nuevos y determinan la zona intermedia o región de transición entre la zona
estática (muy cercana al sistema) y la zona de radiación (muy alejada del sistema). Cabe
anotar que esta clase de términos de carácter intermedio no están presentes en la expresión
para la inducción magnética (7.5.3).
Para ilustrar con más detalle los resultados, se supone que el momento dipolar del
sistema varía con el tiempo según una dependencia sinusoidal

p(t´)=p0exp(-it´)= p0exp[i(kr-t)].

En este caso las expresiones (7.5.3) y (7.5.6) toman la forma

p0 k 3  1 i 
Er    2cosexp ikr  t 
4 0  kr kr2 
 3

p0 k 3  1 1
senexp ikr  t 
i
E    
4 0  kr3 kr2 kr 

 0 p0 k 3  i 1
 kr2  kr senexp ikr  t 
B    
4 0  

En la Figura (7.5.1) se muestran las líneas de fuerza eléctricas de un dipolo oscilante

para diferentes valores de la combinación adimensional t. En dicha figura el dipolo está
situado en el centro y orientado en la dirección vertical. Nótese la disminución de la
longitud de onda con la distancia.

354
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

t=0 t=/4

t=3/4
t=/2

355
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

§ 7.6. POTENCIALES DE LIENARD-WIECHERT.

En las secciones anteriores se ha discutido la radiación de diversos sistemas calculando de

manera aproximada los potenciales retardados en la zona de radiación. Existe, sin embargo,
un sistema físico para el cual estos potenciales pueden ser calculados de manera exacta.
Dicho sistema consiste en una carga puntual que se mueve a lo largo de una trayectoria
arbitraria. Las siguientes tres secciones estarán dedicadas al estudio de este sistema.
Considérese entonces una carga puntual q que se mueve a lo largo de una trayectoria
descrita por el radio-vector r0(t´), tal como se indica en la Figura (7.6.1).

Figura (7.6.1)

Examínese con detalle el potencial escalar en el punto P. Para ello, se re-escribe la

expresión para el potencial retardado  de la siguiente forma

1  r´, t dV ´ 1   r´,t´ t´t dV ´ 

(r , t ) 
4 0 
v
R

4 0  
v
R
dt´ .

(7.6.1)

donde [t]=t-R/c, R  r  r´ . En la segunda igualdad se ha utilizado la identidad

  t´t dt´  1.
Para la carga puntual que se mueve a lo largo de una trayectoria r0(t´) la densidad
volumétrica de carga es

(r´,t´)=q(r´- r0(t´)).

Reemplazando esta expresión en (7.6.1) y recordando que t   t 

1
r  r´ , se obtiene
c

356
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

q   r´r0 t´  1  
( r , t ) 
4 0  dt´   t´t  r  r´ dV ´ .
 v r  r´  c  

La integral con respecto al volumen se calcula con base en la fórmula

 1 
  t´t  r  r´ 
f r´ r´r0 t´dV ´ f r0 t´ , donde f r´    . Se tiene así
c

v
r  r´

   
  t´t  r  r0 t´ 
1 1
  t´t  r  r´ 
 c  r´r t´dV ´  c .
v
r  r´
0
r  r0 t´

Por lo tanto, el potencial escalar es

 
 r  r t´   t´t  c r  r t´ dt´ .
q 1 1
( r , t )  (7.6.2)
4 0
0
0

Con el fin de evaluar la integral resultante, se define ahora una nueva variable

r  r0 t´ .
1
t´´ t´t  (7.6.3)
c
Si se diferencia esta expresión (recordando que dt=0 ya que t es el tiempo fijo de
observación), se obtiene

1 d 
dt´´ dt´  r  r0 t´ dt´ (7.6.4)
c  dt´ 

La derivada temporal que aparece en esta expresión es

   dx
r  r0 t´   r  r0 t´  0i .
d
dt´  x 0i  dt´

dx 0i
Téngase en cuenta que es la i-ésima componente de la velocidad v de la
dt´
partícula en el instante t´. Además, de la definición r  r0 t´ = xi  x0i t´xi  x0i t´ se
deduce que

 r  r0 t´i
r  r0 t´   .
x0i r  r0 t´

357
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

Por lo tanto

v  r  r0 t´
r  r0 t´  
d
.
dt´ r  r0 t´

Introduciendo esta expresión en (7.6.4), se obtiene

 1 v  r  r0 t´ 
dt´´ dt´1  
 c r  r t´ 
 0 

y por lo tanto, la ecuación (7.6.2) para el potencial retardado de la carga puntual

considerada toma la forma

 
 t´´dt´´ ,  t´ t´´t  r  r0 t´  .
q 1 1
( r , t ) 
4 0   1 v  r  r0 t´   c 
r  r0 t´ 1  
 c r  r0 t´ 


La integral que aparece en esta expresión se calcula con base en la fórmula

 f t´´ t´´dt´´ f 0 , f t´´ 

1
donde siendo
 1 v  r  r0 t´ 
r  r0 t´ 1  
 c r  r0 t´ 

t´ t´´t  r  r0 t´ . Así,
1
c

q 1
( r , t )  .
4 0  1 v  r  r0 t´ 
r  r0 t´ 1  
 c r  r0 t´ 
 t´´ 0

Evaluar esta expresión en t´´=0 es equivalente a tomar t´= t- r  r0 t´ /c. Se obtiene
así

q 1
( r , t )  . (7.6.5)
4 0
r  r0 t   v  r  r0 t 
1
c

Para el caso del potencial vectorial, se obtiene de manera análoga

358
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

0q v
A(r , t )  . (7.6.6)
4
r  r0 t   v  r  r0 t 
1
c

Si se introduce la notación

(r,[t])= r  r0 t   v  r  r0 t  ,

1
(7.6.7)
c

las anteriores relaciones se pueden escribir de una forma más compacta

q 1 q 0 v
( r , t )  , A(r , t )  . (7.6.8)
4 0  4 

Los potenciales (7.6.5-6) asociados a una carga puntual que se mueve siguiendo una
trayectoria r0(t) reciben el nombre de potenciales de Lienard-Wiechert.
Interpretación de los potenciales de Lienard-Wiechert. Una interpretación
intuitiva de estas relaciones se puede dar de la siguiente forma: se considera una carga dq
distribuida uniformemente con densidad  en un volumen dV´. Sea R el radio vector de dV´
al punto de observación P. Se rodea ahora el punto P con una esfera de radio R´
suficientemente grande para incluir la carga dq. Permítase ahora que en el instante t-R´/c la
esfera empiece a colapsar con una velocidad c. Durante este colapso se barre la distribución
de carga considerada. Si la carga dq es estacionaria, la cantidad de carga barrida por la
esfera que colapsa al contraerse dR en un tiempo dt´es []dadR (donde da es el área de la
sección transversal). Pero si la carga se mueve con velocidad v, la cantidad de carga barrida
por la superficie que colapsa es dq=[]dadR-[](vR/R)dadt´. Teniendo en cuenta que
dR=cdt´ y dadR=dV´, entonces

  dV ´dq
.
R 1
R vR
c
Introduciendo esta expresión en la fórmula para el potencial retardado (7.1.6), se obtiene el
resultado (7.6.5).

§ 7.7. CAMPOS DE LIENARD-WIECHERT.

Determínense ahora los campos E y B de una carga puntual que se mueve de manera
arbitraria. Para ello se utilizan las relaciones

A
E  grad   , B = rot A. (7.7.1)
t

359
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

Reemplazando aquí las expresiones para los potenciales de Lienard-Wiechert en la

forma (7.6.8), se obtiene

q   1  1   v 
E   grad    2    
4 0     c t    
q  1 1 1  1 1 v 
=  2 grad   2 2 v  , (7.7.2)
4 0   c  t c 2  t 

1 q v 1 q 1 1 
B rot   2  rotv  grad    v  
c 2 4 0    c 4 0   
1 q 1 1 
 2  rotv  2 grad   v  . (7.7.3)
c 4 0    

En estas expresiones hay que diferenciar con respecto a las coordenadas x, y y z del
punto de observación en el instante t. Sin embargo, en la fórmula (7.6.7), que define la
magnitud , entra el tiempo de retardo t´=t-R/c. Por esta razón es preciso calcular
inicialmente las derivadas con respecto a t´. Para ello se tienen en cuenta los siguientes
resultados:
 De la definición R=r-r0(t´) se deduce que R/t´=-v(t´).
 Al derivar con respecto a t´ la identidad RR =R2 se obtiene R(R/t´)=R(R/t´),
R
es decir, (R/t´)=-v(t´) e R , donde eR  es un vector unitario en la dirección del
R
vector R.
 De la definición de tiempo de retardo t´=t-R/c se sigue que R(t´)=c(t-t´) y por lo
tanto

R R t´  t  t´  t´ 

  c  1  c1   .
t t´ t  t´  t  t 
t´  t´ 
 Como (R/t´)=-v(t´) e R , entonces se tiene v t´  eR   c1   , de donde se
t  t 
obtiene inmediatamente

t´ 1 R
  . (7.7.4)
t 1  1 v t´  e  
R
c

 De manera análoga, tomando el gradiente de t´=t-R/c, se halla que

 R 1 R 
grad t´ grad  t     grad R   grad R t 
1
grad t´
 c c c t´ 

360
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

1 R 
Teniendo en cuenta que (gradR)t´=eR, se encuentra que grad t´   eR  grad t´ .
c t´ 
Agrupando los términos que contienen la magnitud gradt´, se obtiene la relación
 1 R  1
1  grad t´  eR . Teniendo en cuenta adicionalmente que (R/t´)=-v(t´)eR,
 c t´  c
 1 
la última expresión toma la forma 1  v t´  eR grad t´  eR y por lo tanto se
1
 c  c
llega a la siguiente relación para el gradiente del tiempo de retardo:

1 eR 1R
grad t´   . (7.7.5)
c 1  v t´  e 
1 c
R
c

Con la ayuda de (7.7.4) se puede de inmediato calcular las derivadas temporales que
aparecen en la expresión (7.7.2). Si se tiene en cuenta que
R  r t´
 r  r0 t´   0  v , siendo v la velocidad de la partícula y v/t´=a , donde
t´ t´ t´
a es su aceleración, entonces se tiene:

v v t´ R R
  a   a , (7.7.6)
t t´ t  

  t´    1  R  R 1 v 1 R  R
    R  R  v     R  v 
t t´ t  t´  c    t´ c t´ c t´  
 1 1 R R  1  1 
  v  e R  R  a  v  v     v   e R  v   R  a 
 c c    c  c 

Introduciendo la magnitud
~e v .
n (7.7.7)
R
c
se obtiene
 R ~ 1 
  n  v  R  a , (7.7.8)
t  c 

~ mediante la
Cabe destacar que la magnitud  está relacionada con la magnitud n
relación
~  R,
n (7.7.7a)

la cual se evidencia de la siguiente cadena de relaciones:

361
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

R R 1 R 1   1 
r , t   R  v  R 
1 ~ R.
 v  R    v   R   eR  v   R  n
c R c R c   c 
Resta determinar grad. Dado que =(R,t´) entonces

grad   grad  t´  grad t´ . Pero (grad)t´=gradR(R-vR/c)=n-v/c. Por lo tanto,
t´

v  ~   n
~  v  1 R  a grad t´ ,
grad   eR   grad t´ n
c t´  c 

o bien, teniendo en cuenta (7.7.5), se tiene:

~   n
grad   n ~  v  1 R  a  R . (7.7.9)
 c  c

Reemplazando ahora (7.7.8) y (7.7.9) en (7.7.2), se llega a la relación

q  1 1 1  1 1 v 
E  2 grad   2 2 v  =
4 0   c  t c 2  t 
q  1 ~ R ~
=  2 n  n  v   R2 R  a   2R 3 v n~  v  1 R  a   2R 2 a  .
4 0    c c   c   c  c  

Reagrupando términos y separando los términos que contienen la aceleración de la

partícula, se obtiene:

q  1 ~ 
E=  3  n   n n  v   2 3 n
R~ ~ R ~
R  a   a  .
4 0    c  c  

~  R , entonces
Si se tiene en cuenta ahora que   n

~
n  v    n~n~  R  R n~n~  v  
R~ ~
 n  n
 c   c 
 ~ ~   ~  ~ 
 n n  R   n  v     n R n  eR   n  v   
R ~ 1 ~
  c    c 
 ~ ~  v    ~  v   ~  v 2 
 n R n   eR      n
v 
R  eR     eR     n R1  2 
   c     c  c    c 
y

n~R  a   a   R  n~  a .
Se tiene así:

362
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

q R  v2 ~ ~  a  .
n  2 3 R  n
R
E=  1  2 
4 0  3  c  c  

~,   R  1v  R y
Retornando a las definiciones de las magnitudes  y n
c
~  e  v , se obtiene finalmente
n R
c

 v2 
 1 2 
q   v   v  
E r , t   c 1
  R  R  R   R  R   a   . (7.7.10)
4 0   v  R
3
 c   v  R 
3
 c  
R  c  c2  R   
   c  

Todas las magnitudes en (7.7.10) se toman en el tiempo retardado t´=t-R/c.

Se procede ahora de manera idéntica para el cálculo de la inducción magnética B. Si

en (7.7.3) se tiene en cuenta que

 t´ v 1
rotv  v    grad t´a   Ra ,
r r t´ c

se obtiene

1 q  1 1 
B=   2 R  a  2 grad   v  ,
c 4 0  c
2
 

o bien,

1 q  1 1 ~ R ~ R  
B=   2 R  a  2  n  n  v  2 R  a   v  .
c 2 4 0  c   c c   

Reagrupando términos, se tiene

1 q  1 ~  
 3  n  v  n  v R  v   2 3 eR  a R  v  ceR  a  .
1 ~ R
B= 
c 4 0   
2
c  c 

~  e  v se ve que
~  R y la definición n
Con base en la relación   n R
c

363
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

 ~   v   v2 
2
 n  v  n  v R  v   1  2
1 ~ ~
 R  v  c1  2  R  n
 c   c   c 

eR  a R  v  ceR  a   ceR  a R  n~  eR  a  cR  eR  a n~   a 

 R 
 cR  eR  a n  n  eR a   ceR  R  n  a .
~ ~ ~

Se obtiene por lo tanto


 q  R  v 2  ~ ~  a 
1
B= eR    
3 
1  
2 
n 
R
R  n  .
 4 0    c 
 c 
2 3
c 

Nótese que la expresión entre llaves es el campo eléctrico (7.7.10). Se llega así al
siguiente resultado:
1
B  eR  E . (7.7.11)
c

Se ve que para una carga puntual que se mueve arbitrariamente la inducción

magnética siempre es perpendicular al campo eléctrico.
El campo eléctrico en (7.7.10) consta de dos partes de naturaleza diferente. El
primer miembro depende sólo de la velocidad de la partícula pero no de su aceleración. A
grandes distancias este término es inversamente proporcional a R2. El segundo término
depende de la aceleración de la partícula a y a grandes R cambia como 1/R. Es precisamente
este término el que está asociado a la radiación de ondas electromagnéticas.

§ 7.8. CASOS PARTICULARES DE LOS CAMPOS DE LIENARD-WIECHERT.

Se propone ahora examinar con detalle los campos de Lienard-Wiechert

 2 
 1 v 
q   v   v  
E r , t   c 2 1
3 
R  R  R   R  R   a   ,
4 0   v  R   c   v  R 
3
 c  
R   c R 
2
 
 
 c   c  

1
B  eR  E
c

en algunos casos particulares de interés.

364
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

Campos de una carga con movimiento uniforme. Supóngase inicialmente que la

partícula se mueve a lo largo del eje x con velocidad uniforme v tal como lo indica la Figura
(7.8.1). En este caso la aceleración de la partícula es cero y los campos de Lienard-
Wiechert en el punto de observación P toman la forma

 2 
 1 v 
q    
E r , t   c 2 v
 3 
R  R , (7.8.1)
4 0 
  R  v  R  
c 

 c  

1
B  eR  E . (7.8.2)
c

Se recuerda que en estas expresiones el radio vector R está definido en el instante

t´=t-R/c. Durante el intervalo de tiempo t-t´ que la información ha gastado en llegar del
punto 1 a P, la partícula se ha desplazado a lo largo de x de la posición 1 a la posición 2
(denominada posición actual). El respectivo radio vector r(t) trazado desde la posición
actual al punto P se denomina radio vector actual. El objetivo es expresar los campos
(7.8.1-2) en términos de la posición actual.
De acuerdo con la figura (7.8.1), tiene lugar la relación vectorial R=r+(t-t´)v, donde
el tiempo de retardo t-t´=R/c. Por lo tanto el radio vector actual de la partícula está
relacionado con R de la siguiente forma

v
r  R R. (7.8.3)
c

Este es precisamente el vector que entra en (7.8.1). Para expresar el denominador en

términos de r, examínese la expresión

v R
2

  R  2 v  R  2 v  R  .
R 1
R 
2 2

 c  c c

Pero, según (7.8.3)

R 1
r 2  R2  2 v  R  2 v2R2 .
c c

365
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

Figura (7.8.1)

Restando estas dos expresiones y teniendo en cuenta que vR=Rvcos, donde  es

el ángulo entre la dirección de movimiento de la partícula y el radio vector actual r, se
obtiene

v R
2
 r 2v 2
R    r  2 sen  .
2 2

 c  c

Téngase en cuenta que las proyecciones transversales de los vectores R y r son

iguales, es decir, Ry=ry. Por otro lado, Ry es igual a la magnitud del producto vectorial Rex;
así mismo, ry es igual a la magnitud del producto vectorial rex. Recordando que entre los
productos escalar y vectorial de dos vectores arbitrarios a y b existe la relación
(ab)2+(ab)2=a2b2, para el cuadrado de la magnitud del producto vectorial Rex se puede
escribir (Rex)2+(Rex)2=R2, o bien, Ry2=R2-(Rex)2; de manera análoga, ry2=r2-(rex)2. La
igualdad Ry2=ry2 toma entonces la forma R2-(Rex)2=r2-(rex)2; dado que v=vex, esta última
relación implica que R2v2-(Rv)2= r2v2-(rv)2, de donde se obtiene entonces
1/ 2
 v  R   v2 2 
 R    r 1  2 sen   . (7.8.4)
 c   c 

366
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

Reemplazando (7.8.3-4) en (7.8.1) se llega a la siguiente expresión:

v2
1 2
E r , t  
q c r
. (7.8.5)
4 0  v 2 
3/ 2
r3
1  2 sen 2 
 c 
En el sistema inercial de referencia con respecto al cual la partícula está en reposo
(v=0) esta expresión se reduce al campo coulombiano, tal como era de esperar. En los
demás sistemas inerciales el campo eléctrico presenta una dependencia con el ángulo 
cuyo origen está relacionado con la existencia de una dirección "privilegiada" determinada
por la velocidad de la partícula. Nótese que a medida que v/c crece desde 0 a 1, la magnitud
del campo eléctrico aumenta en la dirección transversal y decrece a lo largo de la dirección
del movimiento. Este efecto se puede visualizar en términos de diagramas polares para la
magnitud del campo eléctrico, tal como se ilustra en la Figura (7.8.2).

Figura (7.8.2)

367
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

Resta ahora expresar el campo magnético B en términos del radio-vector actual r;

reemplazando (7.8.5) en (7.8.2), se tiene:

  2 
  1 v 
1 1  q  c 2  v  
B  eR  E  eR    3 
R  R   
c c  4 0   R  v  R   c  

 
 c  
v2 v2
1 2 1
1 q c   v  1 q c2  eR  v 
e   R  R  
3  R 3 
 R .
c 4 0  vR   c  c 4 0  vR  c 
R   R  
 c   c 

v v r v
De la relación r  R  R se sigue que R  r  R y por lo tanto eR   . De
c c R c
r v
aquí se deduce que eR  v  . Recordando adicionalmente el resultado (7.8.4), se
R
obtiene finalmente la siguiente expresión para la inducción magnética de una carga que se
mueve uniformemente:

v2
1
0q vr
Br , t   c2 . (7.8.6)
4  v2 
3/ 2
r3
1  2 sen 2 
 c 

Habiendo examinado el comportamiento de los campos de Lienard Wiechert para

una partícula que se mueve con velocidad uniforme, se discute ahora el caso cuando la
partícula se mueve con aceleración. Recuérdese que el campo eléctrico en este caso consta
de dos partes de naturaleza diferente: el primer miembro depende sólo de la velocidad de la
partícula pero no de su aceleración y a grandes distancias este término es inversamente
proporcional al cuadrado de la distancia entre la partícula y el punto de observación, como
lo evidencia el resultado (7.8.5). El segundo término depende de la aceleración de la
partícula a y a grandes distancias cambia como 1/R. Comparando el aporte de estos
términos a grandes distancias, se ve que el término dominante es el que depende de la
aceleración de la partícula y es por esta razón que en la discusión que sigue sólo se
considera dicho término, el cual está asociado a los procesos de radiación.
Radiación de una partícula relativista con velocidad y aceleración colineales.
Así, se examina el término de radiación en los campos de Lienard-Wiechert

368
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

 v  
E rad r , t  
q 1
R   R  R   a 
4c  0 
2
vR
3
 c  
R  
 c 

y se considera el caso particular cuando la velocidad v y la aceleración a son paralelas entre

sí (es decir, va=0). Entonces

E rad r , t  
q
4c  0 
2
1
R  R  a  
q
4c  0 
1

R  a R  R 2 a . 
vR
3 2
vR
3

R   R  
 c   c 

1
Teniendo en cuenta que Brad  Erad  e R , la parte del vector de Poynting que contribuye a
c
la radiación es

0 q 2
Srad 
0
1
Erad  Brad 
1 2
0 c
EradeR =
16 c 
1
R  a R  R 2a 2 eR . 
vR
2 6

R  
 c 
Si  es el ángulo entre la aceleración a y el radio-vector retardado R, entonces

R  a R  R a 2 2
 
 R 2 R 2 a 2  R  a   R 4 a 2 sen 2
2

y se obtiene por lo tanto para el vector de Poynting la expresión

0 q 2 a 2 sen 2
S rad  eR . (7.8.7)
16 2 R 2c  v 
6

1  cos 
 c 

Téngase en cuenta ahora que la cantidad de energía dE radiada en unidad de ángulo

sólido en el intervalo de tiempo dt es dE=(SradeR)R2dt, donde (SradeR) es la componente
radial del vector de Poynting calculada en el instante t y corresponde a la radiación emitida
por la partícula en el instante t´. Por lo tanto la potencia radiada en unidad de ángulo sólido
y que atraviesa una superficie de radio R en el instante t es igual a la pérdida energética en
unidad de tiempo en el instante t´, es decir,

 Srad  eR R 2
dP dW dt
 .
d dt´ dt´

Dado que t´=t-R(t´)/c, entonces

369
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

d  1  1 dRt´ 1 1  dR 
 t´ Rt´  1 
dt 1 d
 1 R  R 1   R
dt´ dt´  c  c dt´ c dt´ c R  dt´ 

dR
Recordando que  v , se tiene
dt´

 1  v  eR   1  cos 
dt 1 v
dt´ c c

y por lo tanto, se obtiene para la potencia radiada en unidad de ángulo sólido la expresión

dP  0 q 2 a 2 sen 2
 . (7.8.8)
d 16 2 c  v 
5

1  cos 
 c 

En el caso no-relativista (v/c<<1)

la ecuación (7.8.8) se reduce a la expresión que describe la distribución angular de la
radiación de una partícula acelerada en la aproximación dipolar (dP/dsen2). Al crecer
la velocidad de la partícula, la intensidad de la radiación aumenta considerablemente en la
dirección del movimiento de la partícula. Esto se ilustra en la Figura (7.8.3). Nótese que a
pesar de que la distribución angular de la radiación se inclina en la dirección del
movimiento, no se produce radiación exactamente en =0. Para partículas aceleradas en un
acelerador lineal, esto quiere decir que hay una estrecha zona de no radiación a lo largo del
camino recorrido por el haz.
Por otro lado, si el haz de partículas es detenido por un bloque de material, tiene
lugar una dispersión ya que la dirección original del movimiento se ha alterado. Esta
radiación emitida durante el proceso de frenado se denomina radiación de frenado
(bremsstrahlung)
La potencia total radiada es
3
 q 2 a 2   v  
2

P 0 1    , (7.8.9)
6c   c  


resultado que se propone demostrar en calidad de ejercicio.

370
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

Figura (7.8.3)

Radiación de una carga en movimiento circular. Considérese una carga puntual

q que se mueve con movimiento circular en el plano yz y centro en el punto (0,, 0) (Figura
(7.8.4)). Se propone determinar los campos en el instante t´ en el punto P, cuando la
partícula pasa en el instante t´ por el punto (0,0,0).
El vector de Poynting asociado a la radiación es

2
1 1 2 0 q 2 1  v  
Srad  Erad  Brad  EradeR  R   R  R   a  eR .
0 0 c 16 c 
2
vR
6
 c  
R  
 c 

De acuerdo con la geometría de la Figura (7.8.4), se tiene

R=R(sencosex+senseney+cosez), v=vez, a=aey y por lo tanto

 v 
2

2 4  v 
2
2  v2 
R   R  R   a  a R 1  cos   sen cos  1  2
2
  ,
 c   c   c 

vR=vRcos .

371
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

Figura (7.8.4)

Así, el respectivo vector de Poynting toma la forma

2  v2 
2
 v 
1  cos   sen cos  1  2 

2

0 q 2  c   c e
Srad 
16 R c
2 2 6 R
 v 
 1  cos  
 c 

y por lo tanto para la potencia radiada en unidad de ángulo sólido dP(t´)/d=SradeR

R2dt´/dt, se obtiene

2  
2
 v  v2
 1  cos    sen 2
cos  1 
 c2 
dP  0 q 2 a 2  c   .
 (7.8.10)
d 16 c 2
 v 
5

1  cos 
 c 

Este patrón de radiación es difícil de visualizar debido a su dependencia de los

ángulos  y . En la Figura (7.8.5) se ilustra su comportamiento en el plano de la órbita
=0. Nótese que la radiación es emitida preferencialmente en la dirección de la velocidad
de la partícula v. Cabe esperar que a medida que la velocidad v se acerca a la velocidad de
la luz, la radiación es emitida en un cono frontal de abertura pequeña.

372
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

Figura (7.8.5)

La potencia total radiada es

2
0q 2a 2   v 2 
P 1     . (7.8.11)
6c  c 
 

Si se compara esta expresión con (7.8.9), correspondiente a la radiación total

emitida por una partícula cuya velocidad es paralela a la aceleración, se ve que a bajas
velocidades ambas expresiones coinciden y se reducen a la fórmula que describe la potencia
de la radiación dipolar. Sin embargo, a velocidades cercanas a la velocidad de la luz estas
potencias difieren considerablemente. Este resultado tiene implicaciones decisivas en el uso
de aceleradores de partículas: en un oscilador lineal la única aceleración está relacionada
con el aumento de la velocidad y una vez ésta se acerca a la velocidad de la luz (vc) la
radiación emitida decrece. En los aceleradores circulares (betatrones, ciclotrones,
sincrotrones) hay una aceleración centrípeta continua a=v2/c2/ ( es el radio de la
órbita) y la respectiva radiación se vuelve dominante cuando vc. Por otro lado, si la
partícula se encuentra sometida a la acción de un campo magnético B0, el radio de la órbita
se puede determinar fácilmente de la ecuación del movimiento

1 v2
m  qvB0 ,
v2 
1 2
c

de donde se obtiene que la aceleración centrípeta es a(qcB0/m)(1-v2/c2)1/2.

373
CAPÍTULO 7. RADIACIÓN DE ONDAS ELECTROMAGNÉTICAS

En este caso la potencia radiada es

1
 q 4c 2 B02   v  
2

P 0 1    .
6cm 2   c  


Se ve que los aceleradores circulares son más eficientes que los lineales, además de
que son más compactos en tamaño. Nótese que en la última expresión la potencia radiada es
inversamente proporcional al cuadrado de la masa de la partícula. Esto indica que las
pérdidas energéticas por radiación son mucho mayores para los electrones en comparación
con las de los protones. Es por esta razón que los aceleradores circulares se eligen para
acelerar protones, mientras que se prefieren los aceleradores lineales para los electrones.

Referencias capítulo 7

Los potenciales retardados y los potenciales y campos de una carga puntual que se mueve arbitrariamente se
discuten, desde diversos puntos de vista, en
Greiner (1998), capítulo 20
Landau y Lifshitz (1987), secciones 62,63
Jackson (1999), secciones 8.1, 8.2, 14.1-14.5
Panofsky and Phillips (1962), capítulos 19-22
Heald and Marion (1995), secciones 8.3-8.8
Thidé (2001), capítulo7
Levich, Vdovin y Miamlin (1974), secciones 23, 25
Brédov, Rumiántsev y Toptigin (1985), secciones 32, 34
Sepúlveda (2009), capítulo15

La radiación de sistemas simples (dipolos y antenas) se discute, por ejemplo, en

Corson y Lorrain (1972), capítulo14
Greiner (1998), capítulo 21
Heald and Marion (1995), secciones 9.1-9.4, 9.6-9.8
Landau y Lifshitz (1987), secciones 66, 67, 71,72,74
Jackson (1999), 9.1-9.5
Levich, Vdovin y Miamlin (1974), secciones 26-28
Becker R. (1981), secciones 64-66, 68, 69
Plonsey and Collin (1961), cap 11
Rossiter (1998), capítulo 8
Thidé (2001), capítulo 8 (rad)
Cooper N. (2001), capítulo5

Problemas capítulo 7
1. Considere un sistema de n partículas con masas m1, m2,... mn y cargas q1, q2,... qn. Demuestre que
cuando la relación carga-masa de estas partículas es la misma, el sistema no puede radiar dentro de la
aproximación dipolar si el centro de masas se mueve con velocidad constante (compárese con el
caso de la radiación dipolar de un sistema de dos partículas).
2. Un cascarón esférico está cargado uniformemente y realiza oscilaciones radiales. Demuestre que este
sistema no puede radiar. (el problema se reduce a demostrar que la inducción magnética de tal
sistema es igual a cero).

374
CURSO DE TEORÍA ELECTROMAGNÉTICA PARA FÍSICOS

3. Dos cargas, +q y –q están separadas una distancia l y rotan con velocidad angular constante 
alrededor de un eje que pasa perpendicularmente por el punto medio de la línea que une las cargas.
Determine las direcciones de los campos de radiación y la intensidad de la radiación (a) en el plano
de rotación, (b) a lo largo del eje de rotación.
4. Un anillo aislante circular de radio r descansa en el plano xy y tiene su centro en el origen de
coordenadas. El disco tiene una carga distribuida con una densidad lineal =0sen, donde  es el
ángulo azimutal. El anillo se hace rotar alrededor de su eje z con una velocidad angular constante .
Dentro de la aproximación dipolar, determine: (a) El momento dipolar de este sistema y el potencial
vectorial. (b) Los campos en la zona de radiación. (c) La potencia radiada.
5. Considere el arreglo de antenas de media onda de la Figura (7.4.3). Suponga que la distancia entre las
antenas es igual a media longitud de onda y las antenas están oscilando en fase. Determine la
potencia radiada en unidad de ángulo sólido a lo largo del eje x (=/2, =0) y a lo largo del eje y
(=/2, =/2). Compare los resultados con el obtenido para una sola antena de media onda.
6. Suponga ahora que en el sistema del problema anterior las antenas están oscilando en oposición de
fase. Determine la potencia radiada en unidad de ángulo sólido a lo largo del eje x (=/2, =0) y a
lo largo del eje y (=/2, =/2). Compare los resultados con el obtenido para una antena
cuadrupolar.
7. Considere el cuadrupolo magnético
constituido por dos espiras paralelas de
radio a separadas una distancia a, tal
como se indica en la Figura. Las
corrientes en las espiras oscilan en
oposición de fase. Determinar los campos
de radiación, la potencia radiada en
unidad de ángulo sólido, la potencia total
radiada y la resistencia de radiación.

8. Mediante la integración sobre una superficie esférica de la expresión (7.8.5) para el campo eléctrico
de una carga que se mueve uniformemente, compruebe que tiene lugar la ley de Gauss Eds=q/0.
9. Con base en la relación (7.8.8) demuestre que el ángulo para el cual es máxima la radiación
producida por una partícula cuyas aceleración y velocidad son colineales está dado por
1  15 2  1 , donde =v/c. Demuestre que cuando 1, 1   . Suponga que
cos max   max 
3 2
una partícula bajo estas condiciones tiene una energía cinética de 100 keV. Dibuje un diagrama polar
del patrón de radiación. ¿Cuál es el ángulo de máxima radiación?
10. Demostrar que la potencia radiada por una partícula cuyas aceleración y velocidad son colineales es
2 3
 q 2a 2   v 

 (fórmula (7.8.9)).
P 0 1  
6c   c  

11. Demostrar que la potencia radiada por una partícula cuyas aceleración y velocidad son
2
  v 2 
perpendiculares es P   0 q a
2 2
1     (fórmula (7.8.11)).
6c  c 
 
12. Demuestre que desde el punto de vista de la electrodinámica clásica un átomo de hidrógeno es
inestable, ya que un electrón orbitando alrededor del protón (con un radio de órbita de 0.5310-8 cm)
radiaría energía con una potencia de 31011 eV/s.

375