R |r − yS | [r + R − 2rR cos(r̂ · ŷS )]

ℓ=0

so

"∞
∂I ∂I rℓ R2 − r 2
−Rr = (2ℓ + 1)P ℓ (r̂ · ŷ S ) = .
Rℓ+1
∂r Electromagnetismo
∂R II |r − yS |3
ℓ=0
Semestre: 2018-1
Therefore, TAREA 3: SOLUCIÓN
Dr. Alejandro Reyes Coronado
2 ! ′
R(R2 − rParra
Ayud. Adolfo Manuel ) ′ ϕ̄(Ω )
Leyva
ϕ(r, Ω) = dΩ
Ayud. Jonathan Alexis4πUrrutia Anguiano |r − yS |3
′
as required because dyS = R2 dΩFecha
. de entrega: jueves, 30 de septiembre de 2018

1. Problema: (15pts) El eje êz corre a lo largo del centro de una lı́nea infinita de un ducto de ventilación con
sección transversal cuadrada. Para un conductor de metal real (no un conductor perfecto), el potencial
7.4 Theφ(x,
electrostático Potential
y) varı́a deinside an Ohmic
forma lineal a lo largoDuct
de las paredes del ducto. Considera que las esquinas
del ducto en (±a, 0) se mantienen a un potencial +V y los vértices del ducto en (0, ±a) se mantienen a
The geometry is
un potencial −V , como se muestra en la figura.

−V

y
V x V

−V
Calcula el potencial dentro del ducto comenzando con la solución de prueba:
The corner values imply that the potential is invariant when x → −x or when y → −y.
This tells us that B = φ(x,
C =y)F= =
A+ 0,Bx
so+ϕ(x, Dx=2 +
Cy +y) AEy
2
+ Dx+ F2 xy.
+ Ey 2 . Moreover, the potential
inside the duct satisfies Laplace’s equation. Therefore,
Solución

2 ∂2 ϕ ∂2 ϕ
Dadas las condiciones de contorno del ∇problema,
ϕ = el2potencial
+ es=invariante
2D + 2Eante= reflexiones
0. de las coordenadas
x → −x o cuando y → −y. Entonces, de la solución ∂y 2 se tiene que las constantes B = C = F = 0
∂x de prueba
(porque no cumplirı́a el potencial con la simetrı́a antes descrita). Además, el potencial cumple la ecuanción de
Laplace,This
dadogives
que noDhay= cargas
−E and dentro
ourdel ducto,
trial es decirbecomes
solution
∂2φ ∂2φ
∇2 φ = + 2 = 2D + 2E 2= 0. 2
ϕ(x,
∂x2 ∂y = A + D(x − y ).
y)
Entonces E = −D, por lo que se tiene que
We get V = A + Da2 from ϕ(a, 0) = V . We get −V = A − Da2 from ϕ(0, a) = −V .
2 2
Combining these two gives Aφ(x, = 0y) and
=A+ = V−/ay2 .). Therefore, the potential inside the duct
DD(x
is
Al incluir las condiciones de frontera φ(a, 0) = A + Da2 = V y φ(0, a) = A − Da2 = −V , se obtiene que A = 0
y D = V /a2 , y por lo tanto la solución es:
101
V 2
φ(x, y) = (x − y 2 ).
a2

FOR ENDORSEMENT PURPOSES ONLY. DO NOT DISTRIBUTE

1
2. Problema: (25pts) Un conjunto de constantes αn parametriza el potencial en un volumen r < a como
∞
X  r n
φext (r, θ) = αn Pn (cos θ).
n=1
R

Sea êz apuntando en la dirección θ = 0 y considera una esfera conductora sólida de radio R < a en el
origen. Muestra que la fuerza ejercida sobre la esfera cuando se conecta a tierra está en la dirección z con
magnitud
X∞
Fz = 4π0 (n + 1)αn αn+1 .
n=1

Hint: Los polinomios de Legendre satisfacen que (n + 1)Pn+1 (x) + nPn−1 (x) = (2n + 1)xPn (x).

Solución

Al existir un potencial para la región r < a entonces existe un campo eléctrico. La esfera conductora de
radio R < a, al estar aterrizada adquiere carga eléctrica en su superficie, por lo que la esfera de radio R va a
sentir una fuerza. El potencial φσ debido a la carga inducida en la esfera cancela al potencial externo φext en
su interior. Como no hay cargas dentro de la esfera, φσ satisface la ecuación de Laplace y es finito dentro de la
esfera, asimismo φs igma y es continuo en r = R, es decir:
X∞ X∞
Bn
P n (cos θ) = φ σ (R, θ) = −φ ext (R, θ) = − αn Pn (cos θ),
n=1
Rn+1 n=1

por lo que Bn = −Rn+1 αn . Entonces el φθ (r, θ) está dado por la expresión

 in P∞

r n
 φσ = − n=1 αn R Pn (cos θ) (r ≤ R)
φσ (r, θ) = .
 out P∞
n+1
φσ = − n=1 αn Rr Pn (cos θ) (R ≤ r)

La densidad de carga inducida está dada por la discontinuidad de la componente normal del campo eléctrico
en la superficie de la esfera
  
out 
  
∂φin  
E ~ in · êr = − ∂φσ 
~ out − E + σ 
= −
∂
φ ext
+ φ out  σ
= ,

∂r r=R+ 
∂r r=R+ ∂r σ σ  ε0
r=R
por lo que despejando σ se obtiene
∞
X αn
σ(θ) = −ε0 (2n + 1) Pn (cos θ).
n=1
R

Por otro lado, la fuerza sobre dicha densidad de carga inducida se obtiene integrando la misma sobre la
superficie de la esfera, es decir Z Z
~ ~ σ2
F = EσdS = êr dS.
S S 2ε0
La independencia en ϕ de la carga inducida y la expresión del vectori normal êr = sin θ cos φ êx + sin θ sin φ êy +
cos θ êz , hace que al integrar en la variable ϕ únicamente la componente en êz contribuya a la fuerza. Como
cos θ = P1 (cos θ), haciendo el cambio de variable x = cos θ, la fuerza puede ser escrita como
∞ X
X ∞ Z 1
Fz = πε0 (2n + 1)(2k + 1)αn αk xPn (x)Pk (x)dx.
n=1 k=1 −1

2
 La integral en Fz puede ser resuelta con la relación entre el polinomio de Legrendre de grado n con los de
grado n + 1 y n − 1, es decir
Z 1 Z 1
xPn (x)Pk (x)dx = Pk (x)[(n + 1)Pn+1 (x) + nPn−1 (x)]dx.
−1 −1

Empleando Además la ortogonalidad de los polinomios de Legendre

Z 1
2δn,k
Pn (x)Pk (x)dx = ,
−1 2n +1
se tiene que
 
X n + 1 2δn+1,k X n + 1 2δn−1,k 
Fz =πε0  (2n + 1)(2k + 1)αn αk + (2n + 1)(2k + 1)αn αk
2n + 1 2k + 1 2n + 1 2k + 1
n,k=1 n,k=1
" #
X X
=πε0 2αn αk (n + 1)δn+1,k + 2αn αk nδn−1,k
n=1 n=1
" #
X X
=2πε0 αn αn+1 (n + 1) + αn αn−1 n .
n=1 n=1

Por hipoteis se tiene que α0 = 0, entonces el primer término no nulo de la seugunda suma es en n = 2, por lo
que podemos redefinirla renombrando el ı́ndice de la suma como n → n + 1, esto hace la segunda suma idéntica
a la primera, y entonces la fuerza es
∞
X
Fz = 4πε0 (n + 1)αn αn+1 .
n=1

3. Problema: (30pts) Dos cascarones esféricos concéntricos de radios a y b (b > a) son divididos en dos
hemisferios por el mismo plano horizontal, el plano XY . El hemisferio de arriba de la esfera interior y
el hemisferio de abajo de la esfera exterior se mantienen a un potencial φ0 . Los otros dos hemisferios se
mantienen aterrizados (a potencial cero).
(i) Calcula el potencial en la región a ≤ r ≤ b como una serie en polinomios de Legendre.
(ii) Incluye términos hasta al menos l = 4 y coteja tu solución con la de los casos lı́mite b → ∞ y a → 0.
(iii) Considera el plano X = 0 (plano Y Z) y grafica la solución que obtuviste en (i) incluyendo suficientes
valores de l en la suma como para que la solución no cambie en las tres primeras cifras significativas.
La gráfica la puedes hacer con código de color. Para el cálculo numérico considera φ0 = 1, a = 2 y
b = 3.

Solución

(i) Tenemos el arreglo que se muestra en la siguiente figura.

Dado que tenemos simetrı́a azimutal, la coordenada ϕ no aparece en este problema. La solución general
a la ecuación de Laplace (∇2 φ = 0) en coordenadas esféricas con esta simetrı́a es de la forma
∞ h
X i
φ(r, θ) = Al rl + Bl r−(l+1) Pl (cos θ).
l=0

3
Entonces, dividiendo el problemas por zonas, como en la figura, tenemos que
 P∞  I l 
 I −(l+1)
φI (r, θ) = Pl=0 Al r + Bl r Pl (cos θ) 0≤r≤a

φ(r, θ) = φII (r, θ) = ∞ A II l
r + B II −(l+1)
r P (cos θ) a≤r≤b .

 Pl=0 l
∞  III l
l
III −(l+1)
l
φIII (r, θ) = l=0 Al r + Bl r Pl (cos θ) b≤r

Como el potencial debe ser finito en el origen y nulo cuando r → ∞ , entonces

BlI = AIII
l =0 ∀l.

Renombrando los coeficientes tenemos que:

 P∞
 φI (r, θ) = Pl=0 Al rl Pl (cos θ)  0≤r≤a
∞
φ(r, θ) = φII (r, θ) = Pl=0 Bl rl + Cl r−(l+1) Pl (cos θ) a ≤ r ≤ b . (1)
 ∞
φIII (r, θ) = l=0 Dl r−(l+1) Pl (cos θ) b≤r

Localizando el eje z sobre el eje perpendicular al plano que separa a los cascarones y que pasa por sus
centros, entonces las condiciones de frontera son
 
φ0 ≤ θ ≤ π/2  0 ≤ θ ≤ π/2
φa0 (θ)|r=a = y φb0 (θ)r=b =
0 π/2 ≤ θ ≤ π φ0 π/2 ≤ θ ≤ π

Empleando las condiciones de contorno con r = a y r = b y los potenciales de la ec. (1) se obtiene
(i)
∞
X ∞ h
X i
φI (r = a, θ) = φa0 (θ) = Al al Pl (cos θ) = Bl al + Cl a−(l+1) Pl (cos θ) = φII (r = a, θ),
l=0 l=0

(ii)
∞ h
X i ∞
X
φII (r = b, θ) = φb0 (θ) = Bl bl + Cl b−(l+1) Pl (cos θ) = Dl b−(l+1) Pl (cos θ) = φIII (r = b, θ).
l=0 l=0

Dado que los los polinomios de Legendre son una base ortonormal, de (i) y (ii) podemos igualr los términos
de la dos suma, por lo que los coeficientes Bl y Cl se relacionan con Al y Dl como
Cl
Al al = Bl al + , (2)
al+1
Dl Cl
= Bl bl + l+1 . (3)
bl+1 b

4
Asimismo, las relaciones de ortogonalidad de los polinomios de Legendre, se obtienen los coeficientes Al y
Dl de (i) y (ii). Sus expresiones en términos de φa0 y φb0 son
Z π Z 1
2
Al al = φa0 (θ)Pl (cos θ) sin θdθ = φ0 Pl (x)dx,
2l + 1 0 0
Z π Z 0
Dl 2
= φb0 (θ)Pl (cos θ) sin θdθ = φ0 Pl (x)dx,
bl+1 2l + 1 0 −1

en donde Además se empleo el cambio de variable x = cos θ. Además como los polinomios de Legendre
tienen paridad definida Pl (−x) = (−1)l Pl (x), los coeficientes Al y Dl se escriben como
Z 1
2l + 1
Al = φ0 Pl (x)dx,
2al 0
Z 1
l+1 2l +1
Dl = b (−1)l φ0 Pl (x)dx.
2 0

La integral en las expresiones de los coeficientes puede resolverse empleando la relación

(2l + 1)Pl (x) = P 0 l+1 (x) − P 0 l−1 (x),

en donde los términos primados se refieren a la derivada respecto al argumento. Dado que l sólo corre en
enternos no negativos, el caso l = 0 se hará de forma independiente. Entonces para l 6= 0
Z 1 1
1 
Pl (x)dx = [Pl+1 (x) − Pl−1 (x)]  ,
0 2l + 1 0

pero Pl (1) = 1 para todo l. Además para l par, l ± 1 es un número impar por lo que al evaluar los
polinomios de Legengre de estos grados en x = 0 el resultado será nulo. Por este motivo sabemos que lo
únicos términos que contribuirán al potencial son l impar para l > 0. Es decir

Al = Dl = 0, ∀ l 6= 0 par, (4)
φ0
Al = − l [Pl+1 (0) − Pl−1 (0)] , ∀ l impar, (5)
2a
φ0 l+1
Dl = b [Pl+1 (0) − Pl−1 (0)] , ∀ l impar. (6)
2

Para l = 0, P0 (x) = 1 por lo que el resultado de la integral en las expresiones de los coeficientes A0 y B0
es 1. Por lo tato
φ0 φ0
A0 = , y D0 = b.
2 2
Resolviendo entonces para B0 y C0 con las ecs. (2) y (3) se obtiene que

φ0
B0 = , y C0 = 0. (7)
2

Si hacemos α = φ20 [Pl+1 (0) − Pl−1 (0)], entonces Al = −α/al y Dl = αbl+1 . Sustituyendo Al y Dl en las
ecs. (2) y (3), y resolviendo para Bl llegamos a las expresiones

α Cl α Cl
Bl = − − 2l+1 = l − 2l+1 .
al a b b

5
EL procedimiento para la resolción de Cl es
   
α Cl α Cl 1 1 1 1
− l − 2+1 = l − 2+1 ⇒ Cl − 2l+1 + 2l+1 = α l + l ,
a a b b a b b a
a2l+1 − b2l+1 al + bl
⇒ Cl = α ,
(ab)2l+1 (ab)l
(al + bl )(ab)l+1
⇒ Cl = α 2l+1 . (8)
a − b2l+1
Asimismo, para Bl
α Cl
Bl = − − 2l+1
al a
 
1 1 (al + bl )(ab)l+1
= −α l + l+1 l 2l+1
a a a a − b2l+1
 l l l+1

1 1 (a + b )b
= −α l + l 2l+1
a a a − b2l+1
al+1 + bl+1
= −α . (9)
(a2l+1 − b2l+1 )

Sustituyendo las ecs. (7), (8) y (9) en la ec. (1) se obtiene la solución del potencial en la región a ≤ r ≤ b,
dada por
( ∞  )
φ0 X Pl−1 (0) − Pl+1 (0)
l a2l+1 b+1 + al+1 b2l+1
l+1 l+1
φII (r, θ) = 1+ Pl (cos θ) − a +b r + . (10)
2 b2l+1 − a2l+1 rl+1
l=1

Recordando que los términos para l para no contribuyen, reetiquetemos l → 2l − 1 para únicamente
considerar términos impares en la suma, Además, factoricemos b2l+1 del numerador dentro de la suma y
reescribamos todo en términos de a/b. Esto lleva al resultado
( ∞
φ0 X P2(l−1) (0) − P2l (0)
φII (r, θ) = 1+ P2l−1 (cos θ)
4l−1 ×
2 1 − ab
l=1
   a 2l   r 2l−1   a 2l−1   a 2l 
× − 1+ + 1+ . (11)
b b b r

La expresin en la ec. (10) es últil para hacer las aproximaciones para b → ∞ y a → 0, mientras que la ec.
(11) se recomiendo emplear para el cálculo numérico del potencial ya que no calcula los l pares, siendo
más eficiente. Asimismo, la expresin en (11) hace evidente que todo dentro de la suma es adimensional.
(ii)

Ahora consideraremos los términos hasta l = 4, recordando que los términos impares no contribuyen.
Entonces
  2   3 2 
φ0 a + b2 a b + a2 b3 1
φII (r, θ) ≈ 1 + P1 (cos θ)[P0 (0) − P2 (0)] r− +
2 a3 − b3 a3 − b3 r2
 4   7 4  
a + b4 3 a b + a4 b7 1
+ P3 (cos θ)[P2 (0) − P4 (0)] r − .
a7 − b7 a7 − b7 r4

Las polinomios de Legendre
para l = 0, · · · , 4 son P0 (cos θ)
= 1, P1 (cos θ) = cos θ,

P2 (cos θ) = 21 3 cos2 θ − 1 , P3 (cos θ) = 12 5 cos3 θ − 3 cos θ y P4 (cos θ) = 81 35 cos4 θ − 30 cos2 θ + 3 ,

6
 entonces
  2   3 2  
φ0 3 a + b2 a b + a2 b3 1
φII (r, θ) ≈ 1 + cos θ r−
2 2 a3 − b3 a3 − b3 r2
 4   7 4  
7 3

a +b 4
3 a b + a4 b7 1
− 5 cos θ − 3 cos θ r − .
16 a 7 − b7 a7 − b7 r4
Para el caso b −→ ∞:
 
φ0 3 a2 7
a4
φII (r, θ)b→∞ = 1 + (cos θ) 2 − 5 cos3 θ − 3 cos θ 4 ,
2 2 r 16 r
y para el caso a −→ 0:
 
φ0 3 r 7
r3
φII (r, θ)a→0 = 1 − (cos θ) + 5 cos3 θ − 3 cos θ 3 .
2 2 b 16 b
(iii)

Para el cálculo numérico se empleó la ec. (11) con l = 15, lo que es equivalente a tener P29 (x) como el
polinomio de Legendre de mayor grado. A continuación se presenta la gráfica del potencial obtenida en
Mathematica.

ϕ(x,y) [V]
1.0

0.8

0.6

0.4

0.2

4. Problema: (30pts) La superficie de una esfera conductora de radio a, centrada en el origen, se divide
en cuatro secciones (“gajos”) iguales por planos que se intersectan en el eje z. Los “gajos” se mantienen
alternadamente a potenciales φ0 y −φ0 .

7
(i) Calcula el potencial como una representación en serie (en base esférica) para todos los puntos interiores
a la esfera.
(ii) Calcula explı́citamente los términos hasta l = 2.
(iii) Calcula como una representación en serie (en base esférica) la densidad superficial de carga en la
superficie de la esfera.
(iv) Grafica en el plano Z = 0 (plano XY ) el potencial que calculaste en (i) considerando φ0 = 1 y a = 1,
tomando suficientes sumandos como para que la solución no cambie en la tercera cifra decimal. Puedes
graficar con código de color.
Solución
El sistema desde una vista superior está representado en la siguiente figura.

(i) Puesto que la carga interior es cero, se tiene que ∇2 φ = 0. La solución en coordenadas esféricas,
considerando que φ → 0 cuando r → 0, es
∞ X
X l
φ(r < a, θ, ϕ) = Al,m rl Yl,m (θ, ϕ).
l=0 m=−l

Al aplicar las condiciones de frontera entre el interior y la superficie se tiene

∞ X
X l
φ(r = a, θ, ϕ) = Al,m al Yl,m (θ, ϕ),
l=0 m=−l

−φ0 0 < ϕ ≤ π/2, π < φ ≤ 3φ/2
donde φ(r = a, θ, ϕ) = . (12)
φ0 π/2 < ϕ ≤ π, 3π/2 < ϕ ≤ 2π
Para calcular de forma explı́cita los coeficientes Al,m , usaremos la ortogonalidad de los armónicos
esféricos. Multiplicando por Yl∗0 ,m0 e integrando en todo el ángulo sólido se tiene
Z 2π Z 1
∗
φ(r = a, θ, ϕ)Yl,m d(cos θ)dϕ =
0 −1

∞ X
X l Z 2π Z 1
0
∗
Al,m a l
Yl,m (θ, ϕ)Yl∗0 ,m0 (θ, ϕ)d(cos θ)dϕ = Al0 ,m0 al ,
l=0 m=−l 0 −1

en particular
Z Z s
2π 1
1 (2l + 1)(l − m)! m
Al,m = l φ(r = a, θ, ϕ) Pl (x) e−imϕ dx dϕ con x = cos θ
a 0 −1 4π(l + m)!
s Z Z
1 2π
1 (2l + 1)(l − m)!
=⇒ Al,m = l Plm (x) dx φ(r = a, θ, ϕ) e−imϕ dϕ.
a 4π(l + m)! −1 0

8
 La segunda integral la separamos por regiones
Z 2π Z π/2 Z π Z 3π/2 Z 2π
φ(r = a, θ, ϕ) e−imϕ dϕ = (−φ0 ) e−imϕ dϕ+ φ0 e−imϕ dϕ+ (−φ0 ) e−imϕ dϕ+ φ0 e−imϕ dϕ
0 0 π/2 π 3π/2

π/2 π 3π/2 2π
iφ0 −imϕ  iφ0 −imϕ  iφ0 −imϕ  iφ0 −imϕ 
=− e  + e  − e  + e  ,
m 0 m π/2 m π m 3π/2

evaluando y usando que e−2imπ = cos(2mπ) − i sin(2mπ) = 1, se tiene

Z 2π
2iφ0  
φ(r = a, θ, ϕ) e−imϕ dϕ = 1 + e−imπ − e−imπ/2 − e−3imπ/2 .
0 m
Ası́, para m = 0, dicha integral es cero. Para m impar, se tiene
Z 2π    
2iφ0 i, m = 1, 5, 9, ... i, m = 3, 7, 11, ...
φ(r = a, θ, ϕ) e−imϕ dϕ = 1−1± ± = 0;
0 m −i, m = 3, 7, 11, ... −i, m = 1, 5, 9, ...
para m par
Z 2π 
2iφ0   4, m = 2, 6, 10, ...
φ(r = a, θ, ϕ) e−imϕ dϕ = 1 + 1 − (−1)k − (−1)k =
0 m 0, m = 4, 8, 12, ...

donde m = 2k . Ası́ pues, sólo sobreviven los múltiplos de 2 y que no son múltiplos de 4

m = 2, 6, 10, 14, 18, ...

Por lo tanto
s  Z 1
1 (2l + 1)(l − m)! 8iφ0
Al,m = l Plm (x) dx. para m = 4n + 2, con n ∈ N.
a 4π(l + m)! m −1

(ii) Ahora al tomar hasta l = 2 se tiene lo siguiente.

Para l = 0, φ(ϕ) = 0.
Para l = 1, como se tienen m = ±1, 0, entonces φ(ϕ) = 0.
Para l = 2, sólo sobreviven m = ±2. Y como P2 (x) = (3x2 − 1)/2, entonces

d2
P22 (x) = (−1)2 (1 − x2 ) P2 (x) = 3(1 − x2 ),
dx2
Z 1
P22 (x) dx = 4.
−1

Además, sabemos que

(1 − m)! m
Pl−m (x) = (−1)m P (x)
(1 + m)! l
1
⇒ P2−2 (x) = P22 (x)
4!
r r
4iφ0 5 4iφ0 5
⇒ A2,2 = 2 y A2,−2 = − 2 .
a 6π a 6π
Con lo anterior, el potencial está dado por
r
4iφ0 5 2
φ(~r) = 2 r [Y2,2 (θ, ϕ) − Y2,−2 (θ, ϕ)] ,
a 6π

9
 y como

1/2
15
Y2,±2 (θ, ϕ) = sin2 θe±2iϕ
32
r  1/2
4iφ0 5 2 15  
φ(~r) = 2 r sin2 θ e2iϕ − e−2iϕ ,
a 6π 32
y usando que sin2 θ = 1 − cos2 θ y 2i sin(2ϕ) = e2iϕ − e−2iϕ , finalmente se obtiene

5φ0 2
φ(~r ) = r sin(2ϕ)(cos2 θ − 1).
πa2
(iii) La carga inducida está dada por

~ ~ ∂φ 
σin = E · n̂ = E · r̂ = Er |r=a = − 
∂r r=a
 
∂  ∞
∂ X X
l 
 l 
σin = − φ(r, θ, ϕ) =− Al,m r Yl,m (θ, ϕ)
∂r r=a ∂r 
l=0 m=−l r=a

∞ X
X l 

σin = Al,m lrl−1 Yl,m (θ, ϕ)

l=0 m=−l r=a

(iv) Podemos meter directamente la expresión del potencial, en términos de los coeficientes Alm en
Mathematica indicando una n especı́fica para la suma. En particular, para l=8 se obtiene

7. Problema TORITO: (20pts) Un anillo de radio R tiene una densidad lineal uniforme λ. Una partı́cula
con carga positiva +q y masa m está inicialmente localizada en el centro del anillo y se le da un pequeño
empujón (lo que coloquialmente llamamos un “garnuchazo”). Si la partı́cula está restringida a moverse
sobre el plano que contiene al anillo, muestra que el movimiento corresponderá a un movimiento armónico
simple (para pequeñas oscilaciones), y calcula su frecuencia. Hint: Calcula la energı́a potencial de la
partı́cula cuando está a un radio r (pequeño), integrando sobre el anillo y luego calculando

10
el negativo de la derivada para calcular la fuerza. Necesitarás usar la ley de los cosenos y la
serie de Tylor de:
1  32
√ ≈1− + .
1+ 2 8
Solución
Consideremos un elemento de carga del anillo infinitesimal dado por Rdθ, como se muestra en la figura.

La energı́a potencial debido al movimiento de la carga q como una función de r está dada por
q(λR)dθ
dU = ,
4π0 l
usando la ley de cosenos para el triángulo mostrado en la figura e integrando sobre todo el anillo obtenemos
que Z π Z π
qλ dθ qλ dθ
U= √ = q .
2π0 0 R1 R2 + r2 − 2Rr cos θ 2π0 0 1+ r2
− 2 r cos θ
R2 R
2
Ahora, utilicemos el hint del problema, sea  = Rr 2 − 2 Rr cos θ la raı́z de la integral, que se puede aproximar
como    2
1 1 r2 r 3 r2 r
q ≈1− − 2 cos θ + − 2 cos θ .
2
1 + r − 2 r cos θ 2 R2 R 8 R2 R
R2 R

Conservando términos hasta segundo orden en r de la expresión anterior, resulta que

Z π   
qλ r2 2 qλ r2
U≈ 1+ (3 cos θ − 1) dθ = 1 + .
2π0 0 2R2 2π0 4R2

El término cuadrático en la energı́a potencial es ya una prueba de que el movimiento de la carga es

oscilatorio, al menos para r  1. Para calcular la frecuencia de la oscilación, recordemos que, F = mr̈ =
−∇U . Para nuestro caso tenemos
qλ qλ qλ
F = mr̈ = −∇U = − 2
r ⇒ mr̈ = − 2
r ⇒ r̈ + r = 0,
4π0 R 4π0 R 4π0 mR2
que resulta en la ecuación de un oscilador armónico con frecuencia de oscilación

r
qλ
ω= .
4m0 R2

11