CAPÍTULO 1

CAMPOS ELECTROSTÁTICOS EN EL VACÍO

Los conceptos de carga eléctrica y campo eléctrico son fundamentales para la
descripción de la interacción entre partículas cargadas en reposo. En el presente capítulo se
discuten las leyes fundamentales que gobiernan la interacción eléctrica entre cargas
estacionarias situadas en el vacío, prestando especial interés a las características
diferenciales del campo eléctrico (divergencia y rotacional), las cuales se determinarán a
partir de la ley experimental que establece la interacción entre distribuciones estacionarias
de carga (ley de Coulomb). Se analiza el comportamiento del campo eléctrico a grandes
distancias de un sistema de cargas estacionarias (expansión multipolar) y se demuestra que
la existencia de un campo eléctrico en una región del espacio lleva a la aparición de una
densidad de energía eléctrica.

§ 1.1. LA CARGA ELÉCTRICA. LEY DE COULOMB.

La evidencia experimental muestra la existencia de un atributo fundamental de las

partículas elementales denominado carga. Se sabe que existen dos clases de cargas:
positivas y negativas1. Dos objetos que tienen exceso de uno de estos tipos de carga ejercen
entre ellos una fuerza repulsiva cuando se acercan entre sí. Dos objetos que tienen exceso
de cargas opuestas, uno cargado positivamente y el otro negativamente, se atraen
mutuamente cuando están cercanos.
Se ha establecido experimentalmente que la carga eléctrica tiene las siguientes
propiedades:
Conservación: En todo sistema cerrado la carga eléctrica neta es constante. Por lo tanto,
la suma algebraica de las cargas fundamentales es la misma en todo proceso físico.
Cuantización: La carga de todo sistema se puede representar como un número entero de
cierta unidad de carga fundamental. Experimentalmente se ha determinado que la menor
partícula estable que posee dicha carga fundamental negativa es el electrón. Así mismo, la
menor partícula estable con carga positiva fundamental es el protón2. Para ambas partículas,
el valor absoluto de la carga (en el sistema de unidades SI) es

1,6×10-19 C,

donde la unidad de carga C (denominada coulombio) se definirá más adelante.

1
Por razones históricas, se denominan cargas positivas a las cargas cuyo signo coincide con las cargas que
surgen en el vidrio frotado por seda. Así mismo, se denominan cargas negativas a las cargas cuyo signo
coincide con las cargas que aparecen en la ebonita, frotada por cuero.
2
El protón no es una partícula elemental, ya que está constituido por entidades más elementales denominadas
quarks, los cuales tienen una carga que es fracción de la carga fundamental. Sin embargo, los quarks no se
han observado en estado libre, por lo que la carga elemental positiva se asigna al protón.
CAPÍTULO 1. CAMPOS ELECTROSTÁTICOS EN EL VACÍO

No se puede dar una definición cuantitativa del concepto de carga hasta que se
establezca la ley de interacción entre partículas con dicha característica. Para simplificar la
discusión de la interacción entre cargas, se considerarán inicialmente las propiedades y las
interacciones que tienen lugar entre cargas eléctricas inmóviles con respecto a un sistema
inercial común de referencia.
Supóngase que se tienen dos cuerpos con cargas q1 y q2 y caracterizados por los
radios vectores r1 y r2, tal como se ilustra en la Figura (1.1.1). Si las dimensiones lineales
de los cuerpos son mucho menores que las demás distancias características involucradas en
el problema, se puede considerar que las cargas son puntuales. Experimentalmente se ha
establecido que entre cargas puntuales en reposo respecto a un sistema inercial de
referencia actúa una fuerza de interacción central, proporcional al producto de las cargas,
e inversamente proporcional al cuadrado de la distancia relativa entre ellas3, es decir,

Figura (1.1.1)

1 q1 q 2 r12
F12 = , (1.1.1)
4πε 0 r122 r12

donde r12=r1-r2 es el radio vector dirigido desde la partícula 2 hacia la partícula 1, tal como
se muestra en la Figura (1.1.1). Nótese que la elección del origen de coordenadas O no
afecta la relación (1.1.1) ya que el radio vector r12 no depende de la posición del punto O.
La expresión (1.1.1) constituye el enunciado matemático de la ley de Coulomb. El
valor de la constante de proporcionalidad en esta fórmula es igual a

1 N × m2
= 9× 109 . (1.1.2)
4πε 0 C2

3
Desde el punto de vista de la teoría cuántica de campos, en la interacción coulombiana se intercambian
fotones. Un análisis elemental basado en la relación de indeterminación para la energía muestra que el largo
alcance de la interacción coulombiana reflejado en la dependencia 1/r2 en la magnitud de la fuerza de
interacción está íntimamente ligado a que la masa de los fotones es nula (Greiner, 1998).

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

La cantidad ε0 se denomina permitividad dieléctrica del espacio vacío.

Si se tienen más de dos cargas, la fuerza que actúa sobre la carga i-ésima es igual a
la suma vectorial de las fuerzas que ejercen las demás cargas sobre ella, es decir,

1 N q j rij
Fi = qi ∑ , rij= ri- rj. (1.1.3)
4πε 0 j ≠i rij2 rij

En la anterior fórmula está implícita una hipótesis no trivial: la fuerza que ejerce la j-ésima
partícula sobre la i-ésima es independiente de la presencia de las otras cargas que están
presentes en el sistema. Esta aseveración constituye la esencia del llamado principio de
superposición, el cual juega un papel fundamental en la electrodinámica.
El anterior análisis se puede extender fácilmente al caso cuando la distribución de
cargas se puede considerar continua. La posibilidad de prescindir del carácter discreto (o
discontinuo) de la carga de un sistema es factible cuando se considera un número enorme
de cargas elementales. En este caso, se pueden definir tres tipos de distribución continua de
cargas:
Distribución volumétrica: Si las cargas se
distribuyen en un volumen V determinado del
espacio (Figura (1.1.2)), dicha distribución se puede
caracterizar introduciendo la densidad volumétrica
de carga ρ(x´,y´,z´) , la cual se define como la razón
entre la carga dq contenida en un volumen
infinitesimal dV´ y dicho volumen, es decir,

dq
ρ= . (1.1.4)
dV ´

Se entiende que el volumen dV´ debe ser lo

suficientemente pequeño como para considerar que
la densidad dentro de sus límites es homogénea, pero
lo suficientemente grande como para permitir
despreciar el carácter discreto que introduce la Figura (1.1.2)
estructura atómica.
La función ρ(x´,y´,z´) es un campo escalar con la ayuda del cual se puede hallar la
carga total contenida en el sistema mediante la relación

q = ∫ ρ ( x´, y´, z´)dV ´ .

V

41
CAPÍTULO 1. CAMPOS ELECTROSTÁTICOS EN EL VACÍO

Distribución superficial: En completa analogía

con las consideraciones anteriores, si las cargas se
distribuyen sobre una determinada superficie S
(Figura (1.1.3)), dicha distribución se puede
caracterizar introduciendo la densidad superficial de
carga σ(x´,y´), la cual se define como la razón

dq
σ= , (1.1.5)
dS´

donde dq es la carga contenida en la superficie

físicamente infinitesimal dS´.
La carga total contenida en el sistema asociada a esta
distribución es
q = ∫ σ ( x´, y´)dS´ .
S

Figura (1.1.3)
Distribución lineal: Si las cargas se distribuyen
a lo largo de una línea L (Figura (1.1.4)), la
respectiva distribución se caracteriza introduciendo
la densidad lineal de carga, la cual se define como

dq
λ= , (1.1.6)
dl´

donde dq es la carga contenida en el elemento de

longitud físicamente infinitesimal dl´. La carga total
contenida en el sistema asociada a esta distribución
es
q = ∫ λ (l´)dl´ .
L
Figura (1.1.4)
La expresión de la Ley de Coulomb en el caso de distribución continua de cargas es
muy sencilla de establecer. Se divide el sistema en elementos infinitesimales de carga dq’;
en virtud del principio de superposición, la fuerza que dicho sistema ejerce sobre una carga
q es

q dq ' rqq'
Fq = ∑r
qq ' rqq '
2
4πε 0
,

donde rqq´ es el radio vector dirigido desde la partícula con carga dq´ hacia la partícula con
carga q. Si en el sistema hay cargas distribuidas volumétricamente entonces dq´=ρdV´ y por
lo tanto se tiene:

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

q (r − r ') ρ (r ' )dV ' .

Fq =
4πε 0 ∫ r − r'
V
3
(1.1.7a)

Así mismo, si en el sistema las cargas están distribuidas superficialmente, entonces

dq´=σds´ y por lo tanto la fuerza que ejerce esta distribución de cargas sobre la carga q es

q (r − r ') σ (r ' )ds' .

Fq =
4πε 0 ∫ r − r'
S
3
(1.1.7b)

Finalmente, si la distribución de cargas es lineal con densidad λ, entonces

q (r − r ') λ (r ' )dl ' .

Fq =
4πε 0 ∫ r − r'
l
3
(1.1.7c)

Estas expresiones generalizan para una distribución continua de cargas los

resultados experimentales establecidos por Coulomb para la interacción electrostática entre
cargas puntuales.

§ 1.2. CONCEPTO DE FUNCIÓN δ DE DIRAC.

En el punto anterior se indicó la forma general de calcular la interacción entre

distribuciones continuas de cargas y para ello se introdujo el concepto de densidad de
carga. Una pregunta que se puede plantear es cómo definir la densidad (ya sea lineal,
superficial o volumétrica) de una carga puntual. La dificultad que plantea esta pregunta
radica en el hecho que una carga puntual (por definición) ocupa una región del espacio de
medida cero, y por lo tanto no se puede utilizar ninguna de las definiciones (1.1.4-1.1.6).
Sin embargo, se puede llegar a la solución del problema planteado examinando una
distribución continua en una cierta región y después hacer tender la medida de dicha región
a cero.
Considérese una carga de 1C distribuida uniformemente a lo largo de un intervalo
de longitud 2ε. La densidad lineal de esta carga es entonces λ=1/(2ε). Se propone examinar
el comportamiento de esta magnitud cuando la longitud del intervalo tiende a cero (esto
corresponde al paso al límite para una carga puntual). De la Figura (1.2.1) se ve que a
medida que disminuye la longitud del intervalo, la densidad lineal de carga se va haciendo
cada vez más grande, pero de tal manera que el área bajo la curva λ(x,ε) (que corresponde a
la carga total del sistema) permanece constante. Por lo tanto, en el límite ε→0, la
respectiva densidad lineal de carga λ(x) cumple las siguientes propiedades:

0, si x ≠ 0
λ ( x) =  . (1.2.1)
∞, si x = 0

43
CAPÍTULO 1. CAMPOS ELECTROSTÁTICOS EN EL VACÍO

Figura (1.2.1)

Pero como la carga total es invariante, debe cumplirse la relación

∫ λ ( x)dx = 1 . (1.2.2)

Toda función que satisfaga relaciones del tipo (1.2.1-1.2.2) se denomina función δ-
de Dirac. En una dimensión, la definición formal es entonces:

 0, si x ≠ 0
δ ( x) = 
∞, si x = 0
, pero tal que ∫ δ ( x)dx = 1 . (1.2.3)

Los límites de integración en la expresión ∫ δ ( x)dx = 1 son arbitrarios, siempre y cuando el

punto x=0 esté dentro de dicho intervalo de integración.
Con base en la anterior definición, se puede entonces escribir para la densidad lineal
de una carga puntual q que se encuentra en el origen de coordenadas la siguiente relación:

λ ( x ) = qδ ( x ) .

La generalización de esta expresión al caso cuando la carga se encuentra en una posición x0

diferente del origen de coordenadas, es obvia:

λ ( x ) = qδ ( x − x0 ) .

Otras representaciones de la función δ de Dirac. Cabe anotar que al concepto de

función δ de Dirac se puede llegar examinando otro tipo de densidades lineales, diferentes a
(1.2.1). En las Figuras (1.2.2-1.2.5) se indican algunas funciones que se aproximan a la
función δ de Dirac cuando el parámetro ε tiende a cero.

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

Figura (1.2.2) Figura (1.2.3)

Figura (1.2.4) Figura (1.2.5)

Otros ejemplos de representaciones de la función δ de Dirac son los siguientes:

 x
exp 
1 1  x 1 ε 
δ ( x) = lim sin   , δ ( x ) = − lim .
π ε → 0 x ε  ε →0 ε  x 
2

1 + exp ε 
  

Es fácil comprobar que cualquiera de estas representaciones posee las propiedades

especificadas en (1.2.3).

45
CAPÍTULO 1. CAMPOS ELECTROSTÁTICOS EN EL VACÍO

Ejemplo 1.2.1. Verifíquese que lim  ε 1 sin 2  x   = δ ( x ) (Figura (1.2.5)).

ε →0 π x 2 ε 
  
Si x=0, λ(x,ε)=1/(πε) y diverge cuando ε→0; por lo tanto, el comportamiento de la función
considerada en x=0 es semejante al de δ(0). Así mismo, si x≠0, la función λ(x,ε) oscila fuertemente
alrededor de cero con amplitud decreciente y cuando ε→0, se obtiene λ(x,0)=0, que es el
+∞
comportamiento de δ(x) cuando x≠0. Finalmente, hay que verificar que lim λ ( x, ε )dx = 1 ; se tiene:
∫ ε →0
−∞

∞ ∞ 2 ∞
ε 1 2  x 1 ε  2x  x  1 1
∫−∞ x 2 sin  ε  dx = π lim
ε →0 ∫ x
  sin   d   = ∫ 2 sin ( y )dy .
2
lim
ε →0 π
− ∞   ε   ε  π −∞ y

La última integral se puede consultar en una tabla de integrales o se puede determinar mediante
programas como MATHEMATICA: el respectivo valor de dicha integral es igual a π. Por lo tanto,
+∞
queda verificado que lim λ ( x, ε ) dx = 1 , lo que corresponde a la última propiedad de la función δ de
∫ ε →0
−∞
Dirac.

Un procedimiento similar al expuesto en el anterior ejemplo se puede aplicar para verificar

que las demás funciones λ(x,ε) aquí introducidas efectivamente son representaciones de la
función δ de Dirac.
Otras propiedades de la función δ de Dirac. Se examinan ahora algunas
propiedades de la función δ de Dirac, las cuales se pueden demostrar fácilmente a partir de
la definición (1.2.3):

δ ( x ) = δ (− x ) , (1.2.3a)
1
δ (ax) = δ ( x) (1.2.3b)
a
 f (a) si a ∈[x1 , x 2 ]
x2

∫ f ( x)δ ( x − a)dx = 
x1
0 si a ∉[x1 , x 2 ]
, (1.2.3c)

1
δ (x 2 − a 2 ) = {δ (x + a ) + δ (x − a )}, (1.2.3d)
2a
d d
δ (− x ) = − δ ( x ) , (1.2.3e)
dx dx

d
x δ (x ) = −δ (x ) . (1.2.3f)
dx

donde a es una constante arbitraria y f(x) es una función arbitraria definida en el punto x=a.
La primera de estas relaciones es obvia, y refleja simplemente el hecho que la función δ es
par.

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

1 1
Para demostrar la segunda relación, se escribe: ∫ δ (ax )dx = ∫ a δ (ax )d (ax ) = ∫ a δ (x´)d (x´) ,
donde se ha hecho x´=ax. En la última integral se puede quitar la prima sobre la variable x.
1
Se obtiene por lo tanto ∫ δ (ax )dx = ∫ δ ( x )dx , de donde se deduce de inmediato (1.2.3b).
a
La propiedad (1.2.3c) se demuestra fácilmente si se tiene en cuenta que la función
δ(x-a) es diferente de cero sólo en x=a, y en el entorno de este punto f(x) es suave y
aproximadamente igual a f(a), por lo que se tiene ∫ f ( x)δ ( x − a)dx =
f (a )∫ δ ( x − a)dx =f(a). Es necesario enfatizar que los límites de integración en la expresión
(1.2.3b) son arbitrarios, siempre y cuando el punto x=a esté dentro de dicho intervalo de
integración. De lo contrario, la integral respectiva se anula.
Para demostrar la relación (1.2.3d) se considera la integral ∫ f ( x )δ ( x 2 − a 2 )dx y se
hace en ella el cambio de variable y=x2-a2. Se tiene entonces
dx
∫ f ( x )δ ( x 2 − a 2 )dx = ∫ f ( x ) δ ( y )dy . Teniendo en cuenta que dx/dy=1/(2x), y que
dy
x = ± y + a 2 , se obtiene

∫ f ( x)δ ( x − a 2 )dx =
2

 1   1 
∫ 2 y + a2
 f y a  ( y )dy ∫  2 − f (− y + a 2 )δ ( y )dy .
2
( + ) δ +
 
  ( y + a2 ) 

Usando ahora la propiedad (1.2.3c) en las integrales de la izquierda en la expresión

1
anterior, se obtiene ∫ f ( x )δ ( x 2 − a 2 )dx = [ f ( a ) − f (− a )]. Sea g(x0)=f(x0), donde
2a
x0=±a. De acuerdo con (1.2.3c) g(x0)=∫g(x)δ(x-x0)dx, es decir,

1 1
∫ f ( x)δ ( x − a 2 )dx = ∫ f ( x)δ (x − a )dx + 2 a ∫ f ( x)δ (x + a )dx .
2

2a

Dado que la función f(x) es arbitraria, de aquí se obtiene de inmediato la relación (1.2.3d).
Las fórmulas (1.2.3e-f) involucran la derivada de la función δ(x) con respecto a su
argumento. La primera de estas fórmulas es consecuencia directa de la paridad de la
función δ . Para entender la forma como actúa la derivada de δ, considérese la integral
d
∫ f ( x) dx [δ ( x)]dx y efectúese en ella una integración por partes. Se obtiene
d  df ( x )   df ( x) 
∫ f ( x) dx [δ ( x)]dx = [ f ( x)δ ( x)] x=0 -∫
 dx  δ ( x)dx =-   .
 dx  x =0

47
CAPÍTULO 1. CAMPOS ELECTROSTÁTICOS EN EL VACÍO

Tomando en particular f(x)=x, se ve que

d
∫ x dx [δ ( x)]dx =−1≡− ∫ δ ( x)dx ,
de donde se obtiene de inmediato (1.2.3f).
El procedimiento seguido para demostrar la relación (1.2.3d) se puede aplicar para
demostrar la siguiente fórmula:

δ (x − a i )
δ ( g ( x )) = ∑ , (1.2.3g)
i g´(a i )

donde g(x) es una función univaluada de x; a1, a2,... , an son los ceros de la función g(x).
Nótese que las relaciones (1.2.3b) y (1.2.3d) son casos particulares de (1.2.3g),
correspondiente a las elecciones particulares g(x)=ax y g(x)=x2-a2.

Ejemplo 1.2.2. Los siguientes ejemplos ilustran las propiedades enunciadas:

(i)δ(±x/5)=5δ(x)- (propiedades (1.2.3a,b))
(ii) δ(x2-25)=[δ(x-5)+δ(x+5)]/10 - (propiedad (1.2.3d))
2
(iii) ya que el punto x=-1 no pertenece al intervalo [0,2].
∫ x δ (x + 1)dx = 0 . La integral es nula
2

0
2
(iv) = 1 . La integral es diferente de cero, ya que el punto x=1 sí pertenece
∫ x δ (x − 1)dx = x
2 2
x =1
0
al intervalo [0,2].

Extensión al caso tridimensional. Al concepto de función δ de Dirac se llegó

considerando una distribución unidimensional de carga. Para la densidad volumétrica de
una carga q se puede escribir, en completa analogía con la situación lineal, la siguiente
expresión

ρ(r)=qδ(r-r0),

donde se entiende que δ(r-r0) es el producto de funciones

δ(r-r0)≡ δ(x-x0)δ(y-y0)δ(z-z0).

Las propiedades (1.2.3a-f) se extienden fácilmente al caso tridimensional considerado. Por

ejemplo, la propiedad (1.2.3c) toma la forma ∫ f (r )δ (r − a )dV = f (a ) , donde el punto con
V
coordenadas (ax, az, az) está dentro del volumen de integración V.

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

Ejemplo 1.2.3. Un sistema consiste en dos cargas puntuales de cargas

–q y +q separadas entre sí una distancia l, tal como lo indica la
Figura-ej(1.2.1). La densidad volumétrica de carga en el punto M´
con coordenadas (x´,y´,z´) es
ρ(x´,y´,z´)=qδ (x´)δ (y´)δ (z´-l/2)-qδ (x´)δ (y´)δ (z´+l/2),
donde el primer término de la derecha corresponde a la densidad
volumétrica de la carga puntual positiva, y el segundo- a la
respectiva densidad volumétrica de la carga negativa. La expresión
obtenida se puede re-escribir como
ρ(x´,y´,z´)=qδ (x´)δ (y´)[δ (z´-l/2)-δ (z´+l/2)]. Figura-ej(1.2.1)

Ejemplo (1.2.4). Al evaluar

∫ r δ (r − a )dV , donde V se extiende a todo el espacio y a es un vector
2

∫ ∫
constante, se tiene: r 2δ (r − a )dV ≡ (r • r )δ (r − a )dV = (a • a ) ≡ a 2 .
V V

Divergencia del campo vectorial r/r3. En las discusiones sucesivas y en muchos

problemas de la física es preciso conocer la divergencia del campo vectorial
A=r/r3≡(1/r2)er. Si se usa la expresión para la divergencia en coordenadas esféricas y si se
tiene en cuenta que este campo vectorial es netamente radial, se tendría div(r/r3)=
1 d  2  1 
r   = 0 ; sin embargo este resultado no es correcto, ya que este campo es
r 2 dr   r 2  
singular en el origen de coordenadas.
Nótese que este campo vectorial surge en los problemas analizados en el parágrafo
anterior, ya que de acuerdo con (1.1.1) la fuerza de interacción entre dos cargas q1 y q2
qq r
situadas a una distancia r es F = 1 2 3 .
4πε 0 r
Considérese el flujo de este campo vectorial a través de una superficie cerrada
arbitraria S. Se tiene

2π π
r
∫ r 3 • dS = ∫0 dϕ ∫0 sinθdθ = (2π )[− cosθ ]0 = 4π .
π

Por otro lado al aplicar el teorema de Gauss a la parte izquierda de la anterior relación se
obtiene

r r
∫r 3
• dS = ∫ div
V'
r3
dV ' .

Por lo tanto, si se igualan estos dos resultados se llega a la siguiente relación:

r
∫ V'
div
r3
dV ' = 4π .

49
CAPÍTULO 1. CAMPOS ELECTROSTÁTICOS EN EL VACÍO

Si se tiene en cuenta que la integral sobre el volumen V’ de la función delta de Dirac es

igual a 1, se puede re-escribir esta última expresión de la siguiente forma:

r
∫
V'
div
r3
dV ' = ∫ 4πδ (r )dV ' .
V'

Igualando las expresiones sub-integrales, se obtiene

 r 
div  3  = 4πδ (r ) . (1.2.4)
r 

Este resultado se puede desarrollar un poco más, si se tiene en cuenta que grad(1/r)=-r/r3.
Por lo tanto, tiene lugar la siguiente relación para el laplaciano del inverso de r:

1
∇ 2   = −4πδ(r). (1.2.5)
r

Estos resultados serán utilizados más adelante cuando se analicen las características
diferenciales de los campos eléctrico y magnético.
Ejemplo 1.2.4. Evaluar ∇ 2 f (r ) , donde f (r ) = 1 exp − r  y λ es una constante. El cálculo mediante
r  λ
la expresión del laplaciano en coordenadas esféricas daría

1 d  2 d 1  r   1 d   r r  r   1  r .
r  exp −   = 2  − exp −  − exp −   = rλ2 exp − λ 
r dr  dr  r
2
 λ   r dr   λ λ  λ    

Si el parámetro λ se hace tender a infinito, la parte derecha de la última igualdad tendería a cero. No
obstante, lim ∇ 2  1 exp − r   = ∇ 2  1  = −4πδ (r ) , lo que muestra que en el cálculo realizado de
 
λ →∞
 r  λ   r
∇ f (r ) se ha perdido un término proporcional a la función delta de Dirac. Esto implica que de la
2

función f(r) es preciso aislar la parte singular 1/r; para ello, se escribe f (r ) = 1 exp − r  − 1 + 1 y
r   λ   r
por lo tanto

1 d  2 d  1   r    r
 exp −  − 1  + ∇ 2 = 2 exp −  − 4πδ (r ) .
1 1
∇ 2 f (r ) = r
r 2 dr  dr  r   λ    r rλ  λ

Solución particular de la ecuación de Poisson mediante el uso de la función δ de

Dirac. En el §0.6, al discutir la determinación de un campo vectorial con base en sus
características diferenciales, se mostró que en la solución de este problema es preciso
resolver la ecuación de Poisson (0.6.5)

∇2φ=f(r). ( 1.2.6)

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

Se puede encontrar una solución particular de esta ecuación mediante los conceptos aquí
desarrollados. En primer lugar, de las propiedades de la función δ de Dirac se deduce que la
función f(r) se puede representar en la forma f(r)=∫f(r´)δ(r-r´)dV´, donde la integral se
extiende a todo el volumen de la región donde está definida la función f. Por lo tanto, la
ecuación de Poisson se puede escribir como

∇2φ=∫ f(r´)δ(r-r´)dV. ( 1.2.7)

De acuerdo con los resultados obtenidos anteriormente (fórmula (1.2.5)),

1 2 1 
δ ( r − r ′) = − ∇  ,
4π  r − r′ 

donde el operador diferencial actúa sobre las variables no primadas. Por esta razón, se tiene
1   1 
∇ 2φ (r ) = − ∫ f (r ′)∇ 2   dV ′ ,

4π V   r − r ′ 

Dado que f(r´) se comporta como constante respecto a la actuación del operador de
Laplace, se puede escribir la anterior expresión de la siguiente forma:

 1  1  
∇ 2 φ (r ) + ∫ f (r ′)  dV ′ = 0 .

 4π V  r − r′  

Así, con exactitud de una constante aditiva, la solución de la ecuación de Poisson tiene
entonces la forma

1 f (r ′)
φ (r ) = −
4π ∫ r − r ′ dV ′ , (1.2.8)

la cual se reportó en (0.6.6) al estudiar la determinación de un campo vectorial mediante

sus características diferenciales (divergencia y rotacional).

51
CAPÍTULO 1. CAMPOS ELECTROSTÁTICOS EN EL VACÍO

§ 1.3. CAMPO ELÉCTRICO Y SUS CARACTERÍSTICAS DIFERENCIALES.

La interacción entre partículas cargadas inmóviles se realiza mediante el llamado

campo electrostático, el cual es un caso particular de una entidad física más fundamental, el
campo electromagnético, que caracteriza la interacción entre partículas cargadas que se
mueven arbitrariamente.
El campo eléctrico E(r) creado en un punto arbitrario r por un sistema estacionario
de cargas (al cual se denomina campo eléctrico estacionario o campo electrostático) se
define como el límite de la razón entre la fuerza Fq que actúa sobre una carga q (llamada
carga de prueba) y la magnitud de dicha carga, cuando ésta tiende a cero:

Fq ( r )
E (r ) = lim . (1.3.1)
q→0 q

El paso al límite garantiza que la introducción de la carga de prueba no altera la

distribución de las cargas que generan el campo. En el caso de una distribución volumétrica
de cargas en un volumen V’, de la definición (1.3.1) y de la expresión (1.1.7a) se ve que el
campo eléctrico tiene la forma

1 (r − r ') ρ (r ' )dV ' .

E (r ) = ∫ r − r'
4πε 0 V ' 3
(1.3.2a)

En la Figura (1.3.1) se especifica el sentido de cada uno de los términos que aparecen en
esta expresión: el elemento de carga dq=ρdV´, situado a una distancia R = r − r ' del punto
1 dq
M , genera en dicho punto un campo dE de magnitud dE = y orientado a lo largo
4πε 0 R 2
del vector unitario R/R.

Figura (1.3.1)

Si la distribución de cargas es superficial en una superficie S, de (1.3.1) y (1.1.7b) se puede

escribir:

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

1 (r − r ') σ (r ' )dS ' .

E (r ) =
4πε 0 ∫ r − r'
S'
3
(1.3.2b)

La expresión para una distribución lineal se escribe de manera inmediata:

1 (r − r ') λ (r ' )dl ' .

E (r ) =
4πε 0 ∫ r − r'
L'
3
(1.3.2c)

Ejemplo 1.3.1. Una carga q está distribuida

uniformemente a lo largo de un segmento de
longitud 2L. Hallar el campo eléctrico en un punto
arbitrario M.
Solución. La densidad lineal de carga es λ=q/(2L).
Dada la simetría del sistema, se puede elegir el
sistema de coordenadas con el eje z a lo largo del
segmento cargado y origen en su punto medio, tal
como lo muestra la Figura (1.3.1ej). Se toma un
elemento de longitud dl´=dz´ y radio vector r´=z´ez.
Se trata de determinar el campo dE generado por
este elemnto en el punto M definido por el radio
Figura (1.3.1ej).
vector r=xex+yey+zez.
Así, r-r´= xex+yey+(z-z´)ez, r − r ' = x 2 + y 2 + ( z − z´)2 ≡ ρ 2 + ( z − z´)2 donde ρ2=x2+y2. Por lo
tanto

1 (r − r ') λ (r ' )dl ' = 1 q  xe x + ye y + ( z − z´)e z 

dE (r ) =   dz´≡
4πε 0 r − r ' 3
[
4πε 0 2 L  ρ 2 + ( z − z´)2 3 / 2  ]
1 q  ρe ρ + (z − z´)e z  .
  dz´
[ ]
4πε 0 2L  ρ 2 + (z − z´)2 3 / 2 

El campo total se obtiene integrando con respecto a z´ entre –L y +L. Se obtiene:

1 q
E (r ) =
4πε 0 2 L
[I1(ρ , z )ρeρ + I2 (ρ , z )ez ],
donde4

L
1  (z − L ) − (z + L ) 
L /2
dz´ 1 (z − z´)
I1 (ρ , z ) = ∫ [ρ 2
+ ( z − z´) ]
2 3/ 2
=
ρ 2
ρ + ( z − z´)
2 2
= 2

ρ  ρ 2 + ( z − L )2 ρ 2 + ( z + L ) 
2
−L / 2
−L 
L
L/2
(z − z´)dz´ 1 1 1
I 2 (ρ , z ) = ∫
− L / 2 [ρ + ( z − z´) ]
2 2 3/ 2
=−
ρ + ( z − z´)
2 2
=
ρ + (z + L )
2 2
−
ρ + ( z − L )2
2
−L

4
Las integrales se calculan a partir de las fórmulas ∫dw/(w2+a2)3/2=(1/a2)[w/(w2+a2)1/2] y ∫wdw/(w2+a2)3/2=
-1/(w2+a2)1/2 con w=(z´-z) y a=ρ.

53
CAPÍTULO 1. CAMPOS ELECTROSTÁTICOS EN EL VACÍO

Si el punto M se sitúa en el plano xy, entonces z=0, I2(x,y,z)=0, I ( x, y , z ) = 2L y

1
(x 2
+y 2
) x +y +L
2 2 2

por lo tanto E (r ) = q  eρ 
 .
4πε 0  ρ ρ 2 + L2 

Si el punto M está suficientemente alejado del sistema (x,y>>L) se tiene E (r ) = q  e ρ  , lo que

4πε 0 ρ2 
 
corresponde al campo de una carga puntual q concentrada en el origen de coordenadas.

Habiendo definido el campo eléctrico mediante la relación general (1.2.3a), se trata

ahora de determinar las características diferenciales de este campo, es decir se propone
investigar el comportamiento de su rotacional y su divergencia.

Rotacional y circulación del campo electrostático. Se demostrará inicialmente

que el rotacional del campo eléctrico asociado a una distribución estacionaria de cargas
es igual a cero. Para simplificar, considérese una distribución volumétrica de cargas. Si se
tiene en cuenta que el operador diferencial actúa sólo sobre las coordenadas del vector r,
entonces

1    r − r '  

rotE (r ) = ∫  rot
 3

 ρ ( r ' ) dV '.
4πε 0 V '  r − r '  
   

El rotacional dentro del signo de integral se puede calcular aplicando la relación (0.5.13):
rot(φA)=φ(rotA)+(gradφ)×A, donde φ=1/r-r´3 y A=r-r´. Pero en este caso
rotA=rot(r-r´)=0 y gradφ=grad[1/r-r´ ]=-3[(r-r´)/r-r´ ]. Por lo tanto
3 5

 r − r'   
rot  = −3 r − r'  × (r − r' ) .
 r − r' 3   r − r' 5 
   

El producto vectorial en el miembro derecho de esta relación es idénticamente igual a cero,

y por lo tanto, el rotacional en la expresión subintegral considerada se anula, lo cual
demuestra que en efecto tiene lugar la relación

rotE(r)=0. (1.3.3)

Se ve entonces que el campo electrostático es un campo potencial, es decir, puede

representarse como el gradiente de cierta función escalar. Las consecuencias de este hecho
se examinarán más adelante, al estudiar la teoría del potencial.
Una consecuencia inmediata de la ecuación rotE=0 se obtiene si se integra esta
relación sobre una superficie S limitada por un contorno cerrado C: ∫ rotE • ds = 0 .
S
Aplicando el Teorema de Stokes (0.4.3), se obtiene

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

∫ E • dl = 0 ,
C

es decir, la circulación del campo electrostático a lo largo de un contorno cerrado

(denominada también fuerza electromotriz) es igual a cero.
Divergencia y flujo del campo electrostático. Se determinará ahora la divergencia
del campo eléctrico. Se tiene

1    r − r '  

divE (r ) = ∫  div  ρ ( r ' ) dV '.
4πε 0 V '   r − r ' 
3

   

r 
Recordando (véase (1.2.4)) que div  3  = 4πδ (r ) , se puede escribir
r 

1   1
divE (r ) =  ∫ δ (r − r ')ρ (r ' )dV ' = ρ (r ) .
ε 0 V '  ε0

Se obtiene entonces la siguiente ecuación para la divergencia del campo

electrostático E:
1
divE (r ) = ρ (r ) . (1.3.4)
ε0

Si se tiene en cuenta que la divergencia es una medida de la densidad de las fuentes que
generan el campo vectorial considerado, se ve que la relación (1.3.4) muestra que la fuente
del campo eléctrico es la densidad volumétrica de las cargas.
1
Las ecuaciones rotE(r)=0 y divE (r ) = ρ (r ) son las ecuaciones fundamentales de
ε0
la electrostática. De acuerdo con el resultado demostrado en el numeral §0.6, el
conocimiento de estas características diferenciales del campo electrostático permite en
principio determinar el campo E en todos los puntos de la región considerada.
Si se considera el flujo del campo eléctrico a través de una superficie cerrada, al
aplicar el Teorema de Gauss se tiene

1
∫ E • dS = ∫ divEdV = ε ∫ ρdV ;
S V o V

es decir, el flujo del campo electrostático a través de una superficie cerrada es

proporcional a la carga total encerrada por dicha superficie. Este resultado se conoce con
el nombre de Ley de Gauss. Es preciso indicar que la ley de Gauss se aplica cuando el
sistema considerado posee simetría plana, cilíndrica o esférica. En los siguientes ejemplos
convencionales se ilustra la metodología.

55
CAPÍTULO 1. CAMPOS ELECTROSTÁTICOS EN EL VACÍO

Ejemplo 1.3.2. Hallar el campo creado por

un sistema de cargas distribuidas
uniformemente en un plano infinito con
densidad superficial σ
Dado que el sistema es homogéneo en el
plano infinito, el campo eléctrico está
dirigido netamente en la dirección normal a
la superficie, tal como lo muestra la Figura
(1.3.2ej). En calidad se superficie
gaussiana se toma un cilindro de altura H y
área transversal A. Dado que el flujo del Figura (1.3.2ej)
campo eléctrico a través de la superficie
lateral del cilindro es igual a cero, sólo resta determinar los flujos a través de la base y la tapa del
cilindro. Dado que Q=σA es la carga contenida dentro de este cilindro, la ley de Gauss toma la
forma E×(2A)=(σA)/ε0, de donde se obtiene que E=(σ/2ε0)n. Nótese que este campo es homogéneo,
lo cual está relacionado con la aproximación de considerar infinitamente extendido el plano cargado.

Ejemplo 1.3.3. Un cilindro de radio R y altura h tiene

una carga Q distribuida con densidad volumétrica
ρ(r)=kr, donde r es la distancia al eje del cilindro5
(Figura (1.3.3ej)). Hállense los campos eléctricos dentro
y fuera del cilindro.
Solución. Nótese que la constante de proporcionalidad k
está relacionada con la carga total del cilindro; en efecto,
de la definición Q = ∫ ρdV se sigue
V
R h 2π
R 3 ; así, 3Q .
Q = k ∫ (r´)r´dr´∫ dz´ ∫ dϕ´ = 2πkh k=
0 0 0 3 2πhR 3 Figura (1.3.3ej)

Si se considera que la altura del cilindro es mucho mayor que su radio (h>>R), el campo eléctrico es
puramente radial.
Para hallar el campo eléctrico E1 dentro del cilindro, tómese una superficie gaussiana consistente en
un cilindro coaxial de radio 0<r<R y altura h, tal como lo muestra la Fig. La carga q(r) encerrada por
r h 2π
esta superficie gaussiana es q (r ) = k (r´)r´dr´ dz´ dϕ´ = 2πkh r . Dado que el campo eléctrico es
3

∫0 ∫0 ∫0 3
netamente radial, el flujo a través de la superficie cerrada se reduce al flujo del campo a través de la
superficie lateral del cilindro, el cual es igual a 2πrhE1. La ley de Gauss toma la forma
2πrhE1=q(r)/ε0 toma la forma 2πrhE1 = 1 2πkh r , de donde E1 = 1 k r . Recordando la
3 2

ε0 3 ε0 3
expresión obtenida más arriba para k, se obtiene E 1 = (Q / h) r e r .
2

2πε 0 R 3
Para hallar el campo eléctrico E2 fuera del cilindro, tómese una superficie gaussiana consistente en
un cilindro coaxial de radio r>R y altura h. La carga encerrada por esta superficie gaussiana es la
carga total Q contenida en el sistema; así, la ley de Gauss toma la forma 2πrhE2=Q/ε0. , de donde se
obtiene que E (r ) = (Q / h ) 1 e .
2 r
2πε 0 r

5
La distancia al eje del cilindro se designa r para evitar confusiones con el símbolo ρ que en el presente caso
designa la densidad volumétrica de carga.

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

Ejemplo 1.3.4. En los núcleos ligeros la distribución de carga se considera esférica y está dada por
ρ(r)=ρ0(1-r2/a2) si r≤a y ρ(r)=0 si r>a. Determínese el campo eléctrico E en el interior (r≤a) y en el
exterior (r>a) del núcleo6.
Solución. La carga total del núcleo es Q = ρ (r´)dV ´ , donde dV´=r´2dr´dΩ´, siendo
∫
V

dΩ´=senθ´dθ´dϕ´ el elemento de ángulo sólido en coordenadas esféricas. Dado que la distribución de

carga es isotrópica (no depende de θ´ y ϕ´), la integración sobre las variables angulares da 4π y por
a
lo tanto se tiene Q = 4πρ 1 − r´ r´2 dr´= 4πa 3 ρ  1 − 1  = 2 ρ V , donde V es el volumen del
2

0∫
a 2 
0 0
0  3 5 5
núcleo.
Debido a la simetría esférica de la distribución de carga, el campo eléctrico, tanto dentro como fuera
del núcleo es netamente radial. Si se toma una superficie gaussiana de radio r, el flujo del campo
eléctrico a través de ella es 4πR2E.
Si el radio r es menor que el radio del núcleo (r<a), la carga q(r) comprendida dentro de la superficie
r
gaussiana es q(r ) = 4πρ 1 − r´ r´2 dr´= 4πρ  r − r  . Si E1 es el campo eléctrico en el interior
2 3 5

0∫ 2  0 2 
0
a   3 5a 
del núcleo, la ley de Gauss toma la forma 4πr 2 E = 4πρ  r − r  , de donde
3 5

1 0 2 
 3 5a 
 1 r2 
E1 = ρ 0  − 2 re r .
 3 5a 
Fuera del núcleo, la carga encerrada por la superficie gaussiana de radio r>a es la carga total Q, por
lo que el campo E2 en esta región es E = Q 1 e , es decir, E = 2 ρ 0 a 1 e .
3

2 r 2 r
4πε 0 r 2 5ε 0 r 2

§ 1.4. POTENCIAL ELÉCTRICO. ECUACIÓN DE POISSON.

En la sección anterior se ha establecido que el campo eléctrico tiene carácter potencial, lo

que se refleja en la ecuación rotE=0. Pero se sabe (véase (0.5.6)) que el rotacional de todo
gradiente es idénticamente igual a cero. Por lo tanto, el campo electrostático puede ser
representado como el gradiente de una cierta función escalar:

E(r)=-gradΦ(r). (1.4.1)

La función escalar Φ(r) recibe el nombre de potencial escalar.

Figura (1.4.1)

6
Para el caso del calcio ρ0=5.0x1025 Coulombios/m3 y a=4,5x10-15 m.

57
CAPÍTULO 1. CAMPOS ELECTROSTÁTICOS EN EL VACÍO

De la relación (1.4.1) se obtiene inmediatamente un resultado fundamental: en un campo

electrostático, el trabajo realizado para llevar una carga de prueba q, a velocidad
constante, desde un punto A hasta un punto B es independiente del recorrido. En efecto,
considérese un contorno cerrado C que contiene los puntos A y B (fig. 1.4.1). El trabajo
total realizado es

W = ∫ F • dl = q ∫ E • dl = −q ∫ [gradΦ (r )] • d l ≡ −q ∫ dΦ = 0 .
C C C C

Pero la integral cerrada es igual a ∫

C
= ∫
C1
+ ∫
C2
, lo cual demuestra que el trabajo realizado

para llevar la carga q desde A hasta B es el mismo tanto si se sigue el camino C1 como si se
sigue C2.
De acuerdo con la ecuación (1.4.1),

E(r)•dl=-gradΦ(r)•dl=-[(∂Φ/∂x)dx+(∂Φ/∂y)dy+(∂Φ/∂z)dz]≡-dΦ(r).

Entonces la diferencia de potencial entre dos puntos A y B está determinada por la integral
de línea del campo eléctrico:
B

ΦA − ΦB = ∫ E ( r ) • dl .
A

Es necesario advertir que el campo eléctrico E(r) sólo determina la diferencia de

potencial entre dos puntos. Si se desea conocer el potencial en un punto dado, es preciso
definir un valor del potencial en un punto fijo tomado como referencia. Por ejemplo, se
puede suponer que el potencial en el infinito es igual a cero. Entonces, el potencial en el
punto 1 es igual a
∞

Φ1 = ∫ E ( r ) • dl .
1

En este caso, el trabajo necesario para traer una carga q desde el infinito hasta el punto 1 es
qΦ1
Considérese una carga puntual q en reposo situada en el origen de coordenadas. En
este caso el campo eléctrico es

1 q
E= r
4πε 0 r 3

y por lo tanto el correspondiente potencial es

∞ ∞ ∞ ∞
q r • dr q dr q  1
Φ ( r ) = ∫ E • dr = ∫ = ∫ = − r  ,
r 4πε 0 r r 3
4πε 0 r r 2
4πε 0 0

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

es decir

q 1
Φ (r ) = . (1.4.2)
4πε 0 r

Nótese que el potencial asociado a una carga puntual en reposo no depende de las variables
angulares θ y ϕ, lo cual es reflejo de la isotropía inherente a este sistema físico.
El anterior resultado se puede generalizar fácilmente a un sistema de cargas
distribuidas arbitrariamente. Para ello, téngase en cuenta la relación (1.1.4). Siguiendo un
procedimiento análogo al que permitió llegar a (1.4.2), se obtiene

1 ρ (r ' )
Φ (r ) =
4πε 0 V ∫ r − r ' dV ' ,
(distribución volumétrica de cargas) (1.4.3a)

1 σ (r ' )
Φ (r ) = ∫
4πε 0 S r − r '
ds' , (distribución superficial de cargas) (1.4.3b)

1 λ (r ' )
Φ(r ) =
4πε 0 ∫ r − r ' dl' .
l
(distribución lineal de cargas) (1.4.3c)

Sin embargo, hallar el potencial de un sistema de cargas mediante las relaciones

(1.4.3) se puede constituir en un problema complejo. Por esa razón es conveniente resolver
la ecuación diferencial que satisface la función escalar Φ(r), conociendo el
comportamiento de dicha función o de sus derivadas espaciales en las fronteras del sistema.
1
Si en la ecuación divE (r ) = ρ (r ) se reemplaza E(r)=-gradΦ(r) y se recuerda la
ε0
definición del operador de Laplace, se ve que el potencial electrostático satisface la
ecuación a derivadas parciales
1
∇ 2 Φ (r ) = − ρ (r ) . (1.4.4)
ε0

Esta ecuación recibe el nombre de ecuación de Poisson. Un caso particular de ésta es la

ecuación de Laplace

∇2Φ(r)=0, (1.4.5)

la cual es válida en las regiones del espacio donde no hay fuentes.

Desde el punto de vista matemático, la determinación del potencial electrostático
para una distribución dada de cargas se reduce a resolver la ecuación de Poisson (o la de
Laplace), imponiendo ciertas condiciones de frontera que definen el comportamiento de los
campos y los potenciales en ciertas regiones del espacio donde se rompe la homogeneidad
([Link]., en la frontera de división de dos medios) o donde los potenciales deben comportarse
adecuadamente. Este problema se abordará con algo de detalle un poco más adelante, en el

59
CAPÍTULO 1. CAMPOS ELECTROSTÁTICOS EN EL VACÍO

Capítulo 3. Sin embargo, para distribuciones simétricas de carga la solución de las

respectivas ecuaciones de Laplace y Poisson es relativamente simple.

Ejemplo (1.4.1). Condensador esférico. Considérense dos

esferas concéntricas de radios r1 y r2 (Figura (1.4.1ej)), las
cuales están sometidos a una diferencia de potencial V.
Debido a la diferencia de potencial, sobre las superficies
de las esferas se distribuyen las cargas +q y –q. Se toma
una superficie gaussiana consistente en una esfera
concéntrica de radio r1<r< r2. Debido a la simetría del
sistema, el campo es netamente radial y por lo tanto el
flujo del campo eléctrico a través de la superficie
gaussiana es 4πr2Er. De acuerdo con la ley de Gauss, este
flujo es igual a q/ε0; por lo tanto Eρ=[1/(4πε0)]q/r2. La
diferencia de potencial entre las placas del condensador es
r2 r
q 2 dr q  1 1 .
∆Φ = ∫ E r dr = ∫ =  − 
r1 4πε 0 r1 r 2
4πε 0  r1 r2 
La relación C=q/∆Φ entre la carga acumulada en una de
las placas del condensador y diferencia de potencial Figura (1.4.1ej)
aplicada se denomina capacitancia. Así, para el
sistema considerado la capacitancia es C = 4πε  r1r2  .
0 
 r1 + r2 

Ejemplo (1.4.2). Condensador cilíndrico. Considérense

dos cilindros coaxiales de altura h y radios r1 y r2
(Figura (1.4.2ej)), los cuales están sometidos a una
diferencia de potencial ∆Φ.
Debido a la diferencia de potencial aplicada, sobre las
superficies de los cilindros se distribuyen las cargas +q y
–q. Se toma una superficie gaussiana consistente en un
cilindro coaxial de radio r1<ρ< r2 y altura h. Si h>>r2, el
campo se puede considerar netamente radial y por lo
tanto el flujo del campo eléctrico a través de la
superficie que encierra al volumen del cilindro
gaussiano es 2πρhEρ. De acuerdo con la ley de Gauss,
este flujo es igual a q/ε0. Por lo tanto Eρ=[1/(2πhε0)]q/ρ.
La diferencia de potencial entre las placas del Figura (1.4.2ej)
condensador es
r2 r
1 q 2 dρ 1 q  r2  . De aquí que la capacitancia de este condensador sea
∆Φ = ∫ E ρ dρ = ∫ = ln 
r1 2πε 0 h r1 ρ 2πε 0 h  r1 
2πε 0 1 .
C=
h r 
ln  2 
 r1 

Ejemplo (1.4.3). Potencial y distribución de carga en un núcleo atómico. Considérese que la carga
Zq en un núcleo atómico de radio R está distribuida uniformemente. Esto significa que la densidad de
carga es ρ (r ) = 3Zq . Sea Φ1 el potencial electrostático en la región 0<r<R. Así mismo, sea Φ2 el
4πR 3
potencial electrostático en la región r>R. Dado que la distribución es isotrópica, las funciones Φ1 y

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

Φ2 satisfacen las ecuaciones 1 d  2 dΦ 1  1 3Zq , y 1 d  2 dΦ 2 

r =− r  = 0 . De la última
r 2 dr  dr  ε 0 4πR 3 r 2 dr  dr 
ecuación se deduce que r 2 dΦ 2 = A1 y el potencial fuera de la esfera es Φ 2 = − A2 + B2 ; dado que
dr r
para r→∞ el potencial se anula, entonces B2=0. La constante A2 se halla teniendo en cuenta que fuera
de la esfera considerada el campo eléctrico corresponde al campo creado por la carga Zq, es decir,
dΦ 2 A Zq 1 , de donde Zq . Así, el potencial fuera del núcleo es
E2 = − = − 21 = A1 = −
dr r 4πε 0 r 2 4πε 0
Zq 1 .
Φ2 =
4πε 0 r
Examínese ahora la solución de la ecuación de Poisson al interior del núcleo.
d  2 dΦ 1  1 3Zq 2 . Integrando, se obtiene  2 dΦ1  1 3Zq r 3
r =− r  r  = − + A1 , y por lo tanto
dr  dr  ε 0 4πR 3  dr  ε 0 4πR 3 3
 dΦ 1  1 3Zq r A1 .
 =− +
 dr  ε 0 4πR 3 3 r 2
Así. el campo eléctrico dentro del núcleo es E =  1 3Zq r − A1 e ; dado que este campo debe
1  ε 4πR 3 3 r 2  1
 0 
ser regular en el centro de la esfera, A1=0. Por lo tanto  dΦ1  = − 1 Zq r ; una nueva integración
 dr  ε 0 4πR 3
da Φ1 = − 1 Zq r 2 + B1 .
ε 0 8πR 3

Dado que el potencial es continuo en la superficie de la esfera, es decir, Φ1(R)= Φ2(R), se obtiene
1 Zq Zq 1 , de donde 3Zq 1 y así, Zq
− + B1 = B1 = Φ1 = (
3R 2 − r 2 )
ε 0 8πR 4πε 0 R 8πε 0 R 8πε 0 R 3

Ejemplo (1.4.4). Distribución de carga en un átomo. Se sabe que el potencial de un átomo de número
atómico Z tiene la forma Φ(r ) = Zq 1 1 + r  exp − r  , donde λ es una longitud característica
4πε 0 r  λ   λ
que caracteriza el apantallamiento de la carga electrónica. Se busca determinar la densidad de carga
que genera dicho potencial.
De la ecuación de Poisson se deduce que ρ (r ) = −ε 0 ∇ 2 Φ (r ) . Dado que el potencial posee simetría
esférica, la parte derecha de esta expresión no contiene la parte angular del Laplaciano, por lo que
1 d  dΦ  . Sin embargo, esta expresión no es válida para distribuciones singulares
ρ (r ) = −ε 0 2  r 2 
r dr  dr 
de carga. En efecto, si λ→∞, el potencial toma la forma Φ(r ) = Zq 1 , que corresponde a la carga
4πε 0 r
puntual Zq del núcleo atómico, situada en el origen de coordenadas; la correspondiente densidad de
carga es ρ(r)=(Zq)δ(r). Por esta razón, del potencial propuesto es preciso separar la parte singular
1/r:
Zq  1  r 1  r  Zq 1 Zq  1   r   1 r 
Φ(r ) =  exp  −  + exp −   ≡ +  exp −  − 1 + exp −  
4πε 0  r  λ  λ  λ  4πε 0 r 4πε  r
0    λ   λ  λ  
La acción del operador de Laplace sobre esta función es:

Zq 2  1  Zq 1 d  2 d  1   r   1  r 
∇ 2 Φ(r ) = ∇  + r  exp −  − 1 + exp  −  
 r  4πε 0 r dr  dr  r   λ   λ
2
4πε 0  λ 

61
CAPÍTULO 1. CAMPOS ELECTROSTÁTICOS EN EL VACÍO

Zq Zq  1  r  r 
=− δ (r ) −  rλ2 1 − λ  exp − λ  
ε0 4πε 0     
y por lo tanto, la densidad volumétrica de carga es

 1 1 r   r 
ρ (r ) = Zq δ (r ) − 2 
− 1 exp  −  = ρ n (r ) + ρ e (r ) .
 4π rλ  λ   λ 

Se tiene así una distribución de carga correspondiente a una carga puntual (en el centro) y una carga
distribuida con densidad ρe(r). Nótese que para r>λ esta densidad de carga es negativa, como es de
esperarse.

§ 1.5. MOMENTO DIPOLAR Y MOMENTO CUADRUPOLAR PARA UN SISTEMA

DE CARGAS DISCRETAS.

Examínese el comportamiento del potencial de un sistema de cargas en un punto situado a

grandes distancias de dicho sistema. Para simplificar la discusión, se analizará inicialmente
el ejemplo sencillo de un sistema de dos cargas (dipolo), el cual permitirá introducir los
principales conceptos.
Dipolo eléctrico. Considérese un sistema de cargas +q y -q separados una distancia
d tal como lo muestra la Figura (1.5.1) y sea P un punto de observación. Tómese un sistema
de coordenadas con origen en el centro del dipolo, con el eje Oz en la dirección de éste y
con el plano xOz conteniendo al punto P. Sean además r1 y r2 los radio-vectores que parten
desde los puntos donde están las cargas +q y -q, respectivamente, hasta el punto P y sea r el
radio vector que une al centro del dipolo con P. Nótese que la posición de la carga positiva
se caracteriza por el vector (d/2)ez, mientras que la de la carga negativa- por el vector
-(d/2)ez.

Figura (1.5.1)

Determínese el potencial en el punto P, bajo la condición que las distancias de las

cargas hasta dicho punto sean mucho mayores que la separación mutua entre éstas.
Analizando la geometría de la Figura (1.5.1), esta condición se puede expresar como
d/r<<1.

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

En una primera aproximación se puede considerar que las cargas que componen al
sistema están reunidas en un punto (esta aproximación se denomina monopolar). Sin
embargo, dado que la suma algebraica de dichas cargas para el caso considerado es igual a
cero, es preciso ir más allá de la aproximación monopolar y tener en cuenta la separación
espacial entre las cargas. La corrección de primer orden de la relación d/r recibe el nombre
de aproximación dipolar.
De acuerdo con el principio de superposición, el potencial en el punto P es igual a
la suma de los potenciales asociados a las cargas +q y -q:

q 1 q 1
Φ (r ) = Φ − (r1 ) + Φ + (r2 ) = − + . (1.5.1)
4πε o r1 4πε o r2

d d
En la Figura (1.5.1) se ve que − e z + r1 = r y e z + r2 = r , de donde se halla que
2 2

2
 d   d  d 1d 
r1 = (r1 • r1 ) =  r + e z  •  r + e z  = r 1 + (er • e z ) +   ,
 2   2  r 4 r 

siendo er=r/r un vector unitario en la dirección radial r. La correspondiente expresión para

la magnitud r2 se puede obtener de la expresión para r1 cambiando d por –d. Así, se tiene:

−1 / 2 −1 / 2
1 1  d 1  d   1 1  d 1  d  
2 2

= 1 + (er • e z ) +   , = 1 − (e r • e z ) +  
r1 r  r 4  r   r2 r  r 4  r  

Las anteriores expresiones se pueden simplificar si se tiene en cuenta que d/r es una
magnitud pequeña y por lo tanto en la expresión subradical se puede despreciar el término
2
1d 
el término infinitesimal de segundo orden   en comparación con el término de
4 r 
d
primer orden (er • e z ) ; así,
r
−1 / 2
1 1  d 1d  
2 −1 / 2
1 d 
= 1 + (e r • e z ) +    ≈ 1 + (e r • e z )
r1 r  r 4  r   r r 

En la última expresión se puede hacer una expansión en serie de Taylor, reteniendo en ella
los términos hasta primer orden en d/r; se obtiene

−1 / 2
1 1 d  1 1 d
= 1 + (e r • e z ) ≈ 1 − (er • e z ) . (1.5.2)
r1 r  r  r 2 r 

Análogamente,

63
CAPÍTULO 1. CAMPOS ELECTROSTÁTICOS EN EL VACÍO

−1 / 2
1 1  d 1  d  
2
1 1 d
= 1 − (er • e z ) +   ≈ 1 + (er • ez ) . (1.5.3)
r2 r  r 4  r   r 2 r 

Nótese que (1.5.3) se puede obtener cambiando d por -d en la expresión para 1/r1, lo cual
d d
se justifica si se observa que las expresiones − e z + r1 = r y e z + r2 = r difieren
2 2
únicamente en el signo de d.
Al reemplazar estas fórmulas en la expresión exacta (1.5.1) para el potencial de un
dipolo, se obtiene

Φ (r ) ≈
1 qd
(er • ez ) = 1 (er •2 p) , (1.5.4)
4πε o r 2
4πε o r

donde se ha definido la magnitud de carácter vectorial

p=qdez (1.5.5)

denominada momento dipolar del sistema. Esta magnitud física es una característica
intrínseca del sistema dipolar, ya que depende de la separación relativa entre las cargas d y
del valor q de las mismas. Designando mediante θ el ángulo que forman los vectores p y r,
el potencial del dipolo considerado toma la forma

p cosθ
Φ (r, θ ) = . (1.5.6)
4πε o r 2

Se ve que el potencial debido a un dipolo disminuye inversamente proporcional al

cuadrado de la distancia al sistema. Este resultado contrasta con el comportamiento del
potencial para una carga puntual, el cual es inversamente proporcional a r. Nótese además
que el potencial de un dipolo depende del ángulo polar θ, lo cual contrasta con el potencial
de una carga puntual, el cual es isotrópico (no depende de las variables angulares). La
anisotropía del potencial dipolar se debe a que el vector momento dipolar p introduce una
dirección privilegiada en el espacio que rompe la isotropía. El sistema, sin embargo, sigue
conservando la simetría de rotación alrededor del eje z, razón por la cual el potencial (1.5.6)
no depende de la variable azimutal ϕ.
Con base en la expresión (1.5.6), se pueden obtener las expresiones para las
componentes del correspondiente campo. Debido a la simetría del sistema, el cálculo es
más cómodo realizarlo en coordenadas esféricas. Recordando la expresión (0.2.11) para el
gradiente de un campo escalar en coordenadas esféricas, se obtiene

∂Φ 1 2p 1 ∂Φ 1 p
Er = − ≈ cosθ , Eθ = − ≈ senθ , Eϕ ≈ 0 .
∂r 4πε o r 3 r ∂θ 4πε o r 3

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

Se ve que la intensidad del campo eléctrico de un dipolo disminuye con el cubo de la

distancia.
Cuadrupolo eléctrico. En algunas situaciones, la aproximación dipolar puede
resultar insuficiente. Para ilustrar esta situación, considérese un sistema de tres cargas +q, -
2q, +q tal como lo muestra la Figura (1.5.2). La carga total de este sistema es cero, por lo
que a grandes distancias la aproximación coulombiana no da ningún aporte. El momento
dipolar de este sistema también es igual a cero, ya que se compone de la suma de dos
momentos dipolares opuestos.
Así, a grandes distancias del sistema de cargas considerado la aproximación dipolar
es evidentemente insuficiente. Por lo tanto, para este sistema de cargas, denominado
cuadrupolo lineal, cabe esperar que el potencial a grandes distancias varíe según 1/r3 . Para
comprobar esto, téngase en cuenta que de acuerdo con la geometría de la Figura (1.5.2),
dez+rb=r, -dez+ra=r. Por lo tanto,

−1 / 2
1 
2
1 1 d d 
= = 1 + 2 (e r • e z ) +    .
ra (r + de z ) • (r + de z ) r  r r 


Figura (1.5.2)

Efectúese una expansión en serie hasta los términos de segundo orden en d/r. Para ello,
2
1 1 3 d d 
téngase en cuenta que ≈ 1 − α + α 2 , donde α= 2 (e r • e z ) +   <<1 y
1+α 2 8 r r
2
2 d 
α ≈ 4  (e r • e z )2 . Se obtiene entonces
r

1 1  d 
2
1d 
ra r  r 2 r 
(
≈ 1 − (er • e z ) +   3(er • e z ) − 1  .
2

) (1.5.7)


Análogamente,

65
CAPÍTULO 1. CAMPOS ELECTROSTÁTICOS EN EL VACÍO

1 1  d 
2
1d 
rb r  r 2 r 
(
≈ 1 + (er • e z ) +   3(er • e z ) − 1  .
2

) (1.5.8)


(Esta última relación se puede obtener fácilmente de la expresión para 1/ra si se observa
que las expresiones dez+rb=r, -dez+ra=r difieren únicamente en un signo de d).
De acuerdo con el principio de superposición, el potencial en el punto P es igual a
la suma de los potenciales asociados a las cargas +q, -2q y +q . Entonces:

Φ (r ) =
q 1 2 1
 − +  ≈
(
 qd 2 3(er • e z )2 − 1
=
) ( 2
)
q 3d 2 (er • e z ) − d 2
. (1.5.9)
4πε o  ra r ra  4πε 0 r 3
4πε 0 r 3

Esta expresión corrobora la aseveración formulada anteriormente: a grandes distancias

(r>>d) el potencial de un cuadrupolo es inversamente proporcional al cubo de la distancia
al punto de observación P.
A la anterior expresión para el potencial cuadrupolar se le puede dar una
formulación más elegante, si se tiene en cuenta que

d2=≡ d2(er•er)≡d2eiei≡(d2δik)eiek.

d2(er•ez)2=(d•er)2≡(d•er)(d•er)≡(diei)(dkek)=didkeiek,

donde ei=xi/r es la i-ésima componente cartesiana del vector unitario er. Entonces

q(3d i d k − d 2δ ik ) ei ek 1 Qik xi x k
Φ (r ) = ≡ , (1.5.10)
4πε 0 r 3
4πε 0 2r 5

donde se ha definido la magnitud

(
Qik = 2q 3d i d k − d 2δ ik ) (1.5.11)

denominada tensor cuadrupolar. Esta magnitud es una característica intrínseca del sistema
de cargas, y para la aproximación cuadrupolar desempeña un papel análogo al que juega el
momento dipolar para la respectiva aproximación.

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

§ 1.6. EXPANSIÓN MULTIPOLAR DEL POTENCIAL ELECTROSTÁTICO.

Se plantea extender la discusión de la sección anterior al caso de una distribución arbitraria

de carga. Para simplificar el análisis, supóngase que las cargas están distribuidas de forma
volumétrica. Las extensiones a otras formas de distribución (superficial y lineal) es
inmediata. La geometría del problema está indicada en la Figura (1.6.1).

Figura (1.6.1)

La carga está localizada en un volumen V´ con una densidad volumétrica ρ(r´). Se

propone hallar el potencial en el punto P, el cual está a gran distancia del sistema. Tómese
un origen arbitrario en el punto O. El potencial en el punto P es entonces

1 ρ (r ')
Φ(r ) =
4πε 0 ∫ r − r ' dV´ .
V

Como el punto de observación está suficientemente alejado del sistema, entonces r´/r<<1
(siendo r´≠0), y por lo tanto

1  1  − 2r • r' + r' 2  3  − 2r • r' + r' 2  

2
1 1
= ≈ 1−  +   
r − r' r 2 − 2r • r' + r' 2 r 2 r2  8 r2  

1 r • r ' 3(r • r ') − r ' 2 r 2
2
≈ + 3 + .
r r 2r 5

Reemplazando este desarrollo en la expresión para el potencial del sistema de cargas, se

obtiene

1  1 r • r ' 3(r • r ')2 − r '2 r 2 

4πε 0 V∫´ r
Φ (r ) ≈  + 3 +  ρ (r ' )dV ' = Φ coul (r ) + Φ dip (r ) + Φ cuadr (r )
r 2r 5  ,


67
CAPÍTULO 1. CAMPOS ELECTROSTÁTICOS EN EL VACÍO

donde

1 1 1 q
Φ coul (r ) = ∫
4πε 0 r V ´
ρ (r ' )dV ' =
4πε 0 r
, (1.6.1)

1  r • r'  1 r 1 r• p
Φ dip (r ) = ∫ 
4πε 0 V ´ r  3 
ρ (r ' )dV ' = • ∫ r ' ρ (r ' )dV ´ ≡
4πε 0 r V ´
3
4πε 0 r 3
, (1.6.2)

1  3(r • r´) − r '2 r 2 

2
1
4πε 0 V∫´ 2 
Φ cuadr (r ) =   ρ (r ' )dV ' =
r5 

1 xi xk 1 1 xi xk Qik
=
4πε 0 r 2 V ´
5 ∫ ( )
3 x'i x 'k −r '2 δ ik ρ (r ' )dV ' ≡
4πε 0 2r 5
. (1.6.3)

En las anteriores expresiones se han introducido las definiciones que se compendian

en la siguiente tabla:

q = ∫ ρ (r ' )dV ´ p = ∫ r´ρ (r ' )dV ´ ( )

Qik = ∫ 3 x 'i x'k −r '2 δ ik ρ (r ' ) dV '
V V V

Carga total del Momento dipolar de un Tensor cuadrupolar de un sistema

sistema sistema de cargas de cargas distribuido (1.6.4)
distribuido continuamente
continuamente

Se ve que el término coulombiano en el desarrollo del potencial corresponde al potencial en

el punto P que crearía una carga q situada en el origen; el siguiente término, Φdip(r), es el
potencial que resultaría si un dipolo puntual igual al momento dipolar de la distribución de
cargas se localiza en el origen de coordenadas; análoga interpretación tiene el tercer
término correspondiente a la aproximación cuadrupolar. En la Figura (1.6.2) se
esquematizan estos resultados.

Figura (1.6.2)

Las expresiones (1.6.2), (1.6.3) junto con las definiciones (1.6.4) generalizan a una
distribución arbitraria de cargas las fórmulas (1.5.4), (1.5.5), (1.5.10) y (1.5.11).
Es preciso recalcar que el momento dipolar y el tensor cuadrupolar dependen del
origen de coordenadas O elegido y, por lo tanto, la expansión multipolar del potencial, en
general depende de la elección de dicho origen de coordenadas. Para ilustrar esta
aseveración, considérese una carga q situada en el origen de coordenadas; el potencial a
grandes distancias es netamente monopolar, ya que el momento dipolar y otros términos de

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

la expansión son nulos. Sin embargo, si el origen de coordenadas se elige en un punto que
no coincide con la posición de la carga puntual, el momento dipolar es diferente de cero; en
efecto, eligiendo el eje z a lo largo de la línea que une el origen de coordenadas y el punto
donde se encuentra la carga puntual, y considerando que la carga q está a una distancia l del
origen, se tiene px=py=0, pz=q∫z´δ(z´-l)dz´=ql; es decir, a grandes distancias el sistema tiene
aportes que van más allá de la aproximación monopolar.
No obstante, si la suma total de las cargas de un sistema cerrado es cero, el momento
dipolar no depende del origen de coordenadas. En efecto, si se desplaza el origen de
coordenadas O al punto O´ (Figura (1.6.3)) y si R es el radio-vector que une estos dos
puntos y r´´= r´+R, es evidente que la densidad volumétrica de carga no depende de esta
transformación (es decir, ρ(r´)=ρ(r´´)); por lo tanto el momento dipolar con respecto al
nuevo origen es

p' = ∫ r ' ' ρ (r ' ' )dV ´ = ∫ (r '+ R )ρ (r ' ' )dV ´ = ∫ r ' ρ (r ' )dV ´ + R ∫ ρ (r ' )dV ´ = p + qR .
V´ V´ V´ V´

Figura (1.6.3)

Se ve que si la carga total q es cero, entonces el momento dipolar no cambia al pasar de un

origen de coordenadas a otro.
El tensor cuadrupolar tiene algunas propiedades generales:
(i) El tensor cuadrupolar es simétrico con respecto a sus índices (Qik=Qki). Esta propiedad
es consecuencia inmediata de la definición.
(ii) La suma de los elementos diagonales del tensor cuadrupolar es igual a cero
(Qxx+Qyy+Qzz≡Qii=0). En efecto,

69
CAPÍTULO 1. CAMPOS ELECTROSTÁTICOS EN EL VACÍO

∑Q ii ( )
= ∑ ∫ 3x' i x' k −r' 2 δ ik ρ ( r' )dV ' =
V´

∫ ∑ [3x' x' −r ' δ ]ρ (r ' )dV ' ≡ ∫ [3r ' −r ' 3]ρ (r ' )dV ' = 0 ,
2 2 2
i i ii
V´ V´
donde se ha tenido en cuenta que la suma de los elementos diagonales del símbolo de
Krönecker es igual a 3.
Ejemplo 1.6.1. Dos cargas 2q y –q están separadas una distancia l, tal
como lo ilustra la Figura (1.6.1ej). Se trata de determinar los términos
monopolar, dipolar y cuadrupolar.
Solución. La densidad volumétrica de carga de este sistema es
ρ ( x´, y´, z´) = [2qδ ( z´−l ) − qδ ( z´)]δ ( x´)δ ( y´) . De aquí podemos
determinar los diversos términos:
• Q = ρdV ´ = [2 qδ ( z´−l ) − qδ ( z´)]dz´ δ ( x´)dx´ δ ( y´)dy´ =
∫
V
∫ ∫ ∫
2q-q=q. Por lo tanto, el término monopolar del potencial a grandes
distancias es Φ coul (r ) = 1 q Figura (1.6.1ej)
4πε 0 r
• p = x´ρdV ´ = [2 qδ ( z´−l ) − qδ ( z´)]dz´ x´δ ( x´)dx´ δ ( y´)dy´ = 0 , ya
x ∫ ∫ ∫ ∫ que ∫ x´δ (x´)dx´ = 0 .
V
Análogamente, py=0. Por otro lado,
p z = ∫ z´ρdV ´ = ∫ [2qz´δ ( z´−l ) − qz´δ ( z´)]dz´∫ δ (x´)dx´∫ δ ( y´)dy´ = 2ql .
V

El término dipolar del potencial toma la siguiente forma: Φ dip (r ) = 1 r • p = ql z .

4πε 0 r 3 2πε 0 r 3
• En cuanto a las componentes del tensor cuadrupolar, es fácil verificar que éstas son iguales a
cero. En efecto, las componentes no diagonales tales como Qxy y Qxz se anulan, ya que son
proporcionales a la integral ∫x´δ( x´)d x´, la cual es igual a cero; en efecto,
[ ][ ]{
Qxy = ∫ (3 x ' y '−r '2 δ xy )ρ ( r ' ) dV ' = 3 ∫ x 'δ (x´)dx´ ∫ y ' δ ( y´)dy´ ∫ [2qδ ( z´−l ) − qδ ( z´)]dz´ = 0 }
V
Por razones análogas se anulan las demás componentes no diagonales. Para las componentes
diagonales, es suficiente calcular Qxx, ya que por simetría ésta es igual a Qyy y por lo tanto
Qzz=-2Qxx. Así
Qxx = ∫ (3 x' 2 −r ' 2 δ xx )ρ ( r ' ) dV ' = ∫ (2 x' 2 − y '2 − z ' 2 )ρ ( r ' ) dV '
V V
Pero
∫ x' ρ (r ' )dV ' =
2
[∫ x' δ (x´)dx´][∫ δ ( y´)dy´]{∫ [2qδ (z´−l ) − qδ (z´)]dz´}= 0 y ∫ y'
2 2
ρ (r ' )dV ' = 0 .
V V
Resta calcular
∫ z' ρ (r ' )dV ' ; se tiene:
2

∫ z ' ρ (r ' )dV ' =

2
[∫ δ (x´)dx´][∫ δ (y´)dy´]{∫ [2qz' δ (z´−l ) − qz' δ (z´)]dz´}= 2ql
2 2 2 .
V
De esta forma, Qzz=Qyy=-Qzz/2=-2ql2. El término cuadrupolar del potencial
1 x Qxx + y Qyy + z Qzz
2 2 2
1 xi xk Qik ql 2 2 z 2 − x 2 − y 2
Φcuadup = = =
4πε 0 2r 5 4πε 0 2r 5 2πε 0 2r 5
Introduciendo coordenadas esféricas, con centro en el punto donde se encuentra la carga
negativa, el potencial a grandes distancias toma así la forma:
q  1 2l cos θ l 2 (3 cos 2 θ − 1) .
Φ(r ) = + +
4πε 0  r r2 r3



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

§ 1.7. ENERGÍA ELECTROSTÁTICA DE UN SISTEMA DE CARGAS.

Considérese un sistema estacionario de N cargas puntuales q1, q2,.... qN, las cuales
interactúan entre sí mediante la interacción coulombiana. Supóngase que inicialmente las
cargas están infinitamente alejadas una de la otra. La energía electrostática de este sistema
será por lo tanto nula. Se realiza ahora el siguiente proceso:
• Se trae la carga q1 desde el infinito hasta un punto 1 caracterizado por un radio-vector r1.
Como las demás cargas q2, q3,... qN permanecen infinitamente alejadas, esta operación
no cambia la configuración del sistema, así que el trabajo realizado será W1=0.
• Se trae ahora la carga q2 desde el infinito hasta un punto 2 caracterizado por un radio-
vector r2. En este caso sobre la carga q2 actúa el campo eléctrico E creado por la carga q1
en el punto donde se encuentra la carga q2. Si la distancia relativa entre las cargas q1 y q2
es r12, entonces el trabajo necesario para traer desde el infinito hasta el punto 2 la carga
q2 es realizado por la fuerza eléctrica F=q2E y es igual a

2 2 2 2
 1 q1 
W2 = ∫ F • dl =q2 ∫ E • dl = − q2 ∫ ( gradφ ) • dl = − q2 ∫ dφ = − q2φ 2 = − q2  .
∞ ∞ ∞ ∞  4πε 0 r12 

• Al traer la tercera carga q3 desde el infinito hasta un punto 3, sobre ella actuarán los
campos producidos por las cargas q1 y q2. El trabajo efectuado será entonces

3 3 3
W3 = q3 ∫ ( E1 + E 2 ) • dl = − q3 ∫ ( gradφ1 ) • dl − q3 ∫ ( gradφ2 ) • dl =
∞ ∞ ∞
3 3
 1 q1 1 q2 
= − q3 ∫ dφ1 − q3 ∫ dφ2 = − q3 (φ1 + φ2 ) 3 = − q3  + ,
∞ ∞  4πε 0 r13 4πε 0 r23 

donde r13= r1-r3, r23= r2-r3.

• Al traer la carga qj, sobre ella actuarán los campos de las (j-1) cargas q1, q2, ..., qj-1 y por
lo tanto el trabajo Wj efectuado será

 1 q1 1 q2 1 q ( j −1)  j −1
1 q j qi
W j = −q j  + + ..... +  = −∑ .
 4πε 0 r1 j 4πε 0 r2 j 4πε 0 r j ( j −1)  i =1 4πε 0 r ji

Así, el trabajo total realizado para conseguir esta configuración de N cargas

puntuales es igual a

N N j −1
1 q j qi
W = ∑ W j = −∑∑ .
j =1 j =1 i =1 4πε 0 r ji

Teniendo en cuenta que de acuerdo con el teorema trabajo-energía, el trabajo realizado es

igual a menos el cambio que experimenta la energía potencial del sistema, y como en el

71
CAPÍTULO 1. CAMPOS ELECTROSTÁTICOS EN EL VACÍO

infinito dicha energía es igual a cero, entonces la energía potencial de origen electrostático
de un sistema de cargas es

N j −1
1 q j qi
U = ∑∑ . (1.7.1)
j =1 i =1 4πε 0 r ji

Por otro lado, la energía del sistema de cargas se puede obtener al llevar cada una de
las cargas del sistema desde su posición inicial hasta el infinito. Los pasos para este proceso
son los siguientes:
• Se toma la carga q1 y se la lleva desde el punto 1 caracterizado por un radio-vector r1
hasta el infinito. En este caso sobre la carga considerada actúa el campo de las cargas q2,
..., qN y por lo tanto el trabajo efectuado será

 1 q2 1 qN 
W1´= q1  + ..... + .
 4πε 0 r21 4πε 0 rN 1 

• Al operar de manera análoga con la carga q2 y teniendo en cuenta que sobre ella la carga
q1 ya no actuará, el trabajo efectuado será

 1 q3 1 qN 
W2 ´= q 2  + ..... + .
 4πε 0 r32 4πε 0 rN 2 

• En general, para llevar la carga qj desde el punto j-ésimo con radio vector rj hasta el
infinito, se debe realizar un trabajo igual a

 1 q ( j +1) 1 qN  N
1 q j qi
W j ´= q j  + ..... +  ∑
= .
 4πε 0 r( j +1), j 4πε 0 rNj  i = j +14πε 0 rij

El trabajo total realizado para disociar este sistema de cargas es entonces

N N N
1 q j qi
W = ∑W j ´= −∑ ∑ 4πε .
j =1 j =1 i = j +1 0 rij

Recordando que el cambio de la energía potencial es -U (ya que la energía potencial en el

infinito es cero), se obtiene para la energía del sistema la siguiente expresión:

N N
1 q j qi
U =∑ ∑ 4πε . (1.7.2)
j =1 i = j +1 0 rij

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

Sumando las expresiones (1.7.1), (1.7.2) y teniendo en cuenta que rij=rji, se obtiene

1 N N 1 q j qi
U= ∑∑
2 j =1 i =1 4πε 0 r ji
, (1.7.3)

o bien,
1 N  N 1 q i  1 N
U= ∑qj ∑ = ∑ q jφ j
2 j =1  i =1 4πε 0 r ji  2 j =1
(1.7.4)

N
1 qi
donde φ j = ∑ es el potencial creado por el sistema de cargas en el j- ésimo punto
i =1 4πε 0 r ji

donde se encuentra la carga qj.

Esta expresión para la energía de un sistema de cargas se puede generalizar al caso
cuando la distribución es continua. Supóngase que el sistema se caracteriza por una
distribución volumétrica con densidad ρ. Se tiene entonces

1 N 1
U= ∑
2 j =1
δq jφ j = ∫ ρ (r )φ (r )dV .
2V
(1.7.5)

Esta expresión se puede transformar un poco más. Para ello se tiene en cuenta que la
densidad volumétrica de cargas se puede expresar a través de la divergencia del campo
eléctrico mediante la relación ρ=div(ε0E)=divD. Se obtiene entonces

ε0
U=
2 ∫ (divE )φ (r )dV .
V

Por otro lado, de la identidad (0.5.11), se obtiene div(φE)=(gradφ)•E+(divE) φ. Teniendo

en cuenta la relación entre el campo eléctrico y el potencial, se obtiene de aquí que
(divE)φ=div(φE)+E•E y por lo tanto, la expresión para la energía electrostática U queda
así:

ε0  
U=  ∫ (divφE )dV + ∫ E dV  .
2

2 V V 

La primera integral, de acuerdo con el teorema de Gauss (0.3.3), es igual a la integral de φE

sobre la superficie que limita el volumen V ; dado que la integración se extiende a todo el
espacio, y en el infinito el potencial se considera nulo, entonces esta integral desaparece. Se
obtiene así la siguiente expresión para la energía electrostática:

ε E2 
U = ∫  0 dV .
 (1.7.6)
V
2 

73
CAPÍTULO 1. CAMPOS ELECTROSTÁTICOS EN EL VACÍO

Esta fórmula expresa la energía electrostática a través del campo eléctrico y lleva a
la concepción de que dicha energía está distribuida en el espacio con una densidad
volumétrica

ε0E2
w= . (1.7.7)
2

Estos resultados se pueden aplicar al caso de una partícula elemental, por ejemplo, a
un electrón, Supóngase que el electrón es una esfera de radio R en cuya superficie está
distribuida la carga q. El campo al interior de la esfera es nulo y en cualquier punto fuera de
q 1
la ella es E = . Por lo tanto, la energía total eléctrica es, de acuerdo con (1.7.6),
4πε 0 r 2

∞
ε0 q2 4πr 2 q2
U= ∫ E dV = ∫R r 4 dr
2
= .
2 V
32πε 0 8πε 0 R

1
En otras palabras, la energía electrostática del electrón de radio R es U = qφ ,
2
donde φ es el potencial creado por el mismo electrón en cualquier punto de su superficie.
Para el caso de una carga puntual se tiene R=0 y por lo tanto U=∞, es decir, de acuerdo con
el formalismo desarrollado, una carga puntual debe poseer una energía electrostática propia
infinita. De esta forma, la concepción de una partícula como objeto puntual, que no tiene
una extensión espacial, lleva a un resultado carente de sentido físico. Esta situación indica
que los resultados obtenidos tienen un límite de aplicación que no se puede extender a
distancias pequeñas. El orden de estas distancias se puede determinar con base en las
siguientes consideraciones: si la partícula elemental (en este caso, el electrón) tiene un radio
1 q2
característico R0, entonces su energía electrostática propia es de la forma U ∝ ,
4πε 0 R0
donde el factor de proporcionalidad depende de la distribución de cargas dentro de la esfera
de radio R0. Pero esta energía debe corresponder a la energía en reposo de la partícula mc2,
1 q2
es decir, ~ mc 2 , de donde se encuentra el tamaño típico de la región dentro de la
4πε 0 R0
cual no son aplicables los resultados de la teoría clásica del campo electrostático:

1 q2
R0 = (1.7.8)
4πε 0 mc 2

La magnitud R0 definida en esta relación recibe el nombre de radio clásico. Para un

electrón se tiene q=1.6×10-19C, m=0.911×10-31 kg, y por lo tanto R0=2.8×10-15 m. Es
necesario aclarar que a R0 no se le puede adjudicar el sentido estricto de radio de una
partícula elemental, ya que la teoría clásica pierde su sentido a distancias

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

considerablemente mayores determinadas por la longitud de Compton λc=h/mc, la cual es

igual a 2.43×10-12 m para el caso del electrón.
Ejemplo 1.7.1. Anteriormente se consideró el caso de los núcleos ligeros, en los cuales la distribución
de carga tiene la forma ρ(r)=ρ0(1-r2/a2) si r≤a y ρ(r)=0 si r>a. Se estableció que si Q es la carga total
del núcleo, Q = 8 πa 3 ρ0 . Se comprobó que el campo eléctrico en el interior del núcleo (r≤a) es
5
1 r 2
2 ρ a 3 1 . Con base
E1 = ρ 0  − 2 re r , mientras que en el exterior (r>a) tiene la forma E 2 = 0 er
 3 5a  5ε 0 r 2
en estos resultados se puede determinar la energía acumulada por este sistema; teniendo en cuenta
que los campos son netamente radiales, el elemento de volumen en la expresión U = ε 0 E 2 dV
2 V∫
tiene la forma dV=4πr2dr y así

ε0 ε0 a ∞
 ε  2ρ 
2
a ∞
 a3  2 
2

U= ∫ E dV =
2
4π ∫ E12 r 2 dr + ∫ E 22 r 2 dr  = 0 4π  0  ( )
∫ r r dr + ∫  2  r dr  =
2 2
 
 5ε 0   0 ar  
2 V 2 0 a  2
48π 2 5 1 3 Q2
ρ0 a =
125ε 0 4πε 0 5 a
Pero
2
a a
  r2
a
r 4 
a
 r 4 r6 r 8   4 
∫0 E r
2 2
1 dr ≡ [
∫0 1
E r ] dr
2
= ∫0  0  3 5a 2 
 ρ  −   dr = ρ 0 ∫ 
2
 − 2 +  dr = ρ 02 a 4 
4  
0  9 15a 2
25a   525 

Compárese este resultado con la energía electrostática que tendría el sistema si su carga Q estuviera
distribuida superficialmente. En este caso E1=0, 1 Q y
E2 = er
4πε 0 r 2
∞
ε0 1 Q2 .1
U= 4π ∫ E22 r 2 dr =
2 a
4πε 0 2 a
¿Cómo cambia el resultado si la carga Q está distribuida uniformemente en el volumen de la esfera?

§ 1.8. SISTEMA DE CARGAS EN UN CAMPO EXTERNO.

Considérese un sistema de cargas que se encuentra en un campo eléctrico externo, tal como
lo ilustra la Figura (1.8.1). Sea Φ(r) el potencial correspondiente a este campo externo. Si el
sistema se caracteriza por tener una densidad de carga ρ(r), entonces la energía potencial
del sistema es

U = ∫ ρ (r´ )Φ (r´ )dV´ . (1.8.1)

V

Es necesario recalcar que esta energía no es la energía necesaria para formar el

sistema de cargas (en el sentido de la discusión del numeral anterior). La expresión (1.8.1)
corresponde a la energía de interacción entre el sistema de cargas con densidad ρ(r´) y el
campo externo con potencial Φ(r´).

75
CAPÍTULO 1. CAMPOS ELECTROSTÁTICOS EN EL VACÍO

Figura (1.8.1)

Se escoge un sistema de coordenadas con su origen en un punto O y sea r´ el radio

vector correspondiente al elemento de carga ρ(r´)dV´ (Figura (1.8.1)). Considérese un
punto O´ dentro del sistema de cargas y sea r el radio vector desde el origen O hasta O´. Se
supone que el campo externo y su potencial asociado cambian suavemente a lo largo del
sistema de cargas. Entonces se puede desarrollar la función Φ(r´) en serie de potencias de
r´-r:

∂Φ ∂ 2Φ
Φ (r´) = Φ (r ) + ( x´i - xi ) + ( x´i - xi )( x´k - xk ) +...
∂xi ∂xi ∂xk

Reemplazando en (1.8.1), se obtiene:

∂Φ 1 ∂2Φ
U = Φ (r ) ∫ ρ (r´ )dV´ +
∂xi V∫ ∫ ρ (r´ )( x´ - x )( x´
ρ (r´ )( x´ i - x i ) dV´ + i i k - x k )dV´ . (1.8.2)
V
2 ∂x i ∂x k V

La primera integral en esta expresión es la carga total q del sistema en el punto O´

con radio vector r. La segunda integral es el momento dipolar p del sistema en el punto O´.
La tercera integral se puede transformar si se tiene en cuenta que el momento cuadrupolar
del sistema en r es

(
Qik = ∫ 3( x`i − x i )( x`k − x k ) − r '−r δ ik ρ (r ' ) dV ' .
2
)
V´

y por lo tanto

1 δ ik
∫ ρ (r ' )( x` − x )( x` − x )dV ' = 3 Q ∫ r '− r
2
i i k k ik − ρ (r ' )dV ' .
V´
3 V´

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

Así, (1.8.2) toma la forma

1 ∂E i
U = qΦ (r ) − E • p − Qik . (1.8.3)
6 ∂x k

Aquí se ha tenido en cuenta que Ei=-∂Φ/∂xi y que δik∂2Φ/∂xi∂xk ≡∇2Φ=0.

Nótese que la carga total del sistema interactúa con el potencial externo, mientras
que el momento dipolar del sistema interactúa con el campo externo y el momento
cuadrupolar del sistema lo hace con el “gradiente” del campo externo (entendiendo como
tal la variación espacial que experimentan las componentes del campo eléctrico).
Conociendo la energía de interacción del sistema de cargas con el campo externo, se
puede determinar tanto la fuerza como el torque que actúan sobre este sistema.
Para el caso de la fuerza, se tiene

∂U ∂E i 1 ∂ 2 Ei
Fl = − = −qEi + pi + Qik . (1.8.4)
∂xl ∂x l 6 ∂x l ∂x k

Se ve en particular que un dieléctrico neutro (q=0) en un campo eléctrico externo uniforme

(∂Ei/∂xl=0) no experimenta la acción de ninguna fuerza de origen eléctrico sobre él.
Los resultados obtenidos se pueden ilustrar con algunos ejemplos concretos.
Considérese inicialmente un dipolo p en el potencial externo de una carga puntual
situada a una distancia r del centro del dipolo (Figura (1.8.2)).

Figura (1.8.2)

En este caso la energía de interacción es U=−E•p, donde E=(1/4πε0)r/r3 es el campo de la

carga puntual. Tomando el momento dipolar p a lo largo del eje z y teniendo en cuenta que
r•p=−rpcosθ (de acuerdo con la geometría de la Figura (1.6.3)), se obtiene
U=(qp/4πε0)cosθ/r2. De aquí se obtiene para las componentes de la fuerza en coordenadas
esféricas las siguientes expresiones:

Fr=−∂U/∂r=(qp/2πε0)cosθ/r3,
Fθ=−(1/r)∂U/∂θ=(qp/2πε0)cosθ/r3.

Se ve que el dipolo en el campo externo considerado está sometido a la acción de una

fuerza radial Fr y experimenta un torque τ=−∂U/∂θ que trata de alinear al dipolo con el
campo externo en el punto donde está situado el dipolo.

77
 Referencias capítulo 1

Las leyes del campo electrostático están expuestas con detalle en

Becker R. (1981), secciones 16-19,24, 31
Brédov, Rumiántsev y Toptigin (1985), secciones 21, 22
Corson y Lorrain (1972), secciones 2.1-2.6; 2.12, 2.14, 2.15
Griffiths D.J. (1999), capítulo 2 (exceptuando la sección 2.5)
Greiner W. (1998), (parte 1)
Jackson (1999), secciones 1.1-1.7, 1.10, 1.11
Landau y Lifshitz (1987), secciones 36,37, 40-42
Levich, Vdovin y Miamlin (1974), secciones 14-18
Panofsky and Phillips (1962), capítulo 1
Plonsey and Collin (1961), capítulo 2
Reitz , Milford y Christy (1996), capítulo 2, capítulo 6 secciones 6.1, 6.2, 6.3
Sepúlveda (1986), capítulo 1

Los complementos 4 de
Born and Wolf (1999)
Brédov, Rumiántsev y Toptigin (1985)
contienen material de utilidad relacionado con la función δ de Dirac y sus diversas representaciones.

Una presentación detallada de la expansión multipolar (al nivel del presente curso) se encuentra en
Heald and Marion (1995), secciones 2.1, 2.2, 2.3, 2.4, 2.5

En el libro de
Sepúlveda (2009), capítulo 5
se pueden seguir con detalle las expansiones multipolares de la energía, la fuerza, y el torque de un
sistema de cargas en un campo externo.

Problemas Capítulo 1
1. Encuentre las unidades de medida en el sistema SI de las siguientes magnitudes físicas: (a) permitividad
dieléctrica del vacío εo; (b) campo eléctrico E; (c) potencial eléctrico Φ; (d) cualquiera de las
componentes cartesianas del momento dipolar p; (e) cualquiera de las componentes cartesianas del tensor
cuadrupolar Qij; (f) densidad de energía electromagnética w.
2. Verificar que las demás funciones λ(x,ε) definidas en las Figuras (1.2.2-1.2.4) efectivamente son
representaciones de la función δ de Dirac.
3. Hállese la carga total de un sistema cuya densidad volumétrica de carga tiene la forma
 1 1 r   r 
ρ (r ) = q δ (r ) − 2 
− 1 exp  −  , donde 0≤r<∞ y λ es una constante.
 4π rλ  λ   λ 
4. ¿Cuál es la densidad superficial en el punto M´ con coordenadas (r´,θ´,ϕ´) de una carga q distribuida
uniformemente sobre un cascarón esférico de radio R?
5. ¿Cuál es la densidad lineal en el punto M´ con coordenadas (ρ´,ϕ´,z´) de una carga distribuida
uniformemente a lo largo de un anillo infinitamente delgado de radio R que descansa en el plano xy y
cuyo centro coincide con el origen de coordenadas?
6. Evaluar las siguientes integrales:
(i) (r 2 + r • a + a 2 )δ (r − a )dV , donde V se extiende a todo el espacio y a es un vector constante.
∫ V
(ii)
∫ r − a δ (5r )dV , donde V es un cubo de lado 2, centrado en el origen de coordenadas y a=4ey+3ez.
V
CAPÍTULO 1. CAMPOS ELECTROSTÁTICOS EN EL VACÍO

(iii)
∫ (r )
+ (r • a )r 2 + a 4 δ (r − a )dV , donde V corresponde al volumen de una esfera de radio R=7 con
4

V
centro en el origen de coordenadas y a=5ex+3ey+2ez. ¿Cuál es el resultado si el radio de la esfera es R=6?.
Usar (1.2.4) para evaluar la integral div r  exp(− r )dV , donde V es el volumen de una esfera de
7.
∫   r 3 
V 
radio R con centro en el origen de coordenadas
8. Se tienen dos campos vectoriales E1(x,y,z)=kxyzer y E2(x,y,z)=k[(2xz+y2)ex+2xyex+x2ez], donde k es una
constante. ¿Cuál de ellos corresponde a un campo eléctrico?
9. ¿Cuál es la densidad de carga que es fuente de los siguientes campos eléctricos?:
(a)E(r)=kr2er, donde k es una constante
(b)E(r)=[aer+asenθsenϕeθ]/r2, donde a y b son constantes.
10. Con base en las fórmulas
(1.3.2b) ó (1.3.2c) hállese el
campo eléctrico en un punto
M situado a una distancia z de
los siguientes sistemas (Fig):
(a) del centro de una espira
cuadrada de lado a sobre la
que está distribuida (a) (b) (c)
uniformemente una carga q;
(b) del centro de una espira circular de radio a sobre la que está distribuida uniformemente una carga q;
(c) del centro de un disco de radio a sobre la que está distribuida uniformemente una carga q;
11. Una placa conductora infinita delgada tiene una densidad superficial de carga σ. Demostrar que la mitad
de la intensidad del campo eléctrico E en un punto situado a z metros de la superficie de la placa se debe
a la carga situada en la placa dentro de un círculo de radio (3)1/2z.
12. Con la ayuda de la ley de Gauss, determínense los campos eléctricos dentro y fuera de las siguientes
configuraciones:
(a) Un cascarón esférico de radio interno a y radio externo b uniformemente cargado (con carga
total q).
(b) Un cilindro de radio R y altura H (H>>R) cargado con densidad volumétrica de carga
ρ(r)=ρ0(1-r/λ), donde r es la distancia al eje del cilindro; ρ0 y λ son constantes.
13. Hallar el potencial dentro y fuera de los sistemas indicados en el punto anterior.
14. Considérense dos líneas paralelas con cargas de polaridad opuesta de densidad lineal λ. Sea M un punto
situado a una distancia r1 de la línea negativa y r2 de la línea positiva. Demuestre que el potencial de este
sistema en P es (λ/2πεo)ln(r1/ r2).
15. Calcular el momento dipolar de una capa esférica con una densidad superficial de carga σ=σ0cosθ, donde
θ es el ángulo polar.
16. Una carga Q está distribuida uniformemente a lo largo de una línea recta de longitud a. (a) Determine el
potencial Φ eligiendo el origen de coordenadas en el centro de la carga. (b) Calcule los términos
monopolar, dipolar y cuadrupolar en el desarrollo en multipolos del potencial Φ. (c) ¿Para qué valores de
la distancia r al centro de la carga es el término cuadrupolar menor que el 1% del término monopolar?
17. Determine las componentes del tensor
cuadrupolar del sistema de cuatro cargas
discutido en el texto (Figura (1.5.2)).
18. Calcular las componentes del tensor cuadrupolar,
el potencial y el campo cuadrupolares de un
cuadrupolo lineal consistente en dos dipolos
idénticos con momentos dipolares p=2qδ
separados por una distancia 2l (Figura P.1.1).
19. Dos cargas puntuales 3q y –q situadas en el eje y
están separadas una distancia l. La carga positiva
se encuentra en el origen de coordenadas.
Determine el potencial de este sistema a grandes
distancias, conservando los términos hasta la
Figura (P.1.1) aproximación cuadrupolar.

40
CAPÍTULO 1. CAMPOS ELECTROSTÁTICOS EN EL VACÍO

20. Determinar las componentes del tensor cuadrupolar de un elipsoide uniformemente cargado (véase la
discusión en Landau L.D, Lifshits E.M., Teoría Clásica de Campos, problema al numeral 41)
21. Determine el potencial a grandes distancias del origen del
sistema de cargas ilustrado en la Figura. Exprese su respuesta
en coordenadas esféricas.
22. Considérese una esfera imaginaria de radio a centrada en un
dipolo de momento p. Integrando la densidad de energía desde
r=a hasta r=∞, determine la energía eléctrica asociada con la
región del espacio exterior a la esfera.
23. ¿Cuál sería el radio del electrón suponiendo que su carga está
distribuida en el interior de una esfera de radio a de acuerdo
con la distribución ρ(r)=ρ0(r/a)? ¿Cuál sería ese radio si
ρ(r)=ρ0(1-r2/a2), (r<a)?
24. Calcular la energía potencial eléctrica de una esfera de radio R con una carga total Q distribuida
uniformemente sobre su volumen. (b) Calcular la energía potencial gravitatoria de una esfera de radio R´
y masa total M y aplique este resultado al caso de la luna. (c) Supóngase que se pudieran reunir protones
hasta formar una esfera de densidad igual a la del agua ¿cuál sería el radio de esta esfera si su energía
potencial fuera suficiente para destruir la luna?.
25. Se tienen dos esferas concéntricas de radios a y b; en el centro de este sistema se tiene una carga puntual
Q, mientras que en la esfera interna se tiene una carga q y en la externa una carga -q distribuidas
uniformemente. ¿Cuál es la energía electrostática de este sistema?

41