Principios de Electrodinámica Clásica

Taller 8

Universidad Distrital Francisco José de Caldas

Principios de Electrodinámica Clásica

Esteban Palomino Gamboa∗

Abril del 2024

PAGS 146-155

• 3.4.1. Approximate potential at large distances

• 3.4.2 The Monopole and Dipole Terms

• 3.4.3 Origin of Coordinates in Multiple Expansions

• 3.4.4 The Electric Field of a Dipole

En este taller, para reducir el tamaño de las ecuaciones y sus resultados, se presentan algunas
ecuaciones que aparecerán frecuentemente en el texto:

• Fuerza sobre una carga puntual q:

F⃗ = q E
⃗ (1)

• Potencial eléctrico.

1 q
V (r) =
4πϵ0 r (2)

• Contribución del dipolo al potencial

1 p⃗ · r̂
Vdip (r) = (3)
4πϵ0 r2

• Campo eléctrico.

⃗ dip (r, θ) = 1 p  
E 2Cos(θ)r̂ + Sen(θ) θ̂ (4)
4πϵ0 r3

• Diferencial dτ .

dτ = r2 Sen(θ)drdθdϕ (5)
∗
20181135003 - epalominog@udistrital.edu.co

1
 Universidad Distrital Francisco José de Caldas

Ejemplos
3.10
Un dipolo eléctrico (físico) consta de dos cargas iguales y opuestas (+q) separadas
por una distancia d. Encuentre el potencial aproximado en puntos alejados del
dipolo.

Figura 1

Para encontrar el potencial aproximado en puntos lejanos como se muestra en la figura 1 se

utiliza el principio de superposición y la ecuación 2

 
q 1 1
V (r) =
4πεo r −r
+ −
(6)

donde las distancias r yr

+ − viene dada por la ley de cosenos como

 2  2 !
r2
+
2
=r +
d
2
d
− 2r Cos(θ) = r2
2
d
1 − Cos(θ) +
r
d
2r
(7)

 2    2 !
r 2
− = r2 + −
d
2
− 2r −
d
2
Cos(θ) = r2
d
1 + Cos(θ) +
r
d
2r
(8)

como se quiere estudiar el caso cuando r >> d, el termino (d/2r)2 → 0, por consiguiente

   1/2  −1/2
r 2
+
d
= r2 1 − Cos(θ)
r
→ r +
d
= r 1 − Cos(θ)
r
→
1
r
+
=
1
r
d
1 − Cos(θ)
r
(9)

   1/2  −1/2
r 2
−
d
= r2 1 + Cos(θ)
r
→ r −
d
= r 1 + Cos(θ)
r
→
1
r
−
=
1
r
d
1 + Cos(θ)
r
(10)
los dos primeros términos en la serie de Maclaurin para (1−x)−1/2 y (1+ x)−1/2 dan la siguiente
aproximación

2
 Principios de Electrodinámica Clásica

x
(1 − x)−1/2 ≈ 1 + (11)
2

x
(1 + x)−1/2 ≈ 1 − (12)
2

por consiguiente, si x = dr Cos(θ), se tiene que

 −1/2  
1 1 d 1 d
r
+
=
r
1 − Cos(θ)
r
≈
r
1 + Cos(θ)
2r
(13)

 −1/2  
1 1 d 1 d
r
−
=
r
1 + Cos(θ)
r
≈
r
1−
2r
Cos(θ) (14)

al reemplazar las ecuaciones anteriores en 6

    
q 1 d 1 d
V (r) ≈ 1 + Cos(θ) − 1 − Cos(θ)
4πεo r 2r r 2r
 
1 q d d 1 qd
≈ 1 + Cos(θ) − 1 + Cos(θ) ≈ Cos(θ) (15)
4πϵ0 r 2r 2r 4πϵ0 r2

El potencial de un dipolo cuando r es muy grande decrece mas rápido que para una sola carga, de
igual forma, a medida que se juntan dipolos, el decrecimiento es mayor, donde para el cuádruplo
es V ∼ r13 , para un óctuplo seria V ∼ r14 y así sucesivamente Para desarrollar una expansión
sistemática del potencial de una distribución de carga localizada arbitraria (Fig.2)en potencias
de 1r , el potencial en r esta dado por

Figura 2

Z
1 1
V (r) =
4πϵ0 r ρ(r )dτ
′ ′
(16)

donde r por ley de cosenos es

"  ′ 2  ′ #
r2 2 ′ 2 ′
= r + (r ) − 2rr Cos(θ ) = r ′ 2
1+
r
r
−2
r
r
′
Cos(θ ) (17)

si tomáramos el valor de ε como

 ′ 2  ′  ′ ′ 
r r ′ r r ′
ε≡ −2 Cos(θ ) ≡ − 2Cos(θ ) (18)
r r r r

3
 Universidad Distrital Francisco José de Caldas

entonces

r = r √1 + ε (19)

Para puntos muy fuera de la distribución de carga, ε es mucho menor que 1, y esto sugiere la
siguiente expansión binomial
 
1 1 1 1 3 5
r = (1 + ε)−1/2 =
r r
1 − ε + ε2 − ε3 + ........
2 8 16
(20)

remplazando ε para tener la expresión en términos de r, r′ , y θ

   ′     2
1 1 1 r′ r 3 r′ 2 r′
=
r r
1−
2 r r
′
− 2Cos(θ ) +
8 r r
′
− 2Cos(θ )

 ′ 3  ′ 3    ′  ′ 2  
5 r r ′ 1 r ′ r 3Cos2 (θ′ ) − 1
− − 2Cos(θ ) + ... = 1+ (Cos(θ )) +
16 r r r r r 2

 ′ 3   
r Cos3 (θ′ ) − 3Cos(θ′ )
+ + ... (21)
r 2

en la ecuación anterior, se puede observar como los términos entre paréntesis corresponden a
los polinomios de Legendre, por ende

∞  ′ n ∞
1 1X r X (r′ )n
r =
r
n=0
r
′
pn (Cos(θ )) =
rn+1
n=0
pn (Cos(θ′ )) (22)

donde θ′ es el ángulo entre r y r′ . Al sustituir la expresión anterior en 16 se obtiene que el

potencial para de una distribución de carga localizada arbitraria es:

∞ Z
1 X 1
V (r) = (r′ )n pn (Cos(θ′ ))ρ(r′ )dτ ′ (23)
4πϵ0 rn+1
n=0

Ejercicios
3.26
Una muestra de radio R, centrada en el origen, lleva densidad de carga.

R
ρ(r, θ) = k (R − 2r)Sen(θ) (24)
r2
donde k es una constante y r, θ son las coordenadas esféricas habituales. Encuentre
el potencial aproximado para puntos en el eje z, lejos de la esfera.
La expansión múltiple dada por la ecuación 23 al definirla explícitamente se tiene que
 Z Z
1 1 1
V (r) = ρ(r )dτ + 2 r′ Cos(θ)ρ(r′ )dτ ′
′ ′
4πϵ0 r r

4
 Principios de Electrodinámica Clásica

Z   
1 3 1
+ (r)´2 Cos2 (θ′ ) − ρ(r′ )dτ ′ + .... (25)
r3 2 2

en la ecuación anterior se ven las contribuciones del monopolo, dipolo, y cuadrupolo, de tal
forma que al evaluar estas por separadas se tiene que
• Contribución del monopolo:
Z Z
1 1 R
ρ(r, θ)dτ = k (R − 2r)Sen(θ)dτ (26)
r r r2

donde dτ viene dado por la ecuación 5, por consiguiente

Z
1 R
Vmon (r) = k (R − 2r)Sen(θ)r2 Sen(θ)drdθdϕ
r r2

Z R Z π Z 2π R
kR π 2 kR
= (R − 2r)dr Sen (θ) 2
dϕ = (Rr − r2 )
r 0 0 0 r 0

π 2 kR 2
= (R − R2 ) = 0 (27)
r

• Contribución del dipolo:

Z Z  
1 1 R
Vdip (r) = 2 rCos(θ)ρ(r, θ)dτ = 2 rCos(θ) k 2 (R − 2r)Sen(θ) dτ (28)
r r r

donde dτ viene dado por la ecuación 5, por consiguiente

Z  
1 R
Vdip (r) = rCos(θ) k 2 (R − 2r)Sen(θ) r2 Sen(θ)drdθdϕ
r2 r

Z R Z π Z 2π
kR
= 2 (Rr − 2r )dr 2
Cos(θ)Sen (θ)dθ 2
dϕ (29)
r 0 0 0

donde u = Sen(θ), por lo tanto du = Cos(θ)dθ, entonces

  RZ π
2πkR r2 2
Vdip (r) = R − r3 u2 du
r2 2 3 0 0

   π  
2πkR R3 2 3 Sen3 (θ) 5πkR4 Sen3 (π) Sen3 (0)
= − R = − =0 (30)
r2 2 3 3 0 3r2 3 3

• Contribución del cuadrupolo:

Z  
1 3 1
Vcua (r) = 3 r 2
Cos2 (θ) − ρ(r, θ)dτ
r 2 2

Z   
1 3 1 R
= 3 r 2
Cos2 (θ) − k 2 (R − 2r)Sen(θ) dτ (31)
r 2 2 r

5
 Universidad Distrital Francisco José de Caldas

donde dτ viene dado por la ecuación 5, por consiguiente

Z   
1 3 1 R
Vcua (r) = 3 r 2
Cos2 (θ) − k 2 (R − 2r)Sen(θ) r2 Sen(θ)drdθdϕ
r 2 2 r

Z R Z π   Z 2π
kR 3 1
= 3 (Rr − 2r )dr
2 3
Cos2 (θ) − Sen (θ)dθ2
dϕ
r 0 0 2 2 0

  RZ π
2πkR Rr3 r4 1

= −2 3Cos2 (θ) − 1 Sen2 (θ)dθ
r3 3 4 0 0 2

  Z π Z π 
2πkR R4 1
= − 3Cos2 (θ)Sen2 (θ)dθ − Sen2 (θ)dθ (32)
r3 6 2 0 0

1−Cos(4θ)
donde Cos2 (θ)Sen2 (θ) = 8 , por ende

Z   
πkR5 1 π
1 − Cos(4θ) π
Vcua (r) = − dθ −
3r3 2 0 8 2

πkR5  π  π 2 kR5
=− − = (33)
3r3 16 48r3

por consiguiente, el potencial aproximado a lo largo del eje z es

π 2 kR5
V (r) ∼
= (34)
48r3

3.27
Cuatro partículas (una de carga q, una de carga 3q y dos de carga −2q) se colocan
como se muestra en la Fig. 3, cada una a una distancia a del origen. Encuentre una
fórmula aproximada simple para el potencial, válida en puntos alejados del origen.
(Exprese su respuesta en coordenadas esféricas.)

Figura 3

La carga total Q = 0, por consiguiente el aporte del monopolo al potencial es

6
 Principios de Electrodinámica Clásica

1 Q
Vmon (r) = =0 (35)
4πϵ0 r

por el contrario, el aporte al potencial total del dipolo se da mediante la ecuación 3, donde
Vdip (r) =∼
= V (r). El momentum para la ecuación 3 se expresa de la siguiente forma (note que
al tener cuatro cargas, n = 4 para la sumatorio)

X
4
p⃗ = qi⃗ri = 3qaẑ + qa(−ẑ) + (−2qaŷ) + (−2qa(−ŷ))
i=1

= 3qaẑ − qaẑ − 2qaŷ + 2qaŷ = 2qaẑ (36)

por consiguiente

1 2qaẑ · r̂
V (r) ∼
= (37)
4πϵ0 r2

donde el ángulo entre los vectores unitarios ẑ y r̂ es ẑ · r̂ = Cos(θ), al reemplazar, se obtiene

que la aproximación al potencial válida en puntos alejados del origen es

1 2qaCos(θ)
V (r, θ) ∼
= (38)
4πϵ0 r2

3.28
En ej. 3.9 derivamos el potencial exacto para una capa esférica de radio R, que
tiene una carga superficial σ = kCos(θ)

Figura 4

(a) Calcule el momentum dipolar de esta distribución de carga.

Al tener una distribución de carga σ como se observa en la Fig. 4, su momentum esta dado por

Z Z
′ ′ ′
p= r σ(r )da = zσda (39)

7
 Universidad Distrital Francisco José de Caldas

donde da = R2 Sen(θ)dθdϕ, y z = RCos(θ) por consiguiente al reemplazar dichos términos (y

σ) en la ecuación anterior se tiene que

Z Z π Z 2π
2 3 2
p= (RCos(θ))(kCos(θ))R Sen(θ)dθdϕ = kR Cos (θ)Sen(θ)dθ dϕ (40)
0 0

donde u = Cos(θ) por ende du = Sen(θ)dθ. Reemplazando

Z π   π
u3 2
p = 2πkR 3 2
u du = 2πkR 3
= πkR3 (Cos(π) − Cos(0))
0 3 0 3

4 4
= πkR3 → p⃗ = πkR3 ẑ (41)
3 3
(b) Encuentre el potencial aproximado, en puntos alejados de la esfera, y compare
las respuestas exactas (3.87). ¿Qué puedes concluir sobre los multipolos superiores?
El contribución del dipolo al potencial se da por la ecuación 3, donde al utilizar el resultado de
la ecuación 41 se obtiene que


1 4πkR3 ẑ · r̂
Vdip (r) = (42)
4πϵ0 3r2

donde ẑ · r̂ = Cos(θ), por consiguiente

kR3 Cos(θ)
Vdip (r) = (43)
3ϵ0 r2

3.29
Para el dipolo del Ej. 3.10, expanda 1/r± para ordenar (d/r)3 y use esto para deter-
minar los términos cuádruple y óctuple en el potencial.
El potencial de un dipolo físico se define en el Ej 3.10 mediante la ecuación 6 (teniendo como
r r
referencia la Fig. 1), donde se expande + y − hasta el orden (d/r)3 . El valor obtenido para
r+ según la ecuación 7 es

 2 !1/2
r +
d
= r 1 − Cos(θ) +
r
d
2r
(44)

Para simplificar la expresión, se hace el reemplazo

  
d d
ϵ= − 2Cos(θ) (45)
2r 2r

por consiguiente se tiene que

1 1 1
r +
=
(1 + ϵ) 1/2
= (1 + ϵ)−1/2
r

"          #
1 1 1 3 (1) −5/2 1 3 5 (1) −7/2
= (1)−1/2 + − (1)−3/2 ϵ + − − ϵ2 + − − − ϵ3
r 2 2 2 2! 2 2 2 3!

8
 Principios de Electrodinámica Clásica

       
1 1 3 2 5
= 1− ϵ+ ϵ − ϵ3 (46)
r 2 8 16

Reescribiendo la ecuación anterior en términos de r, θ hasta el segundo orden en el que hay una
contribución del cuadrupolo obtenemos
    2      4    3
1 1 1 d d 3 d 3 d
r
+
=
r
1−
2 2r
+
2
Cos(θ) +
8 2r
−
8 r
Cos(θ)

   2  "    2 #
3 d 1 d d
+ Cos2 (θ) = 1− Cos(θ) + (3Cos2 (θ) − 1) + ... (47)
2 r r 2 2r

Repitiendo el mismo procedimiento para r , el cual se obtiene mediante la ecuación 8 dando

−
como resultado
 2 !1/2
r −
d
= r 1 + Cos(θ) +
r
d
2r
(48)

donde   
d d
ϵ= + 2Cos(θ) (49)
2r 2r

Siguiendo el mismo proceso que 1/ r , se obtiene que 1/r

+ − expresado en términos de r, θ es.

    2      4     3
1 1 1 d d 3 d 3 d
= 1− − Cos(θ) + + Cos(θ)
η− r 2 2r 2 8 2r 8 r

   2  "    2 #
3 d 1 d d
+ 2
Cos (θ) = 1+ Cos(θ) + (3Cos (θ) − 1) + ...
2
(50)
2 r r 2 2r

por consiguiente, el potencial del cuadrupolo (considerando los términos de segundo grado) es:

 2
1 q d

V (r, θ) = 3Cos2 (θ) − 1 − (3Cos2 (θ) − 1) = 0 (51)
4πϵ0 r 2r

como r = r , el potencial de un cuadrupolo es (V = 0).

+ −

El desarrollo anterior se sintetiza en la primera parte de la ecuación 22, donde pn es el momento

dipolar y n indica el orden de la aproximación. Para determinar el potencial del octópolo n = 3;
r r
así para + y − se tiene que

∞  ′ 3
1 1X r
r +
=
r
n=0
r
p3 (Cos(θ′ )) (52)

∞  ′ 3
1 1X r
r
−
=
r
n=0
r
p3 (Cos(π − θ′ )) (53)

En el desarrollo del cuadrupolo el momento dipolar es una función impar, por lo que pn (−x) =
−pn (x) donde x = Cos(θ) y su signo en la sumatoria alterna como (−1)n . Así, − seria r
9
 Universidad Distrital Francisco José de Caldas

∞  ′ 3
1 1X r
r−
=
r
n=0
r
(−1)3 p3 (Cos(θ′ )) (54)

donde r′ = d/2, por consiguiente la diferencia de las distancias r y r− está dada como
+

"∞   ∞  ′ n
#
1 1 1 X r′ n X r
r −r
+ −
=
r
n=0
r
pn (Cos(θ′ )) −
r
(−1)n p3 (Cos(θ′ ))
n=0
"   ′ 2  ′ 3
1 r′ ′ r ′ r
= p1 (Cos(θ )) + p2 (Cos(θ )) + p3 (Cos(θ′ ))
r r r r

 ′  ′ 2  ′ 3 #
r ′ r ′ r ′
+ p1 (Cos(θ )) − p2 (Coss(θ )) + p3 (Cos(θ )) (55)
r r r

El primer termino (n = 1) corresponde al dipolo físico, (n = 2) corresponde al cuadrupolo (el

cual se anula), y (n = 3) corresponde a el término del octopolo y el que se usará para determinar
el potencial.

 ′ 3 !
1 2q r
Voct (r, θ) = p3 (Cos(θ′ )) (56)
4πϵ0 r r

De la expansión binomial al factorizan los términos de tercer orden y se obtiene:

1 qd3
Voct (r, θ) = (5cos3 (θ) − 3 cos(θ)) (57)
4πϵ0 8r4

3.30
Dos cargos de puntos, 3q y −q, están separados por una distancia a. Para cada uno
de los arreglos en Fig.5, encuentre:

Figura 5

(i) El momento monopolar.

El momento de un monopolo se puede definir como
Z
pmon = ρ(r)dτ (58)

10
 Principios de Electrodinámica Clásica

donde la densidad de carga ρ(r) es cero en todos los puntos excepto en el origen, donde es
infinita (una delta de Dirac). Entonces, la integral se reduce a:

pmon = Q (59)

donde la carga total Q en esta ocasión seria la suma de las cargas, por consiguiente el potencial
del monopolo de cada configuración de la figura 5 son:
• Figura 5a

pmon = 3q − q = 2q (60)

• Figura 5b

pmon = 3q − q = 2q (61)

• Figura 5c

pmon = 3q − q = 2q (62)

(ii) El momento dipolar.

El momento dipolar se da como

X
n
p⃗ = qi⃗ri (63)
i=1

por consiguiente el potencial del dipolo de cada configuración es

• Figura 5a

p⃗ = −q(0) + 3qaẑ = 3qaẑ (64)

• Figura 5b

p⃗ = −qa(−ẑ) + 3q(0) = qaẑ (65)

• Figura 5c

p⃗ = −q(0) + 3qaŷ = 3qaŷ (66)

(iii) El potencial aproximado (en coordenadas esféricas) en r grande (incluye las

contribuciones monopolo y dipolo).
El potencial aproximado para una distribución de cargas discreta es la suma del potencial
monopolar Vmon y el potencial dipolar Vdip , de tal forma que
 
1 Q p⃗ · r̂
V (r) ∼
= Vmon + Vdip ∼
= + 2 (67)
4πϵ0 r r

por consiguiente el potencial aproximado de cada configuración es

11
 Universidad Distrital Francisco José de Caldas

• Figura 5a

 
1 2q 3qaẑ · r̂
V (r) ∼
= + (68)
4πϵ0 r r2

Recordando que ẑ · r̂ = Cos(θ), entonces

 
1 2q 3qaCos(θ)
V (r, θ) ∼
= + (69)
4πϵ0 r r2

• Figura 5b

 
1 2q qaẑ · r̂
V (r) ∼
= + (70)
4πϵ0 r r2

Recordando que ẑ · r̂ = Cos(θ), entonces

 
∼ 1
V (r, θ) =
2q qaCos(θ)
+ (71)
4πϵ0 r r2

• Figura 5c

 
1 2q 3qaŷ · r̂
V (r) ∼
= + (72)
4πϵ0 r r2

Recordando que ŷ · r̂ = Cos(θ), entonces

 
1 2q 3qaCos(θ)
V (r, θ) ∼
= + (73)
4πϵ0 r r2

3.31
Un dipolo ”puro” p está situado en el origen, apuntando en la dirección z.
(a) ¿Cuál es la fuerza sobre una carga puntual q en (a, 0, 0)(coordenadas cartesia-
nas)?

Figura 6

12
 Principios de Electrodinámica Clásica

Según la Fig.6, r = a, θ = π/2, ϕ = 0, r̂ = x̂, y θ̂ = −ẑ. Utilizando la ecuación 4 se obtiene que

el campo eléctrico del dipolo es

1 p  π  π   1 p
⃗ dip (r, θ) =
E 3
2Cos x̂ + Sen (−ẑ) =− ẑ (74)
4πϵ0 a 2 2 4πϵ0 a3

por consiguiente, si describimos la fuerza como en la ecuación 1 teniendo en cuenta el campo

eléctrico dipolar encontrado para (a, 0, 0), se obtiene que las fuerza sobre la carga puntual q es

1 qp
F⃗ = − ẑ (75)
4πϵ0 a3

(b) ¿Cuál es la fuerza sobre q en (0, 0, a)?

Figura 7

Según la Fig.7, r = a, θ = 0, ϕ = 0, r̂ = ẑ, y θ̂ = x̂. Utilizando la ecuación 4 se obtiene que el

campo eléctrico del dipolo es

⃗ dip (r, θ) = 1 p 1 2p
E (2Cos(0)ẑ + Sen(0)x̂) = ẑ (76)
4πϵ0 a3 4πϵ0 a3

por consiguiente, si describimos la fuerza como en la ecuación 1 teniendo en cuenta el campo

eléctrico dipolar encontrado para (0, 0, a), se obtiene que las fuerza sobre la carga puntual q es

1 2qp
F⃗ = ẑ (77)
4πϵ0 a3

(c) ¿Cuánto trabajo se necesita para mover q de (a, 0, 0) a (0, 0, a)?

13
 Universidad Distrital Francisco José de Caldas

Figura 8

El trabajo necesario para mover q de (a, 0, 0) a (0, 0, a) es

W = q V(0,0,a) − V(a,0,0) (78)

donde el potencial de un dipolo en función de r, θ es

1 pCos(θ)
Vdip (r, θ) = (79)
4πϵ0 r2

donde θa,0,0 = π/2, θ0,0,a = 0, y r = a se obtiene que el trabajo necesario para mover la carga q
es:

1 qp 1 qp
Vdip (r, θ) = 2
(Cos(0) − Cos(π/2)) = (80)
4πϵ0 a 4πϵ0 a2

3.32
Tres cargas puntuales están ubicadas como se muestra en la Fig. 9, cada una a una
distancia a del origen. Encuentre el campo eléctrico aproximado en puntos alejados
del origen. Exprese su respuesta en coordenadas esféricas e incluya los dos órdenes
más bajos en la expansión multipolar.

Figura 9

14
 Principios de Electrodinámica Clásica

Se debe tener en cuenta el campo del monopolo donde la carga total es Q = q − q − q = −q,
por consiguiente:

1 Q 1 q
⃗ mon (r) =
E 2
r̂ = − r̂ (81)
4πϵ0 r 4πϵ0 r2

El campo eléctrico dipolar se puede describir como

1 1
⃗ dip =
E p · r̂)r̂ − p⃗]
[3(⃗ (82)
4πϵ0 r3

donde el momentum para n = 3 (tres cargas) es

X
3
p⃗ = qi r̂i = qaẑ + (−qaŷ)(−qa(−ŷ)) = qaẑ − qaŷ + qaŷ = qaẑ (83)
i=1

por ende

1 qa
⃗ dip (r) =
E [3(qaẑ · r̂)r̂ − qaẑ] (84)
4πϵ0 r3

el ángulo entre ẑ y ẑ es la magnitud del producto de los dos unitarios el cual es 1, y el ángulo
entre estos dos es θ, por consiguiente

1 qa   1 qa  
⃗ dip (r, θ) =
E 3Cos(θ)r̂ − Cos(θ)r̂ + Sen(θ) θ̂ = 2Cos(θ)r̂ + Sen(θ) θ̂ (85)
4πϵ0 r3 4πϵ0 r3

⃗ mon y
por consiguiente, el campo eléctrico aproximado en puntos alejados del origen al sumar E
⃗
Edip

 
⃗ θ) = − 1 q r̂ + 1 qa 2Cos(θ)r̂ + Sen(θ)θ̂
E(r,
4πϵ0 r2 4πϵ0 r3
 
1 q 2aCos(θ) aSen(θ)
= −r̂ + r̂ + θ̂ (86)
4πϵ0 r2 r r

3.33
Demuestre que el campo eléctrico de un dipolo (”puro”) (ecuación 3.103) se puede
escribir en la forma libre de coordenadas

1 1
⃗ dip (⃗r) =
E p · r̂)r̂ − p⃗)]
[3((⃗ (87)
4πϵ0 r3

El potencial de un dipolo en la forma libre se expresa como

1 p⃗ · r̂
Vdip (r) = (88)
4πϵ0 r2

Como el campo eléctrico se obtiene a partir del gradiente del potencial, se tiene que

15
 Universidad Distrital Francisco José de Caldas

 
1 p⃗ · r̂
Edip (r) = −∇V = − ∇ (89)
4πϵ0 r2

por propiedades se tiene que

  
1 1 1
Edip (r) = − 2
∇(⃗
p · r̂) + (⃗
p · r̂)∇ 2 (90)
4πϵ0 r r

Solucionando cada uno de los gradientes por separado, se tiene que aplicando la identidad

∇(p · r̂) = p × (∇ × r̂) + r̂ × (∇ × p) + (p · ∇)r̂ + (r̂ · ∇)p (91)

teniendo en cuenta que ∇ × r̂ = 0 debido a que la derivada de una constante es cero, ∇ × p = 0
porque sus derivadas en cualquier dirección son cero ya que el momentum dipolar es fijo, y
(r̂ · ∇)p = 0, por lo tanto el termino p
(p · ∇)r̂ al ser diferente de cero, y al solucionarlo teniendo
en cuenta que r̂ = ⃗r/|r| donde |r| = r2 + y 2 + z 2 , por consiguiente el termino expandido es

   
∂ x ∂ x ∂ x ∂ y  ∂ y  ∂ y 
(p · ∇)r̂ = px + py + pz x̂ + px + py + pz ŷ
∂x r ∂y r ∂z r ∂x r ∂y r ∂z r

 
∂ z  ∂ z  ∂ z 
+ px + py + pz ẑ (92)
∂x r ∂y r ∂z r

La derivada de |r| es r3 , de tal forma que dr = (r2 + y 2 + z 2 )3/2 , por ende se puede reescribir
la ecuación anterior como

    xy   xz    xy     yz 
1 x2 1 y2
(p · ∇)r̂ = px − 3 − py − pz x̂ − px + py − 3 − pz ŷ
r r r3 r3 r3 r r r3

  xz   yz   
1 z2
− px − py + pz − ŷ (93)
r3 r3 r r3

Al momento de agrupar los términos por el momentum en px , py y pz se tiene que

px xpx py ypy pz
∇(p · r̂) = (p · ∇)r̂ = x̂ − 3 (xx̂ + y ŷ + z ẑ) + ŷ − 3 (xx̂ + y ŷ + z ẑ) + ẑ
r r r r r
zpz 1 1 1
− 3 (xx̂ + y ŷ + z ẑ) = [px x̂ + py ŷ + pz ẑ] − 3 [px x + py y + pz z]⃗r = [⃗
p − (r̂ · p⃗)r̂] (94)
r r r r

En la ecuación 90 falta por resolver el siguiente gradiente el cual solo depende del vector, por
lo tanto

       
1 ∂ 1 ∂ 1 ∂ 1
∇ = x̂ + ŷ + ẑ
r2 ∂x x + y2 + z2
2 ∂y x + y2 + z2
2 ∂z x + y2 + z2
2

   
2x 2y 2z 2 xx̂ + y ŷ + z ẑ 2 1 2
= − 4 x̂ 4 ŷ − 4 ẑ = − 3 = − 3 (⃗rr) −
→∇ =− r̂ (95)
r r r r r r r2 r3

16
 Principios de Electrodinámica Clásica

Reemplazando los resultados obtenidos en las ecuaciones 94 y 95, se encuentra la expresión para
la energía de un dipolo ”puro”libre de coordenadas
     
1 1 1 2
Edip (r) = − p − (r̂ · p⃗)r̂] + (⃗
[⃗ p · r̂) − 3 r̂
4πϵ0 r2 r r
   
1 1 1 2 1 1
=− 3
p⃗ − 3 (r̂ · p⃗)r̂ + − 3 (p · r̂)r̂ = − p − (r̂ · p⃗)r̂ − 2(p · r̂)r̂]
[⃗
4πϵ0 r r r 4πϵ0 r3

1 1
= p · r̂)r̂ − p⃗)]
[3((⃗ (96)
4πϵ0 r3

17