Erl.

ferianz{]

Revista Mexicana de Fsica 34 No. (1988) 6.5-669

El desarrollo multipolar

por fuera y por dentro

E. Ley-Koo y Araceli Gngora-T. Instituto de Fsica, Universidad Nacional Autnoma de Mxico, apartado postal 20-364, 01000 Mxico, D. F.
(recibido el 17 de febrero de 1988; aceptado el lo. de agosto de 1988)

Resumen. Se destacan la validez y la unidad del desarrollo multipolar no slo en el exterior sino tambin en el interior de una regin donde se encuentren distribuidas las fuentes de campo electrosttico o magnetosttico. Se ilustran la unidad del desarrollo y las caractersticas de cada uno de sus trminos al reconocer qu fuentes de una multipola.ridad definida distribuidas sobre una superficie esfrica dan lugar a campos con la multipolaridad correspondiente tanto fuera como dentro de la esfera. Se dan ejemplos de trampas de partculas cargadas y de partculas neutras con momento magntico, en las cuales el carcter multipolar de los campos en el interior es crucial para la dinmica de esas partculas. PACS: 41.1O.Dq; 41.70.+1

Introduccin El estudio del desarrollo multipolar en electrosttica y magnetosttica en los niveles introductorio [1), intermedio [2] y avanzado [3] comnmente se restrige al exterior de una regin donde se encuentran confinadas las fuentes de los campos. Algunos autores hacen mencin al desarrollo multipolar en el interior de la regin de las fuentes (4), pero no llevan a cabo un anlisis detallado de la situacin. En este trabajo se presenta un estudio del desarrollo multipolar con nfasis en su validez general y en la unidad de su aplicacin tanto fuera como dentro de la regin de las fuentes en los casos electrosttico y magnetosttico. La situacin electrosttica se analiza en la seccin 2, partiendo del ejemplo familiar de una carga elctrica puntual y su campo coulombiano. Como punto de referencia, se presenta el desarrollo multipolar del potencial coulombiano para fuera, y se da la interpretacin comn de sus trminos y de la generacin de los mismos para distancias grandes, asocindolos a distribuciones de cargas puntuales que se introducen al tratar de centrar en el origen de coordenadas las fuentes multipolares sucesivas. En contraste, al considerar el desarrollo correspondiente del mismo potencial para dentro, se reconoce la interpretacin alternativa, vlida tanto adentro como afuera, en trminos de distribuciones de carga con multipolaridad definida sobre una superficie esfrica centrada en el origen y que contiene al punto donde se encuentra la carga original; la necesidad de los trminos sucesivos est asociada con la representacin de la carga puntual corno una superposicin de cargas superficiales. A continuacin, el principio de superposicin permite escribir el desarrollo multipo-

646

E. Ley.l\oo yA. Gngom.T.

lar para cualquier distribucin de carga elctrica en general, destacndose que cada componente de la fuente con una multipolaridad definida produce una componente del campo con la misma multipolaridad, tanto por fuera como por dentro. La siluaci(n magnetosttica se estudia en la seccin 3, siguiendo dos caminos alternativos; el primero es general y se ba..qa en el uso del potencial vectorial, mientras que el segundo usa pOl<-'ncialcs escalares y slo es vlido en regiones (Ionde no haya fuent.es. El desarrollo multipolar del potencial vectorial magnetosttico es anlogo al del potencial electrosttico ('011 la diferencia natural de las fuentes, que son el vector de densidad ele corriente elctrica ('stacionaria y la densidad de carga e1ct.rica, respect.ivamente. Adems, micntras el campo de densidad elctrica se obtiene calculando el gradiente del potrnfal clectrstatico, el campo de induccin magntica se ohtiene como el rotacional <1('1 potencial vectorial magnetosttico. En aquellas regiones del espacio donde no hay corrientes elctricas el campo es irrotacional, y por lo tanto, se puede obtener (01110 el gradiente de un potcncial escalar magnelosttiC'O. Seguirnos a Bronzan [5] para construir este tipo de potencial cn las regiones externa e interna a la regin donde estn localizadas las fucntcs, escribindolos bajo la forma de d<'SarrollosTnllltipolares e identificarHlo las densidades de carga de magn!'tizacin que juegan el papel de las fuentes. El uso de ('sl.os potcnciai('s ('Scalares en la situacin magnetosttica tiene la ventaja de poder aprovechar la analoga e\ectrosttic;" f('colloci('uclo desde luego las diferencia..". Se muest.ra la equivalellcia de los dos camillos p"ra ohtpner el campo de induccin magntica para fuentC's de multipolaridacl definida restringidas sobre la superfi(ic de lIna esfera. En la seccin .1se discuten algunos dispositivos que utilizan campos electrostti. cos o magnetostticos en su interior para controlar el movimiento de partculas con carga elctrica, poniendo nfa..qis en ('1car('tcr multipolar de los campos involucrados. Tambi'n se hace una presentacin (it.1inter('S actual de desarrollar trampas para confinar partculas rlctricamentc ncutras pero COII motn('ntos dipolares magnticos, y la importancia de la multipolaridad de los ('ampos magnticos que se usan.

2. Electrosttica Una carga elctrica puntual q en rqlOso es ftlcnte del campo de Coulomb ca. raclerizado por un potencial invr'rsamrnt.e proporcional a la distancia de la carga y una intensidad de campo radial e inversamente proporcional al cuadrado dc la distancia,

,;(R) =

R'

(1 )

E(r) = -\1';= JI"

qR

(2 )

El dC$(Irmll" multipolrlr por fucm y por dcntm

647

donde R = r - r' es el vedar de posicin de los puntos del espacio desde la posicin r' de la carga. El potencial eledrosttico se fH1t'(leescribir en la forma explcita q
q

~(e) ;

-1

e - e

-'1;

Jr2 + r'2 _ 2n'(f.

i-')'

(3)

la cllal puede desarrollarse en series de potencias tic la coordenada radial

161

l'

.,

OC>

"

'() q qre.e </lr;;;;-+---+ 2

r r

Lqrp(""') ... ;;;; --lr.r,

1=0 r +1 '

(3e)

vlida para afuera, donde r >

r', y
'1

.: q ~ (e); r'

qrr. r + --,+ ... ; r'w

L -PI(e"e), qr ._,
00

1=0

r"+1

(3i)

vlida para adentro, donde r < r'. Las funciones P,(f. r') son los polinomios de Legendre y su argumento r . .' ;;;; cos I slo depende del ngulo I entre los dos vectores de posicin. El inverso de la distancia entre el punto del campo y el punto de la fuente juega el papel de funcin generadora de los polinomios de Legendre [6]. Las series de las Ecs. (3c, i) constituyen el desarrollo multipolar del potencial coulombiano en el exterior y en el interior de una esfera centrada en el origen y que pasa por la posicin de la carga, rcspectivamen1<l. A continuacin, vamos a dar una interpretacin posible de los trminos del desarrollo muItipolar. Empezamos COII el desarrollo para el exterior, Ec. (3e), cuyo primer trmino domina a grandes dist.ancias y tiene la forma de la Ec. (1) con la carga q situada en el origen de coordenadas. La Fig. 1 lIlucstra como sc podra tener esle trmino monopolar al inlroducir de a) a b) ulla carga q y otra -q en el origen, adems de la carga original q en r/; la carga original y su pOlencial se estn reemplazando por la carga q cenlrada eTlel origen con su campo coulombiano, y el par de cargas de signos opuestos que constiluyen un dipolo elctrico p = qr', con la carga -q en el origen y la carga original en r'. Notamos que el dipolo no esl centrado en el origen, pero se puede agregar un par de dipolos opuestos entre s y centrarlos en el origen, pil."iandode h) a eL teniendo entonces la carga centrada en el origen, el dipolo cenlrado en el origen y un par de dipolos opuestos que constituyen un cuadrupolo. El segundo ti'rmino del d(';~;arrollo de la Ec. (3e) est asociado al dipolo centrado y resulta de superponer los pot.enciales cottlombianos de cada .una de sus dos cargas, de modo que a (Iislancias grandes sus contribuciones mono polares se cancelan entre s y solo sobrevive el potencial proporcional a la proyeccin del dipolo a lo largo de la dir('('cio del vector ele posicin del plinto de observacin, p. r, e inversamente proporcional al cua<lrado de la distancia. Una vez ms, se nota que el

648

E. Ley-Koo y A. Gngom- T.

0 00 00 0<l>

-, r o

J
b)

r
j

0f 0<l> o 0 0~
00

O)

el

d)

FIGURA 1. Distribuciones de cargas equivalentes: a) carga puntual desplazada del origen, b) carga puntual centrada y dipolo desplazado, e) tTlonopolo y dipolo centrados y cuadrupolo desplazado y d), monopolo, dipolo y cuadrupolo centrad08 y octupolo desplazado.

cuadrupolo en e) no est centrado, y el costo de centrarlo, como se muestra en d), es que resulta acompaiado de dos cuadrupolos OPIJ(,'Stos que constituyen un octupolo que . su vez lampoco est centrado. En general, para centrar un multipolo de orden 2', es necesario que se acompae de un par de ffiultipolos del mismo orden y opuestos entre s, los cuales constituyen un rnultipolo del siguiente orden, es decir, 2/+1, Para el multipolo de orden 21 con magnilud qr/I, las conlribuciones de mullipolaridad ms baja son nulas tan lo para la distrioucin de carga como para el pOlencial, y su propia conlribucin al polencial es inversamenle proporcional a la potencial (l + 1) de la dislancia y esl modulada angularmenle por el correspondienle polinomio de Legendre de orden J, exhibiendo J conos nodales cuyo eje es la direccin del vector de posicin de la carga original. El nmero infinito de lrminos en el desarrollo mullipolar se puede enlender sobre la base de (tiC cada nuevo mullipolo que se introduce para cenlrar el anlerior esl a su vez descenlrado. Consideremos ahora el desarrollo mullipolar del polencial electrosltico por dentro, Ec. (3i). El primer lrmino q/r' no depende de la posicin en el inlerior de la esfera, por lo que puede inlerprelarse como debido a la dislribucin mono polar uniforme de la carga sobre la superficic de la esfera. El siguienle lrmino es lineal en la coordenada cartesiana en la direccin del vector de posicin de la carga original y corresponde a un campo uniforme E = - 'V4J = -qr' /r'2, cuyas lneas se originan en cargas positivas en el hemisferio donde se encuentra la carga original y terminan en las cargas negativas del olro hermisferio; las cargas y el pOlencial lienen la dis. tribucin cosenoidal asociada al polinomio de Legendre de orden uno con un nodo en plano ecualorial y magnilud mxima o mnima en los polos. En general, el trmino J-simo e~ un polencial de mullipolaridad 2 caracterizado por su variacin ' proporcional a la potencia 1 de la coordenada radial, modulada por el polinomio de Legendre de orden 1 y se puede interprelar como debido a una distribucin de

el

El df'lmrrollo mu/tipolar por fuero y por dentro

649

a)

b)

el

d)
e)

FIGURA 2. Distribuciones de cargas superficiales: a) monopolar uniforme b) dipolar, p(f.e), cuadrupolar, P,(f. e), y d) odupolar, P3(f. f').

carga sobre la superficie esfrica dada por el mismo polinomio de Legendre. La aparicin de los trminos sucesivos en el desarrollo multipolar se puede entender como resultado de la sustitucin de la carga puntual original por distribuciones de carga con multipolaridad definida sobre la superficie de la esfera, empezando con el trmino monopolar de distribucin uniforme; sta tiene que ser compensada con carga positiva adicional en el polo d")nde se debe reconstruir la carga original y con carga negativa en el resto de la esfera y de manera creciente hacia el otro polo. La Fig. 2 muestra esquer:nticamente las distribuciones de carga sobre la superficie esfrica para los multipolos ms bajos. Es bien sabido que el primer trmino del desarrollo multipolar para afuera, Ec. (3e), puede asociarse a la carga puntual en el origen o a la carga uniformemente distribuida sobre una esfera centrada en el origen. De hecho, los trminos sucesivos del desarrollo tambin se pueden interpretar en base a las distribuciones de carga de multipolaridad definida sobre la superficie esfrica mencionadas en el prrafo anterior, reconociendo la unidad de origen del desarrollo multipolar por fuera y por dentro, Ecs. (3e, i). Efectivamente, a partir del potencial podemos calcular el campo de inlensidad elctrica usando la Ec. (2),

E' =

-Y",

,..,~, = -

L'"
1=0

qr

" [.r-a

. a . 1 a] ar +'-rao +fI--r sen O a,,

(.

--r'+1

P,(r.r')

~r'+2

'"

qr"

[. .., r(l+l)P,(r.r)-

1 'aO+flscnOa,, P,(r.r)

a .

a)

.., ] ,

(40)

q [E; =-'V'' ; =- L'" -,r-+'-+fI--r' + I

1=0

a - a . 1 a], .., ar rao r senO a,, rP,(r.r)

650

E. Ley-Koo y A. Gngom-T.

= ="-L~I
1=0

qr'-t [_ _ _, -rll',(r-r)-

(_ {)

8-+~--M _O~

_ 1

{))

_ _,]

I',(r.r)

(4i)

Notamos quc para todos los puntos sohre la superficie esfrica r = r' las componentes tangenciales del campo elctrico. es decir en las direcciones; y ~, son continuas al pasar del interior al exterior. Pero las componentes normales, es decir en la direccin radial r, muestran una discontinuidad

r. (E' - Ei) =

L
1=0

~2(21 + 1)1',(T - r')

(5 )

que, de acuerdo con la ley de Gauss [1-4]' corresponde a una densidad de carga superficial

(6)

Aqu se identifican las densidades de carga superficial asociadas a cada multipolo

O"l(r

_)

---z--p r. r),
r
471"

q 2/+1

( __ ,

( 61)

confirmando la afirmacin de que siguen las variaciones de los polinomios de Legendre. En las Ecs. (4c, i) quedan por calcular las derivadas angulares de los polinomios de Legendre, y en la Ec. (6) queda por realizar la suma sobre l. Si la carga se toma sobre el eje z entonces " . k = cos O y en las Ecs. (4c, i) los trminos en ~ son nulos, reflejando la simetra bajo rotaciones alrededor del eje polar. Esa misma simetra se refleja en la independencia de 0"1(0) en Ec. (61) con respecto a 1>. Como se sabe que 1',(1) = 1 Y que y'(2/ + 1)/2I',(c050) constituyen una base completa de funciones orto normales al integrar sobre d(cosO) en todas direcciones de () = Oa 11", se reconoce que la suma en la Ec_ (6) es

u(O)

2,,"'2

= yf2i+l L -;--2-I',(I)Y
1=0

f2i+l

-;--2-I',(cosO)

211"r<28(cosO - 1)_

(6' )

Esta ecuacin muestra que la carga q est concentrada en () = O Y en consecuencia, el campo elctrico es discontinuo solamente ah. Si la carga no se toma sobre el eje z, la fsica no cambia, aunque s las expresiones

El dC8arrollo mu1tipolar por fuero y por dentro

651

en trminos de las coordenadas angulares. La conexin se logra a travs de una rotacin para llevar la carga de la posicin sobre el eje z a (r',O',4/). Entonces,

jo. jo/ =

cos O cos O/

+ sen O sen O/ cos(4J-

4J1),

(7)

y los polinomios de Legendre generan a su vez los armnicos esfricos a travs del

teorema de la adicin [61

(71)
Los armnicos esfricos son eigenfunciones del cuadrado del operador generador de rolaciones

1=

.
-Ir

XV

= -,

.(

1 -, sen O &1> + ~&1>

& . &)

(8)

y de su componente a lo largo del eje z

(8z) con eigenvalores 1(1+ 1) Y m, respectivamente. Entonces, la Ik doble suma sobre los ndices de muItipolaridad 1 y m, (6) involucra la

,,(0,1

= =

~2 r

LL
1=0 m=-l

00

Y,;"(OI,1>')lm(O,1
cos O')( 1>-1>'),

~2(

cos 0-

(6")

que nuevamente se puede reducir usando la complet.ez de la base ortonormal que forman los armnicos esfricos, y que muestra que la carga q est concentrada en (7'/. O/, 4J'). Se pueden identificar tambin sus componentes multipolares de densidad , -':lrga superficial

(6Im) que ahora siguen las variaciones angulares de los armnicos esfricos.

652

E. Ley-Koo y A. Gngom-T.

,--

,-- (
1

al
z
, ,,

/
y ,

"

i \, \
'. " 1 '/ I ./"-

_..L_

-----r.:-

.. "

.-

I Y

bl

,, ,
x

el
FIGURA 3 DistrihucionC's cuadrupolarf"S,

d)

e)

f)
)'1_1,

de cargas superficiales: dipolarl.'S, a) YII + )'1-1 ): b) 'VII e) )':!1 + }'2-1, d) )'21 - Y2-1, e) )'n + Y2-2 Y f), Y22 + Y2-2.

La conexin explcita entre los armnicos esfricos y los polinomios de Legendrc 16)

\"m(O,

,M ~

(_)m

(21

+ 1)(1-

411"

(1 +

m )! m ' scn m ).

J' ( O) im' d( COS O) m ,eos e

(9)

para m POSitivO y }-m(O,4 = (- t'};n(O, ,p), mlH'stra que la parte real y la parte imaginaria del armnico esfrico tienen 1 superficies Ilodalcs cada una: 1 - m son conos nodales asociados con las races de la m-sima derivada del polinomio de Legendre, y m son planos nodales asociados con las races de cos mq, o sen m,p, respeclivamente. La Fig. 3 ilustra las distribuciones de cargas puntuales y la carga superficial que acompaan a los multipolos de las Figs. 1 y 2, notndose los efectos de recentrado y rotacin. Al usar la expresin del teorema de la adicin para los polinomios de Legcndrc, Ec. (7/), en las Ecs. (3) del desarrollo multipolar, el ltimo toma la forma siguiente
1 I

00

Ir,-r'1 ~ ~ ~

'"

'"

1=0 m=-l

4~ r < \" (")\' 21+ 1 rl+1 1m r

>

1m

(')

r.

(31m)

que involucra la doble suma sobre los ndices de multipolaridad. Se reconoce tambin

El delmrrollo muitil'()/nr por fuero y por dentro

653

la correspondencia trmino por trmino entre las Ecs. (6) y (31m) y que implica que las fuentes de multipolaridad definida dan lugar a potenciales con la misma multipolaridad. El uso de la misma Ec. (7/) para los polinomios de Legendre en las Ecs. (4e, i) permite completar el clculo pendicnt.e de las derivadas angulares
. I

. 8 ( , 80

8), (' ., _ '\' + ~sen O 81> E, r . r ) - ~

m=-J

h "(' , (. 8 2/ + 11'm O .1 , 80

+ ~sen O 81>

8 )

( ) lm 0.1>

i\tese de las Ecs. (-1) y (8) que la parte angular del gradiente es equivalente a -ir x 1; esta relacin y otras anlogas son de utilidad para propsit.os de clculo y de comparacin, como se ilustra ms adelante. En el caso general de una distrihucin de carga eJ('ctrica con densidad p(r), el campo de intensidad elctrico obedece la ley de Gauss

v .E
y es conservativo o irrotacional

= '~p,

(10)

vx

E = O.

(11 )

Esta tltima propiedad permite cxprL'Sarel campo de intensidad elctrica en trminos del potencial electrosttico

E=-V1>,
y la Ec. (10) se convierte en la ecuacin de Poisson para el potencial

( 12)

V'1>= -'~p. Para la carga puntual, la densidad de carga toma la forma

(13)

654

E. Ley-lioo y A. Gngom. T.

p(r)
y el potencial y la intensidad y (2), respectivamente.

ql>(r - r').

(14)
de las ecs. (1)

de campo toman la forma coulombiana

El carcter lineal de las Ecs. (13) y (10) en el potencia! y en la intensidad de

campo elctrico, respectivamente, permite construir las soluciones correspondi~ntes aplicando el principio de superposicin a las soluciones coulombianas de las Ec5. (I) y (2) para cada elemento de carga p(r') d1r':

'( ) = >!' r E( ) =

r
para el potencial electrosttico:

J J

p(r') d'r'

Ir _ .-'1
Ir-.-'I'

(15 )

(r - r')p(r') d'r'

(16)

El uso de la Ec. (31m) en la Ec. (15) permite escribir el desarrollo multipolar

1>( ) = "
r

L 21+
I,m

{(1~ p(r')};"W)d'r')
r

r/1+ 1

r'~ (') r
1m

+ ([

p(r')r"Y,;"W)d'r')

Y~~~t)}.

(151m)

vlido para cualquier punto del espacio. En cada trmino de la doble sumatoria sobre los ndices de multipolaridaJ hay dos trminos, el primero debido a la carga externa a una esfera ccnlrda en el origen y que pasa por el punto de inters, y el segundo a la carga en el interior de la misma esfera. Si la carga est localizada dentro de una esfera de radio a centrada en el origen, el potencial en los puntos externos a la esfera se reduce a

(15e)

en lrminos

de los momenlos

mullipolares

de la dislribucin

de carga inlerna

Q'Im--

f." (') "}"

o
pr r

(")d' lmr

, r.

(17e)

El desarrollo multipolar por Juera y por dentro

655

Anlogamente, si la carga est localizada afuera de una esfera de radio b centrada en el origen, el potencial cn los puntos en el interior de la esfera es

<P ;() r

= '" ~ 21h + 1 Q;lmr

I.m

'r1m(') r ,

(15i)

en trminos de los momcntos multipolares de la distribucin de carga externa

(17i) En c1yaso ?C la carga puntual, los momcntos multipolarcs se reducen a Qm = qr"}~(.r) y Q'm = q}:T1(f')/r'/+l, respectivamente. Las expresiones para la intensidad de campo elctrico se pueden obtener cal. culando el gradiente de las Ecs. (15c,i) o usando las Ecs. (4e,i), reemplazando los momentos multipolares de la carga puntual por los de la distribucin de nuestro inters, Ecs. (17e,i):

E'(r)

L 2/:
I.m

1 ~~; [i-(l + 1)Y,m(i-)

+ ii- x

lY,m(i")),

(18e)

( 18i)

Para destacar la unidad del desarrollo multipolar, consideremos la situacin en que la carga est distribuida sobre la superficie de una esfera r = a = b. La densidad de carga

(19)

indica a travs de los coeficientes N).II la intensidad de cada una de sus componentes de ffillltipolaridad definida. Los momentos rrultipolares de esta distribucin de carga se obtienen de las Ecs. (17c,i)

(20)
Q/m = t+i"Ntm.
a ; q

(21)

656

E. Ley.Koo y A. Gngora-T.

Entonces, las Ecs. (18e, i) describen el campo elctrico en todos los puntos por fuera y por dentro de la esfera, respectivamente. Notamos, una vez ms, que al atravesar la superficie de la esfera, r = a las componentes tangenciales de cada contribucin multipolar son continuas, y las componentes normales muestran la discontinuidad

debida a la densidad de carga superficial a travs de cada tina de sus componentes multipolares.

3. Magnetosttica El campo de induccin magntica B es un campo solenoidal

V'. B = O,
porque no existen cargas o polos magnticos aislados. La ley de Ampere

(23)

V' x B = -J,
e

4rr

(24 )

reconoce que dicho campo se origina en corrientes e1tricas, donde J(r) es el vector de densidad de corriente y, adems, esas corrientes son estacionarias en el caso magnetosttico,

V' . J

= O.

(25)

Dada una distribucin de corriente J(r), el mtodo general para determinar el campo de induccin magntica asociado en cualquier punto del espacio introduce un potencial vectorial A, de modo que

B = V' x A,

(26)

garantizando que B es solenoidal, como lo expresa la Ec. (23). Entonces, la ley de Ampere se convierte en la ecuacin para el potencial

El dC!l(Jrrolio multipolar por fuero y por dentro

657

\7 X (\7 X A) o bien,

11r e

-J,

\7(\7.

A) - \72 A = 11r J.
e

(27)

El potencial \'cclorial A no est definido unvocamente, pues siempre es posible agregarle el gradiente de una funcin escalar sin que sc modifiquc el campo de induccin magntica de la Ec. (26). Las transformaciorH's del potencial de este tipo sc conocen corno trilllsforll1ilciollCS de norma, y la libertad asociada a ellas puede aprovecharse para simplificar la Ec. (27) que debe satisfacer el potencial. Ulla norma, comnmente usada, lIamilda norma trilnsvcrsal, impone la condicin de que el potencial sea solelloicial

\'.

O,

(28)

de modo que la Ec. (27) se redllce a la ("('Ilacin de Poisson

(29)
La solucin de esta ecuacin para el potencial vectorial para el potencial escalar, Ec. (I.i). que satisface la ecuacin la sustitucin de las fucntes correspondientes
t..'S anloga de Poisson

a la solucin Ec. (13), con

A(r) = ~

J(r')d''r'

Ir- r'1
escribir el (Icsarrollo

(30)
multipolar

Una vez ms, el uso de la EL (31m) permite para el potcncial vectorial magntico:

A( ) = "~{
r ~
I,m

21 + 1

(~1O
e
r

J(r')YI:'W)d r'1+ I

r')

'y

1m r

(')

(301m)

658

E. Ley-Koo

A. Gngom- T.

vlido para cualquier punto {!f'1("Spacio. Si las corrientes estn localizadas dcntro de una esfera r = a, el primer trmino dentro del corchete se anula y el desa.rrollo multipolar por fuera se reduce a

(30e)

en trminos de los momentos multipolar<'S de la distribucin interior

de corrientes en el

ni '1' 1m ::;::-

11'
O

J( r ') r "}' 1m (")d'" r

r .

(31e)

Si las corrientes ('Stn localizadas fuera de una ("S(cTa r = b, el segundo trmino dentro del corchete de la Ec. (301m) se anula y el primero contribuye al desarrollo multipolar por denLro:

(30i)

en trminos de los momentos Illultipolarcs

de la distribucin de corrientes externas

M;

1m -

~;.oo
e
b

J( ')
r

Y,;"W) ,-1'+ 1 d" r , .

(3li)

Es til comparar las Ecs. (30c, i) del desarrollo multipolar para el pOlencial veclorial magntico con sus contrapartes Ecs. (I5c, i) para el potencial electrosttico. Se conoce que todos son desarrollos en las bases de funciones armnicas'y difieren en sus coeficientes, los cuales son los momentos tTlultipolares de las fuentes respectivas, es decir, las densidades de corriente Ecs. (:lIc, i) y las densidades de carga elctrica Ecs. (17e,i), respectivamente. La. forma genera) del campo de induccin magntica se obtiene de las Ecs. (26) y (30),

= ~\7 x

J(r'),J"r'

Ir-r'l

- _~
= ~

jJ(r')

x \7

(_1_) d"r' Ir-r'1

J(r') x (r - r')d.'r'

Ir - r'l"

El desarrollo multipolar por fuera y por dentro

659

Los desarrollos multipolares para el campo de la induccin magntica siguen de las Ecs. (26) y (30e, ). Antes de proceder a su clculo, vamos a suponer que las corrientes estn restringidas a la superficie de la esfera r = a. La densidad de corriente

J(r') = 8(r' - a)
a

:L K~.( -il')Y~.("),

~,.

(33)

incluye componentes (-il}l'll.l) de multipolaridad definida, cada una con intensidad l(ll.l. Para justificar la forma de cada componente se puede sealar que el operador vectorall, Ec. (8), garantiza que la corriente se mantiene sobre la superficie esfrica y que es estacionaria, Ec. (25). Vamos aprobar primero, que cada una de estas corrientes de multipolaridad definida produce un potencial vectorial con la misma multipolaridad. Algunas relaciones tiles para los clculos que siguen incluyen otras representaciones del operador 1 y los efectos de sus componentes sobre los armnicos esfricos:
t

(8') donde

Ix = 1, :1:i1.

" = J2 1 (":1: .") ex I lJ.

(8:1:)

I,Y/m = mlm, Ix lm

(9z)

q(/, m)lmx'

)(/:: m)(l:I: m

+ 1)lmx.

(9:1:)

Los momentos multipolares de cada componente multipolar de la distribucin de corriente, Ec. (33), se obtiene de las Ecs. (31et i), involucrando la misma integracin angular:

660

E. Ley./{oo y A. Gngora.T.

(31e') Mi fm

Mm
a 2f+J'

(31 i') se obti('lwtl al realizar las sumas de las Ec. (:JOe, i)

Los potenciales correspondientes

AA.I

L- 2/
I.m

4r.

+ 1 -;- l \A,jIU

-i ,

'+1,
U),f

(30e' )

(30i')

La comparacin de las Ecs. (30c, i') con los trminos de la Ec. (33), muestra que los potenciales por fuera y por dentro c.omparten la misma dependencia vectorial y angular con las corrientes dc la misma lTlultipolaridad. El desarrollo lTlultipolar se reducc a la siguiente forma:

A'( ) =
r

L- (2'\
A"
~

"

+ I)e

l'

'A"a

A+'(_ '1)
I

(l'A"(f))
rA+!

(30e")

t '

4r.

A (r) = L- (2'\

A."

+ I)e

g ),l

'

),,.

aA (-"Hr

lA,,).

(30i")

El desarrollo multipolar ],or juera Y por dentro

661

En el clculo del rotacional de estos potenciales vectoriales para obtener los campos de induccin magtica vamos a distinguir entre las componentes radiales y transversales B = rr. B + B - rr' n = rr. B - f x (f x B). La componente radial depende del operador

-ir. v x I = ~(-ir x V) . I = ~I. I = ~2. r r r

Para componentes trans\'{'rsait.s se necesita

(34)

irx Ir x (V x I)J''= ir x {r x IV x (-ir x V'')]} =rx =rx =rx {rx [rV2''+V''.Vr-V.rV''-r.VV,,J) {rx [v,,.r-3v,,-r:rV'']} {rx [-2v,,-r%rv,,]}

= -i~ x {r x (-iV) (1 + r :r) } ,,= -i~ x I (1 + r :r) ,,.

Al desarrollar el triple producto vectorial, pasando del primer rengln al segundo, se toma en cuenta el carcter de \7 como vector y como operador de derivada. Del segundo rengln al tercero desaparece el trmino del laplaciano debido a que su coeficiente i- x r se anula; y se han sustituido el gradiente y la divergencia de r. Del tercer regln al cuarto se reducen los terminos con el gradiante de "p. Del cuarto rengln al quinto se usa el conmutador ralar, \7] = -\7. En en quinto rengln se reconoce la presencia del operador 1. La aplicacin de los operadores de las Ecs. (34) y (35) a los desarrollados mul. tipolares del potencial vectorial, Ecs. (30e, i") conduce a los desarrollos correspondientes de la induccin magnt.ica

" 41r B;(r) __ L.,(2.\+I)c

f{,.r'-I [0 (

'.

a'

r.\ .\+1

Y,.(0) r -irxlY,.

(0)(

r I+A

')]

(36i)

Ntese la continuidad de las componentes normales de la induccin magntica al pa.

662

E. y-Koo yA. Gngom-T.

5ar de un lado a otro de la superficie esfrica y la discontinuidad en los componentes tangenciales

( /" ) K).Pir x IY).p(r)(~ 2+1co

+ 1 + ~),

(37)

que es proporcional y transversal a la componente multipolar correspondiente de la corriente superficial, Ec. (33). En aquellas regiones donde no hay corrientes, la Ec. (24) indica que el campo de la induccin es rrotacional J, por lo tanto, derivable de un potencial escalar magnetico

B;

-'V4>m.

(38)

Al sustituir esta expresin en la Ec. (23), se encuentra que el potencial escalar magntico debe satisfacer la ecuacin de Laplace

'V'4>m; O.

(39)

En la referencia [5] se discute como construir este potencial escalar magnetica

y su desarrollo multipolar por fuera de la regin de las fuentes. Aqu seguimos este

mtodo construyendo los potenciales y sus desarrollos multipolares por fuera y por dentro de las fuentes. De la Ec. (38) se puede escribir

(40)

donde la trayectoria de integracin de a a r es arbitraria, pero no debe pasar por dond haya corrientes, y la posicin inicial a se puede escoger a nuestra conveniencia. Si se usa la expresin integral del campo de induccin magntica en trminos de la corriente, Ec. (32), el potencial escalar se puede reescribir:

o/m r

~ ( ); _~1.' e J
c a d .

d'r'J(r') x (e - r')

le _

r'I'

(40')

Vamos a escoger la trayectoria de integracin como una recta en la direccin del vector de posicin del punto de nuestro inters ( = r( y a reescribir el integrando resultante

El desarrollo multipolar por fuera y por dentro

663

f. J(r')

(U - r')

f.r'

J(r')

1 (U - r'). r'

J(r')

IU - r'1'
= ~V' ( . 1 ~ l~r-r'l ) . r'
X

IU - r'13

IU - r'1'
_ 1 V'. (~' X J(r'). ~ l~r-r'1

J(r') = ~V'. ~

(r' ~ J(r')) l~r-r'l

Al sustituir en la expresin para el potencial escalar, Ec. (40'), la integral volumtrica del primer trmino con la divergencia se transforma en una integral de superficie mediante el teorema de Gauss y su valor es cero, pues la corriente est confinada al interior del volumen; para la integral del segundo trmino intercambiamos el orden de integracin obteniendo

1 <ilm(r)= -

J ' ('
V. r

J(r)

') '1.'
d3r

~ ~r-

l' d~ r'1"

(40")

El factor de la divergencia en el integrando juega el papel de una densidad de carga de magnetizacin anlogo a la densidad de carga de polarizacin -V. p en electrosttica (aj. A continuacin, distinguimos entre los casos externo e interno, tomando a en el infinito y a en el origen, respectivamente. En ambos casos usamos la Ec. (31m) para oblener

<iI:"(r) =

~
e

J V' .
L 21
I,m

(r'

X J(r'))

d r'

L 2tr+ 1r'IY,;"(r')V,m(f)
l,m

loo ~
'

r ~!2

= - "

~~

+1c

d3r'r"y'

(r')V' . (r'
1m

J(r'))

v,m(r') (l + 1 )r +1

(40.)

en trminos de los multipolos de la carga de magnetizacin interna

Mm
y anlogamente

= __c(_1_) 1+1

V', (r'

J(r'))r"Y;"(r')d3r',

(41e)

664

E. Ley-Koo yA.

Gngom-T.

~:'(r)

= -1
e

V' . (r' x J(r'd'r'

L 2/'1+ ~ Yo' ';;'+(i-' 1 1'm(i-) f.' d{{I-l 1 r o

l,m

=L
1m

2t:

1 MI,"r'Y'm(i:).
externa estn dados por

(40i)

donde los multipolos

de la carga de magnetizacin

H'

1m

=.!el
=

V' . ('

r X

J( ')} };"(i-'} d' ,

r
r,1+1

(41i)
de

Si el factor de la divergencia

en las Ecs. (41e,

i)

se toma como una densidad

carga de magnetizacin Pm

(1 fe)V' .(r' x J(r', de inmediato se notan las diferen-

cias de sus momentos multipolares en comparacin con los del caso electrosttico, Ecs. (l7e,i), dehido a los denominadores extras de -(1 + 1) Y 1, respectivamente. Para explicar el origen y la necesidad de estos denominadores vamos a calcular los momentos rnultipolarcs, Ees. (41 e, i), y potenciales escalares, Ecs. (40e, i), para la

distribucin de corriente de la Ec. (33).

Conviene reconocer que la carga de magnetizacin nos del operador I en la forma se puede reescribir en trmi-

PM =~ [(V' x r') . J - r' . V' x J]

(42) =Entonces, e

1, ') i( .,. i -(r xV .J=--Ir xV).J=--I.J.

e e

al usar la corriente de la Ec. (33) se obtiene

(42')

Los momentos

multipolares

Ecs. (41c, i) de esta distribucin

son

(41c')
y

(4li')

El desarrollo multipolnr por fuero y por dentro

665

En consecuencia, los desarrollos multipolares del potencial escalar, Ecs. (40e,i), se reducen a las formas

(40e')

4>' ( )
mr

= _ '"

f;;:

(2/

+ l)e

K'm(l a'

+ 1) '}' e)
r

1m r .

(40i')

El campo de induccin magntica se obtiene calculando el gradiente de estos potenciales, Ec. (38), encontrndose que es el mismo que ya se calcul mediante el uso del potencial vectorial, Ecs. (36e, i), corno debe de ser. Notamos que los potenciales escalares magnticos y sus componentes de multipolaridad definida, Ecs. (40e, i), muestran una discontinuidad al pasar del interior al exterior de la superficie esfrica, debido a la presencia de los factores diferentes 1 y -(1+ 1), que ya se seal en conexin con los momentos de la carga de magnetizacin. Esto se puede contrastar con el carcter continuo del potencial electrosttico en la situacin anloga de la esfera con carga elctrica distribuida sobre su superficie. La analoga elctrica tambin permite f('conocer que la discontinuidad de los potenciales escalares magnticos se puede interpretar romo asociado con una capa dipolar sobre la superficie esfrica. La necesidad de esa capa dipolar est impuesta por la diferencia entre el carcter solenoidal del campo de induccin magntica y el carcter irrotacional del campo de intensidad electrosttico, que se traduce en las diferentes condiciones de frontera en la superficie esfrica donde estn las fuentes para cada tipo de campo. La comparacin directa de las Ecs. (18e, i) y (36e, i) para los desarrollos multipolares de los campos de intensidad elctrica y de induccin magntica, tambin permite reconocer la necesidad de los factores extra de 1 por fuera y de -(1+ 1) por dentro, para ir de las componentes radiales discontinuas de E, Ec. (22), a las componentes radiales continuas de D, Ecs. (36e, i) y al mismo tiempo, de las componentes tangenciales continuas de E, Ecs. (l8e, i), a las componentes tangenciales discontinuas de D, Ec. (37). La Fig. 4 muestra las distribuciones de carga de magnetizacin de multipolaridades ms bajas rcescaladas por fuera y por dentro, y la Fig. 5 muestra las distribuciones de corriente superficial correspondientes. La comparacin con las figuras de los casos respectivos de electrosttica permite destacar que en el caso magntico no hay contribucin monopolar puesto que el trmino correspondiente con>. :;;;: O en las fuentes Ecs. (42) es nulo desde el principio; si las lneas de campo de Em Y Dlm por fuera coinciden en direccin y sentido, las lneas correspondientes de E Y Bim por (Ientro tienen la misma direccin y sentidos opuestos, y viceversa, asegurando que las lneas de B son cerradas y las lneas de E empiezan en cargas positivas y terminan en cargas negativas.
I

'm

666

E. Ley-Koo y A. Gngom- T.

N N S S S S N N N

N S N N N S S S N N S S N S N

N S N N S S N S S S

al

bl
superficiales:

el
a) dipolar YIO, b) ~uadrupolar

FIGURA 4 Distribuciones de cargas de magnetizacin Y20 y e), octupolar )'30.

------

al

bl

el

FIGURA 5 Distribuciones de corrientes supt'fficiales: a) dipolar, -IYIO, b) cuadrupolar, -ilY20 y e) octupolar, -jI}'30.

4. Discusin Se ha presentado un estudio del desarrollo multipolar en electrosttica y magnetosttica, destacando su unidad y validez general como lo ilustran las Ecs. (151m) y (301m). En el caso de una carga puntual, se ha reconocido la conveniencia de
interpretar el desarrollo multipolar tanto por fuera como por dentro, Ecs. (30e, i),

en trminos de distribuciones superficiales de carga con multipolaridad definida, Figs. 2 y 3, Y Ecs. (6/) y (6/m), en vez de la interpretacin lIsuallimitada al exterior en trminos de distribuciones de cargas puntuales como las de la Fig. l. Con respecto a esta ltima, se debe mencionar de paso que adems del centrado de los multipolos sucesivos, es necesario transformarlos a sus versiones puntuales para que cada uno de ellos corresponda a un ~olo trmino de la Ec. (3e)j de otro modo, los multipolos finitos de la Fig. 1 contribuyen no solo al trmino del potencial de la misma mul.

El de. mrrollo multipolar por Juera y por dentro

667

tipolaridad sino tambin a otros de multipolaridad ms alta. El uso de fuentes con multipolaridad definida distribuirlas sobre la superficie de una esfera, Ecs. (19), (33) Y (42'), Y Figs. 2-5, permite destacar que cada una de ellas da lugar a contribuciones a los potenciales, tanto por fuera romo por dentro, con la misma multipolaridad ilustrando as la unidad dd desarrollo multipolar, Ecs. (ISe, i) y (40c, i). Tambin se ilustran las discontinuidades en los campos de intensidad elctrica, Ec. (22), Y de induccin magntica, Ec. (37), asociadas a la..,fuentes superficiales respectivas. Aunque en este trabajo se ha utilizado una geometra esfrica para la distribucin de las fuentes, otras geometras pueden ser ms prcticas para la produccin de campos de las diferentes multipolaridadcs. En todo caso, nuestros resultados permiten determinar algunas de esas geometras alternativas. A manera de ilustracin, consideramos los campos dipolar y cuadrupolar en la situacin electrosttica, ana. lizando los trminos con 1 = 1 Y 2 de la Ec. (3i). En el primer caso, las superficies equipotencialcs son planos paralelos.

ePdr)

= --,,r

qrcosO

q ---;z,
r-

(3d)

sugiriendo la geometra de dos piaras cquipotcnciales paralelas para producir el campo uniforme entre ellas. En el segundo caso, las superficies equipotenciales

<p,( r)

,
qr , (3 cos' O _ 1) 2r'

= ...!L"
2,.

(2z' _ x' _ y')

(3c)

se identifican como hiperholoi(lcs de revolucin de dos hojas cuando la cantidad dentro del parntesis es positiva, y de una hoja cuando esa cantidad es negativa. Correspondientemente, se sugieren electrodos ("(]uipotenciales con las polaridades apropiadas y con la forma de esos hiperboloides para generar el campo cuadrupolar. En el caso magnetosttico, el campo de induccin magntica uniforme se puede producir con un emhobinado sobre una esfera romo el de la Fig. Sal, o entre los polos paralelos de un imn sugeridos por la analoga electrosttica de la Fig. 4a), O con un embobinado uniforme sobre un cilindro recto. A continuacin, describimos algunos dispositivos para controlar el movimiento de partculas cargadas en su interior, poniendo nfasis en el carcter multipolar de los campos utilizados. En un selector de partculas cargadas de acuerdo con sus velocidades, se usan un campo elctrico uniforme E = -iE Y un campo magntico uniforme B = kB, perpendiculares entr s; las partculas que se mueven a travs de la regin de estos campos en la direccin perpendicular a ambos v = jv, estn sujetas a la fuerza F = q(-iE + ivil/c), la cual se anula para aquellas partculas con velocidad v = cE/ /J. En un selector de partculas cargadas de acuerdo con sus cantidades de movimiento; se usa un campo magntico uniforme en el cual las partculas describen trayectorias circulares: (q/c)v x B' = -mv2R/ R, de modo que la cantidad de movimiento de una partcula es proporcional al radio de su trayectoria mv = qlJ' !l/c. En un espcctrmetro de masas, se combinan las acciones sucesivas

668

E. Ley-Koo yA. Gngom-T.

,,
\

,,
\

, ,,

al

bl

el
nf"ulras: a)"cuadrupolar", b) "octupolar" y y "octupolar", y transversal "cuadrupolar".

FIGURA 6 Trampas magnetostticas para partculas e) hbrida con campos axiales "uniforme"

de seleccin de \'e1ocidadcs y seleccin de cantidad de movimiento para determinar la razn de carga a masa de las partcula.. bajo estudio q/m = Ec2/ BB' R. En una trampa de Pcnning, se logra el confinamiento de partculas cargadas con movimiento armnico bajo la accin combinada de un campo cuadrupolar elctrico, Ec. (3c), y un campo magntico uniforme [7]. Este tipo de trampa se ha utilizado desde hace poco ms de una dcada para confinar electrones individuales [8] y, posteriormente, tambin iones atmicos [9] para estudios espectroscpicos de la ms alta precisin. En todos estos dispositivos se puede apreciar la importancia de la pureza en la multipolaridad de los campos involucrados para definir la dinmica de las partculas cargadas y mejorar la precisin de las mediciones. Recientemente, se han desarrollado trampas para confinar partculas neutras, como neutrones ultrafros [10] y tomos enfriados con lser [11]. El funcionamiento de estas trampas se basa en la accin de sus campos magnticos no uniformes sobre los momentos magnticos de las partculas, F = J1. VD. La Fig. 6 muestra tres tipos de trampas magnticas con las que se ha logrado el confinamiento de partculas neutras [12]: a) la trampa "cuadrupolar" est formada por dos espiras coaxiales idnticas con corrientes en sentidos opuestoS b) la trampa "octupolar" consiste de tres espiras coaxiales colocadas sobre la superficie de una esfera en las posiciones 0= 1r/4, 1r/2 Y 31r/4, con la corriente en la espira ecuatorial igual y en el sentido opuesto a la corriente en las otras dos espiras: y c) la trampa hbrida tiene dos espiras coaxiales y cuatro conductores rectores y paralelos igualmente espaciados, en la cual las corrientes iguales y en el mismo sentido en las espiras contribuyen con campos

El desarrollo multipolar por fuera y por dentro

669

axiales "uniforme" y "octupolar", y las corrientes iguales y con sentidos alternados en los conductorcs rectos contribuycn con un campo transversal "cuadrupolar'. En esta descripcin de las trampas, hemos entrecomillado las multipolaridades de los campos para destacar que esas son las contribuciones al orden ms bajo, puesto que las espiras y conductores rectos individuales contribuyen con campos de mul. tipolaridades ms altas. La presencia de los ltimos indudablemente complica la dinmica de las partculas confinadas. Concluimos sealando que el diseo de este tipo dc trampas y el estudio de la dinmica de las partculas neutras en su interior son problemas abiertos a la investigacin en la actualidad, y que el papel de la rnultipolaridad de los campos es importante como lo seala la experiencia con el caso de partculas cargadas.

Referencias l. 2. 3. 4. 5. 6. 7. 8_ 9. D. lIalliday and R. Resnick, Physics for Students o/ Science and Engineering, Wiley, New York (1960), p. 722. J. R. Reitz, F.J. 1.1ilford and R.W. Christy, Foundations of Electromagnetic Theory, Addison-lVesley, Reading (1979), p. 41. J. D. Jaekson, Classical Elcctrodynamics, Caps. 4 y 5, Wiley, New York (1962). L. Eyges, The Classical Elcclromagnetic He/d, Dover, New York (1972), p. 30. .1. n. Ilronzan, Am. J. l'hys. 39 (1971) 1357. G. Arfken, it Mathematieal Methods for Physicists, Second Edition, Aeademic Press, New York, (1970). E. Ley-Koo yA. Gngora, T., Rev_ Mex. Fis. 27 (1980) 69. R.S. Van Dyek, Jr., P. B. Sehwinberg y !l. G. Dehmelt, Phys. Rev. Letl_ 38 (1977) 310. D. A. Chureh, R. A. Kenefiek, IV. S_ Buens, C.S.a.R. !lolmes, S. !luid!, S. Berry, M. llreing, S. EIston, J.P. Rozet, J. A. Sellin, D. Taylor and B. Thomas. Phys. Rev. Letl. 51, (1983) 1636. K. J_ Kugler, IV. Paul and U. Trinks, Phys. Lett. 72B (1978) 422. A.L. Migdall, J.V. Prod,l", IV.D. Philiips, T.Il. Bergeman and 1!.J. Meteal!, Phys. Rev. Letl. 54, (1985) 2596. T. Bergeman, G. Erez, n.J. 1.fetcalf,"1'1agnetostatic Trapping Fields foc Neutral Atoms" SUNY, Stony Brook, Preprint (1986). Abstract. \Ve emphasize the validity and unity of the multipole ex. pansion not ollly outside but also insirle a region where the sources of eleclrostatie or magnetostatic ficld are distributed. \Ve illuslrate the unity of the expansion and the eharacteristies of each of its terms by recognizing that sources of a definite multipolarity distribuled on a spherical surface give rise to fields with the corresponding multipola. rity both outsidc and nside the sphere. \Ve give examples of traps for chargcd particles and for neutral particles with a magnetie moment, in which the multipole character oC the fields in their interior is crucial foc the dynamics of those particles.

ID.
11. 12.