Geser GeneralPhysics 2018

Lectures: Problemas de fı́sica general
por
Lic. Geser, Federico Alejandro
UNIVERSIDAD NACIONAL DE CÓRDOBA
30 de octubre de 2018
FAMAF
c - UNC 2018
Índice
1. Problemas de cinemática y dinámica 1
1.0.1. Cuerpo en el manto de un cono . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
2. Problemas de electrostática básica 4

2.1. Comparación: fuerza eléctrica vs. fuerza gravitatoria . . . . . . . . . . . . . . . . . . . 4
2.2. Cristales y fuerzas electrostáticas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
2.3. Cuadrado de cargas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
2.4. Tubo de rayos catódicos: funcionamiento básico . . . . . . . . . . . . . . . . . . . . . . 8
2.5. Fuerzas moleculares: Cloruro de Cesio . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2.6. Hilo cargado uniformemente . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.7. Disco cargado uniformemente . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.8. Anillo cargado uniformemente . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.9. Esfera cargada uniformemente . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.10. Hilo cargado uniformemente: Gauss . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.11. Cilindro de carga no-uniforme . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
2.12. Tubo metálico cargado . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.13. Placa infinita cargada . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
2.14. Campo eléctrico “Rutherfordiano” . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
2.15. Concepto de material conductor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
2.16. Potential: carga puntual . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
2.17. Potencial: cargas colineales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
2.18. Potencial: anillo cargado uniformemente . . . . . . . . . . . . . . . . . . . . . . . . . . 29
2.19. Potencial: superficies equipotenciales en conductores . . . . . . . . . . . . . . . . . . . 30
2.20. Par de placas infinitas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
3. Problemas de magnetostática básica 33

3.1. Protones en movimiento en un campo magnético . . . . . . . . . . . . . . . . . . . . . 33
3.2. Electrones en movimiento en un campo magnético . . . . . . . . . . . . . . . . . . . . 33
3.3. Conductor en un campo magnético . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
3.4. Comparación: fuerza magnética vs. fuerza gravitatoria . . . . . . . . . . . . . . . . . . 35
4. Problemas de óptica 37
4.1. Reflexión y refracción en un prisma . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
5. Matemática general 40
5.1. Cambio de sistemas de coordenadas . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
5.2. Ecuaciones diferenciales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
5.2.1. Ecuaciones diferenciales de primer orden homogéneas . . . . . . . . . . . . . . 41
5.2.2. Ecuaciones diferenciales de primer orden inhomogéneas . . . . . . . . . . . . . 41
5.2.3. Ecuaciones diferenciales de segundo orden . . . . . . . . . . . . . . . . . . . . . 42
5.2.4. Ecuaciones diferenciales en derivadas parciales . . . . . . . . . . . . . . . . . . 43
5.3. Tensores . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
5.3.1. La métrica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
5.3.2. Simetrı́a y antisimetrı́a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
5.3.3. Transformación de tensores en coordenadas curvilineas . . . . . . . . . . . . . . 44
A. Scattering cuántico 45
A.1. Determinación de la amplitud de dispersión . . . . . . . . . . . . . . . . . . . . 47
A.2. Aproximación de Born . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
A.3. Método de las ondas parciales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
A.4. Teorema óptico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
A.5. Espı́n en dispersión . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
B. Interacciones nucleares (fuerte y débil) 53
C. Teorı́a de perturbaciones 53
C.1. Perturbaciones constantes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
C.1.1. Regla de oro de Fermi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
C.1.2. Sección eficaz de transición . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
D. Aproximación dipolar 56
E. Cálculos usuales 56
E.1. Colisión frontal y máximo intercambio energético . . . . . . . . . . . . . . . . . 56
E.2. Energı́a de formación de una esfera cargada . . . . . . . . . . . . . . . . . . . . 58
E.3. Modelado de un dieléctrico: modelo de Drude . . . . . . . . . . . . . . . . . . . 59
E.4. Consideraciones sobre (ω) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
F. Ecuación de Dirac 59
G. Potencial, energı́a potencial, trabajo: campo electrostático 59
ii
Notación:
NA → número de Avogadro, NA = 6, 022 × 1023 mol−1 .
c → velocidad de la luz, c = 299,792,458 m/s.
me → masa del electrón, me = 9, 109 × 10−31 kg.
e → carga del electrón, e = −1,602 × 10−19 C.

e2
r0 → radio clásico del electrón, r0 = m e c2
= 2, 817 × 10−15 m.
h, ~ → constante de Planck, h = 6,626 × 10−34 Js , ~ = 1,054 × 10−34 Js.
G → constante de gravitación universal, G = 6,7 × 10−11 N m2 /kg 2 .
k → constante de Coulomb, k = 9 × 109 N m2 /C 2 .
mn , mp → masa neutron y proton respectivamente, mn = 1,674×10−27 kg, mp = 1,672×10−27 kg.
z ∗ → si z ∈ C entonces z ∗ es el complejo conjugado.
iii
Tabla de fórmulas usuales: Sistemas de Coordenadas
Sistema cartesiano Sistema cilı́ndrico Sistema esférico
~r(t) = xî + y ĵ + z k̂ ~r(t) = ρρ̂ + z k̂ ~r(t) = rr̂

~v (t) = ẋî + ẏ ĵ + ż k̂ ~v (t) = ρ̇ρ̂ + ρϕ̇ϕ̂ + ż k̂ ~v (t) = ṙr̂ + rϕ̇ sin(θ)ϕ̂
2
~a(t) = ẍî + ÿ ĵ + z̈ k̂ ~a(t) = (ρ̈ − ρϕ̇ )ρ̂ + +rθ̇θ̂
(2ρ̇ϕ̇ + ρϕ̈)ϕ̂ + z̈ k̂ ~a(t) = ar t̂ + aϕ ϕ̂ + aθ θ̂
Momento angular:
ar = r̈ − rθ̇2 − rϕ̇ sin2 (θ)
Momento angular:
Mx = m(y ż − z ẏ) aθ = rθ̈ + 2ṙθ̇ −
My = m(z ẋ − xż) Mx = m(ρż − z ρ̇) sin(ϕ) − rϕ̇ sin(θ) cos(θ)
Mz = m(xẏ − y ẋ) mρz ϕ̇ cos(ϕ) aϕ = 2ṙϕ̇ sin(θ) +
My = −m(ρż − z ρ̇) cos(ϕ) − 2rθ̇ϕ̇ cos(θ) +
Cartesianas en cilı́ndricas:
mρz ϕ̇ sin(ϕ) rϕ̈ sin(θ)
x = ρ cos(ϕ) Mz = mρ2 ϕ̇
Momento angular:
y = ρ sin(ϕ) M2 = m2 ρ2 ϕ̇2 (ρ2 + z 2 ) +
z = z m2 (ρż − z ρ̇)2 Mx = −mr2 (θ̇ sin(ϕ) +
î = cos(ϕ)ρ̂ − sin(ϕ)ϕ̂ ϕ̇ sin(θ) cos(θ) cos(ϕ))
Cilı́ndricas en cartesianas:
ĵ = sin(ϕ)ρ̂ + cos(ϕ)ϕ̂ My = −mr2 (θ̇ cos(ϕ) −
p
k̂ = k̂ ρ = x2 + y x ϕ̇ sin(θ) cos(θ) sin(ϕ))
y
ϕ = arctan Mz = mr2 ϕ̇ sin2 (θ)
Cartesianas en esféricas: x
z = z M2 = m2 r4 (θ̇2 + ϕ̇2 sin2 (θ))
x = r sin(θ) cos(ϕ) ρ̂ = cos(ϕ)î + sin(ϕ)ĵ
Esféricas en cartesianas:
y = r sin(θ) sin(ϕ) ϕ̂ = − sin(ϕ)î + cos(ϕ)ĵ p
z = r cos(θ) k̂ = k̂ r = x2 + y 2 + z 2
y
î = sin(θ) cos(ϕ)r̂ + ϕ = arctan
Elementos de longitud x
cos(θ) cos(ϕ)θ̂ − p !
x2 + y 2
sin(ϕ)ϕ̂ dlρ = dρ ; dlϕ = ρdϕ ; dlz = dz θ = arctan
z
ĵ = sin(θ) sin(ϕ)r̂ +
Elementos de área r̂ = sin(θ) cos(ϕ)î +
cos(θ) sin(ϕ)θ̂ +
cos(ϕ)ϕ̂ dSρ = ρdϕdz sin(θ) sin(ϕ)ĵ + cos(θ)k̂
k̂ = cos(θ)r̂ − sin(θ)θ̂ dSϕ = dρdz ϕ̂ = − sin(ϕ)î + cos(ϕ)ĵ
dSz = ρdρdϕ θ̂ = cos(θ) cos(ϕ)î +
= cos(θ) sin(ϕ)ĵ − sin(θ)k̂
Elemento de volumen
Elementos de longitud
dV = ρdρdϕdz
dlr = dr ; dlθ = rdθ
Relaciones entre versores dlϕ = r sin(θ)dϕ
î × ĵ = k̂ Elementos de superficie
ρ̂ × ϕ̂ = k̂ dSr = r2 sin(θ)dθdϕ
r̂ × θ̂ = ϕ̂ dSθ = r sin(θ)drdϕ
r̂˙ = θ̇θ̂ dSϕ = rdrdθ
˙
θ̂ = −θ̇r̂
Elemento de volumen
dV = r2 sin(θ)drdθdϕ
iv
Tabla de fórmulas usuales: Cálculo Vectorial
Identidades vectoriales
~ ) = ∇(∇ · V
∇ × (∇ × V ~ ) − (∇ · ∇)V
~ = ∇(∇ · V
~ ) − ∇2 V
~
~) = 0
∇ · (∇V
∇×∇×∇×V ~ = −∇2 (∇ × V ~)
∇ × (∇φ) = 0
~ ×W
∇ · (V ~) = W~ · (∇ × V
~)−V ~ · (∇ × W
~)
~ ×W
∇ × (V ~ ) = (W~ · ∇)V
~ − (V
~ · ∇)W~ +V~ (∇ · W
~ )−W
~ (∇ · V
~)
∇(V~ ·W
~ ) = (W~ · ∇)V
~ + (V
~ · ∇)W~ +W~ × (∇ × V~)+V~ × (∇ × W
~)
ψ(∇2 φ) = ∇ · (ψ∇φ) − ∇ψ · ∇φ
Sistema cartesiano
∂φ ∂φ ∂φ
∇φ = î + ĵ + k̂
∂x ∂y ∂z
~ ∂Vx ∂Vy ∂Vz
∇·V = + +
∂x ∂y ∂z

~ ∂Vz ∂Vy ∂Vx ∂Vz ∂Vy ∂Vx
∇×V = − î + − ĵ + − k̂
∂y ∂z ∂z ∂x ∂x ∂y
∂2φ ∂2φ ∂2φ
∇2 φ = + 2 + 2
∂x2 ∂y ∂z
Sistema cilı́ndrico
∂φ 1 ∂φ ∂φ
∇φ = ρ̂ + ϕ̂ + k̂
∂ρ ρ ∂ϕ ∂z
~ 1 ∂ 1 ∂Vϕ ∂Vz
∇·V = (ρVρ ) + +
ρ ∂ρ ρ ∂ϕ ∂z

~ 1 ∂Vz ∂Vϕ ∂Vρ ∂Vz 1 ∂ ∂Vρ
∇×V = − ρ̂ + − ϕ̂ + (ρVϕ ) − k̂
ρ ∂ϕ ∂z ∂z ∂ρ ρ ∂ρ ∂ρ
1 ∂ 2 φ ∂ 2 Vz

1 ∂ ∂φ
∇2 φ = ρ + 2 +
ρ ∂ρ ∂ρ ρ ∂ϕ2 ∂z 2
Sistema esférico
∂φ 1 ∂φ 1 ∂φ
∇φ = r̂ + θ̂ + ϕ̂
∂r r ∂θ r sin(θ) ∂ϕ
∇·V~ = 1 ∂ (rVr ) + 1 ∂
(sin(θ)Vθ ) +
1 ∂Vϕ
r2 ∂r r sin(θ) ∂θ r sin(θ) ∂ϕ

~ 1 ∂ ∂Vθ 1 1 ∂Vr ∂
∇×V = (sin(θ)Vϕ ) − r̂ + − (rAϕ ) θ̂+
r sin(θ) ∂θ ∂ϕ r sin(θ) ∂ϕ ∂r

1 ∂ ∂Vr
(rVθ ) − ϕ̂
r ∂r ∂θ
1 ∂2φ

2 1 ∂ 2 ∂φ 1 ∂ ∂φ
∇ φ= 2 r + sin(θ) +
r ∂r ∂r sin(θ) ∂θ ∂θ sin2 (θ) ∂ϕ2
Identidades de Green
Z Z
f (~x)∇2 g(~x) − g(~x)∇2 f (~x) d3 ~x = (f (~x)∇g(~x) − g(~x)∇f (~x)) · n̂ d2 S

Z V Z∂V
2
3
g(~x)∇ f (~x) + ∇g(~x) · ∇f (~x) d ~x = f (~x)∇g(~x) · n̂ d2 S
V ∂V
v
Polinomios de Legendre:
Pl (−x) = (−1)l Pl (x)
Z 1
P0 (x) = 1 2
Pl (x)Pl0 (x)dx = δll0
P1 (x) = x −1 2l + 1
1
(−1)l (2l)!
Z
1 P2l+1 (x)dx =
P2 (x) = (3x2 − 1) (2l + 1)(2l l!)2
2 0
5 3 (l + 1)Pl+1 (x) = x(2l + 1)Pl (x) − lPl−1 (x)
P3 (x) = x3 − x
2 2 dPl+1 dPl−1
1 (2l + 1)Pl (x) = (x) − (x)
35x4 − 30x2 + 3

P4 (x) = dx dx
8 dPl
(x2 − 1) = lxPl (x) − lPl−1 (x)
(
0 si l es impar dx
Pl (0) = (−1)l/2 (l−1)!!
2l/2 (l/2)!!
si l es par
Polinomios de Legendre asociados:
1 l+m
2 m/2 d
Plm (x) = (−1)m (1 − x ) (x2 − 1)l ∀m ∈ Z
2l l! dxl+m
Z 1
2 (l + m)!
Plm (x)Pl0 m dx = δll0
−1 2l + 1 (l − m)!
Plm (−x) = (−1)l+m Plm (x)
Armónicos esféricos:
∗
s Yl,−m (θ, φ) = (−1)m Yl,m (θ, φ)
2l + 1 (l − m)! Z 2π Z π
Yl,m (θ, φ) = Plm (cos(θ))eimφ
4π (l + m)! dφ sin(θ)dθYl,m (θ, φ)Yl0 ,m0 (θ, φ)
0 0
1
Y0,0 (θ, φ) = √ = δll0 δmm0
4π ∞ X
l
r X
∗
Y1,0 (θ, φ) =
3
cos(θ) Yl,m (θ, φ)Yl,m (θ0 , φ0 ) =
4π l=0 m=−l
r
3 = δ(φ − φ0 )δ(cos(θ) − cos(θ0 ))
Y1,±1 (θ, φ) = ∓ sin(θ)e±iφ
8π
r
2l + 1
Yl,0 (θ, φ) = Pl (cos(θ))
r
5 4π
Y2,0 (θ, φ) = (3 cos2 (θ) − 1)
16π
r
15
Y2,±1 (θ, φ) = ∓ sin(θ) cos(θ)e±iφ
8π
r
15
Y2,±2 (θ, φ) = sin2 (θ)e±2iφ
32π
vi
Tabla de fórmulas usuales: Mecánica
Formulación Hamiltoniana Corchetes de Poisson

X X ∂f ∂g ∂f ∂g

H= pi q̇i − L [f, g] = −
i ∂pk ∂qk ∂qk ∂pk
k
∂H ∂H ∂f
ṗi = − , q̇i = [f, pk ] = − [pk , f ] = −
∂qi ∂pi ∂qk
∂f
Sistema cartesiano [f, qk ] = − [qk , f ] =
∂pk
1 2
H= (p + p2y + p2z ) + U (x, y, z) [pi , pk ] = [qi , qk ] = 0 ; [pi , qk ] = δik
2m x
pi ∂U Identidad de Jacobi
ẋi = , ṗi = −
m ∂xi [f, [g, h]] + [g, [h, f ]] + [h, [f, g]] = 0
Sistema cilı́ndrico
! Teorema de Poisson
1 p2ϕ
H= p2ρ + 2 + p2z + U (ρ, ϕ, z) df dg d [f, g]
2m ρ = =0 ⇒ =0
dt dt dt
pρ pϕ pz Condiciones de transformación canónica
ρ̇ = ; ϕ̇ = ; ż =
m mρ2 m
[Qi , Qk ]p,q = [Pi , Pk ]p,q = 0 ; [Pi , Qk ]p,q = δik
2p2ϕ ∂U ∂U ∂U
ṗρ = 3
− ; ṗϕ = − ; ṗz = −
ρ ∂ρ ∂ϕ ∂z Transformaciones canónicas
Sistema esférico ∂F1 ∂F1
F1 = F1 (q, Q, t) ⇒ p= ; P =−
! ∂q ∂Q
1 p2ϕ p2θ
H= p2r + + + U (r, θ, ϕ) ∂F2 ∂F2
2m r2 sin (θ)2 r2 F2 = F2 (q, P, t) ⇒ p= ; Q=
∂q ∂P
Pr pϕ pθ ∂F3 ∂F3
ṙ = ; ϕ̇ = ; θ̇ = F3 = F3 (p, Q, t) ⇒ q=− ; Q=
m 2 2 mr 2 ∂p ∂Q
mr sin (θ)
∂F4 ∂F4
F4 = F4 (p, P, t) ⇒ q=− ; Q=
d∂p ∂P
2p2ϕ p2 ∂U
ṗr = + 3θ − ∂Fi
r 3 r ∂r K = H+
∂t
2p2ϕ cos(θ) ∂U
ṗθ = −
r2 sin3 (θ) ∂θ
∂U
ṗϕ = −
∂ϕ
Problema de Kepler
Velocidad angular con
ángulos de Euler Vector de Lenz
α ~ = ~v × J~ − αr̂
Si U = − ⇒L
Ω1 = φ̇ sin(θ) sin(Ψ ) + θ̇ cos(Ψ ) r
Ω2 = φ̇ sin(θ) cos(Ψ ) − θ̇ sin(Ψ ) Trayectoria
Ω3 = φ̇ cos(θ) + Ψ̇
J2
r =
m(L cos(ϕ) + α)
P
Energı́a cinética de sistema =
1 + e cos(ϕ)
no-inercial rotante
J2 L
1 P = ; e=
T = ~ × ~r) + 1 m(Ω
mv 2 + m~v · (Ω ~ × ~r)2 αm α
2 2
vii
Tensor de inercia Pequeñas oscilaciones
Z
Iik = ρ(r2 δik − xi xk )dV Si se llega a que
X 1 ˙ ˙T 1
Iik = m(r2 δik − xi xk ) T = ~x T ~x y U = ~x U ~xT
2 2
Teorema de Steiner entonces, las ecuaciones de movimiento son
I = ICM + M d2 ¨ + U~x = 0
T~x
Donde d es la distancia entre ejes paralelos El problema se resuelve encontrando un sistema

diagonalizado
~q¨ = −H~q
Relación entre las velocidades
del centro de masa para dos donde las coordenadas ~q cumplen
sistemas inerciales ligados
a un cuerpo rı́gido ~q = A−1 ~x ; ~x = A~q
~0 =V~CM + Ω
~ × ~a siendo A la matriz de
V CM autovalores de H.
en donde ~a apunta del sistema “viejo” al
“nuevo” (primado).
Para un eje de rotación

~0 =0
instantáneo VCM
viii
Tabla de fórmulas usuales: Electromagnetismo
Ecuaciones de Maxwell Energı́a del campo electromagnético: vacı́o

Z
Caso general 1 ~ ·E ~ +B ~ ·B

~ d3 ~x
= E
8π R3
~
∇·D = ρ Z
d
~ = P− ~s · n̂d2 S
∇·B = 0 dt ∂V
~
∂B
Z Z
~
∇×E = − = − ~ ~ 3
j · Ed ~x − ~s · n̂d2 S
∂t V ∂V
∂D~ c ~ ~ el vector de Poynting.
~
∇×H = J~ + ~s = E×B
∂t 4π
Caso vacı́o Energı́a de la ecuación de onda

ρ
~
∇·E = 1 ∂2φ
0 0 = 2 ∂t2
− ∇2 φ
c
~
∇·B = 0
Z
1
(∂t φ)2 + ∇φ · ∇φ d3 ~x

~ ⇒ =
~ ∂B R3 2
∇×E = −
∂t 1 ∂2φ
~ 0 = − ∇2 φ + m2 φ
~ ∂E c2 ∂t2
∇×B = µ0 J~ + µ0 0 Z
∂t 1
(∂t φ)2 + ∇φ · ∇φ + m2 φ2 d3 ~x

⇒ =
Caso dependencia temporal armónica R3 2
~ x, t)
E(~ = ~ x)eiωt
E(~ Tensor de energı́a-momento
~ x, t)
B(~ ~ x)eiωt
= B(~
~ 1 2
∇·D = ρ Tij = 0 Ei Ej − δij E +
2
~
∇·B = 0
1 1 2
~
∇×E = −iω B~ Bi Bj − δij B
µ0 2
~
∇×H ~ + J~
= iω D δij ~ ~ ~ ·B
~

Tij = Ei Dj + Hi Bj − E·D+H
2
Electrostática
ρ
Se resuelve la ecuación de Poisson ∇2 φ = 0 , y en vacı́o la de Laplace ∇2 φ = 0, en donde se tiene que
~ = −∇φ.
E
Solución de Laplace con simetrı́a azimutal

∞
X
φ(~x) = Al rl + Bl r−(l+1) Pl (cos(θ))
l=0
Solución de Laplace sin simetrı́a azimutal

∞ X
X l
φ(~x) = Alm rl + Blm r−(l+1) Yl m (θ, ϕ)
l=0 m=−l
Condiciones de contorno:
Neumann → n̂ · ∇φ|S = 0 o f (~x) con ~x ∈ S.
Dirichlet → φ|S = 0 o f (~x) con ~x ∈ S.
Cauchy → combinación de las anteriores.
ix
Expansión multipolar del potencial
∞ X
l
Y m (θ, ϕ)
Z
X 1
φ(~x) = 4π Qlm l l+1 ; Qlm = ρ(~x)rl (Yl m (θ, ϕ))∗ d3 ~x
2l + 1 r
l=0 m=−l V
Dipolo eléctrico:

p~ · ~r ~ = −∇ p~ · ~r
φ= ; E
r3 r3
ε = p~ · E ~ → Energı́a del dipolo
Expansión de la función de Green

∞ l
1 X r<
= P (cos(γ))
l+1 l
|~x − ~x 0 | r
l=0 >
∞ X
l l
1 X 1 r< ∗
= −4π Y m (θ, ϕ) (Yl m (θ, ϕ))
|~x − ~x 0 | 2l + 1 r> l
l+1
l=0 m=−l
en donde r> = max {|~x|, |~x 0 |}, r< = min {|~x|, |~x 0 |}, γ el ángulo entre ~x y ~x 0 .
Transformaciones de Lorentz
Campos generales
2
~ − γ β~ β

~0
E = γ E~ +β ~×B ~·E
~
γ+1
2
~ − γ β~ β

~0
B = γ B~ − β~ × E ~·B
~
γ+1
Campos alineados con la dirección del movimiento relativo entre marcos de referencia
~k 0
E = ~k
E

~⊥ 0
E = γ E ~ ⊥ − cβ~ × B
~
~k 0
B = ~k
B
!
~
~⊥
B 0
= γ ~ ⊥ + β~ × E
B
c
x
Magnetostática:
I Z
Ley de Amper: ~ · d~l = 4π
B J~ · n̂da
C c S
~
Ley de Biot-Savart: dB~ = I l × ~r
c r2
1
Z ~ x 0)
J(~
~
Potencial vector: A(~x) = d3 ~x0
c R3 |~x − ~x 0 |
Magnetización: Ĥ = B̂ − 4π M ~
!
~ × ~r
M M~ × ~r
~=
Dipolo magnético: A ; B ~ =∇×
r3 r3
Z
∇0 · M ~ (~x 0 ) I ~ (~x 0 )
n̂ 0 · M
3 0
Potencial con magnetización: φm = − d ~
x + da0
V |~x − ~x 0 | S |~ x − ~x 0 |
Condiciones de contorno dieléctricos:

~1 −D
0 = (D ~ 2 ) · n̂
S
= (ε1 ∇φ1 − ε2 ∇φ2 ) · n̂|S

~1 − E
4πσ = (E ~ 2 ) · n̂ carga inducida sobre la interfaz
S
~1 −D
σ = (D ~ 2 ) · n̂ carga libre
S
)
0 ~1 − E
= (E ~ 2 ) × n̂
S si no hay cargas (0.1)
0 = (∇φ1 − ∇φ2 ) × n̂|S
Condiciones de contorno magnetostática:

~1 − B
0 = (B ~ 2 ) · n̂
S
4π
k̂ = (∇φm1 − ∇φm2 ) · n̂|S
c
Condiciones de contorno conductor:

n̂ × E~ = 0 ; n̂ · E
~ = σ
S S

~ ~
n̂ · H = 0 ; n̂ × H = k̂

S S
xi
Tensor de campo electromagnético: Ecuaciones varias:
E
− Ecx − Ecz
 
0 − cy
Ex
0 −Bz By  ~ = E
D ~
F µν = 

c 
Ey
Bz 0 −Bx  ~ = 1B

c H ~
Ez
c −By Bx 0 µ
Ecuación de continuidad: J~ = σ E
~ , σ = conductividad
I
∂ρ P = ~ · d~a
S
∇ · J~ + =0 S
∂t
Densidad de momento lineal y angular:
Z
~ ~ ~ ~ ~
Pe.m. = µ0 0 S = 0 E × B → P = P~e.m. d3 ~x
V
Z
~ ~ ~
Je.m. = ~r × Pe.m. → L = ~r × P~e.m. d3 ~x
V
Maxwell para guias de ondas: Ecuaciones varias guia de onda:
∂E ~t P = e−2γz Pideal , Potencia guia de

~ t = ∇t Ez
+ iωẑ × B
∂z ondas con atenuación γ
~ t ) = iωBz
ẑ · (∇t × E ~k
K̂ef = n̂ × H
∂B~t
− iµωẑ × E ~ t = ∇t Bz dPloss 1 ~ 2
∂z = |Kef |
da 2σδ
~ t ) = −iµωEz
ẑ · (∇t × B 1 ~ 2
= ~k
|Hk | si n̂ ⊥ H
~ t = − ∂Ez 2σδ
∇t · E
∂z
~t = − z ∂B
∇t · B
∂z
Modo TE: Modo TM:
Ez = 0 Bz = 0

∂Bz Ez |S = 0
=0
∂ n̂
S ~ t = ik ∇t Ez
E
ik γ2
Ĥt = ∇t Bz 2
γ2 ~ t = iµω ẑ × ∇t Ez
B
~ t = − i ωẑ × ∇t Bz γ 2
E
γ2 2 2
∇t Ez + γ Ez = 0
∇2t Bz + γ 2 Bz = 0
Siempre γ 2 = µω 2 − k 2 , con µ los coeficientes

del material de la cavidad y k es el de la propa-
gación eikz .
xii
Modo TEM:

∂ ∂
E z = B z = 0 ; ∇t = , ,0
∂x ∂y
~ t = 0 ∇t × B
∇t × E ~t = 0 ⇒ E ~ t = −∇t φ ∆t φ = 0
√
~ t = ± µẑ × E
~t ; E
~ =E ~ t e−γz eikz−iωt ; B
~ =B
~ t e−γz eikz−iωt
B
Potencia asociada a una guia de ondas:
1 ∂P
γ es el coeficiente de atenuación ⇒ γ = −
2P ∂z
2 1/2
ω2
Z
1 ω
P = √ 1 − λ2 |ψ|2 da
2 µ ωλ ω µ A
 2 
2 I  1 ∂ψ
dP 1 ω 2 2
µωλ ∂ n̂

= 2 dl
dz 2σδ ωλ 1 ω 2 ω2
∂A 
µω 2
1 − ωλ2 |n̂ × ∇t ψ| + 2 |ψ|2
ωλ

λ
En modo TE ⇒ se eligen las opciones de abajo de las llaves y ψ = Bz

En modo TM ⇒ se eligen las opciones de abajo de las llaves y ψ = Ez
xiii
Tabla de fórmulas usuales: Termodinámica
Ecuación de difusión
1 ∂ψ
∇2 ψ − = 0.
κ ∂t
xiv
1. Problemas de cinemática y dinámica
1.0.1. Cuerpo en el manto de un cono
El cuerpo de masa m = 120 g que se muestra en la

figura, está apoyado sobre un cono de ángulo α = 60o .
Dicho cono gira alrededor de su eje con frecuencia an-
gular ω = 10 rev/min. El hilo que sostiene el cuerpo
es inextensible, y su longitud es L = 50 cm. Calcular:
1. La velocidad tangencial del cuerpo.
2. La tensión del hilo.
3. La velocidad angular necesaria para reducir la

reacción del plano a cero.
Solución:
Primero que nada, conviene re-expresar los datos del problema en unidades compatibles para que, por
ejemplo, los cálculos de fuerza resulten en unidades de Newtons (N = kg m/s2 ). Entonces:
m = 120 g = 0,12 kg , L = 50 cm = 0,5 m ;

−1
ω = 10 rev/min = 1,047 s , α = 60o = π/3 .
Dado que el cuerpo no está deslizando “hacia arriba” o “hacia abajo” sobre el manto del cono, y utilizando
la suposición implı́cita de que además se mantiene sobre dicha superficie (no se “despega”), se puede visualizar
el problema como un movimiento circular en un plano, ubicado a la altura del cono donde está apoyado el
cuerpo, lo cual está esquematizado en el siguiente gráfico.
Para que exista un movimiento circular de este tipo, el cuerpo debe estar sujeto a una aceleración centrı́peta
~ac , apuntando hacia el interior del cono, y ubicada en el mismo plano donde está efectuando dicho movimiento
circular, es decir, el plano de la imagen. Ahora si, pasamos a la resolución:
1. Queremos conocer la velocidad tangencial. Sabemos que para una frecuencia angular constante, y una
distancia fija al centro de giro, la magnitud de la velocidad tangencial queda de la forma vt = ωR, en donde
1
R es la distancia perpendicular entre el eje del cono y el cuerpo. Dado que el dato que conocemos es el
largo del hilo, que coincide con la hipotenusa del triángulo que se ve en la figura, utilizando trigonometrı́a,
tenemos que R = L sin(α). Reemplazando, se obtiene que:
m
vt = ωL sin(α) = 1,047 s−1 × 0,5 m × 0,866 = 0,45 . (1.1)
s
2. Queremos ahora conocer cuanto vale la tensión (fuerza) que hace el hilo para soportar el cuerpo sin que
este se caiga por efecto de la gravedad. Para poder plantear las ecuaciones, debemos utilizar un “sistema
de coordenadas” que dé sentido a las posibles direcciones donde apuntan dichas fuerzas. Como el cuerpo
se está moviendo, utilizaremos el concento de “diagrama de cuerpo aislado”, en donde, utilizando el
cuerpo de forma instantanea, podemos, también de forma instantanea, poner un sistema de coordenadas,
y escribir las fuerzas. ¿Que queremos decir con instantáneamente? Simplemente intentamos reflejar el
hecho de que en un instante anterior o posterior, el cuerpo no esta más en el lugar donde pusimos el
sistema de coordenadas, lo cual implica que las ecuaciones que vamos a escribir dejan de ser válidas.
El diagrama que usaremos para un determinado in-

stante en este problema, tiene la forma que se ve
en la figura. Elegimos los ejes cartesianos de forma
paralela a las fuerzas T~ de tensión del hilo y N ~ de
normal que aplica la superficie por la tercera ley
de Newton. Podrı́amos en principio, haber elegido
los ejes cartesianos paralelos al peso P~ y la acel-
eración centrı́peta ~ac , y los resultados que se obten-
gan serı́an los mismos, pero dado que se pide la ten-
sión de la cuerda, utilizar estos ejes va a permitir un
cálculo más directo. Es necesario escribir las compo-
nentes d ela aceleración centrı́peta y del peso para
poder escribir la segunda ley de Newton para cada
eje cartesiano.
Nuevamente, utilizando trigonometrı́a, tenemos que:
P~ = P xî + P y ĵ = −P sin(α)î − P cos(α)ĵ ,

~ac = axc î + ayc ĵ = −ac cos(α)î + ac sin(α)ĵ .
Aquı́, P = mg y ac = vr2 /R = ω 2 R. Con todo esto, estamos en condiciones de escribir el par de ecuaciones
P~
correspondientes a los ejes x e y respectivamente, utilizando la segunda ley de Newton F = m~a:

N − m g sin(α) = −m ac sin(α) ,
(1.2)
T − m g cos(α) = m ac cos(α) .
De la ecuación para y, la segunda, podemos despejar T sin problemas, y ya conocemos todos los datos
como para realizar el cálculo. Ojo, hay que reemplazar ac = ω 2 R y R = L sin(α):
T = m(g + ω 2 L sin(α)) cos(α)

0,12 kg × 9,81 + (1,047)2 × 0,5 × 0,866 m/s2 × 0,5

=
= 0,64 N .
3. Queremos ahora conocer con que frecuencia angular ω deberı́a rotar el cono para que la reacción del plano
~ = 0. Primero, ¿que significa que esto suceda? Dado que la normal al plano es la reacción que
satisfaga N
se produce por el efecto de apoyar el cuerpo sobre el cono, que la normal sea cero significa lógicamente que
2
el cono no esta más sobre la superficie. Sin embargo, si el cono estuviera “muy lejos” de la superficie, las
ecuaciones que planteamos en el punto anterior dejan de ser válidas. Solamente podemos predecir lo que
sucede cuando la velocidad es tal que N ~ = 0 pero el cuerpo sigue muy cerca del cono, tal que podamos
seguir utilizando α como ángulo para descomponer las fuerzas. En esta aproximación, y utilizando lo que
queremos, que N = 0, podemos utilizar la ecuación para el eje x del punto anterior. Ojo, nuevamente
reemplazar P = mg, ac = ω 2 L sin(α):
r
g g
− mg sin(α) = −m ω 2 L sin(α) cos(α) ⇒ ω 2 = ⇒ ω=± . (1.3)
L sin(α) L sin(α)
Los dos signos que aparecen para ω solamente indican que si el cono rota en sentido horario o antihorario
únicamente. Reemplazando los datos, se obtiene que:
s
9,81 m/s2
ω= = 4,75 s−1 . (1.4)
0,5 m × 0,866
3
2. Problemas de electrostática básica
2.1. Comparación: fuerza eléctrica vs. fuerza gravitatoria
1. Calcular el cociente q/m entre la carga y la masa de dos partı́culas idénticas que se repelen electrostática-
mente con la misma fuerza con que se atraen gravitatoriamente. Comparar el valor hallado con el cociente
e/me para el electrón.
2. Calcular la fuerza gravitatoria entre dos esferas de cobre de un centimetro de diámetro, separadas una
distancia de un metro. Si se retirara de cada esfera un eletrón por átomo, calcular la fuerza de repulsión
electrostática entre ambas. Datos: ρCu = 9 g/cm3 , ACu = 63,5 g/mol.
Solución:
1. Tenemos que la fuerza electrostática obedece la ley

de Coulomb Fe = k qr1 q22 , mientras que la fuerza
de gravedad puede escribirse según la expresión de
Newton Fg = G mr1 m 2
2
. Si asumimos que ambas
tienen la misma magnitud, la esfera sobre la que
dibujamos las fuerzas estará en equilibrio.
Esto se esto puede expresarse de la forma:

X
F = Fe − Fg = 0 ⇒ Fe = Fg . (2.1)
Como ambos cuerpos tienen las mismas propiedades, tenemos que q1 = q2 y m1 = m2 . Si reemplazamos
las fórmulas de las fuerzas, encontramos que:
v
2
6,7 × 10−11 Nkgm2
r u
q G u C
= = t
N m 2 = 8,628 × 10−11 . (2.2)
m k 9
p × 10 C 2 kg
Comparandolo con el mismo valor para el electrón, tenemos:
e 1,6 × 10−19 C C
= = 1,756 × 1011 . (2.3)
me 9,11 × 10−31 kg
La diferencia de veintidos órdenes de magnitud demuestra lo débil de la fuerza gravitatoria respecto de
la fuerza electrostática.
2. Vamos a asumir ahora que las cargas del punto anterior estan repartidas en esferas con volumen (i.e, no
son puntuales). Les asignamos un diámetro d = 1 cm = 0,01 m. Sabemos además que son de cobre, con
lo que podemos inmediatamente calcular su masa:
πd3
m = ρCu V = ρCu . (2.4)
6
Ambas esferas tienen la misma masa, con o que puedo escribir la magnitud de la fuerza de atracción
gravitatoria entre ellas de la forma:
2 2
m2 π 2 d6 2
π (0,01 m)6

−11 N m kg
Fg = G 2 = Gρ2Cu = 6,7 × 10 9000 2 = 1,487 × 10−15 N .
r 36r2 kg 2 m 36(1 m)2
Ahora, en la situación en la que están las esferas, son neutras. Si sacaramos un electrón de cada átomo,
quedarı́an cargadas positivamente. Calculemos la carga neta. La cantidad de átomos contenida debe ser
proporcional al número de moles en la masa de cobre de la esfera e inversamente proporcional a la cantidad
de materia contenida en dicho mol. Matemáticamente esto es:
NA m ρCu NA πd3
n= = = 4,467 × 1022 . (2.5)
ACu 6ACu
4
De esta forma, la carga que tiene cada esfera será q = +ne = 7147,2 C. Podemos entonces calcular la
fuerza electrostática entre ambas:
q2 2
9 N m (7147,2 C)
2
Fe = k = 9 × 10 = 4,596 × 1017 N .
r2 C2 (1 m)2
Podemos ver que aunque los cuerpos en cuestión no sean puntuales, la gravedad está muchisimos ordenes
de magnitud por debajo de la fuerza electrostática.
2.2. Cristales y fuerzas electrostáticas

La sal común N aCl es un cristal iónico. Para forma el enlace entre el Cl y el N a, el átomo de sodio cede
un electrón al átomo de cloro. Si en dicho cristal la distancia entre los átomos es de 2,82 Å, calcular la fuerza
de atracción entre los iones.
Solución:
El cristal N aCl conocido como Clouro de Sodio, es la sal común de consumo humano. Al entregar el N a un
electrón al Cl, queda cargado positivamente con carga qN a = +e, mientras que la neutralidad del Cl desaparece
para tener carga qCl = −e. Esto es lo que usualmente se describe como que el N a queda en estado iónico +1 y
el Cl en −1. Entonces la fuerza electrostática entre dichos iones puede escribirse con la ley de Coulomb:
qCl qN a N m2 (−(1, 6 × 10−19 C)2 )

F =k = 9 × 109 = −2,897 × 10−9 N .
d 2 C 2 (2,82 × 10−10 m)2
2.3. Cuadrado de cargas

Cuatro cargas iguales q están situadas en los puntos de coordenadas cartesianas ~x1 = (a, 0, 0), ~x2 = (−a, 0, 0),
~x3 = (0, a, 0) y ~x4 = (0, −a, 0).
1. Determinar la fuerza ejercida sobre la carga situada en el punto ~x3 por el resto de las cargas.
2. Determinar el campo eléctrico creado por esa distribución de cargas en el eje z.
3. Verificar el valor del campo en puntos del eje z alejados del sistema de cargas.
5
Solución:
1. Comenzemos por determinar el lado del cuadra-

do formado por las cargas como vértices. Para ello
hacemos uso del triángulo rectángulo que se forma
con los ejes coordenados, del cual dicho lado es la
hipotenusa h. Entonces podemos aplicar el teorema
de Pitágoras:
h2 = a2 + a2 , (2.6)
con lo que:
√
h= 2a . (2.7)
Podemos entonces conocer la magnitud de la fuerza sobre la partı́cula en ~x3 . La notación que usaremos es
F~13 para la fuerza que siente la partı́cula en ~x3 por efecto de la partı́cula en ~x1 . Ahora, como se ve en la
figura, dicha fuerza puede descomponerse en el plano xy, y la componente x se anulará con la componente
x de la fuerza F~23 que tiene la misma magnitud y dirección, pero sentido opuesto. Entonces, la fuerza
~ 3 será:
sobre la partı́cula en X
y y
F = F13 sin(θ) + F23 sin(θ) + F43 , (2.8)
donde θ puede obtenerse de consideraciones geométricas en el triángulo de antes:
cat. op. a 1
sin(θ) = =√ =√ . (2.9)
hip. 2a 2
De esta forma:
q2 1 q2 1 q2
F = k √ √ +k √ √ +k
( 2a) 2 2 ( 2a)2 2 (2a)2
q2 √
= k 2 1+2 2 .
4a
Si incluimos la dirección del vector a la magnitud encontrada, tenemos que:
q2 √
F~ = k 2 1 + 2 2 ĵ . (2.10)
4a
6
2. Quiero el campo eléctrico que dichas cargas gener-
an sobre puntos del eje z. Vamos a suponer que las
cargas son todas positivas, lo que por convención
sitúa las lineas de campo en dirección externa a las
cargas puntuales. Conviene realizar el tratamiento
por pares de cargas en cada eje. Ahora, como todas
las cargas son iguales, con resolver para un par es
suficiente, extendiendose el resultado al doble por
presencia de las otras cargas. He aquı́ un uso básico
de la superposición en la teorı́a del electromag-
netismo.
Descomponiendo el campo que genera la carga en ~x1 en sus componentes en el plano xz, haciendo uso
del ángulo θ que se ve en la figura, tenemos que:
~1
E = |E~ 1 |(− cos(θ)î + sin(θ)k̂)
kq a z
= 2
− î + k̂
d d d
kq
= (−aî + z k̂) .
(a2 + z 2 )3/2
El mismo analisis para el campo generado por la carga en la posición ~x2 da:
~2 = kq
E (aî + z k̂) .
(a2 + z 2 )3/2
Entonces, la suma de este par de cargas queda:
~ 1−2 = E
~1 + E
~2 = 2kq
E z k̂ . (2.11)
(a+ z 2 )3/2
Por simetrı́a, la suma del par de cargas sobre el eje y en el plano yz dará exactamente lo mismo, con lo
que el campo total sobre puntos del eje z resulta:
~ = 4kqz
E k̂ . (2.12)
(a2 + z 2 )3/2
3. Si queremos analizar el valor del campo en puntos muy alejados, queremos verificar que pasa cuando
z >> 0. Podemos calcular en primera aproximación, el lı́mite z → ∞, que queda:
~ 4kqz
lı́m E = lı́m k̂
z→∞ z→∞ (a2 + z 2 )3/2
4kq
= lı́m 3/2 k̂ (2.13)
z→∞ 2 a2
z z2 + 1
4kq
∼ k̂ → 0 .
z2
(2.14)
Vemos que se obtiene el tı́pico decaimiento cuadrático del campo eléctrico de una carga puntual, de valor
total 4q. Esto significa que en el infinito, las dimensiones del sistema de carga son despreciables, pudiendo
asumirse que lo que generó el campo es una carga puntual.
7
2.4. Tubo de rayos catódicos: funcionamiento básico
En la siguiente figura, se observa como un electrón es lanzado con velocidad inicial v0 en la dirección del eje
central ubicado entre dos láminas. Vamos a asumir que las láminas generan un campo eléctrico E ~ perpendicular
a dicho eje.
1. Describir el movimiento del electrón desde el momento en que ingresa a la región entre las láminas.
2. Si el campo eléctrico es uniforme de intensidad 1,2 × 104 N/C dirigida hacia abajo, y la velocidad inicial
del electrón es 2 × 107 m/s, determinar la distancia por encima del eje a la cual alcanzará la pantalla
fluorescente S.
Solución:
1. Para encontrar la trayectoria del electrón, conviene separar el problema en dos regiones: aquella entre las
láminas, donde el campo neto no es cero, y aquella entre la salida de las láminas y la pantalla fluorescente.
Ubicamos el sistema de coordenadas con el cero donde comienzan las láminas, y el eje x apuntando
en la dirección de la velocidad inicial. Tenemos entonces las condiciones iniciales x(t = 0 s) = 0 m y
v(t = 0 s) = v0 . La fuerza sobre el electrón no es otra que la fuerza de Coulomb, dada por F~ = −eE. ~ De
esta forma, la ecuación que se debe resolver para la trayectoria es:
d2 ~x
me~a = me = F~ = −eE
~. (2.15)
dt2
Solo tenemos movimiento en el plano xy, con lo que las ecuaciones quedan:
d2 x
me = 0
dt2
d2 y
me 2 = −eE .
dt
Esto se integra en el tiempo, dando:
dx
= vx (t = 0 s)
dt
dy e
= − Et + vy (t = 0 s) .
dt me
Usando el dato de velocidad inicial sale que vx (t = 0 s) = v0 y vy (t = 0 s) = 0. Volvemos a integrar en el

tiempo, para obtener:
x(t) = v0 t + x 0
e
y(t) = − Et2 + y0 .
me
Al imponer las condiciones iniciales nuevamente, se obtiene:
8

e 2
~x(t) = v0 tî + − Et ĵ . (2.16)
me
Ahora, si x1 corresponde con la posición donde se acaban las láminas, podemos por la ecuación en x
conocer cual es el tiempo t1 para el cual llega, i.e.
x1
t1 = . (2.17)
v0
Con este tiempo podemos conocer la altura sobre el eje al tiempo t1 cuando sale de las placas, dada por:
2
e e x1
y1 = − Et2 = − E , (2.18)
me 1 me v0
y las velocidades en ambas direcciones:
vx1 = v0 (2.19)
e e x1
vy1 = − Et1 = − E . (2.20)
me me v0
Lo conveniente entonces, ahora que el electrón llegó a la salida de las placas y tenemos toda la información
sobre la altura y la velocidad con la que esto sucede, es cambiar el sistema de coordenadas, corriendo el
cero del eje x justo al punto donde terminan las placas, y volviendo a poner t = 0 s en este punto. Esto
se traduce en la condición inicial de la trayectoria en x, siendo ahora x(t = 0 s) = x1 = 0 como antes.
Entonces, lo que resta el resolver un problema simple de movimiento rectilineo uniforme, pues ahora:
F~ = 0 . (2.21)
Entonces, las ecuaciones de la trayectoria son simples:
x(t) = at + b ; y(t) = ct + d . (2.22)

Las condiciones iniciales nos dicen que:
x(t) = vx1 t + x1 = v0 t (2.23)

2
e x1 e x1
y(t) = vy1 t + y1 = − E t− E (2.24)
me v0 me v0
(2.25)
Si el electrón llega a la pantalla a distancia x2 respecto de las placas, podemos obtener de la ecuación
para t2 , i.e., el tiempo respecto de la salida de las láminas. Se tiene que:
x2
,
t2 = (2.26)
v0
y con esto podemos dar la altura sobre el eje y a la cual el electrón golpeará la pantalla fluorescente:
2
e x1 e x1
y2 = − E t2 − E
me v0 me v0
2
e x1 x2 e x1
= − E 2 − E
me v0 me v0
eE
x1 x2 + x21 .

= −
me vo2
Entonces, a modo de conclusión, la presencia de la fuerza eléctrica producida por el campo en las láminas,
produce una trayectoria parabólica para el electrón, hasta el momento en que este sale de las placas, donde
la fuerza neta sobre el mismo es nula. Pero en este punto la velocidad no es nula, con lo que el electrón
describirá un movimiento rectilineo uniforme hasta chocar con la pantalla fluorescente.
9
2. Queremos ahora dar numéricamente la altura a la cual golpea la pantalla. Esto significa, aplicar la fórmula
2.27 con los datos numéricos dados. Tenemos que:
E = −1,2 × 104 N/C ,

v0 = 2 × 107 m/s ,
me = 9,109 × 10−31 kg ,
e = 1,602 × 10−19 C ,
x1 = 4 cm = 0,04 m ,
x2 = 12 cm = 0,12 m ,
obteniendose finalmente con todo esto, que:
y2 = 3,4 cm . (2.27)
2.5. Fuerzas moleculares: Cloruro de Cesio
En el compuesto CsCl o Cloruro de Cesio, los áto-

mos de Cs se encuentran en los ángulos de un cubo,
con un átomo de Cl en el centro de él. La longitud de
la arista del cubo es de 0,4 nm. A los átomos de Cs les
falta un electrón, mientras que el átomo de Cl tiene un
electrón en exceso.
1. Calcular la intensidad de la fuerza eléctrica neta que actúa sobre el átomo de Cl.
2. Si el átomo de Cs marcado con una flecha en la figura no estuviera presente (defecto del cristal), calcular
ahora la fuerza eléctrica neta sobre el átomo de Cl.
Solución:
1. Dado que todas las cargas en los vértices corresponden a átomos Cs, cuyo exceso es q = +e, y además
por la simetrı́a, todas son equidistantes al átomo de Cl, cuya carga es q = −e, las magnitudes de las
fuerzas son siempre las mismas. Además, por cada carga en un vértice que genere una fuerza sobre el Cl,
existe el vértice opuesto que genera una fuerza en sentido contrario. La fuerza neta es finalmente cero.
2. Si falta la carga del Cs señalada con flecha en la

figura, todas las demás cargas siguen teniendo una
que es diagonalmente opuesta, con lo que la fuerza
neta se puede calcular poniendo un sistema de co-
ordenadas centrado en el Cl. La figura muestra co-
mo la posición del Cs respecto del Cl se escribe por
medio del vector ~xCs , siendo todas sus componentes
de largo a/2 pero en direcciones distintas, es decir:
a a a
~xCs = , ,− . (2.28)
2 2 2
La fuerza de Coulomb escrita de forma vectorial está dirigida en la dirección del vector posición del Cs,
con lo que:
10
|qCs qCl |
F~ = k x̂Cs
|~xCs |2
e2
= k ~xCs
|~xCs |3
ke2 d
= ( 1, 1, −1) .
3d2 3/2 2

4
Con esto finalmente se llega a:
4ke2
F~ = 3/2 2 (î + ĵ − k̂) . (2.29)
3 d
La magnitud evaluada numéricamente resulta F = 3,324 × 10−27 N .
2.6. Hilo cargado uniformemente

Un hilo rectilineo de longitud 2L está cargado uniformemente con una densidad de carga λ. Calcular el
campo eléctrico en cualquier punto P del plano perpendicular al hilo y que pasa por su centro. Estudiar además
el valor del campo si L → ∞. En el caso del hilo de longitud finita, calcular el valor del campo en puntos muy
alejados del hilo e interpretar el resultado.
11
Solución:
Elegimos un sistema de coordenadas de forma que

el eje z sea paralelo al hilo, y z = 0 esté a una dis-
tancia L del extremo del mismo. Queremos poder dar
el valor del campo eléctrico en los puntos del plano xy
que corta al hilo en z = 0. Haremos uso de la simetrı́a
cilı́ndrica del problema, poniendo por consiguiente en
dicho plano, coordenadas polares (ρ, φ). Para estudiar
el efecto de la distribución, elegimos una pequeña can-
tidad de carga dq = λdz sobre el hilo, y estudiamos
el campo que genera en uno de los puntos del plano
xy. Como se ve en la figura, es posible descomponer el
campo en una componente paralela al plano xy y otra
perpendicular al mismo:
dE⊥ = dE sin(θ) (2.30)

dEk = dE cos(θ) (2.31)
Reescribamos la primera componente, perpendicular al plano:

z z
dE⊥ = dE sin(θ) = dE = dE p . (2.32)
r z + ρ2
2
Ahora, por la ley de Coulomb, la carga dq = λdz genera una intensidad de campo a la distancia r dado por:
dq λdz
dE = k =k 2 , (2.33)
r2 z + ρ2
con lo que:
zdz
dE⊥ = kλ . (2.34)
(z 2 + ρ2 )3/2
Sumar todas las contribuciones de la linea a la componente perpendicular es equivalente a integrar:
Z Z L
zdz
E⊥ = dE⊥ = kλ = 0, (2.35)
−L (z 2 + ρ2 )3/2
porque el integrando es impar. Puede realizarse la cuenta de todas maneras y verificar que esto sucede. Una
forma más rápida de obtener este resultado es haciendo uso de la simetrı́a de cargas. Esto quiere decir, para el
dq que se encuentra en la posición z, siempre existe un dq en −z que genera la misma componente perpendicular
al plano pero con sentido opuesto, cancelandose entre si.
Vamos ahora a calcular la componente paralela al plano. La intensidad del campo es la misma, con lo que:
dz ρ dz
dEk = dE cos(θ) = kλ = kλρ 2 . (2.36)
r2 r (z + ρ2 )3/2
Nuevamente, sumar todas las contribuciones de carga de la linea al campo, corresponde a integrar sobre la
misma:
Z Z L
dz
Ek = dEk = kλρ . (2.37)
−L (z 2 + ρ2 )3/2
Realizar esta integral resulta en:
L
z 2kλL
Ek = kλρ 2 2 2 1/2
= . (2.38)
ρ (z + ρ )
−L ρ(L + ρ2 )1/2
2
Finalmente, el campo eléctrico escrito vectorialmente queda de la forma:
12
~ = 2kλL
E ρ̂ . (2.39)
ρ(L2 + ρ2 )1/2
Vamos a estudiar ahora los casos lı́mites de este campo. Si el hilo fuera infinito, tenemos que:
2kλL −−−−−→ 2kλ
L→∞ , (2.40)
ρ(L2 + ρ2 )1/2 ρ
es decir se vuelve independiente de la longitud del hilo y decae radialmente con potencia uno.
Por otro lado, si el hilo es finito y estudiamos el comportamiento lejos del mismo, se puede realizar la
aproximación para x << 1 en forma de serie de Taylor:
p x2
1 + x2 ∼ 1 + − ··· (2.41)
2
De esta forma:
2kλL 2kλL
=
ρ(L2 + ρ2 )1/2
r 2
ρ2 1 + Lρ
−−−−−→ 2kλL
ρ >> L
L2
ρ2 1 + 2ρ 2
4kλL
∼ .
(2ρ2 + L2 )
Puede verse la tendencia al decaimiento cuadrático esperado para una carga puntual, que es como puede
tratarse el hilo cuando se lo observa desde una distancia apreciablemente mayor a la longitud del mismo.
13
2.7. Disco cargado uniformemente
Un disco conductor de radio a está uniformemente cargado (σ = Q/πa2 ). Calcular el campo eléctrico en
los puntos del eje perpendicular que pasa por el centro del disco. Estudiar el valor del campo para puntos muy
próximos y muy alejados del disco.
Solución:

el eje z pasa por el centro del disco, situando al mismo
en el plano xy. Queremos poder dar el campo eléctri-
co sobre puntos del eje z. Para estudiar el efecto de
la distribución de cargas, elegimos una pequeña canti-
dad dq = σda sobre el disco, y estudiamos el campo
que genera sobre un punto del eje z. Como se ve en
la figura, es posible descomponer el campo dE ~ en las
direcciones paralela al eje z y perpendicular al mismo,
o sea sobre el plano xy.
De nuevo, se tiene que dEk = dE cos(θ), y dE⊥ = dE sin(θ), y es posible escribir las funciones trigonométricas
haciendo uso de su definición y la geometrı́a de la figura:
ρ ρ z z
sin(θ) = =p , cos(θ) = = p . (2.42)
r 2
ρ +z 2 r ρ + z2
2
Ahora, la magnitud del campo generado es:

dq σρ dρ dϕ ρ dρ dϕ
dE = k =k = kσ 2 , (2.43)
r2 r2 ρ + z2
donde φ es el ángulo azimutal sobre el disco. De esta forma, la componente paralela al eje z queda:
zρ dρ dϕ
dEk = kσ . (2.44)
(ρ2 + z 2 )3/2
Para conocer la contribución del disco en su totalidad se debe integrar en el mismo:
Z
Ek = dEk
Z
zρ dρ dϕ
= kσ 2 + z 2 )3/2
disco (ρ
Z 2π Z a
ρ dρ dϕ
= kσz
0 (ρ + z 2 )3/2
2
Z a0
ρ dρ
= 2πkσz 2 2 3/2
0 (ρ + z )
!a
1
= 2πkσz − p

2
ρ + z 0 2

z
= 2πkσ 1 − √ .
z 2 + a2
Nuevamente se pueden usar argumentos de simetrı́a para demostrar que la componente transversal al eje
z es cero, pues otro dq en el lado opuesto del disco genera la componente de misma magnitud dirección pero
sentido contrario. Es útil también realizar la integral correspondiente y verificar que es cero. De esta forma el
campo total sobre puntos del eje z es:
14

~ z
E = 2πkσ 1 − √ k̂ . (2.45)
z 2 + a2
Vamos a chequear lo lı́mites del campo encontrado. Si nos acercamos al disco, es equivalente a pedir que
z ∼ 0. Esto se traduceen que:
E(z ∼ 0) ∼ 2πkσ = cte. (2.46)

Como se verá más adelante, que sea constante es una conclusión natural para el campo generado por una
superficie infinita de distribución uniforme de carga. Dado que muy cerca del centro del disco, este se extiende
hacia todas las direcciones simulando ser infinito, obtenemos este resultado.
Por otro lado, al alejarse del disco se está pidiendo que z >> a, lo cual lleva a que:

z
E(z >> a) ∼ 2πkσ 1 − √ = 0, (2.47)
z2
lo cual es consistente como siempre con la idea de que lejos de la distribución del disco, este parece una carga
puntual, decayendo el campo a cero como se espera.
Nota: es interesante notar que, si d es la distancia a la distribución de cargas, el campo de una carga
puntual decae como d−2 , el campo de un hilo infinito decae como d−1 mientras que el campo de una
superficie infinita decae como d0 , es decir, no decae. Serı́a entonces esperable una relación del tipo:
si n es la dimensión de la distribución de cargas infinita (cantidad de parámetros necesarios para su
descripción), entonces el decaimiento del campo será d−(n−2) . Si nos abstraemos aún más, y decimos
que la distribución de cargas está inmersa en un espacio de dimensión m con m > n, entonces el campo
que generan decae de la forma d−(n−m+1) . En el caso de electrostática que estamos tratando, m = 3.,
pues el espacio es el clásico euclideo R3 . Habrı́a que chequear la demostrabilidad de esta propiedad de
los campos vectoriales con fuentes desde un punto de vista matemático.
15
2.8. Anillo cargado uniformemente
Un anillo de radio a tiene repartida uniformemente en su contorno una carga Q. Calcular el campo eléctrico
en un punto P situado en el eje del anillo a una distancia z de su centro. ¿Por qué se anula el campo en el
centro del anillo? ¿Cuánto vale el campo en puntos muy alejados del anillo?
Solución:

el eje z pasa por el centro del anillo, situando al mismo
en el plano xy. Queremos poder dar el campo eléctri-
co sobre puntos del eje z. Para estudiar el efecto de
la distribución de cargas, elegimos una pequeña canti-
dad dq = λdl sobre el anillo, y estudiamos el campo
que genera sobre un punto del eje z. Como se ve en
la figura, es posible descomponer el campo dE ~ en las
direcciones paralela al eje z y perpendicular al mismo,
o sea sobre el plano xy, como hicimos en el caso del
disco cargado.
dEk = dE cos(θ) , (2.48)
dE⊥ = dE sin(θ) . (2.49)
Como antes, también se tiene que cos(θ) = z/d y sin(θ) = a/d, y el diferencial de carga que produce el
diferencial de campo electrico es dq = λdl = λadϕ, quedando:
dq kλa
dE = k = 2 dϕ . (2.50)
d2 d
Nuevamente por simetrı́a se aprecia que la componente perpendicular al eje z se anula por la existencia de
una carga diametralmente opuesta sobre el anillo que la cancela. De esta forma:
kλa z z
dEk = dϕ = kλa 2 dϕ . (2.51)
d2 d (a + z 2 )3/2
Integrar esto en todo el anillo resulta simplemente en:
~k = 2πλaz
E k̂ . (2.52)
(z 2 + a2 )3/2
Notar que la carga total del anillo es Q = 2πaλ, con lo que finalmente:
~k = kQz
E k̂ . (2.53)
(z 2 + a2 )3/2
Ahora, los casos lı́mites. Como era de esperarse, en el centro del anillo con z = 0 el campo de anula debido
a que todas las componentes del campo electrico apuntan allı́ en la misma dirección por la uniformidad en la
distribución de Q sobre el anillo. Por otro lado, si z >> a,
kQz kQ
Ek → 3/2 ∼ z 2 , (2.54)
a2
z3 z2 +1
es decir, decae cuadráticamente como una carga puntual.
16
2.9. Esfera cargada uniformemente
Una distribución volumétrica uniforme de carga ρ se encuentra dentro de una esfera de radio a. En el espacio
que rodea la esfera no hay cargas. Calcular el campo eléctrico en función de la distancia r al centro de la esfera.
Analizar la continuidad en r = a.
Solución:
Para resolver un problema de este tipo, es útil hac-

er uso de la ley de Gauss, la cual relaciona el valor del
campo eléctrico sobre una superficie cerrada con las
cargas electrostáticas que dicha superficie encierra:
Z I
ρenc (~x)d3 ~x = 0 ~ · dS
E ~, (2.55)
V ∂V
donde ∂V es el borde del volumen V limitado por la

superficie gaussiana que sea de interés para el proble-
ma a resolver, y dS ~ = n̂dS es el vector que apunta en
la dirección normal saliente del volumen V . La elección
de dicha superficie suele ser intrı́nseca del problema y
sus simetrı́as. En este caso, se eligen las superficies S1
y S2 que se muestran en la figura.
Vamos a aplicar la ley de Gauss primero para la esfera exterior S2 . La carga encerrada es la total, dada por
q = 43 πa3 ρ. La otra integral, que es el flujo del campo eléctrico, consiste en evaluar el campo en la superficie,
calcular el producto correspondiente y escribir el diferencial de superficie correcto. Usamos la simetrı́a esférica
que implica que el campo no puede depender de los ángulos polar y azimutal porque la distribución de cargas
no depende de ellos, y la geometrı́a tampoco. Finalmente igualamos a la carga encerrada.
I I
0 ~ · dS
E ~2 = 0 ~ 2 ) · r̂dS
E(r
S2 S2
I
= 0 E(r2 )dS
S2
Z 2π Z π
= 0 E(r2 )r22 sin(θ)dθdϕ
0 0
= 0 4πr22 E(r2 )
4 3
= πa ρ ,
3
a partir de lo que despejando se obtiene:
ρ a3
E(r2 ) = . (2.56)
30 r22
Esta es la expresión correcta que contempla los datos que brinda el enunciado. Puede reescribirse la misma
multiplicando y dividiendo por 4π y formando la carga total q con la densidad de carga ρ, de forma que:
1 q
E(r2 ) = , (2.57)
4π0 r22
que no es otra cosa que el campo eléctrico generado por una carga puntual ubicada en el centro de la esfera.
Por otro lado, vamos a calcular el campo en la superficie interna S1 . En este caso, la carga encerrada es:
17
Z Z
ρenc dV = ρ dV
V V
Z r1 Z 2π Z π
= ρ r2 sin(θ)drdθdϕ
0 0 0
4 3
= πr ρ .
3 1
La integral del flujo es análoga pero ahora se evalúa al campo eléctrico en la superficie interna, con lo que:
4 3 ρ
0 4πr12 E(r1 ) = πr ρ −→ E(r1 ) = r1 . (2.58)
3 1 30
Se tiene entonces como expresión general:
( ρ
30 r si r < a
E(r) = ρ a3 (2.59)
30 r 2 si r > a .
Ahora, tomando los diferentes lı́mites acercandose a la superficie r = a tanto desde la izquierda como desde
la derecha, se obtiene que:
ρ
lı́m+ E(r) = a = lı́m− E(r) , (2.60)
r→a 30 r→a
lo cual da sustentabilidad a la continuidad del campo eléctrico sobre la superficie de la esfera. Un esquema del
campo eléctrico en función de la distancia al origen puede verse en la siguiente figura:
ρa/3ε0
Campo electrico E(r) [N/C]
0 a
Distancia radial r [cm]
18
2.10. Hilo cargado uniformemente: Gauss
Considerar un hilo recto muy largo (infinito) con una distribución lineal uniforme de carga λ. Utilizar la ley
de Gauss para calcular el campo eléctrico generado por dicha distribución.
Solución:
Nuevamente queremos calcular el campo de un hi-

lo de carga infinito, pero en este caso utilizando la ley
de Gauss. En este caso, conviene elegir una superficie
gaussiana como la cilı́ndrica mostrada en la figura. Si
este cilindro tiene altura h y radio r, se puede calcular
la carga encerrada haciendo uso de la densidad lineal
de carga λ, de forma que q = λh. Para integrar el flu-
jo en la superficie del cilindro, se debe dividir en sus
distintas componentes, las tapas y el manto:
I Z Z Z
~ · dS
E ~= ~ · dS
E ~1 + ~ · dS
E ~2 + ~ · dS
E ~3 .
S S1 S2 S3
Sin embargo, sobre las tapas, los elementos de volumen son dS ~1−2 = ±k̂dS1−2 , y dado que el campo tiene
simetrı́a cilı́ndrica como consecuencia de la distribución uniforme, solo tiene componente r̂, anulando el producto
interno, i.e.:
Z Z
~ ~
E · dS1−2 = ± EdS1−2 r̂ · k̂ ≡ 0 .
S1−2 S1−2
En el manto, nuevamente por simetrı́a de la distribuciónde cargas, el campo solo puede depender de la
distancia al hilo, con lo que:
Z Z
~ · dS
E ~3 = E dS3 r̂ · r̂
S3 S3
Z h Z 2π
= E(r)rdzdϕ
0 0
= 2π r h E(r)
1
= λh ,
0
con lo que:
~ = λ r̂
E . (2.61)
2π0 r
Dicho valor de campo coincide con el encontrado anteriormente en (2.40).
19
2.11. Cilindro de carga no-uniforme
Un largo cilindro posee una densidad de carga que es proporcional a la distancia r al eje del mismo; ρ(r) = kr,
con k constante. Utilizando la ley de Gauss escribir el campo eléctrico interno.
Solución:
Nuevamente tenemos simetrı́a cilı́ndrica en la dis-

tribuciónde cargas. En este caso tenemos un cilindro
de densidad de carga no uniforme, con lo que conviene
elegir una superficie gaussiana que sea otro cilindro in-
terior como el punteado de la figura. Si este cilindro
tiene altura h y radio r, se puede calcular la carga
encerrada haciendo uso de la densidad volumétrica de
carga ρ, de forma que:
Z Z h Z 2π Z r
ρdV = kr0 r0 dr0 dzdϕ
V 0 0 0
2
= πh k r3 .
3
Para integrar el flujo en la superficie del cilindro,

se debe dividir en sus distintas componentes como en
el problema anterior.
Entonces, debemos despejar el campo de la ecuación:

2
0 2π r h E(r) = πh k r3 , (2.62)
3
con lo que:
~ int = k r2 r̂ .
E (2.63)
30
k 2
Si se evalua el campo en el borde, se obtiene E = 30 R .
Por otro lado, si se usara la ley de Gauss en el exterior, la integral del flujo del campo no cambia su expresión,
pero la carga si, quedando la ecuación:
2
0 2π r h E(r) = πh k R3 , (2.64)
3
y el campo afuera finalmente queda:
3
~ ext = k R r̂ ,
E (2.65)
30 r
que, en primera instancia es el campo de un hilo infinito, y en segunda instancia al evaluarlo en la superficie, se
pega con el campo interno, lo cual asegura la continuidad del campo eléctrico.
20
2.12. Tubo metálico cargado
Un tubo metálico de paredes delgadas y de radio R,tiene en su superficie una carga por unidad de longitud
λ. Encontrar una expresión para el campo eléctrico en función de la distancia r perpendicular al eje del tubo
dentro y fuera. Realizar una gráfica cualitativa del mismo.
Solución:
Conviene resolver el problema nuevamente utilizando la ley de Gauss. En este caso de dispone de un tubo
cargado con densidad de carga lineal λ en su superficie. Si elegimos una superficie cerrada S de cualquier forma
dentro del tubo, la carga encerrada es nula, con lo que el flujo del campo eléctrico es cero. Esto se satisface
~ int = 0.
únicamente si E
Si elegimos una superficie gaussiana con forma de cilindro de radio mayor al del tubo y de altura h, se puede
integrar la carga contenida, que será simplemente q = λh y el flujo del campo eléctrico es análogo al de los
últimos dos problemas, con lo que:
0 2πr hE(r) = λh , (2.66)

y el campo finalmente queda:

~ = 0 si r < R
E λ r̂ (2.67)
2π0 r si r > R .
Un esquema del mismo puede verse en la siguiente figura:
λ/(2πε0R)
Campo electrico E(r) [N/C]
R 0 R
Distancia radial r [cm]
21
2.13. Placa infinita cargada
Una región del espacio con forma de placa infinita de ancho 2L está repleta de carga uniformemente dis-
tribuida con una densidad ρ. Encontrar el campo eléctrico en cada punto del espacio.
Solución:
Nuevamente, la forma más simple de abordar este

problema es utilizando la ley de Gauss. Como la pla-
ca es infinita, al observar el comportamiento el cam-
po afuera, toda componente paralela a la superficie
contará con su contraparte apuntando en la dirección
opuesta, con lo cual el campo eléctrico siempre será
perpendicular a la superficie de la placa. Se utilizará
además un sistema de coordenadas centrado en la pla-
ca, de forma que cada una de sus caras esté ubicada en
−L y en L. En el interior de la placa, conviene elegir
una superficie tipo pastilla ∂V como la que observamos
en la imagen, cuya cara presenta un área A. Es con-
veniente situar una cara de dicha superficie con centro
en el origen del sistema de coordenadas. Entonces, la
carga encerrada se calcula fácilmente como q = Axρ.
La ley de Gauss teniendo en cuenta que el campo es
perpendicular a la superficie, queda:
I Z Z
~ · d~a =
E Ex da + (−Ex )(−da)
∂V ∂Vx+ ∂Vx−
Z
= 2 Ex da
A
= 2AEx ,
en donde ∂Vx− es la cara en el origen, cuya normal apunta en el sentido de las x negativas, mientras que ∂Vx+
es la cara cuya normal apunta en la dirección positiva de las x, ubicada en la posición x. Esto se iguala a la
carga encerrada:
Axρ
2AEx =
0
xρ
Ex = ,
20
siempre que |x| < L.
Para calcular en el exterior, utilicemos la misma superficie tipo pastilla pero con la cara ∂Vx+ ubicada en
una posición x > L. Entonces, la carga encerrada tiene ahora un lı́mite, porque más allá de x = L no hay
cargas. De esta forma, q = ALρ. Las integrales del campo eléctrico en la superficie ∂V es de la misma manera
que antes, de forma que la ley de Gauss dice ahora que:
ALρ
2AEx =
0
Lρ
Ex = ,
20
22
siempre que x > L.
En los valores de x < −L, el módulo del campo eléctrico será el mismo que para x > L, pero su sentido
será el opuesto, dado que estamos asumiendo que las cargas de la placa son positivas (el campo va hacia afuera
de las mismas). De esta manera, escribir el campo en todo punto del espacio requiere notación vectorial:

Lρ
 − 20 î x < −L

~ xρ
|x| < L
E(x) = 20 î (2.68)
 Lρ î

x > L
20
Notar que, como en los casos anteriores, existe continuidad en las caras de la placa. Además, dada su
infinitud en las direcciones z e y, el campo no puede depender de ninguna de estas variables.
23
2.14. Campo eléctrico “Rutherfordiano”
En un trabajo publicado en 1911, Ernest Rutherford señalaba, a fin de hacerse una idea de las fuerzas
necesarias para desviar una partı́cula α (alfa) un gran ángulo: “... consideremos un átomo que contenga una
carga puntual positiva Ze en su centro y que esté rodeado de una distribución de electricidad negativa −Ze
uniformemente distribuida en una esfera de radio R. El campo eléctrico E a una distancia r del centro en un
punto dentro del átomo es:”

Ze 1 r
E= − . (2.69)
4π0 r2 R3
Verificar la ecuación anterior.
Solución:
Por la simetrı́a esférica del problema, es útil en este

caso utilizar la ley de Gauss. Además, se utilizará el
principio de superposición. Sabemos que la carga pun-
tual en el centro del átomo, que simula ser el núcleo
atómico, es:
Ze
En = k , (2.70)
r2
con r la distancia radial a la carga como siempre. Aho-
ra, para calcular el campo de la distribución uniforme
de carga −Ze, escribimos en primera instancia la den-
sidad de carga:
Q Ze 3Ze
ρ= =−4 3 =− . (2.71)
V 3 πR
4πR3
Pensando en una esfera de Gauss de radio r como la punteada en la figura, es posible expresar la ley de
Gauss de la forma:
Z Z
~ 1
Ee · d~a = ρdV (2.72)
∂V 0 V
donde V es la esfera de radio r y ∂V su superficie. Entonces:
Z Z
3Ze
ρdV = − dV
V 4πR3 V
3Ze 4 3
= − πr
4πR3 3
r 3
= −Ze ,
R
Z Z π Z 2π
~ e · d~a
E = Ee (r)r2 sin(θ)dθdϕ
∂V 0 0
= Ee (r)r2 4π ,
e igualando ambas expresiones se obtiene:

Ze r 3 Ze r r
Ee (r)r2 4π = − ⇒ Ee (r) = − = −Zek 3 . (2.73)
0 R 4π0 R3 R
Si se aplica superposición se obtiene:

1 r
E = En + Ee ⇒ E(r) = Zek − 3 . (2.74)
r2 R
24
2.15. Concepto de material conductor
Considerar una carga q colocada en el centro de una cavidad esférica metálica de radio a. Utilizando la ley
de Gauss, encontrar el campo eléctrico en cualquier punto tal que 0 < r < a y en cualquier punto en el interior
del material conductor.
Solución:
Vamos a pensar que luego de colocada la carga en el interior de la cavidad metálica, pasó el tiempo suficiente
como para que el sistema haya llegado al equilibrio electrostático. Sabemos que las caracterı́sticas del material
conductor implican que los portadores de carga pueden desplazarse libremente en su interior, pero que al llegar
a una superficie limitante del mismo, no pueden escapar.
En estas condiciones, es normal suponer que sobre

la cara interna de la cavidad se reunan cargas opuestas
a q, atraidas por la fuerza de Coulomb ejercida por q,
y formando una distribución uniforme σ. Esto debe ser
ası́ porque sino no habrı́a más equilibrio electrostático.
Tenemos entonces que la carga total en la superficie de
la cavidad es −q = 4πR2 σ. Además, podemos elegir
diferentes superficies de Gauss Si y Se para el interior
y exterior de la cavidad respectivamente, donde “ex-
terior” en este caso es correspondiente con “interior al
conductor”. En el interior no existen más cargas encer-
radas por Si que la carga puntual q con lo que:
q
E=k . (2.75)
r2
En el exterior de la cavidad, la carga encerrada es:
Z
qenc = q+ σdSe
Se
Z
= q+ σR2 dΩ
Se
= q + 4πR2 σ
= q − q = 0.
Esto simplemente indica que en el interior del conductor (o fuera de la cavidad esférica) E = 0. Esta es la
condición que electrostáticamente define a un material conductor perfecto.
25
2.16. Potential: carga puntual
Para una carga puntual q, calcular la diferencia de potencial entre los puntos A y B en los siguientes casos:
Solución:
Existen dos formas de resolver este problema, una que llamaré formal y otra informal. Mientras conceptual-
mente se entienda que es lo que se está haciendo, cualquiera de las dos formas es totalmente válida.
La solución formal consiste lógicamente en utilizar la definición de diferencia de potencial entre dos puntos
del espacio, es decir:
ZB
φAB = − ~ · d~l ,
E (2.76)
A
en donde d~l el vector diferencial tangente en cada punto a la linea que va desde A hasta B. Sabemos que el
campo eléctrico generado por una carga puntual es:
~ r) = k q r̂ ,
E(~ (2.77)
r2
pero por ejemplo, en el primer caso, los puntos A y B y la carga q son colineales, con lo cual utilizar un eje,
digamos x, como parámetro, con origen en la carga, es totalmente válido. Entonces, el campo sobre este eje
puede escribirse fácilmente:
~ q
E(x) = k 2 î . (2.78)
x
Por otro lado, como recorremos el camino desde A hacia B sobre la recta x en todos los puntos tenemos
que d~l = dxî. Entonces:
ZB Z2d Z2d
~ · d~l = − q dx
φAB = − E k 2 dxî · î = −kq . (2.79)
x x2
A −d −d
Esta integral no está definida dado que pasa por la singularidad x = 0, por lo cual el área bajo la curva
es “infinita”. Entonces, si ya se tiene el conocimiento suficiente de la teorı́a electrostática, se puede argumentar
que como las lineas equipotenciales de una carga puntual son circulares. El valor del potencial en A no cambia
al moverlo entre la primera y la segunda situación, y dado que el punto B permanece en el mismo lugar, la
diferencia de potencial debe ser la misma. Se puede entonces proceder a realizar el cálculo formal para la segunda
situación. Es interesante pensar en realizar la integral de lı́nea primero entre A y el origen, y luego entre el
origen y B. Es evidente que este camino no es posible, porque nuevamente en algún punto nos encontraremos
con la necesidad de evaluar 1/x|x=0 lo cual está indeterminado matemáticamente.
Conviene entonces utilizar directamente la lı́nea Vectorialmente, podemos escribir los puntos A y B
que une A con B, considerando alguna parametrización de la forma:
adecuada que llamaremos γ(t). Esta lı́nea tiene la par-
ticularidad de que “rodea” la singularidad x = 0, por
lo cual la integral puede hacerse sin problemas.
26
A = (0, d, 0)
B = (2d, 0, 0) .
Entonces, una parametrización γ(t) con t ∈ [0, 1] se puede construir como:
γ(t) = tA + (1 − t)B = (2d t, (1 − t)d, 0) ⇒ d~l = dγ = (2d dt, −d dt, 0) = (2d, −d, 0)dt .
Si escribo el campo de forma conveniente, evaluado sobre γ(t):
~ q n o
E(γ(t)) = k x(t)î + y(t)ĵ + z(t) k̂
r3 (t)
q n o
= k 2d tî + (1 − t)d ĵ
(4d2 t2 + (1 − t)2 d2 )3/2
q n o
= k 2d tî + (1 − t)d ĵ ,
(5d2 t2 − 2d2 t + d2 )3/2
se puede calcular la diferencia de potencial sin mayores problemas:
n o n o
ZB Z1 2d tî + (1 − t)dĵ · 2dî − dĵ
φAB = − ~ · d~l = −kq
E dt
(5d2 t2 − 2d2 t + d2 )3/2
A 0
Z1 Z1
4d2 t + (t − 1)d2 5d2 t − d2
= −kq dt = −kq dt
(5d2 t2 − 2d2 t + d2 )3/2 (5d2 t2 − 2d2 t + d2 )3/2
0 0
Z1 Z4
kq 5t − 1 kq du
= − dt = −
d (5t2 − 2t + 1)3/2 2d u3/2
0 1
4
kq 1 1 kq 1
= − = −1
2d − 32 + 1 u3/2−1 1

d 41/2
kq
= − .
2d
en donde definı́ u = 5t2 − 2t + 1. Este es el resultado de la diferencia de potencial entre los puntos A y B en
la segunda situación, y puede extenderse a la primera situación utilizando las lı́neas equipotenciales de la carga
puntual.
La solución informal, consiste en utilizar la definición de diferencia de potencial entre dos puntos para
calcular directamente el potencial de una carga puntual en cualquier punto. Sabemos por la teorı́a electrostática
que este potencial se conoce asignandole el valor φ(~r) → 0 si |~r| → ∞, y queda:
q
φ(~r) = k , (2.80)
r
donde r es la distancia entre la carga y el punto de observación. Entonces, el potencial en los puntos A y B
queda:
27
q q
φA = k
; φB = k , (2.81)
d 2d
y como son lı́neas equipotenciales, tanto en la primera situación como en la segunda, se tiene que:

kq 1 kq
φAB = φB − φA = −1 =− , (2.82)
d 2 2d
obteniendose el mismo resultado para las dos situaciones por estar los puntos A y B sobre lı́neas equipotenciales.
2.17. Potencial: cargas colineales
Para una configuración de cargas como la que se

muestra en la figura, mostrar que el potencial en el eje
r, V (r), está dado por:

q 2qa
V (r) ≈ k + 2 ,
r r
siempre que r >> a. Luego, calcular el potencial para

puntos que no pertenezcan al eje, pero que sigan satis-
faciendo la condición |~r| >> a (aproximación dipolar).
Solución:
La clave está en utilizar la superposición. La carga que está en el origen de la coordenada r tiene un potencial
dado por:
q
V0 (r) = k . (2.83)
r
Si desplazamos esta carga hacia r = a o r = −a, y además ponemos el signo correspondiente para la carga
en r = −a, estos potenciales pueden escribirse simplemente como:
q q
; V−a (r) = −k
V+a (r) = k . (2.84)
r−a r+a
Notar que estamos tratando r como una coordenada unidimensional en el eje que une las cargas, i.e., no
representa la magnitud de un vector de tres dimensiones ~r como usualmente. La superposición dicta que el
potencial total sobre el eje queda:
q q q
V (r) = V0 (r) + V−a (r) + V+a (r) = k − k +k
r r+a r−a

q r + a − (r − a) q 2a
= k + kq = k + kq 2
r r 2 − a2 r r − a2
q 2qa
≈ k +k 2 ,
r r
en donde en el último paso se utilizó que r >> a.
Si ahora queremos evaluar el potencial en puntos fuera del eje, es necesario cambiar el nombre del mismo
por algo más adecuado. Si las cargas están ubicadas colinealmente en el eje z, el sistema tiene simetrı́a de
rotación con este eje, o simetrı́a azimutal. Esperamos entonces que el potencial resultante sea una función
V (r, θ). Siguiendo la idea anterior, el potencial V0 no cambia, solo que ahora r si representa la magnitud de |~r|.
En los otros dos casos, conviene notar que el desplazamiento es solo sobre el eje z, con lo que:
28
q q
V+a (~r) = k p ; V−a (~r) = −k p .
2 2
x + y + (z − a) 2 x + y + (z + a)2
2 2
Trabajo una de estas expresiones, la otra es análoga. Tenemos que:
q q
V+a (~r) = kp = kp
x2 y2
+ + (z − a)2 x + y + z 2 − 2az + a2
2 2
q q
≈ k√ = kp
r2 − 2za r2 − 2ar cos(θ)
en donde se utilizó la definición de la coordenada z en esféricas. El otro potencial queda:

q
V−a (~r) = −k p .
r2 + 2ar cos(θ)
La superposición es igual que en el caso anterior:
q q q
V (~r) = V0 (r) + V−a (~r) + V+a (~r) = k − k p + kp
r 2
r + 2ar cos(θ) 2
r − 2ar cos(θ)
(p p )
q r2 + 2ar cos(θ) − r2 + 2ar cos(θ)
= k + kq p
r r4 − 4a2 r2 cos2 (θ)
(p p )
q q r2 + 2ar cos(θ) − r2 + 2ar cos(θ)
= k +k p
r r r2 − 4a2 cos2 (θ)
q q p 2
n p o
≈ k +k 2 r + 2ar cos(θ) − r2 − 2ar cos(θ) .
r r
Notar como la simetrı́a azimutal aparece por si sola, y la dependendia en el ángulo ingresa a través de
z = r cos(θ). La última expresión podrı́a trabajarse un poco más, pero opté por dejarla ası́ para poder comparar
con el potencial sobre el eje que encontramos en la primera parte del problema.
2.18. Potencial: anillo cargado uniformemente

Un anillo de radio a posee una carga q uniformemente distribuida. Encontrar el potencial eléctrico en un
punto del eje normal al plano del anillo y que pasa por el centro del mismo. Obtener el campo eléctrico en este
eje a partir del resultado obtenido para dicho potencial.
Solución:
La figura es en este caso igual al problema ya resuelto del anillo con carga uniforme. Como ahora nos interesa
el potencial en el eje, el aporte de una pequeña cantidad de carga diferencial dq = λdl puede pensarse como
puntual, es decir:
dq dl
dφ = k = kλ √ , (2.85)
d a + z2
2
con lo cual el aporte de todo el anillo corresponde a sumar todos los aportes diferenciales (integrar) a lo largo
del anillo:
I
1 2πa
φ(z) = kλ √ dl = kλ √ . (2.86)
2
a +z 2 a2 + z 2
Por otro lado, sabemos que E(z) = − dφ
dz , con lo que:

2πakλ 1 2πλaz
E(z) = − 2 2 3/2
2z − = 2 , (2.87)
(a + z ) 2 (a + z 2 )3/2
que no es otra cosa que lo que ya encontramos cuando resolvimos el campo eléctrico en el eje del anillo (2.52).
29
2.19. Potencial: superficies equipotenciales en conductores
Sean dos esferas metálicas de radios R1 = 1 cm y R2 = 2 cm alejadas una de la otra inicialmente. Se carga
la esfera más pequeña con q0 = 2 × 10−7 C y luego se conectan ambas esferas con un alambre también metálico,
pero delgado. Calcular la carga final en cada esfera, la densidad de carga y el potencial que adquieren luego de
conectarlas.
Solución:
Mientras tenemos la esfera 1 cargada, el sistema está en equilibrio electrostático, y como la esfera es metálica
(conductora), toda la carga excedente se distribuirá uniformemente en su superficie. Luego, al conectar las esferas
utilizando un hilo conductor delgado, parte de esta carga se desplazará por el mismo hacia la esfera 2. Esto no
puede suceder por tiempo indeterminado hasta vaciar la esfera 1, sino que en algún momento las superficies de
las esferas y el conductor estarán a un mismo potential, impidiendo el movimiento de cargas, es decir, se llega
al equilibrio electrostático. En esta situación, y utilizando la información de que el hilo es delgado, podemos
aproximar los potenciales de ambas esferas como iguales, φ1 = φ2 . Para la superficie de la esfera, el potential es
como el de una carga puntual ubicada en su centro, con lo que:
q1 q2
φ1 = φ2 ⇒ k =k . (2.88)
R1 R2
Además, al llegar al equilibrio, se debe cumplir que q0 = q1 + q2 , dado que no existe una fuente de carga
que haga que la cantidad neta que existe en el sistema cambie (conservación de la carga). Tenemos entonces un
sistema de dos ecuaciones con dos incognitas:
q1 q2
R1 − R2 = 0
(2.89)
q1 + q2 − q0 = 0
−1
R2
Resolviendolo se obtiene que q1 = 1 + R 1
q0 = 6,66 × 10−8 C, y q2 = 1,33 × 10−7 C. Además:
q1 q1 C
σ1 = = = 5,29 × 10−9
S1 4πR12 cm2
q2 q2 C
σ2 = = = 2,64 × 10−10
S2 4πR22 cm2
son las correspondientes densidades de carga, y el potencial de ambas es el mismo:

q1
φ1 = φ 2 = φ = k = 59940 V . (2.90)
R1
2.20. Par de placas infinitas

Se disponen dos planos conductores muy grandes de forma paralela separados una distancia 2d. Uno de ellos
tiene densidad superficial de carga +σ y el otro con −σ en sus caras internas. Calcular el campo eléctrico en
todo el espacio y la diferencia de potencial entre las placas.
Solución:
Vamos primero a utilizar la ley de Gauss para re- problema, el campo eléctrico siempre será perpendic-
solver el campo generado por una sola placa ubicada en ular al plano cargado. De esta forma, sobre el manto
el origen. Para esto, voy a utilizar una superficie tipo de la pastilla no hay flujo de campo eléctrico, pero si
pastilla como se muestra en la figura, con las caras de sobre las caras.
área A equidistantes al plano cargado. La carga encer-
rada queda:
qenc = σA . (2.91)
Por otro lado, como la placa es muy grande, to-
da componente del campo eléctrico que sea paralela a
la placa tiene una contribución de la misma magnitud
en sentido contrario, por lo cual, dada la simetrı́a del
30
Si la placa tiene densidad de carga +σ, en la cara ubicada en los x negativos el campo eléctrico es saliente
hacia la dirección de los x negativos al igual que la normal a dicha cara. Lo mismo sucede en el caso de las x
positivas, con lo que:
Z Z
ΦE~ = ~ · dS
E ~x<0 + E~ · dS
~x>0
Sx<0 Sx>0
Z Z
= EdS(−î) · (−î) + EdS î · î
Sx<0 Sx>0
= 2EA .
Recordar que todos los puntos de las caras de la pastilla estan a la misma distancia del plano cargado, con
lo cual el campo eléctrico vale lo mismo en todos ellos y por eso es posible sacarlo de la integral.
qenc
Utilizando la ley de Gauss ΦE~ = 0 , obtenemos que:
σA σ
2EA = ⇒ E= . (2.92)
0 20
Este es el módulo del campo, el cual vectorialmente puede escribirse como:
(
~ σ x − 2σ0 î x < 0
E= î . = σ (2.93)
20 |x| 20 î x>0
Si ahora volvemos al caso de las dos placas, una está ubicada en la posición x = d mientras que la otra está
en x = −d, con lo que el campo de cada una puede escribirse como:
(
~+ σ x−d − 2σ0 î x<d
E = î = σ
20 |x − d| 20 î x>d
(
σ
~− σ x+d 20 î x < −d
E = − î =
20 |x + d| − 2σ0 î x > −d
Con esto es muy fácil darse cuenta que E ~ T = 0 si |x| > d, es decir, solo existe
~ T = − σ î si −d < x < d y E
0
campo eléctrico en la región entre las placas. Con esta información podemos calcular la diferencia de potencial
entre las mismas:
31
ZB Za
~ · d~l = − σ
∆Φ = − E − dx î · î
0
A −a
Za
σ σ
= dx = 2a .
0 0
−a
Todas estas cantidades son importantes para el estudio de los condensadores, que son elementos de un
circuito eléctrico conformados justamente por un par de placas paralelas con densidades de carga opuestas.
32
3. Problemas de magnetostática básica
3.1. Protones en movimiento en un campo magnético
Un haz de protones se desplaza a 3,2 × 105 m/s a través de un campo magnético de magnitud 0,2 T , dirigido
a lo largo del eje z positivo. El vector velocidad de cada protón está en el plano xz formando un ángulo de
θ = 30o respecto del eje z positivo. Determinar la fuerza que se ejerce sobre el protón.
Solución:
Este problema se aborda haciendo uso de la expresión de la fuerza magnética sobre cargas en movimiento,
derivada experimentalmente y siendo:
F~B~ = q ~v × B
~. (3.1)
En este caso, cada protón presenta q = e = 1,602 × 10−19 C. Como B ~ = B0 k̂ y ~v = vx î + vz k̂, tenemos que:

î ĵ k̂
F~B~ = q vx 0 vz = q(−vx B0 )ĵ ,

(3.2)
0 0 B0
donde las barras verticales indican el cálculo del determinante. Por otro lado, tenemos que vx = v0 sin (θ), con
lo que la fuerza sobre uno de los protones será:
F~B~ = −qv0 sin(θ)B0 ĵ = −5,1264 × 10−14 N ĵ . (3.3)
3.2. Electrones en movimiento en un campo magnético

Un haz de electrones cuya energı́a cinética es K,
sale de una “ventana” de lámina delgada en el extremo
de un tubo acelerador. Existe una placa de metal a una
distancia d de esta ventana, perpendicular a la direc-
ción de movimiento inicial de los electrones. Demostrar
que se puede impedir que el haz choque contra q la placa
si se aplica un campo magnético B tal que B ≥ 2m eK
e2 d2
donde me es la masa del electrón. ¿Como deberı́a estar
orientado B?~
Solución:
Asumiendo que se trata de un problema de energı́as no relativı́stas, tenemos que K = 21 me v02 . Además, la
carga de los electrones es q = −e, y en el momento que salen del tubo, sienten la fuerza del campo magnético
cuya expresión es la misma que antes:
F~B~ = q~v × B
~ = −e~v × B
~. (3.4)
Como la fuerza total sobre los electrones es no nula, la segunda ley de Newton se puede utilizar sobre cada
electrón:
r
2K
FB~ = me a = ev0 B = e B. (3.5)
me
Notar que asumimos que ~v0 ⊥ B, ~ lo cual usaremos más adelante, y además esta ecuación es válida solo para
las magnitudes de la fuerza y la aceleración. Sabemos que la trayectoria que describen las partı́culas cargadas
en presencia de campos magnéticos uniformes tiene forma de cı́rculo, y para movimientos circulares uniformes
de radio R la aceleración es centrı́peta, es decir:
v02 2K
a= = . (3.6)
R me R
33
Podemos entonces encontrar el radio del cı́rculo que describirán los electrones en presencia del campo
magnético:
r r
2K 2K 2K me
=e B ⇒ R= .
R me eB 2K
Para que el haz de electrones no colisiones contra la pantalla, necesariamente R ≤ d, con lo cual:
r
2K me
d ≥
eB 2K
r r
2K me 22 K 2 me
B ≥ =
ed 2K e2 d2 2K
r
2me K
= ,
e2 d2
que no es otra cosa que lo que querı́amos probar.
Por otro lado, sabemos que la fuerza que sienten los electrones apenas salen del tubo debe ser perpendicular
a su dirección de movimiento y hacia el positivo del eje z, basándonos en el dibujo. Con esta dirección asociada
a la fuerza, y la dirección que tiene el vector velocidad de los electrones dentro del tubo, usamos que la relación
entre la fuerza, el campo magnético y la velocidad está dada por un producto vectorial. Por esta razón, el campo
magnético necesariamente es perpendicular tanto a ~v como a F~ , debiendo ser un vector perpendicular al plano
de la hoja. La regla de la mano derecha permite deducir que éste debe apuntar “hacia adentro de la hoja”.
Esto puede verificarse matemáticamente de la siguiente forma: en el momento que los electrones salen del tubo,
podemos elegir F~B~ = FB~ k̂, ~v = v î y con esto:
~ = −e v î × B ĵ ,
FB~ k̂ = −e v î × B (3.7)
y para que la componente de la fuerza apunte hacia el positivo del eje z, si o si B = −|B|. Entonces, el campo
magnético queda:
~ = −|B|ĵ .
B (3.8)
34
3.3. Conductor en un campo magnético
Una barra recta horizontal de cobre, transporta una corriente de 52 A de oeste a este en una región com-
prendida entre los polos de un gran electroimán.
En esta región hay un campo magnético horizontal

hacia el noreste (45o al norte del este), cuya magnitud
es de 12 T .
1. Encontrar la magnitud y dirección de la fuerza

sobre una sección de 1 m de la barra.
2. Conservando la barra en posición horizontal,

¿como se debe orientar para que la magnitud de
la fuerza sea máxima? ¿cuál es dicha magnitud
en este caso?
Solución:
Sabemos que sobre una carga q en movimiento en un campo magnético, se produce una fuerza dada por
F~ = q~v × B.
~ Si la carga está en movimiento dentro de un conductor recto, de forma tal que en un intervalo de
tiempo ∆t recorra la longitud completa del mismo, digamo ∆x, entonces es directo que:
∆~x
F~ = q ~,
×B (3.9)
∆t
donde únicamente calculamos la velocidad promedio de desplazamiento de las cargas en la barra. Ahora, sabemos
q
que en este tiempo, la carga q recorre la barra, generando por conservación de la carga una corriente i = ∆t .
Entonces, si cambiamos el nombre de la longitud de la barra por ~l = ∆~x, donde ~l es un vector cuyo módulo es
la longitud de la barra, y su dirección y sentido estan dados por la circulación de la corriente, tenemos que la
fuerza sobre dicho conductor es de la forma:
F~ = i~l × B
~. (3.10)
Con esta fórmula, estamos en condiciones de resolver el problema.
1. Si elegimos Norte = ĵ, Este = î, tenemos que el campo magnético solo tiene componentes en estas
~ = Bx î + By ĵ. Entonces, como ~l = lî por definición de vector “longitud” de
direcciones, de forma que B
la barra conductora, podemos calcular:

î ĵ k̂
F~ = l

0 0 = ilBy k̂ = 441,23 N k̂ . (3.11)

Bx By 0
2. La magnitud de la fuerza, dado que tenemos un producto vectorial, es F = ilB sin(α), donde α es el
ángulo entre ~l y B.
~ En el inciso anterior, este ángulo era de 45o , y si hubieramos escrito la componente
By , verı́amos que la magnitud de la fuerza es ilBy = ilB sin(45o ). Ahora, Fmax coincide justamente con
α = π2 , donde sin(π/2) = 1. Entonces, el alambre deberı́a estar orientado hacia el sureste (45o al sur del
este) para que el ángulo resultante entre los vectores sea de 90o sin modificar la orientación del campo.
En esta posición la magnitud de la fuerza serı́a F = 624 N .
3.4. Comparación: fuerza magnética vs. fuerza gravitatoria

Un segmento de alambre de cobre, recto y horizontal, porta una corriente i = 28 A. ¿Cuáles son la magnitud
y dirección del campo magnético necesarias para “hacer flotar” el alambre, i.e., para equilibrar su peso, si su
densidad lineal de masa es de 46,6 g/m?
35
Solución:
Los datos que brinda el problema en un sistema de unidades coherentes son λ = 46,6 g/m = 4,66×10−2 kg/m,
i = 28 A. Esta corriente, para un circuito tı́pico (orden del miliamper), es MUY alta. Si la fuerza ejercida por
el campo gravitatorio es:
P~ = −mg ĵ = −λl g ĵ , (3.12)

necesitamos un campo B ~ tal que F~ ~ = λl g ĵ para que el segmento de alambre esté en equilibrio. Si éste está
B
posicionado sobre el eje x, tenemos que:

î ĵ k̂ n o

λl g ĵ = l 0 0 = i −lBx ĵ + lBy k̂ . (3.13)

Bx By Bz
Para satisfacer la ecuación, necesariamente By ≡ 0. Entonces:

λg ~ = −0,016 T k̂
lλ g = −ilBz ⇒ Bz = − = −0,016 T ⇒ B (3.14)
i
Como se puede apreciar, si queremos utilizar campos magnéticos pequeños, conseguibles sin mucha sofisti-
cación (militesla es del orden de un imán de refrigerador), la corriente debe ser muy alta para acercanos al valor
de la fuerza gravitatoria.
36
4. Problemas de óptica
4.1. Reflexión y refracción en un prisma
Un haz delgado y colimado de luz blanca, incide sobre un prisma de vı́drio de ı́ndice de refracción 1.6,
ubicado de dos maneras diferentes como se muestra en la figura.
1. Trazar cualitativamente la marcha del haz principal en cada caso, justificando donde se produce reflexión
total interna.
2. Indicar dónde se produce descomposición del haz en colores, mostrando cualitativamente la trayectoria
de diferentes colores.
Solución:
(i) Para abordar la solución del problema, se deben tener en cuenta las consideraciones relacionadas con el
principio de Fermat y su consecuencia, la ley de Snell. Según el principio de Fermat, la luz recorre el
camino de más corta duración entre dos puntos en un mismo medio. A partir de este concepto se pueden
deducir las leyes de reflexión y refracción (o ley de Snell). Según esta última, en una interfaz entre
dos medios, midiendo los ángulos de manera consistente utilizando la normal a dicha interfaz en el punto
donde incide un rayo de luz, se cumple que:
ni sin(θi ) = nr sin(θr ) , (4.1)

donde ni y nr son los ı́ndices de refreacción del material desde donde incide el rayo y el material hacia
el cual penetra el rayo respectivamente, y θi , θr son los ángulos de incidencia y refracción respectivamente.
Cada uno de los rayos paralelos entre si, que conforman el haz de luz blanca, y que incidan sobre las caras
del prisma, obedecen esta ley. Por ejemplo, en el caso (i), el haz incide con ángulo θi = 0, por lo cual la
ley de Snell dice que θr = 0, y todos los rayos siguen en lı́nea recta. Esto no es del todo cierto, dado que
la ley de reflexión dice que θi = θR = 0, es decir, existen rayos que son reflejados en la misma dirección
en que venı́an, sobre el haz incidente.
Ahora, al llegar a la segunda interfaz, esta presenta una inclinación tal que el ángulo que forma cualquiera
de los rayos con la misma es 45o . Según la ley de reflexión, el ángulo de cada rayo reflejado será idéntico
al ángulo de incidencia de cada correspondiente rayo incidente. Esto puede verse en el gráfico donde se
traza el seguimiento del rayo principal.
37
Si queremos calcular el ángulo de refracción, la ley
de Snell dice que:
nv
sin(θr ) = sin(θi )
na
1,6 π
= sin
1 4
≈ 1,13
Pero esto es imposible, dado que
− 1 < sin(θ) < 1 . (4.2)
Lo que sucede aquı́ es lo que se conoce como re-

flexión total interna.
Bajo esta condición, no existe rayo refractado, produciendose una reflexión de todos los rayos que llegan
a la interfaz. Lo interesante de esta situación es que existen ángulos de incidencia para los cuales sucede,
pero otros para lo cuales no. Evidentemente, se puede estimar un ángulo crı́tico θic a partir del cual no
existirá rayo refractado, es decir, asumiendo que la función seno alcanza su máximo valor permitido:
1 =sin(θr )
nv
= sin(θic )
na
na
sin(θic ) = .
nv
Un par de aclaraciones sobre esta fórmula son necesarias: únicamente es válida cuando el rayo incidente,
lo es desde el medio de mayor ı́ndice de refracción, siendo además un θic una magnitud dependiente única-
mente de la configuración de materiales en la interfaz.
Es importante que θic es el ángulo a partir del cual existe reflexión total interna, es decir que todo rayo con
θi > θic sufrirá reflexión total interna. Para el problema que estamos trabajando, tenemos que θic ≈ 38o ,
con lo que el ángulo de incidencia de 45o es más que suficiente para que estemos en condición de reflexión
total interna.
(ii) En este caso, el haz incidente se encuentra en primera instancia con una interfaz que presenta una incli-
nación respecto de la dirección de incidencia. Por la simetrı́a del triángulo, tenemos que θi = 45o . Sabemos
que los materiales presentan variaciones en su ı́ndice de refracción según los distintos colores (el ı́ndice de
refracción es función de la frecuencia o de la longitud de onda asociada a la luz). Por ejemplo, el vı́drio
Flint presenta:
nrojo
v = 1,644 ; namarillo
v = 1,650 ; nazul
v = 1,665
Por lo tanto, por cada rayo de luz blanca incidente sobre la cara del prisma, habrá más de un rayo
refractado, satisfaciendo:

na
θrrojo = arcsin sin(θi ) = 25o 290
nrojo
v

na
θramarillo = arcsin sin(θi ) = 25o 220
namarillo
v

na
θrazul = arcsin sin(θi ) = 25o 70
nazul
v
38
Como vemos, el color rojo presenta un ángulo de refracción mayor que el amarillo, y estos dos uno mayor
que el azul. Esto significa que el rayo de luz blanca sufrirá una descomposición en colores debido a la
diferente respuesta del medio material a las diferentes longitudes de onda (recordar que cada color esta
asociado a una longitud de onda). Los rayos se refractan ángulos diferentes al atravesar la interfaz, y
consecuentemente se pueden distinguir por separado.
EL RESTO DEL EJERCICIO ESTA COMENTADO EN EL .TEX PERO HAY QUE VERIFICAR SI
ESTÁ BIEN LO QUE DICE DE LA SEGUNDA INTERFAZ Y SI EL GRÁFICO ES COHERENTE.
39
5. Matemática general
5.1. Cambio de sistemas de coordenadas
Supongamos que tenemos un sistema de coordenadas, por ejemplo cartesiano {x, y, z}, y quiero escribir un
vector ~v en otro sistema {x1 , x2 , x3 }. Para esto es necesario conocer las ecuaciones que vinculan ambos sistemas,
es decir:
x = x(x1 , x2 , x2 ) ; y = y(x1 , x2 , x2 ) ; z = z(x1 , x2 , x2 ) . (5.1)
Por ejemplo, con el sistema de coordenadas esféricas:
x = r cos(φ) sin(θ) ; y = r sin(φ) sin(θ) ; z = r cos(θ) . (5.2)
Pensando en el producto interno canónico del sistema cartesiano, un pequeño desplazamiento en cada
coordenada conforma un camino diferencial cuyo largo puede calcularse como:
d~s · d~s = ds2 = dx2 + dy 2 + dz 2 . (5.3)
Al realizar un cambio en el sistema de coordenadas, estamos interesados en que ds2 sea invariante (mantenga
el mismo valor). En el sistema {x1 , x2 , x3 }, este diferencial puede tomar la forma general:
ds2 = h21 dx21 + h22 dx22 + h23 dx23 + · · · + hij dxi dxj + · · · (5.4)
Las funciones hj se conocen como factores de escala. En un sistema curvilineo ortogonal, hi,j = 0 para
los términos cruzados. Entonces, en un sistema curvilineo ortogonal genérico, existen versores mutuamente
ortogonales x̂1 , x̂2 y x̂3 tales que:
~v = vx î + vy ĵ + vz k̂ = v1 x̂1 + v1 x̂1 + v1 x̂1 . (5.5)
¿De que forma se relacionan? Sabemos que:
d~s = dxî + dy ĵ + dz k̂ = h1 dx1 x̂1 + h2 dx2 x̂2 + h3 dx3 x̂3 . (5.6)
Esto permite leer los diferenciales con respecto a las coordenadas del nuevo sistema de la forma:
∂x ∂x ∂x
dx = dx1 + dx2 + dx3 , (5.7)
∂x1 ∂x2 ∂x3
∂y ∂y ∂y
dy = dx1 + dx2 + dx3 , (5.8)
∂x1 ∂x2 ∂x3
∂z ∂z ∂z
dz = dx1 + dx2 + dx3 . (5.9)
∂x1 ∂x2 ∂x3
Entonces:

∂x ∂x ∂x
d~s = dx1 + dx2 + dx3 î + (· · · ) ĵ + (· · · ) k̂ (5.10)
∂x1 ∂x2 ∂x3

∂x ∂y ∂z
= î + ĵ + k̂ dx1 + (· · · ) dx2 + (· · · ) dx3 . (5.11)
∂x1 ∂x1 ∂x1
De esta forma, por comparación, se obtiene que:
∂x ∂y ∂z
h1 x̂1 = î + ĵ + k̂ , (5.12)
∂x1 ∂x1 ∂x1
∂x ∂y ∂z
h2 x̂2 = î + ĵ + k̂ , (5.13)
∂x2 ∂x2 ∂x2
∂x ∂y ∂z
h3 x̂3 = î + ĵ + k̂ . (5.14)
∂x3 ∂x3 ∂x3
(5.15)
r 2 2 2
∂x ∂y ∂z
Entonces, obtenemos que |hj x̂j | = hj = ∂xj + ∂xj + ∂xj . Observamos entonces que para encontrar
d~s en el nuevo sistema de coordenadas curvilineo ortogonal, estamos aplicando una transformación cuya matriz
es:  
∂x ∂x ∂x
∂x1 ∂x2 ∂x3
J˜ =  ∂y ∂y ∂y
, (5.16)
 
∂x1 ∂x2 ∂x3
∂z ∂z ∂z
∂x1 ∂x2 ∂x3
40
donde cada columna es como el gradiente de la coordenada cartesiana correspondiente, por lo que son ortogonales
entre ellos, i.e., J˜ es invertible. Sin embargo, como base de esta transformación lineal, se toma la matriz:
 
1 ∂x 1 ∂x 1 ∂x
h1 ∂x1 h2 ∂x2 h3 ∂x3
1 ∂y 1 ∂y 1 ∂y
J = . (5.17)
 
h1 ∂x1 h2 ∂x2 h3 ∂x3
1 ∂z 1 ∂z 1 ∂z
h1 ∂x1 h2 ∂x2 h3 ∂x3
De esta forma, cualquier vector ~v puede transformarse de la forma:

   
v1 vx
 v2  = J T  vy  . (5.18)
v3 vz
Entonces, podemos obtener las relaciones entre sistemas ortogonales para diferentes operadores diferenciales,
y de esta forma obtener las formulas de las tablas al principio de estas lectures: CHEQUEAR LOS INDICES,
SOBRE TODO DE ∇ × ~v .
1 ∂f 1 ∂f 1 ∂f
∇f = x̂1 + x̂2 + x̂3 ,
h1 ∂x1 h2 ∂x2 h3 ∂x3

1 ∂(h1 h2 h3 v1 ) ∂(h1 h2 h3 v2 ) ∂(h1 h2 h3 v3 )
∇ · ~v = + + ,
h1 h2 h3 ∂x1 ∂x2 ∂x3

1 ∂h3 v3 ∂h2 v2 1 ∂h2 v2 ∂h3 v3 1 ∂h3 v3 ∂h1 v1
∇ × ~v = − x̂1 + − x̂2 + − x̂3 ,
h1 h2 ∂x2 ∂x2 h1 h3 ∂x3 ∂x1 h1 h2 ∂x2 ∂x2

1 ∂ h2 h3 ∂f ∂ h1 h3 ∂f ∂ h1 h2 ∂f
∆f = + + .
h1 h2 h3 ∂x1 h1 ∂x1 ∂x2 h2 ∂x2 ∂x3 h3 ∂x3
5.2. Ecuaciones diferenciales

5.2.1. Ecuaciones diferenciales de primer orden homogéneas
dy dy
Separable: → dx = f (x)g(y) → g(y) = f (x)dx.
dy y
y du u f (u)
Homogénea 1: → dx =f x . Se reemplaza u = x y queda dx + x = x .
Homogénea 2: En el caso más complejo en que:

dy a1 x + b1 y + c1
=f , (5.19)
dx a2 x + b2 y + c2
se plantea que:
a1 x + b1 y + c1 = 0 a b
⇒ det(sist) = 1 1 = a1 b2 − a2 b1 ~ .
a2 x + b2 y + c2 = 0 a2 b2

dy
Si ~ = 0 entonces dx = f a(x+y)+c 1
b(x+y)+c2 , y llamando u = x + y se vuelve separable.
Si ~ 6= 0 entonces existen x0 , y0 tales que las rectas del numerador y denominador en el argumento de f se
cruzan. Entonces proponemos η = x − x0 e ξ = y − y0 , lo cual al introducirlo en la ecuación la transforma en
una ecuación homogénea 1 en ξ/η.
5.2.2. Ecuaciones diferenciales de primer orden inhomogéneas

dy
Factor integrante: → dx + f (x)y = g(x), se multiplican ambos miembros por una función arbitraria u(x),
de forma que: uy + uf (x)y = ug ⇒ u0 = uf (x) ecuación diferencial separable para u. Entonces (uy)0 = ug es
0
integrable.
dy
Ecuación de Bernoulli: → dx + f (x)y = g(x)y n , entonces se hace la sustitución y = z n − 1 y se obtiene una
ecuación de tipo factor integrante.
dy
Ecuación exacta: → F (x, y) + G(x, y) dx = 0, se resuelve con formalismo de diferenciales:
F (x, y)dx + G(x, y)dy = 0 ,
41
∂ψ ∂ψ ∂ψ ∂ψ
entonces, si ∂F ∂G
∂y = ∂x , existe una función ψ tal que dψ = ∂x dx + ∂y dy siendo F (x, y) = ∂x y G(x, y) = ∂y . De
0
esta forma, la ecuación queda (ψ) = 0 y finalmente ψ = constante, pudiendo despejar y.
Ecuación inexacta: → igual que la anterior pero con ∂F ∂G

∂y 6= ∂x . Entonces, proponemos u = u(x, y) tal que
uF dx + uGdy = 0 (u la vuelve exacta). Se debe entonces satisfacer que: ∂uF ∂uG
∂y = ∂x . Generalmente se pide que
u = u(x) o u = u(y), se obtiene el valor de u y luego se integra exactamente.
n
d y dy
Ecuación exacta de orden n: → an (x) dx n + · · · + a1 (x) dx + a0 (x)y = f (x), cumpliendo que
(n)
a0 (x) − a01 (x) + a002 (x) − · · · + (−1)n + an (x) = 0, entonces la ecuación es exacta. Se puede poner:
dn y dn−1 y

d
an (x) n + · · · + a0 (x)y = bn−1 (x) n−1 + · · · + b0 (x)y ,
dx dx dx
y comparar los coeficientes para armar una ecuación diferencial de menor orden.
Ecuación de Clairaut: → (xy 0 − y) = f (y 0 ). Se derivan ambos miembros con respecto a x:
y 0 + xy 00 − y 0 = xy 00 = f 0 (y 0 )y 00 → y 00 (x − f 0 (y 0 )) = 0 ,
y entonces y 00 = 0 dá la solución general, x − f 0 (y 0 ) = 0 dá la solución sigular.
5.2.3. Ecuaciones diferenciales de segundo orden

y0 yy 00 −(y 0 )2
Random: para ecuaciones del tipo yy 00 − (y 0 )2 + f (x)y 2 = 0, sirve elegir z = y y entonces dz
dy = y2 .
Ecuación con solución en serie:→ y 00 + P (x)y 0 + Q(x)y = 0. Si P y Q tienen puntos singulares regulares, se
propone una solución en serie de Frobenius:
∞
X ∞
X
y = xσ an xn = an xn+σ ,
n=0 n=0
previamente habiendo factorizado el comportamiento en infinito (si existe un punto singular regular). Esto se
verifica proponiendo y = z1 y a la nueva ecuación se le busca el punto singular regular en z = 0. La serie se
puede centrar en cualquier punto.
Polinomio indicial: puede tener tres posibles soluciones:

1) 2 raices iguales: solo una solución en forma de serie de Frobenius.
2) Raices difieren en un entero: la raı́z más grande es la que dá la solución de Frobenius.
3) Raices no difieren en entero: las dos brindas soluciones linealmente independiente.
Segunda solución (para 1) y 2) anteriores): si σ1 y σ2 son las raices del polinomio indicial con σ1 > σ2 e
y1 corresponde a la solución con σ1 , la segunda solución se obtiene de la forma:
R z 1 − R u P (v)dv
a) Método del Wronskiano: y2 (z) = y1 (z) y12 (u)
e du.
∂
b) Método de la derivada: y2 (z) = ∂σ [(σ − σ2 )y(z, σ)] σ=σ2 , en donde:
∞
X
y(z, σ) = an (σ)xn+σ ,
n=0
con an (σ) saliendo de la relación de recurrencia sin reemplazar el valor de σ.

∂
Si las raices son iguales, y2 (x) = ∂σ y(z, σ) σ=σ2 .
42
5.2.4. Ecuaciones diferenciales en derivadas parciales
5.3. Tensores
Sea V un espacio vectorial de dual V ∗ y T : V ∗ × · · · × V ∗ × V × · · · × V → C, donde V ∗ aparece r veces
y V aparece s veces, y C es un cuerpo. Entonces, T es un tensor de rango (r, s) que toma r vectores de V ∗
(covectores) y s vectores de V y da como resultado un escalar, pudiendo escribirse utilizando las bases:
r s
{~ei }|i=1 ∈ V ; {~ej }i=j ∈ V ∗ ,
siendo dim(V ) = n, de la forma:

···ir
T = Tji11ji22···j ~e ⊗ · · · ⊗ ~eir ⊗ ~ej1 ⊗ · · · ⊗ ~ejs ,
s i1
(5.20)
···ir
donde Tji11ji22···js
son las componentes del tensor en esta base. Ejemplos simples:
Tensor de rango (1, 0) es un vector ~v = v i~ei (vertical).

Tensor de rango (0, 1) es un vector dual ~v ∗ = vj ~ej (horizontal).
Antiguamente se utilizaba la siguiente denominación:
Base Componentes Clasificación

Espacio V covariante ~ei contravariante xi Vectores
Espacio V ∗ contravariante ~ej covariante xj Covectores
5.3.1. La métrica
La métrica g ij es un tensor dos veces covariante, es decir, toma componentes de vectores (componentes
contravariantes).
La métrica gij es un tensor dos veces contravariante, es decir, tomar componentes de covectores (compo-
nentes covariantes).
La métrica g ij es un tensor mixto.
Además, g ij = (gij )−1 y g ir grj = g ij = δ ij siempre.
El producto interno puede escribirse en términos de la métrica:
~a · ~b = g ij ai bj = δ ij ai bj = gij ai bj . (5.21)
5.3.2. Simetrı́a y antisimetrı́a

A C todo tensor se le puede extraer sus partes simétrica y antisimétrica, definidas de la siguiente forma:
1
C [αβ] = C αβ − C βα → parte antisimétrica,

2
1
C (αβ) C αβ + C βα → parte antisimétrica,

=
2
con las propiedades que C [αβ] = −C [αβ] y C (αβ) = C (βα) .
43
5.3.3. Transformación de tensores en coordenadas curvilineas
∂x0β
Si tenemos xα y x0β relacionados de la forma x0β = Lβα xα , tenemos entonces que ∂xα = Lβα . Ası́, podemos
ver como transforman los tensores de primeros rangos fácilmente:
∂x0α β
Vector contravariante → a0α = a = Lαβ aβ ,
∂xβ
∂xβ
Vector covariante → a0α = aβ = (L−1 )αβ aβ ,
∂x0α
∂xγ ∂xδ
Tensor covariante → Tαβ = Tγδ = (L−1 )αγ (L−1 )βδ Tγδ ,
∂x0α ∂x0β
∂x0α ∂x0β γδ
Tensor contravariante → T αβ = T = Lαγ Lβδ T γδ ,
∂xγ ∂xδ
∂x0α ∂xδ γ
Tensor mixto → T αβ = T = Lαδ (L−1 )βδ T γδ .
∂xγ ∂x0β δ
44
A. Scattering cuántico
Se puede estudiar la materia de varios libros de grado de mecánica cúantica, como por ejemplo de,1 ,2 etc.
Supongamos que disponemos un detector D que mide la cantidad de partı́culas dn dispersadas en un ángulo
sólido dΩ por un centro dispersor ubicado en el origen del sistema de coordenadas O y según el cual el potencial
de interacción entre las partı́culas dispersadas y el haz incidente es de la forma V (~r1 − ~r2 ).
Figura 1: Esquema de dispersión de una onda plana incidente en una onda esférica dispersada un ángulo θ
respecto del vector de onda incidente.
Dado que observaremos al sistema desde su centro de masa, nos queda que V (~r1 − ~r2 ) = V (~r). Designamos
con Fi el número de partı́culas incidentes en el punto O que atraviesan una superficie perpendicular a la dirección
de incidencia. Evidentemente, esperamos que el número dn cumpla que dn ∝ dΩ, pues si agrandamos el área
de detección obtendremos una mayor probabilidad de colectar partı́culas. Asimismo, dn ∝ Fi pues cuantas
más partı́culas incidan mayor es la probabilidad de tener partı́culas dispersadas. Definimos el coeficiente de
proporcionalidad de la forma σ(θ, φ) satisfaciendo:
dn = σ(θ, φ)Fi dΩ . (A.1)

σ tiene unidades de superficie y se conoce como sección eficaz diferencial de dispersión. Podemos entonces
interpretar esta ecuación como que el número de partı́culas que llegan al detector por unidad de tiempo corre-
sponden a aquellas que atraviesan un área σ(θ, φ)dΩ perpendicular a la dirección del haz de incidencia.
Los estados cuánticos dispersivos deben estar asociados a una solución de la ecuación de Schrödinger. Es
esperable que la partı́cula libre incidente tenga un comportamiento de onda plana, y dado que en general la
onda es dispersada en distintas direcciones según las caracterı́sticas del sistema en cuestión, proponemos que:
~ eikr
ψscattering ∼ eik·~r + f (θ, φ) , r→∞, (A.2)
r
donde ~k es el vector de onda asociado a la energı́a del sistema, que satisface k 2 = 2mE}2 , y se asume dispersión
~
elástica (|k| = k) con decaimiento asintótico de la onda esférica dispersada. La función f (θ, φ) se conoce como
amplitud de dispersión. En otras palabras, suponemos la existencia de una solución a la ecuación de Schrödinger
para la dispersión de dos partı́culas en el centro de masa del sistema de forma que asintóticamente (distancia
rel-

1
ativa entre partı́culas r → ∞) tome la forma anterior. Esta aproximación es válida siempre que V (~r) ∼ O |~
r| ,
lo cual deja afuera el potencial coulombiano. Este último es un caso particular entre decaimiento rápido (ex-
ponencial) y decaimiento lento, y debe resolverse aparte. Nos interesa entonces encontrar formas concretas de
calcular f (θ, φ) conociendo información básica sobre el sistema fı́sico en estudio.
Por otro lado, sabemos que a la ecuación de Schrödinger le corresponde una corriente de la forma:
}
J~ = [ψ̄∇ψ − ψ∇ψ̄] , (A.3)
2mi
lo cual se puede interpretar como el flujo diferencial de partı́culas por unidad de área y tiempo. Calculemos la
corriente para el término dispersado de nuestra función de onda. Para ello necesitamos el operador diferencial
∇ escrito en coordenadas esféricas:
45
∂ 1 ∂ 1 ∂
∇ = r̂ + θ̂ + φ̂ . (A.4)
∂r r ∂θ r sin(θ) ∂φ
ikr
Aplicándolo al término ψscattered ∼ f (θ, φ) e r , obtenemos que:
ikeikr eikr eikr ∂f eikr ∂f

∇ψscattered ∼ r̂f − r̂f 2 + θ̂ 2 + φ̂ 2 , (A.5)
r r r ∂θ r sin(θ) ∂φ
y que:
ike−ikr e−ikr e−ikr ∂f ∗ e−ikr ∂f ∗
∇ψ̄scattered − ∼ r̂f ∗ − r̂f ∗ 2 + θ̂ 2 + φ̂ 2 . (A.6)
r r r ∂θ r sin(θ) ∂φ
Vemos que salvo el primer término, al calcular ψ∇ψ̄ y ψ̄∇ψ, todos los términos son de orden O(1/r3 ), y
como estamos viendo r → ∞, nos quedamos solo con este. De esta forma, tenemos:
−ikr
eikr eikr ∗ e−ikr

} ∗e
J~scattered ∼ f r̂f ik +f r̂f ik
2mi r r r r
2
}k|f (θ, φ)|
= r̂ .
mr2
El número de partı́culas asociadas a esta corriente en un ángulo sólido entonces está dado por:
}k
dn = |r̂ · J~scattered |r2 dΩ =
|f (θ, φ)|2 . (A.7)
m
Por otro lado, podemos calcular la cantidad de partı́culas incidentes sobre la superficie perpendicular a la
~
dirección de incidencia utilizando la corriente asociada al estado incidente ψincident ∼ eik·~r . Tenemos que:
~
∇ψincident ∼ i~keik·~r , (A.8)
y que:
~
∇ψ̄incident ∼ −i~ke−ik·~r , (A.9)
con lo cual:
}~k
J~incident = . (A.10)
m
Entonces, Fi que definimos antes se puede calcular de la forma Fi = |k̂ · J~incident | = }k

m. Tenemos entonces
que:
dn = σ(θ, φ)Fi dΩ (A.11)

}k
= σ(θ, φ) dΩ (A.12)
m
}k
= |f (θ, φ)|2 dΩ , (A.13)
m
con lo que obtenemos:
σ(θ, φ) = |f (θ, φ)|2 . (A.14)

Tenemos entonces una forma práctica de calcular la sección eficaz de interacción para distintos procesos,
simplemente conociendo la solución correspondiente de la ecuación de Schrödinger para el sistema en estudio,
y separando la parte correspondiente a la onda dispersada. Existen diversos métodos para conocer la amplitud
de dispersión, como el método de ondas parciales que se puede ver en proximas secciones.
46
A.1. Determinación de la amplitud de dispersión
Si asumimos que nuestra solución de la ecuación de Schrödinger para dispersión se comporta asintóticamente
2
como una onda libre, sabemos que E = }k 2m , lo cual puede ingresarse en la ecuación misma, reacomodar un
poco los términoa resultantes y obtener que nuestra función de onda debe satisfacer la ecuación:
2mV (~r)
{∇2 + k 2 }ψ(~r) = U (~r)ψ(~r) , U (~r) = . (A.15)
}2
Esta es la ecuación de Helmholtz inohomegenea. La solución general de una ecuación inhomogenea se puede
conocer como la suma de una solución particular de la inhomogenea con una solución general de la homogenea.
Una solución de la parte inohomegenea se puede conocer por medio de la función de Green, que es aquella
función dependiente de un par de conjunto de coordenadas que satisface:
{∇2 + k 2 }G(~r − ~r 0 ) = δ (3) (~r − ~r 0 ) . (A.16)

De esta forma, se puede verificar que:
Z
0
ψinh (~r) = G(~r − ~r )U (~r 0 )ψinh (~r 0 )d3~r 0 , (A.17)
0
es solución de la ecuación A.16. Notar que la variable ~r está restringida a aquella región donde U (~r) 6= 0, con
lo que una de las condiciones para que el problema tenga solución es que esta función tenga soporte compacto
(un dominio finito donde es no nula). A.17 puede resolverse de forma iterativa para ψinh como veremos. De esta
forma, la solución dispersiva de nuestro sistema queda:
Z
0
i~
ψ(~r) ∼ e k·~
r
+ G(~r − ~r )U (~r 0 )ψ(~r 0 )d3~r 0 , (A.18)
~
pues eik·~r es una solución de la ecuación homogenea. Comparando esta última con A.2 es esperable que la
información sobre la amplitud de dispersión aparezca en la función de Green. Para dar con una expresión de la
misma, escribimos los desarrollos en Fourier:
Z
1 0
G(~r − ~r 0 ) = 3
G̃(~q)ei~q·(~r−~r ) d3 ~q (A.19)
(2π)
Z
1 0
δ (3) (~r − ~r 0 ) = 3
ei~q·(~r−~r ) d3 ~q . (A.20)
(2π)
Aplicamos el operador diferencial {∇2 + k 2 } a G e igualamos a la delta de Dirac. Reemplazamos momen-

~ = ~r − ~r 0 :
taneamente R
Z
~ 1 ~
{∇2 + k 2 }G(R) = d3 ~qG̃(~q){∇2 + k 2 }ei~q·R
(2π)3
Z
1 ~
= d3 ~qG̃(~q){−q 2 + k 2 }ei~q·R
(2π)3
Z
1 ~
= ei~q·R d3 ~q .
(2π)3
Esto tiene solución únicamente si:

1
G̃(~q) = . (A.21)
k2 − q2
Esto no es otra cosa que la transformada de Fourier de la función de Green que buscamos. De esta forma,
para conocer la función de Green en variables espaciales, que es nuestro interés, debemos aplicar la transformada
inversa.
47
~
e−i~q·R 3
Z
~ 1
G(R) = d ~q
(2π)3 k2 − q2
Z ∞ Z π −iqR cos(θ)
1 e
= 2
q 2 dq sin(θ)dθ
(2π) 0 0 k2 − q2
Z ∞ Z 1 iqRu
1 e
= q 2 dqdu
(2π)2 0 −1 k 2 − q2
Z ∞ −iqR Z ∞ iqR
i e e
= qdq − qdq
R(2π)2 0 k 2 − q 2 0 k2 − q2
Z ∞ −iqR Z 0
e−ivR

i e
= qdq + vdv , (A.22)
R(2π)2 0 k 2 − q 2 2
−∞ k − v
2
donde usamos que u = cos(θ) y v = −q. Ahora, si quisieramos juntar estas integrales en una sola extendiendose
en el intervalo q ∈ (−∞, ∞), debemos tener cuidado con las singularidades en q = ±k. Para ello es recomendable
realizar la cuenta en el plano complejo, utilizando algún contorno que encierre a uno u otro de los polos:
Figura 2: Posibles integrales de contorno que encierran alguno de los dos polos.
Cual elijamos nos dirá que tipo de solución es la que estamos contemplando. Como veremos, una de ellas nos
devuelve una onda que “viaja hacia afuera” mientras que el otro “viaja hacia adentro”. Elegimos el contorno
C1 que encierra el polo q = k.
qeiqR
Z
~ i
G(R) = 2
dq
R(2π) C1 (k − q)(k + q)
Z
i f (q)
= dq
R(2π)2 C1 q − k
i
= (2πif (k))
R(2π)2
i keikR
= 2
(2πi) ,
R(2π) 2k
donde utilizamos el teorema de Cauchy-Goursat, que establece

R fque si C1 encierra un polo simple en z = z0 de
(z)
f y ésta no presenta ninguna singularidad en C1 , entonces C1 z−z 0
dz = 2πif (z0 ). Obtenemos ası́ la función de
Green que buscamos:
0
~ =− 1 eikR 1 eik|~r−~r |
G(R) =− . (A.23)
4π R 4π |~r − ~r 0 |
La expresión asintótica de esta función se calcula aproximando la magnitud de ~r −~r 0 por aquella de ~r menos
la proyección de ~r 0 sobre ~r, pues al suponer que |~r| → ∞ y ~r 0 es no nula únicamente en el soporte compacto de
U , ésta proyección es entonces muy pequeño frente a |~r|.
48
Figura 3: Aproximación asintótica basada en el soporte compacto del potencial para la función de Green de
estados dispersivos. La región encerrada por la linea de puntos indica cual es la región donde es válido el
potencial.
Entonces:
|~r − ~r 0 | −−−−→ r − r̂ · ~r 0 , (A.24)

|~
r |→∞
y además:
0 0
eik|~r−~r | eik(r−r̂·~r ) eikr −ikr̂·~r 0
0
−−−−→ ∼ e . (A.25)
|~r − ~r | |~r|→∞ r − r̂ · r~0 r
Si ponemos esto en la solución A.18, sale que:
eikr
Z
~ 1 0
ψ(~r) ∼ eik·~r + − e−ikr̂·~r U (~r 0 )ψ(~r 0 )d3~r 0 , (A.26)
r 4π
con lo que obtenemos una expresión para la amplitud de dispersión, en vistas de la expresión A.2:
Z
1 2m 0
f (θ, φ) = − e−ikr̂·~r V (~r 0 )ψ(~r 0 )d3~r 0 . (A.27)
4π }2
A.2. Aproximación de Born

Si hubieramos integrado A.22 utilizando C2 en vez de C1 , obtenemos:
0
± 01 e±ik|(~r−~r )|
G (~r − ~r ) = − , (A.28)
4π |~r − ~r 0 |
con lo que la ecuación integral A.18 equivalente a la ecuación de Schrödinger queda:
U (~r 0 ) ±ik|~r−~r 0 | ± 0 3 0
Z
(±) ~ 1
ψ~ (~r) = eik·~r − e ψ~ (~r )d ~r , (A.29)
k 4π |~r − ~r 0 | k
conocida como ecuación de Lippmann-Schwinger. Esta ecuación puede resolverse de forma iterativa, proponiendo
(0)
que en primera instancia la función de onda del integrando es una onda plana. Insertamos la expresión ψ~ (~r) =
k
~
eik·~r , y obtenemos la aproximación a orden uno:
0
U (~r 0 )eik|~r−~r | (0) 0 3 0
Z
(1) (0) 1
ψ~ ≈ ψ~ (~r) − ψ~ (~r )d ~r . (A.30)
k k 4π |~r − ~r 0 | k
Asi siguiendo con las iteraciones, obtenemos una expresión:

0
U (~r 0 )eik|~r−~r | (n) 0 3 0
Z
(n+1) ~ 1
ψ~ (~r) ≈ eik·~r − ψ~ (~r )d ~r . (A.31)
k 4π |~r − ~r 0 | k
49
Podemos, de forma análoga a la deducción que hicimos en A.27, obtenemos una muy buena aproximación
de la amplitud de dispersión:
Z
1 0 (n)
f (n) (θ, φ) ≈ − U (~r 0 )e−ikr̂·~r ψ~ (~r 0 )d3~r 0 . (A.32)
4π k
Esta serie se conoce como serie de Born, y el caso n = 0 de se conoce como aproximación de Born. Podemos
reescribir la expresión resultante como:
Z
1 ~0 0 ~ 0
f B (θ, φ) = − U (~r 0 )e−ik ·~r eik·~r d3~r 0
4π
1 (2π)3/2
Z
~ ~0 0 3 0
= − 3/2
U (~r 0 )ei(k−k )·~r d ~r
4π (2π)
π 1/2 2m
= − F(V (~k − ~k 0 )) , (A.33)
2 }2
donde F simboliza la transformada de Fourier. De esta forma, obtenemos la expresión en aproximación de Born
para la amplitud de dispersión, y por consiguiente para la sección eficaz:
2πm2
dσ B (θ, φ) = |F(V (~k − ~k 0 ))|2 dΩ . (A.34)
}4
A.3. Método de las ondas parciales

Presentamos aqui uno de los métodos usuales para encontrar la amplitud de dispersión en el caso de
potenciales centrales de interacción que decaen más rápido que O(1/r). Sabemos que en estos casos, el momento
angular orbital L ~ se conserva, existiendo estados estacionarios con valores de momento angular bien definidos
que dan lugar a autoestados bien definidos. Estos autoestados forma una auto-base común de los operadores Ĥ,
~ 2 y Lz , y los llamarémos las “ondas parciales”. Como en r → ∞ el potencial decae completamente, debemos
L
obtener que nuestras ondas parciales corresponden a soluciones de la ecuación de Schrödinger simétricamente
esféricas, estacionarias, de vacio y tres dimensionales de la forma:
1/2
2k 2

(0)
φk,l,m (~r) = jl (kr)Yl m (θ, φ) . (A.35)
π
Donde jl son las funciones de Bessel esféricas. En el caso de un V (r) arbitrario con la anterior condición,
podemos escribir la onda parcial de la forma:
1
uk,l (r)Yl m (θ, φ) .
φk,l,m (~r) = (A.36)
r
uk,l satisface la ecuación radial reducida que conocemos:
}2 ∂ 2 l(l + 1)}2 }2 k 2

− + + V (r) u k,l (r) = uk,l (r) , (A.37)
2m ∂r2 2mr 2m
donde m es la masa reducida del sistema de partı́culas interactuantes. Esta ecuación debe satisfacer la condición
uk,l (r = 0) = 0 pues es básicamente una ecuación unidimensional de Schrödinger para el potencial Vef f (r) =
l(l+1)}2
2mr + V (r) que satisface Vef f (r < 0) = ∞. Evaluamos la solución asintótica de esta ecuación. Ésta queda
de la forma:
∂ 2 u∞
k,l
− = k 2 u∞
k,l , (A.38)
∂r2
y las soluciones son bien conocidas:
u∞
k,l (r) ∼ Ae
ikr
+ Be−ikr . (A.39)
Si pensamos en términos de onda reflejada y transmitida, como el r < 0 no puede haber onda transmitida a
través del potencial, necesariamente la onda incidente y la reflejada tienen la misma intensidad, i.e. |A| = |B|.
~ ~
Por otro lado de una ecuación real se esperan soluciones reales, con lo que A = |A|eiδ(k) y B = −|A|e−iδ(k) .
Como vemos no podemos suponer desde un principio que dicha fase sea independiente de todos los parámetros
del problema. Esto nos dice que:
50
u∞ ~
k,l ∼ C sin(kr + δ(k)) . (A.40)
Comparemos esta solución con la expresión asintótica de la solución A.35. Tenemos que:

1 lπ
jl (kr) ∼ sin kr − , (A.41)
kr 2
con lo que conviene expresar δ(~k) = − lπ ~
2 + δl (k) para poder comparar. Tenemos entonces que:

∞ lπ ~
uk,l ∼ C sin kr − + δl (k) (A.42)
2
lo cual puede reescribirse en términos de exponenciales y reemplazarlo en la expresión general de la función de
onda para obtener:
1/2
2k 2 [eikr e−ilπ/2 − e−ikr eilπ/2 ] m

(0)
φk,l,m = Yl (θ, φ) (A.43)
π 2ikr
~ ~
C [eikr e−ilπ/2 eiδl (k) − e−ikr eilπ/2 e−iδl (k) ] m
φ∞
k,l,m = Yl (θ, φ) . (A.44)
r 2i
Podemos elegir C de forma que:
~ ~
[eikr e−ilπ/2 eiδl (k) − e−ikr eilπ/2 eiδl (k) ] m
φ∞
k,l,m = C1 Yl (θ, φ) . (A.45)
2ikr
Además, como las fases globales no cambian la interpretación fı́sica de la función de onda en términos
~
probabilı́sticos, conviene sacar factor común e−iδl (k) y absorber el factor con C1 . Obtenemos:
~
[eikr e−ilπ/2 e2iδl (k) − e−ikr eilπ/2 ] m
φ∞
k,l,m = C1 Yl (θ, φ) . (A.46)
2ikr
Vemos entonces que la diferencia de fase en la parte de la onda que “viaja hacia adelante” con respecto a la
onda en vacio (sin interacción con un potencial) es la que contiene toda la información fı́sica sobre lo que hace
el potencial. Comparar A.46 con A.43.
VER EN COHEN COMO ES EL TEMA CON LA ROTACION. SINO EN EL SAKURAI
Por las consideraciones anteriores, tenemos que nuestro sistema presenta simetrı́a azimutal alrededor del eje
sobre el cual yace el vector de onda incidente ~k. Esto nos indica que f (θ, φ) = f (θ). Podemos proponer entonces:
∞
X
f (θ) = f~k (l)(2l + 1)Pl (cos(θ)) , (A.47)
l=0
en donde:
Z π
1
f~k (l) = f (θ)Pl (cos(θ)) sin(θ)dθ . (A.48)
2 0
Con esto, podemos calcular la sección eficaz diferencial, utilizando la ortogonalidad de los polinomios de
Legendre:
Z π
2
Pl (cos(θ))Pl0 (cos(θ)) sin(θ)dθ = δl, l0 , (A.49)
0 2l + 1
obteniendose:
51
Z
σ = σ~k (θ)dΩ (A.50)
Z
= |f (θ)|2 dΩ (A.51)
XX Z
= (2l + 1)(2l0 + 1)f~k (l)f~k (l0 ) Pl (cos(θ))Pl0 (cos(θ))dΩ (A.52)
l l0
∞
X
= 4π (2l + 1)|f~k (l)|2 (A.53)
l=0
Dado que la simetrı́a azimutal aplica a las soluciones de la ecuación de Schrödinger, es posible escribir:
∞
X
ψ~+ (r, θ) = (2l + 1)C~k R~k,l (r)Pl (cos(θ)) . (A.54)
k
l=0
R~k,l son las funciones que encontramos anteriormente, que asintóticamente se comportan de la forma A.42,
i.e.:

∞
X sin kr − lπ
2 + δl (~k)
ψ~+ (r, θ) ∼ (2l + 1)C~k (l) Pl (cos(θ)) . (A.55)
k kr
l=0
Por otro lado, mezclamos los desarrollos ?? y A.47 (ARREGLAR DONDE APARECEN ESTOS DESAR-
ROLLOS Y POR ENDE EL METODO DE REFERENCIA, I.E. EL DESARROLLO DE LA EXP DEBERIA
APARECER POR PRIMERA VEZ ACA), tenemos que la expresión asintótica de A.2 queda:
∞
eikr

X sin(kr − lπ/2)
ψ~+ (r, θ) ∼ (2l + 1)Pl (cos(θ)) il + f~k (l) . (A.56)
k kr r
l=0
Proponiendo que A.55=A.56 obtenemos una ecuación que relaciona los coeficientes del desarrollo asintótico
de la función de onda con el desarrollo de la amplitud de dispersión:
sin(kr − lπ/2 + δl (~k)) sin(kr − lπ/2) eikr

C~k (l) = il + f~k (l) . (A.57)
kr kr r
Desarrollando la función seno del primer miembro en términos de exponenciales, sumando y restando el
término eikr e−ilπ/2 y acomodando correctamente se puede obtener que:
eikr −ilπ/2 −iδl (~k) 2iδl (~k) ~ sin(kr − lπ/2) sin(kr − lπ/2) eikr
C~k (l) e e (e − 1) + C~k (l)e−iδl k = il + f~k (l) , (A.58)
2ikr kr kr r
lo cual finalmente da un par de ecuaciones a resolver:
( ilπ/2 ~ ~
C~k (l) e 2ik e−iδl (k) (e2iδl (k) − 1) = f~k (l)
~ (A.59)
C~k (l)e−iδl (k) = il
~
Conviene despejar de la segunda C~k (l) = il eiδl (k) y reemplazar en la primera para obtener:
~
e2iδl (k) − 1 1 ~
f~k (l) = = eiδl (k) sin(δl (~k)) , (A.60)
2ik k
Siendo finalmente la sección eficaz de interacción:
∞
4π X
σ= (2l + 1) sin2 (δl (~k)) . (A.61)
k2
l=0
52
A.4. Teorema óptico
Utilizando la expansión A.47, ya sabiendo el resultado A.60 se obtiene:
∞
1X ~
f~k (θ) = (2l + 1)eiδl (k) sin(δl (~k))Pl (cos(θ)) . (A.62)
k
l=0
Si observamos como es la transmisión de la onda incidente en la misma dirección de incidencia, el ángulo

de dispersión es θ ≈ 0. Como Pl (1) = 1, obtenemos:
∞
1X ~
f~k (0) = (2l + 1)eiδl (k) sin(δl (~k))
k
l=0
= <(f~k (0)) + i=(f~k (0)) ,
P∞
donde =(f~k (0)) = 1
k l=0 (2l + 1) sin2 (δl~k), lo cual es básicamente la expresión de la sección eficaz. Tenemos
entonces que:
4π
σ~k = =(f~k (0)) . (A.63)
k
Esta expresión es la conocida como Teorema Óptico, y relaciona lasección eficaz total de dispersión con la
amplitud frontal de la onda dispersada, es decir la parte dispersada “hacia adelante”, conocida como forward
scattering.
A.5. Espı́n en dispersión
B. Interacciones nucleares (fuerte y débil)

Arreglatelas como puedas....
C. Teorı́a de perturbaciones
C.1. Perturbaciones constantes
C.1.1. Regla de oro de Fermi
ESTARIA BUENO COMPLETAR; MUY BUENAS IDEAS EN GALINDO.
Como siempre, tenemos un hamiltoniano de la forma Ĥ = Ĥ0 + V , donde V será tratado como una
perturbación, en este caso independiente del tiempo. Ésto puede tomarse de la forma:

0 si t < 0
V (t) = (C.1)
V si t ≥ 0 (independiente del tiempo).
Aplicamos entonces la teorı́a de perturbaciones dependientes del tiempo a ésta perturbación. Recordemos
que la amplitud de transición de un estado |si a uno |ki (ambos autoestados de Ĥ0 ) al aplicar la perturbación,
a primer orden queda:
i t
Z
0
hk|T (t, 0)|si = Tks (t, 0) = − hk|V |sieiωks t dt0 . (C.2)
} 0
Ek −Es
Aquı́, ωks = } como siempre. Como en t > 0 V no depende del tiempo, podemos calcular:
Z t
i 0
Tks (t, 0) = − hk|V |si eωks t dt0
} 0
hk|V |si iωks t
= − e −1 .
Ek − Es
2
La probabilidad de transición está dada por Ps→k (t) = |Tks (t, 0)| . Si calculamos esto se obtiene:
53
|hk|V |si|2 iωks t
− 1 e−iωks t − 1

Ps→k (t) = 2
e (C.3)
(Ek − Es )
|hk|V |si|2
= 2 (1 − cos(ωks t)) . (C.4)
(Ek − Es )2
Tenemos entonces que la probabilidad de transición es una función oscilatoria con un tiempo caracterı́stico
dado por τ = ω2π
ks
.
Recordemos que para poder tratar V como una perturbación, es necesario que Ek ' Es , es decir, que la
energı́a sin perturbación “se conserve”. En el sentido estricto esto no es cierto, quien se conserva es la energı́a
total pues hĤi es rigurosamente una cantidad conservada, pero la energı́a perturbada Ek será muy próxima a
Es si la perturbación es débil. Si esto se satisface, tenemos entonces que ωks ∼ 0. Comparado con el tiempo
caracterı́stico, ωks τ ∼ 2π, y como argumento del coseno es análogo evaluar ωks t ∼ 0. Entonces:
(ωks t)2 1 2 2
1 − cos(ωks t) ∼ 1 − 1 + = ωks t . (C.5)
2! 2
La probabilidad queda entonces:
1
|hk|V |si|2 t2 .
Ps→k (t) ∼ (C.6)
}2
Esto indica que para tiempos menores a la oscilación caracterı́stica de la probabilidad, la misma crece de
forma cuadrática con el tiempo, lo cual es una muy buena caracterización de la probabilidad de transición.
Supongamos ahora que se tiene una variedad de estados finales perturbados {|ki}k∈f con energı́as Ek ∼ Es ,
que en alguna aproximación pueden suponerse continuos. Esto es válido pues sabemos que los estados finales son
autoestados de Ĥ0 que es el hamiltoniano sin perturbar, y por ende esperamos que |hk|ki|2 < ∞. Se introduce
entonces la densidad de estados como la cantidad de estados finales perturbados nf con energı́as en el intervalo
(Ek , Ek + dEk ):
dnf
ρf (Ek) = . (C.7)
dEk
Con esta herramienta es posible calcular la probabilidad total de transición desde el estado |si cambiando
las sumatorias por integrales de la forma que sigue, a partir de (C.4):
X X |hk|V |si|2
P (t) = Ps→k (t) = 2 (1 − cos(ωks t)) (C.8)
(Ek − Es )2
k k
Z Ek +∆E
2 |hk 0 |V |si|2
= (1 − cos(ωk0 s t))ρf (Ek0 )dEk0 (C.9)
}2 Ek −∆E ωk20 s
−1
2 (Ek +∆E)} |hk 0 |V |si|2
Z
= (1 − cos(ωk0 s t))ρf (ωk0 s )dωk0 s . (C.10)
} (Ek −∆E)}−1 ωk20 s
Evaluar esta integral requiere entender ciertas aproximaciones del modelo de perturbación utilizado. Sabe-
mos que debemos satisfacer la condición Ek ∼ Es . Entonces ρf (Ek ) ∼ ρf (Es ) ≈ cte. Ahora, los estados {|ki}k∈f
son similares a |si porque pertenecen al mismo hamiltoniano con energı́as muy parecidas entre si, con lo que
|hk|V |si|2 ≈ cte también. Podemos verificar que esto es cierto observando el factor (1 − cos(ωt))/ω 2 , con t > 0.
Si graficamos esto en función de ω para distintos tiempos, se tiene que:
54
Figura 4: Comportamiento a lo largo del tiempo del factor del integrando para calcular la probabilidad de
transición.
Como vemos, la función se hace cada vez más estrecha alrededor de ω = 0 a medida que el tiempo pasa.
Esto nos dice que para un tiempo suficientemente largo, pero todavı́a lo suficientemente chico como para que
(C.6) < 1 (recordar que es una probabilidad), podemos evaluar la integral con lı́mites infinitos, pues el integrando
tendrá un soporte muy pequeño (casi δ(ω)). Ahora, como para este tiempo ω ≈ 0 y ω = E − Es tenemos que
E ≈ Es con lo que la aproximación que nos permitió sacar afuera de la integral los otros factores es válida.
Falta entonces evaluar la integral:
∞
(1 − cos(ωt)) 1 ∞ 1 − cos(x)
Z Z
dω = dx
−∞ ω2 t −∞ x2
" ∞ Z ∞ #
cos(x) 1
= t − − − (− sin(x))dx
x −∞ −∞ x
Z ∞
sin(x)
= t dx
−∞ x
= πt ,
quedando finalmente la probabilidad de transición expresada de la forma:
2πt
|hk|V |si|2 ρf (Ek ) .
P (t) = (C.11)
}
La tasa de transición asociada es simplemente ws→k = dP dt , lo cual resulta en:
2π
ws→k = |hk|V |si|2 ρf (Ek ) , (C.12)
}
i.e., tenemos una tasa constante. Esta es la conocida como regla de oro de Fermi.
C.1.2. Sección eficaz de transición

Si se interpreta esta tasa de probabilidad de transición como la probabilidad de que luego de una interacción
se obtenga una partı́cula con estado |ki, es posible relacionarla con la sección eficaz a través de la especificación
de la tasa del flujo incidente de partı́culas Fi de la forma:
ws→k
dσ = . (C.13)
Fi
Ahora, la tasa del flujo de partı́culas incidentes con momento k lo conocemos de (A). Escribimos entonces
la densidad de estados correspondiente a un continuo de estados en el intervalo (Ek , Ek + dEk ). Sabemos que, a
55
más de una constante de normalización del espacio de las fases que tiene que ver con el volumen L3 de una caja,
que se cancela con la misma normalización utilizada para los estados (por eso no lo incluimos), la densidad de
estados satisface:
VER PORQUE LO SATISFACE
d3 p~k
ρf (Ek )dEk = (C.14)
(2π})3
p2k dpk dΩ
= (C.15)
(2π})3
(C.16)
p
Ahora, sabemos que Ek = p2k c2 + m2 c4 , con lo que, usando que p = γmv y E = γmc2 , se puede calcular:
dEk 1 1
= p 2pk c2
dpk 2 p2k c2 + m2 c4
1
= 2pk c2 (C.17)
Ek
1
= γk mvk c2
γk mc2
= vk (C.18)
Haciendo uso de esta relación, tenemos que:
p2k dΩ
ρf (Ek ) = . (C.19)
vk (2π})3
Con esto, se tiene que:
2π p2 dΩ
ws→k = |hk|V |si|2 k , (C.20)
} vk (2π})3
y reemplazando en (C.13) se obtiene:
2π p2 dΩ m
dσ = |hk|V |si|2 k (C.21)
} vk (2π})3 }ks
2π E 2 vk dΩ 1
= |hk|V |si|2 k4 (C.22)
} c (2π})3 vs
donde se uso que pk = Ec2k vk despejando desde (C.17). Finalmente la sección eficaz a partir de la regla de oro
de Fermi queda de la forma:
dσ 1 E 2 vk
= 2 4
|hk|V |si|2 4k . (C.23)
dΩ (2π) } c vs
D. Aproximación dipolar
ACA COMPLETAR CON LA APROXIMACION DIPOLAR PARA EL CALCULO DE LA SECCION
EFICAZ DE UNA INTERACCION DE UNA PARTICULA CARGADA CON UN CAMPO ELECTROMAG-
NETICO, INTRODUCIR EL DIPOLE STRENGH OSCILLATOR Y LA REGLA DE LA SUMA QUE DA
UNO O Z PARA MUCHOS ELECTRONES.
E. Cálculos usuales
E.1. Colisión frontal y máximo intercambio energético
Supongamos que dos partı́culas de masas en reposo m1 y m2 colisionan a energı́as relativı́stas (vi ∼ c,
i = 1, 2). Podemos plantear las condiciones de conservación de energı́a y momento relativı́stas para estimar
56
información sobre lo que sucede con nuestras partı́culas luego de la colisión. Recordemos que este tipo de
problemas suelen tratarse haciendo un cambio entre sistemas de referencia, considerandose uno de ellos el
laboratorio, desde donde v1lab 6= 0 mientras que v2lab = 0, i.e. una de las partı́culas está en reposo. El otro es
el conocido centro de masa, pero no haremos uso del mismo para la resolución del problema. Entonces, para
simplificar notación utilizaremos únicamente v1 . Luego de la colisión ambas partı́culas tendran velocidades no
nulas dadas por vif , i = 1, 2.
Figura 5: Esquema de una colisión entre dos partı́culas.
Si escribimos las energı́as correspondientes tenemos:
Ei = γ1i m1 c2 + m2 c2 , (E.1)
para el estado inicial, y:
Ef = γ1f m1 c2 + γ2f m2 c2 = γ1i m1 c2 − ∆ + m2 c2 + ∆ .

(E.2)
para el estado final, donde:

1 v
γ=p ,β= , (E.3)
1 − β2 c
i.e., el factor de Lorentz ya conocido de la relatividad especial y ∆ es la transferencia de energı́a en el proceso
de colisión.
Por otro lado, podemos escribir el momento total desde el sistema del laboratorio antes y después de la
colisión. Tenemos que:
q
pi1 = (γ1i )2 m21 c2 − m21 c2 (E.4)
pi2 = 0 (E.5)
y para el estado final:
s s
2
(E1f )2 γ1i m1 c2 − ∆
pf1 = − m21 c2 = − m21 c2 (E.6)
c2 c2
s s
2
(E2f )2 (m2 c2 + ∆)
pf2 = − m22 c2 = − m22 c2 , (E.7)
c2 c2
donde utilizamos la conocida expresión relativı́sta E 2 = p2 c2 + m2 c4 . Por conservación del momento, debemos
satisfacer la ecuación:
57
pi = pf → pi1 + pi2 = pf1 − pf2 , (E.8)
donde el signo menos pone de manifiesto que en la suma vectorial (si la hubieramos hecho) elegimos θ = π, i.e.
la colisión es frontal y por ende es la que más energı́a entrega. Ver en la figura 5 que esto implica que cada
partı́cula se mueve en dirección contraria a la otra, i.e. tienen momentos opuestos. De esta forma, ∆ = Wmax
es la máxima transferencia de energı́a en una sola colisión. Por más simplicidad aún llamaremos m1 = M y
m2 = m, y dado que solo γ1i entra en juego, lo llamaremos γ. Tenemos entonces que resolver:
s s
2 2 2
(γM c − Wmax ) 2 c2 −
(mc2 + Wmax ) p
2
− M 2
− m2 c2 = γ 2 M 2 c2 − M 2 c2 , (E.9)
c c
Resolver esta ecuación es bastante largo, pero el resultado de despejar la máxima transferencia de energı́a
queda:
2mc2 γ 2 β 2
Wmax = . (E.10)
m m 2

1+ Mγ + M
m
Si ponemos s = M y η = γβ, obtenemos:
2mc2 η 2
Wmax = p . (E.11)
1 + 2s 1 + η 2 + s2
Notar que podemos tomar varios lı́mites, como ser M >> m, y v << c lo cual deriva en s ∼ 0 o η ∼ β.
El primer lı́mite sirve por ejemplo en el caso de que la partı́cula incidente sea un ión pesado y colisiona con un
electrón, cálculo útil para determinar el stopping power electrónico. El segundo lı́mite es simplemente la versión
no relativı́sta.
E.2. Energı́a de formación de una esfera cargada

Supongamos que queremos traer cargas desde infinito hasta formar una esfera de radio R y carga total Q
con distribución uniforme de cargas de forma que su densidad es
Q Q
ρ= = 4 3
. (E.12)
V 3 πR
Vamos a asumir que ya tenemos formada una esfera de radio r cuya carga es Q(r) = 34 πr3 ρ, y que acercamos
a la misma una pequeña carga dq distribuida uniformemente en un carcarón esférico de ancho dr. La carga del
mismo ya sobre la esfera será dq = 4πr2 drρ. Sabemos que el campo externo a la primera esfera coincide con el
de una carga puntual, de forma que:
Q(r)
E(r, r0 ) = k , (E.13)
r0
donde r es el radio de la esfera que contiene la carga Q(r) y r0 es el radio exterior a dicha esfera, es decir, el
punto campo observado. Una carga dq que se mueve en este campo desde el infinito hasta el punto de interés r,
debe entregar una energı́a que se calcula de la forma:
Z r Z r
Q(r) dqQ(r)
dW = −dq E(r, r0 )dr0 = −dq k 0 dr0 = k . (E.14)
∞ ∞ r r
Si escribimos todo en términos de ρ y de r, tenemos que:
16 2 4
dW = kπ r ρdr . (E.15)
3
El truco está en pensar que esta es la energı́a para formar una primera capa de ancho dr, habiendo una
pequeña esferita antes, con lo que hacer la integral entre 0 y R nos dirá la energı́a que se requiere para formar
la esfera en su totalidad.
Z R
16 2 2 5
W = dW =kπ ρ R . (E.16)
0 15
Reemplazando la expresión de ρ para la distribución uniforme, obtenemos:
3 Q
W = k . (E.17)
5 R
58
E.3. Modelado de un dieléctrico: modelo de Drude
Aqui
VER JACKSON
E.4. Consideraciones sobre (ω)

Sabemos que la expresión general para una respuesta lineal de un material dieléctrico3 puede escribirse
como:
Z ∞
(ω) = 1 + f (τ )eiωτ dτ , (E.18)
0
donde f es una función que depende de las propiedades del medio. En general, es complejo. Sean 1 (ω) =
Re((ω)) y 2 (ω) = Im((ω)). Entonces se satisface que (−ω) = ∗ (ω). Veamos que significa esto en términos
de las partes real e imaginaria de :
(−ω) = 1 (−ω) + i2 (−ω)

Z ∞ Z ∞
= 1+ f (τ ) cos(−ωτ )dτ + i f (τ ) sin(−ωτ )dτ
0 0
Z ∞ Z ∞
= 1+ f (τ ) cos(ωτ )dτ − i f (τ ) sin(ωτ )dτ
0 0
= 1 (ω) − i2 (ω)
= ∗ (ω) .
Se tiene entonces que 1 (−ω) = 1 (ω) y 2 (−ω) = −2 (ω), es decir son par e impar respectivamente. Esto
lleva además a que:
1 1 1 (ω) 1 (ω)
= = 2
−i , (E.19)
(ω) 1 (ω) + i2 (ω) |(ω)| |(ω)|2
con lo que es evidente que Re|1/(ω)| también es par, mientras que Im|1/(ω)| es impar. Esto nos lleva a la
conclusión de que al calcular:
Z ∞
ω
dω , (E.20)
−∞ (ω)
pueda despreciarse el término real por ser impar, y reemplazar media integral de la parte compleja por un dos:
Z ∞ Z ∞
ω 1
dω = 2i ω Im
dω . (E.21)
−∞ (ω) 0 (ω)
F. Ecuación de Dirac
G. Potencial, energı́a potencial, trabajo: campo electrostático
Conviene realizar un resumen de todos los casos posibles y sus significados para tener bien en claro lo que
sucede con los signos de la energı́a potencial (por ende del trabajo) en función de los signos de la carga y la
diferencia de potencial electrostático.
Por convención, decimos que el campo eléctrico de una carga puntual positiva “apunta hacia afuera” de
dicha carga. En el caso de una carga negativa, decimos que “apunta hacia adentro”.
Si elegimos una carga positiva, y la soltamos cerca de otra carga positiva, la carga que elegimos tiende a
alejarse de la segunda por repulsión. Si la misma carga positiva la soltamos cerca de una carga negativa, ahora
tiende a acercarse a la segunda carga por atracción. Entonces, podemos decir que las cargas positivas se mueven
siempre “a favor” del campo electrostático.
59
Si elegimos una carga negativa, y la soltamos cerca de una carga positiva, la carga que elegimos tiende a
acercarse de la segunda por atracción. Si la misma carga negativa la soltamos cerca de otra carga negativa,
ahora tiende a alejarse a la segunda carga por repulsión. Entonces, podemos decir que las cargas positivas se
mueven siempre “en contra” del campo electrostático.
La definición de potencial en el punto ~rf se define a partir del valor del potencial en el punto de partida ~ri
utilizando la integral de camino recorrido (y haciendo uso de que la fuerza electrostática es conservativa):
Zf Zf
V (~rf ) = V (~ri ) − ~ · d~r ⇒ ∆V = V (~rf ) − V (~ri ) = −
E ~ · d~r .
E (G.1)
i i
De aquı́ podemos ver que:

~ · d~r < 0, o lo que es lo mismo,
Si ∆V > 0, el potencial en ~rf es mayor que en ~ri , pero esto significa que E
que el camino recorrido va en contra del campo. Entonces, el campo eléctrico apunta en la dirección
de zonas de de mayor potencial a zonas de menor potencial.
~ · d~r > 0, o lo que es lo mismo,
Si ∆V < 0, el potencial en ~rf es menor que en ~ri , pero esto significa que E
que el camino recorrido va a favor del campo. Entonces, la conclusión es la misma, el campo eléctrico
apunta en la dirección de zonas de de mayor potencial a zonas de menor potencial.
Con todo esto estamos en condiciones de analizar algunos casos dependientes del signo de la carga en juego.
Supongamos que la carga que traigo desde ~ri hacia ~rf es positiva, es decir q > 0. Entonces:
1. Si ∆V > 0, el camino recorrido es de una zona de mayor potencial a una de menor. En el gráfico,
menor potencial a una de mayor. En el gráfico, se ve que la carga, al ser positiva, naturalmente
la carga, al ser positiva, naturalmente se moverı́a se moverı́a utilizando la lı́nea de color rojo (pun-
utilizando la lı́nea de color rojo (punteada), a fa- teada), a favor del campo. Además, en este caso,
vor del campo. Pero en este caso, el camino recor- el camino recorrido casualmente va a favor del
rido va en contra. Esto implica que el agente ex- campo. Entonces, el agente externo no debe hac-
terno deberı́a hacer un trabajo para poder mover er fuerza porque la carga naturalmente sigue este
la carga. En efecto, ∆V > 0, q > 0 y por lo tanto camino. Entonces, la carga “realiza un trabajo
el trabajo es positivo (gana energı́a potencial): sobre el agente externo”. En efecto, ∆V < 0,
q > 0 y por lo tanto el trabajo es negativo (pierde
W = ∆Ep = q∆V > 0 . energı́a potencial):
W = ∆Ep = q∆V < 0 .
2. Si ∆V < 0 el camino recorrido es de una zona de
Supongamos ahora que la carga que traigo hasta ~rf desde ~ri es negativa, es decir, q < 0. Entonces:
60
Si ∆V > 0 el camino recorrido es de una zona de Si ∆V < 0, el camino recorrido es de una zona
menor potencial a una de mayor. En el gráfico se de mayor potencial a una de menor. En el gráfico
ve que la carga, al ser negativa, naturalmente se se ve que la carga, al ser negativa, naturalmente
moverı́a utilizando la lı́nea de color azul (puntea- se moverı́a utilizando la lı́nea de color azul (pun-
da), en contra del campo. Además en este caso, teada), en contra del campo. Pero en este caso,
el camino recorrido casualmente va en contra del el camino recorrido es a favor del campo. Esto
campo. Entonces, el agente externo no debe hac- implica que el agente externo deberı́a hacer un
er fuerza porque la carga naturalmente sigue este trabajo para evitar que la carga siga su camino
camino. Entonces, la carga “realiza un trabajo natural. En efecto, ∆V < 0, q < 0 y por lo tanto
sobre el agente externo”. En efecto, ∆V > 0, el trabajo es positivo (gana energı́a potencial):
q < 0 y por lo tanto el trabajo es negativo (pierde
energı́a potencial):
W = ∆Ep = q∆V > 0 .
W = ∆Ep = q∆V < 0 .
Notar que la lı́nea equipotencial más grande representa el valor de potencial mayor en todos los gráficos,
para que sea más facil entender la dirección asignada al campo eléctrico.
61
Referencias
1
B.D.F.L. Claude Cohen-Tannoudji. Quantum Mechanics Volume 2. Hermann, 2006.
2
Eugene Merzbacher. Quantum Mechanics. John Wiley & Sons, Inc., New York, second edition, 1970.
3
L.D. Landau and E.M. Lifshitz. Electrodynamics of continuous media. Pergamon Press, Oxford, first edition,
1960.
62

Geser GeneralPhysics 2018

Cargado por

Copyright:

Formatos disponibles

Geser GeneralPhysics 2018

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Geser GeneralPhysics 2018

Cargado por

Copyright:

Formatos disponibles

Lectures: Problemas de fı́sica general

Lic. Geser, Federico Alejandro

UNIVERSIDAD NACIONAL DE CÓRDOBA

2. Problemas de electrostática básica 4

3. Problemas de magnetostática básica 33

B. Interacciones nucleares (fuerte y débil) 53

G. Potencial, energı́a potencial, trabajo: campo electrostático 59

NA → número de Avogadro, NA = 6, 022 × 1023 mol−1 .

c → velocidad de la luz, c = 299,792,458 m/s.

me → masa del electrón, me = 9, 109 × 10−31 kg.

e → carga del electrón, e = −1,602 × 10−19 C.

h, ~ → constante de Planck, h = 6,626 × 10−34 Js , ~ = 1,054 × 10−34 Js.

G → constante de gravitación universal, G = 6,7 × 10−11 N m2 /kg 2 .

k → constante de Coulomb, k = 9 × 109 N m2 /C 2 .

mn , mp → masa neutron y proton respectivamente, mn = 1,674×10−27 kg, mp = 1,672×10−27 kg.

z ∗ → si z ∈ C entonces z ∗ es el complejo conjugado.

Sistema cartesiano Sistema cilı́ndrico Sistema esférico

~r(t) = xî + y ĵ + z k̂ ~r(t) = ρρ̂ + z k̂ ~r(t) = rr̂

Relaciones entre versores dlϕ = r sin(θ)dϕ

Polinomios de Legendre asociados:

Formulación Hamiltoniana Corchetes de Poisson

Donde d es la distancia entre ejes paralelos El problema se resuelve encontrando un sistema

Para un eje de rotación

Ecuaciones de Maxwell Energı́a del campo electromagnético: vacı́o

Caso vacı́o Energı́a de la ecuación de onda

Solución de Laplace con simetrı́a azimutal

Solución de Laplace sin simetrı́a azimutal

Dirichlet → φ|S = 0 o f (~x) con ~x ∈ S.

Cauchy → combinación de las anteriores.

Expansión de la función de Green

Condiciones de contorno dieléctricos:

Maxwell para guias de ondas: Ecuaciones varias guia de onda:

∂E ~t P = e−2γz Pideal , Potencia guia de

Modo TE: Modo TM:

Siempre γ 2 = µω 2 − k 2 , con µ los coeficientes

Potencia asociada a una guia de ondas:

En modo TE ⇒ se eligen las opciones de abajo de las llaves y ψ = Bz

El cuerpo de masa m = 120 g que se muestra en la

1. La velocidad tangencial del cuerpo.

2. La tensión del hilo.

3. La velocidad angular necesaria para reducir la

m = 120 g = 0,12 kg , L = 50 cm = 0,5 m ;

El diagrama que usaremos para un determinado in-

Nuevamente, utilizando trigonometrı́a, tenemos que:

P~ = P xî + P y ĵ = −P sin(α)î − P cos(α)ĵ ,

T = m(g + ω 2 L sin(α)) cos(α)

1. Tenemos que la fuerza electrostática obedece la ley

Esto se esto puede expresarse de la forma:

Comparandolo con el mismo valor para el electrón, tenemos:

2.2. Cristales y fuerzas electrostáticas

qCl qN a N m2 (−(1, 6 × 10−19 C)2 )

2.3. Cuadrado de cargas

1. Comenzemos por determinar el lado del cuadra-

Si incluimos la dirección del vector a la magnitud encontrada, tenemos que:

Esto se integra en el tiempo, dando:

Usando el dato de velocidad inicial sale que vx (t = 0 s) = v0 y vy (t = 0 s) = 0. Volvemos a integrar en el

Al imponer las condiciones iniciales nuevamente, se obtiene:

x(t) = at + b ; y(t) = ct + d . (2.22)

x(t) = vx1 t + x1 = v0 t (2.23)

Siempre γ 2 = µω 2 − k 2 , con µ los coeficientes

0 2πr hE(r) = λh , (2.66)