Capı́tulo 1:

Electrostática y Magnetostática

Los primeros experimentos interesantes en electrostática fueron hechos por Cavendish (∼ 1770) y Coulomb
(∼ 1780), y en el caso magnetostático por Biot y Savart (1820) y Ampere (∼ 1820). Maxwell resumió
todo esto en las leyes del electromagnetismo. En un sentido estas leyes definen unas “cantidades extrañas ”
llamadas cargas y corriente que regulan la magnitud de las respectivas fuerzas. En orden a relacionar estos
dos conceptos es necesario tener una teorı́a microscópica para las partı́culas elementales, en particular del
electrón como una entidad discreta.

Una descarga electrostática fractal

1
Índice
1. Electrostática 4

2. Campo eléctrico y potencial eléctrico 5

3. Comentarios sobre la ley de la fuerza 7

3.1. Variación para una ley de fuerza cuadrática . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
3.2. Masa del fotón . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8

4. Condiciones de Borde 8

5. Ecuación de Laplace y funciones ortogonales 11

5.1. Bases rectangulares . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
5.2. Bases polares . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
5.3. Bases esféricas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
5.3.1. Funciones de Legendre . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
5.3.2. Funciones de Legendre asociadas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
5.4. Bases cilı́ndricas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
5.5. Proyectores en general . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25

6. La Ecuación de Poisson y la función de Green 26

6.1. Función de Green en el espacio infinito . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
6.1.1. En 3D . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
6.1.2. En 2D . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
6.1.3. En 1D . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
6.2. Función de Green para un volumen arbitrario . . . . . . . . . . . . . . . . . . . . . . . . . . 29
6.3. Método de las imágenes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
6.4. Expansión de la función de Green en bases ortonormales . . . . . . . . . . . . . . . . . . . . 35
6.4.1. G en coordenadas esféricas para espacio infinito por método directo . . . . . . . . . . 35
6.4.2. G en coordenadas esféricas entre dos esferas conductoras por método directo . . . . . 37
6.4.3. G en coordenadas esféricas fuera de una esfera conductora por método indirecto . . . 38
6.4.4. Expansión en bases de un problema relacionado . . . . . . . . . . . . . . . . . . . . . 41

7. Magnetostática 42

8. Campo magnético y potencial magnético 43

9. Método de solución 45

10.Fuerza de Lorentz 48

11.Solución numérica 49
11.1. Problemas en 1D . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
11.1.1. “Shooting Method” . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50

2
 11.1.2. Elementos finitos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
11.2. Problema en 2-D y 3-D . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
11.2.1. Solución por relajación . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
11.2.2. Solución directa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
11.2.3. Solución por elementos finitos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57

A. Convención de las sumatorias de Einstein 59

B. Funciones Delta 60

3
1. Electrostática
Cavendish (∼ 1770) y Coulomb (∼ 1780) propusieron que la fuerza producido sobre un “paquete” de carga
q1 debido a un “paquete” de carga q2 satisface la fuerza vectorial central
q1 q 2
F1,2 ∼ (r1 − r2 ) = q1 E2 (r1 ).
| r1 − r2 |3

Charles de Coulomb (1736-1806) Henry Cavensidh (1731-1810)

Sacado de la Wikipedia, Frontispiece of The Life of the Hon. Henry
Sacado de la Wikipedia, https://www.photo.rmn.fr/archive/99-003694-
Cavendish
2C6NU0X9TUSI.html

Esta dependencia fue descubierta por los cuidadosos experimentos de Cavendish (∼ 1770) y Coulomb
(∼ 1780), pero requiere unos comentarios (ver más adelante). En un sentido esto define estas “extrañas
cantidades” denominadas cargas, que pueden tener dos signos (positivo o negativo). En el sistema CGS la
constante de proporcionalidad es k = 1, mientras que en MKS es

k = 9 · 109 kgm3 C −2 s−2

Como en general estamos interesados en describir distribuciones macroscópicas de cargas, se acostumbra a

definir una función continua, la densidad de carga ρ, para escribir la dependencia espacial de las cargas.
Esta función continua es un promedio espacial de la densidad de carga
P
qi Nq
ρ∼ i ∼
V V
en una escala intermedia (entro lo microscópico y lo macroscópico) particular, si está escala existe. Esto asu-
me que hay una escala espacial (y temporal) donde tal promedio es significativo, i.e. donde las desviaciones
son pequeñas en todo tiempo.
Esto permite re-escribir la fuerza entre densidades de carga como
Z Z
ρ1 (x1 )ρ2 (x2 )
F1,2 ∼ dx1 dx32
3
(x1 − x2 )
| x1 − x2 |3

4
 Intermedia

Micro Macro

λT
Figura 1: (a) Escalas espaciales y temporales. (b) Interpolar una función más suave.

2. Campo eléctrico y potencial eléctrico

El campo eléctrico está definido como una fuerza por unidad de carga, independiente de la carga de prueba,
que se puede expresar como
x−y 3
Z
E(x) = ρ(y) dy
| x − y |3
donde, hemos asumido la idea de superposición.

Aquı́ es importante notar que esta definición se aplica al campo eléctrico producido fuera de un paquete de
carga descrito en el continuo por ρ(x), osea no podemos tomar el limite y → x porque la integral diverge
a menos que ρ(y) → 0 también, o que hayan ciertas situaciones de simetrı́a. Notemos que en el limite
discreto esto no es relevante ya que en general estamos interesados en calcular el campo eléctrico lejos de las
cargas puntuales. Esto sera importante cuanto calculemos la energı́a de una distribución de carga, donde
la descripción continua tiene ciertas falencias. En el plasma este problema del paso al continuo se resuelve
asumiendo que las fluctuaciones se hacen pequeñas, como veremos después (descripción tipo Vlasov)

Podemos definir fácilmente el potencial escalar

Z
ρ(y) 3
E(x) = −∇x d y = −∇Φ
|x−y |
donde
Z
ρ(y) 3
Φ= dy
|x−y |
Por lo tanto

∇×E=0
El inverso también es cierto. De hecho el ”teorema de Helmholtz” estipula que cualquier campo vectorial
suficientemente suave, se puede escribir como la suma

5
 B = −∇ψ + ∇ × A
Si lı́mr→∞ B → 0 mas rapido que 1/r, entonces

1 R ∇y · B(y) 3
ψ(r) = dy
4π |x − y|
1 R ∇y × B(y) 3
A(r) = dy
4π |x − y|
Por lo tanto esta claro que el campo vectorial esta completamente definido como B = B1 + B2 si somos
capaces de encontrar 4 campos vectoriales tal que

∇·B = ∇ · B1 = ρ
∇ × B = ∇ × B2 = C
De la misma forma vemos que si ∇ × B = 0, entonces

B = −∇Ψ.

Entonces, qué es ∇ · E en el caso de la electrostatica? Hay muchas formas de proceder, pero tomaremos
una ruta simple. Tomemos una carga q y una unidad de superficie dS como en la figura 2.

n E

dS
δΩ

Q
Figura 2: Superficie orientada infinitesimal.

El campo eléctrico satisface

cos θ
E · ndS = q dS = qdΩ
r2
independiente de la forma de la superficie S. Si integramos sobre la superficie, tenemos
I Z
E · n dS = 4π ρ d3 x
S V

6
para la carga contenida. Este resultado puede ser usado para calcular campos eléctricos para casos altamente
simétricos. Usando el teorema de calculo vectorial obtenemos
Z
(∇ · E − 4πρ) d3 x = 0
V
y por lo tanto obtenemos las ecuaciones de Maxwell para la electrostática,

∇·E = 4πρ
∇×E = 0
Dado que ∇ × E = 0 existe un potencial E = −∇Ψ, que satisface

∇2 Ψ = −4πρ
que es la ecuación de Poisson. Esta ecuación es usada para resolver casos más generales, ya sea utilizando
el campo eléctrico o la forma potencial escalar. Representa un problema con condiciones de borde para el
potencial o campo eléctrico, definido en una superficie.

Notemos que es claro que si ponemos una carga unitaria en una posición con condiciones de borde naturales
en infinito, obtenemos (mirar el apéndice B)
1
∇2 = −4πδ (3) (x − y)
|x−y |
este es un ejemplo especı́fico de una función de Green como veremos más adelante. Con este resultado, el
cual es derivado en el apéndice, podemos calcular directamente
R ρ(y) 3
∇ · E = −∇2x dy
|x−y |
R 1
= − dy 3 ρ(y) ∇2x
|x−y |
R
= 4π dy 3 ρ(y) δ (3) (x − y)

= 4πρ(x)

3. Comentarios sobre la ley de la fuerza

En un sentido la ley de fuerza define esta “cantidad extraña” denominada carga, que puede tener dos signos.
En el sistema CGS la constante de proporcionalidad es k = 1, mientras en MKS es k = 9 · 109 .

7
3.1. Variación para una ley de fuerza cuadrática
Notemos que gracias a la exactitud de los instrumentos originales, fue postulado que la fuerza puede ser
escrita también en su forma generalizada
1
E∼
r2+2ε
Con esta forma, podemos calcular el campo eléctrico, o la fuerza que sentirı́a una carga, dentro de una
esfera como
1 ε 1 ε
   
 −   − 
π
(1 − 2z + z 2 ) 2 2 (1 + 2z + z 2 ) 2 2
Z
sin θ
= −
0 (1 + z2 + 2z cos θ)(1+ε)/2 z(−1 + ε) z(−1 + ε)
lo cual implica que E(ε = 0) = 0 dentro de la esfera, y da una forma de encontrar |ε| < 10−16 en forma
experimental.

3.2. Masa del fotón

En mecánica cuántica una fuerza que usa una partı́cula virtual para actuar a una distancia satisface una
ecuación de potencial escalar, la ecuación de Yukawa, dada por

h m c i2 4π
γ
∇2 Ψ − Ψ = ρ
h̄ c
mγ c e−µr
µ= → Ψ∼
h̄ r
Esta ecuación aparece en muchos contextos, desde mecánica cuántica, fı́sica de plasmas, astrofı́sica, etc.
Nos dice que si la masa del fotón tiende a cero, recuperaremos nuestro potencial 1/r. Se encuentra que la
masa del fotón, usando medidas satelitales, es a lo menos mγ < 10−48 gm.

4. Condiciones de Borde
Hasta ahora hemos discutido la situación de un volumen infinito, y para ese caso podemos utilizar la integral
sobre la densidad de carga dada anteriormente. En muchos casos estamos interesados en situaciones en un
espacio restringidos y con condiciones de borde. Por ejemplo con conductores en los bordes se inducen
cargas en las superficies las cuales en principio no las conocemos. En este caso es mucho más fácil tratar de
resolver la ecuación de Poisson como un problema con condiciones de borde

∇2 Ψ = −4πρ
directamente.

8
Hay dos tipos de condiciones de borde que aparecerán naturalmente en la construcción de la solución general
de la ecuación de Poisson, estas son las cargas de la superficie, y las capas del dipolo. Ellos corresponden a
diferentes condiciones de borde.

Tomemos la ecuación de Poisson cerca de un borde y tratemos de deducir la carga inducida σ en la superficie.
Usando la ecuación de Poisson

∇ · E = −∇2 Ψ = 4πρ
y una caja infinitesimal de tapa y fondo dS, y lado dl → 0 (ver Fig. 3). Usando el teorema de Gauss para
la integral de volumen obtenemos
H R
Ω
∇ · E · dx3 = 4π Ω ρ d3 x
H R
δΩ
E · n dS = 4π Ω
ρ d3 x
R
(E2 − E1 ) · n dS ≈ dS lı́mdh→0 ρdh
con lo cual obtenemos una condición de borde para los campos

(E2 − E1 ) · n = 4πσ
en el limite dh → 0, con
Z
σ = lı́m ρdh
dh→0

Por lo tanto, la carga de la superficie genera una discontinuidad en el campo eléctrico normal debido a una
carga de la superficie. Si asumimos un borde conductor (las cargas son libres de moverse en respuesta a un
campo eléctrico), podemos obtener la carga superficial desde el potencial en la superficie
−∇Ψ · n
σ=
4π
Debemos tener cuidado en definir el normal n. El problema de condiciones con conductores, por ejemplo,
es que en general no conocemos σ de antemano, y debe encontrase a partir de la solución completa.

Para el caso de campo tangencial, tomemos un rectángulo perpendicular, como se muestra en la Fig. 3b.
Integrando sobre esta superficie que atraviesa la interfase tenemos
H
(∇ × E) · n dS = 0
H
E · d` = 0
(E2 − E1 ) × nd` ≈ 0
con lo cual obtenemos la segunda condición de borde

(E2 − E1 ) × n = 0

9
 E2 E2
n
dl
∆S dh
dh

E1 n E1
Figura 3: Condición de borde en la superficie con dh → 0.

En muchas situaciones esta condición de borde es equivalente a que el potencial eléctrico es continuo en
interfases, aunque esto no es correcto en todas las situaciones.

10
5. Ecuación de Laplace y funciones ortogonales
Estamos interesados en resolver la ecuación de Laplace

∇2 Ψ(x) = 0 x ∈ Ω
sujeta a la condición de borde

Ψ(x)|δΩ
Este problema generalmente se resuelve utilizando una base ortonormal
Z
Ψ∗n (x)Ψm (x) dx = δn,m
Ω
y completa
X
Ψ∗n (x)Ψn (y) = δ(x − y).
n

Con esto cualquier funcion f (x) se puede expandir en términos de esta base
X
f (x) = An Ψn (x)
n

y podemos encontrar los coeficientes An usando el proyector

Z
dx Ψ∗m (x) dx
Ω
sobre la ecuación de arriba, obteniendo
Z X Z X
∗
dx Ψm (x)f (x) = An dxΨ∗m (x)Ψn (x) = An δn,m = Am .
Ω n Ω n

En el lı́mite infinito, cuando no hay discretizacion de los indices de la base ortonormal, las sumatorias
generalmente convergen a integrales.

las bases ortonormales (completas) son muy útiles si además son también soluciones de la ecuación de
interés. Luego, la expansión es una expansión natural, pero esto no es requerido.

La forma natural de construir bases ortonormales es utilizando una separación de variables, y a través del
teorema de Sturn-Liouville obtenemos una ortogonalización y completitud de las funciones que permiten
expandir cualquier condición de borde general usando el proyector para invertir relaciones como la de arriba.

11
5.1. Bases rectangulares
En una base rectangular los operadores diferenciales para un campo escalar Ψ son

2 ∂ 2Ψ ∂ 2Ψ ∂ 2Ψ
∇Ψ = + +
∂x2 ∂y 2 ∂z 2

∂Ψ ∂Ψ ∂Ψ
∇Ψ = x̂ + ŷ + ẑ
∂x ∂y ∂z
Para un campo vectorial A la divergencia es
∂Ψ ∂Ψ ∂Ψ
∇·A= + + ,
∂x ∂y ∂z
el rotor es

x̂ ŷ ẑ
∂ ∂ ∂
∇×A=
∂x ∂y ∂z
Ax Ay Az
y el Laplaciano es

∇2 A|x = ∇2 Ax
∇2 A|y = ∇2 Ay
∇2 A|z = ∇2 Az

La ecuación de Laplace se puede resolver por separación de variables.

∂ 2Ψ ∂ 2Ψ ∂ 2Ψ
+ + =0 → Ψ ∼ U1 (x)U2 (y)U3 (z)
∂x2 ∂y 2 ∂z 2
para lo cual obtenemos

d2 U1
= −a2 U1 → U1 = Aa eiax + Ba e−iax
dx2
d2 U2
= −b2 U2 → U2 = Ab eiby + Bb e−iby
dy 2
d2 U3
= −c2 U3 → U3 = Ac eicz + Bc e−icz
dz 2
con la restricción

c 2 + a2 + b 2 = 0

12
Por lo tanto tenemos una sumatoria en a y b (con c restringido por la relación c2 + a2 + b2 = 0). En generar
la restricción en a y b dependen de las condiciones de borde de interés.

Aparte de las soluciones exponenciales podemos utilizar soluciones oscilatorias en el intervalo [0, L] (esto
implica condiciones de borde periódicas en x = L y x = 0. Las condiciones de ortonormalidad son

∞ 
1 X  πmx   πmx 
f (x) = Ao + Am sin + Bm cos
2 m=1
L L

con
Z L  πmx   πnx  Z L  πmx   πnx  L
sin sin dx = cos cos dx = δn,m
0 L L 0 L L 2
y
Z L    
2πmx 2πnx
sin cos dx = 0
0 L L

En el lı́mite infinito L → ∞ (aquı́ es mejor mover el intervalo a L/2 ≤ x ≤ L/2 y después tomar el lı́mite),
con
2πm
→k
L
tenemos
1 R∞ 1 R∞
f (x) = √ −∞
A(k)eikx dx ↔ A(k) = √ −∞
f (x)eikx dx
2π 2π
1 R ∞ i(k−k0 )x 1 R ∞ ik(x−x0 )
√ e dx = δ(k − k 0 ) ↔ √ e dk = δ(x − x0 )
2π −∞ 2π −∞

φ =+V
Problema: Calcule el potencial dentro de
φ =0
una caja cubica de tamaño L, que tiene to-
das sus caras a potencial 0, excepto la cara
de arriba que esta a potencial V y la de abajo
a potencial -V.

φ =−V

13
5.2. Bases polares
En base polares tenemos

1 ∂ 2Ψ
 
2 1 ∂ ∂Ψ
∇Ψ = ρ + 2 2
ρ ∂ρ ∂ρ ρ ∂φ

1 ∂Ψ 1 ∂Ψ
∇Ψ = ρ̂ + φ̂
ρ ∂ρ ρ ∂φ
La ecuación de Laplace se puede resolver por separación de variables. Estos son problemas en los cuales la
solución o las fuentes no dependen en una de las variables, e.g z

1 ∂ 2Ψ
 
1 ∂ ∂Ψ
ρ + 2 2 = 0 → Ψ ∼ U1 (ρ)U2 (φ)
ρ ∂ρ ∂ρ ρ ∂φ
Para el caso m 6= 0 tenemos
 
d dU1
ρ ρ = m2 U1 → U1 = Am ρm + Bm ρ−m
dρ dρ
d2 U2
2
= −m2 U2 → U2 = Cm cos mφ + Dm sin mφ
dφ
y para el caso m = 0 tenemos
 
d dU1
ρ ρ =0 → U1 = A0 + B0 ln ρ
dρ dρ
d2 U 2
=0 → U2 = C0 + D0 φ
dφ2
La sumatoria es para todos lo m reales, pero generalmente la restricción en m depende en las condiciones
de borde de interés. Por ejemplo, m es entero si el problema incluye φ ∈ [0, 2π].

Problema: Calcule el potencial dentro de un

cilindro infinito que tiene un potencial Ψ = ψ
V cos φ sobre su superficie.

5.3. Bases esféricas

En una base rectangular los operadores diferenciales para un campo escalar Ψ son

14
 1 ∂ 2 (rΨ) ∂ 2Ψ
 
2 1 ∂ ∂Ψ 1
∇Ψ = + sin θ +
r ∂r2 r2 sin θ ∂θ ∂θ r2 sin2 θ ∂φ2

∂Ψ 1 ∂Ψ 1 ∂Ψ
∇Ψ = r̂ + θ̂ + φ̂
∂r r ∂θ r sin θ ∂φ
Notemos que también tenemos

1 ∂ 2 (rΨ)
 
1 ∂ 2 ∂Ψ
= 2 r
r ∂r2 r ∂r ∂r
Para un campo vectorial A la divergencia es
 
1 ∂ 2
∂ ∂Aφ
∇·A= 2 sin θ r Ar + r (sin θAθ ) + r ,
r sin θ ∂r ∂θ ∂φ
el rotor es

r̂ rθ̂ r sin θφ̂

1 ∂ ∂ ∂
∇×A=
r2 sin θ ∂r ∂θ ∂φ
Ar rAθ r sin θAφ
y el Laplaciano es
   
2 2 2 2 ∂ 2 cos θ 2 ∂
∇ A|r = ∇ Ar − 2 Ar + − 2 − Aθ + − 2 Aφ
r r ∂θ r2 sin θ r sin θ ∂φ
1 2 ∂Ar 2 cos θ ∂Aφ
∇2 A|θ = ∇2 Aθ − Aθ + 2 − 2 2
r2
sin θ r ∂θ r sin θ ∂φ
1 2 ∂Ar 2 cos θ ∂Aθ
∇2 A|φ = ∇2 Aφ − 2 2 Aφ + 2 − 2 2
r sin θ r sin θ ∂φ r sin θ ∂φ

La ecuación de Laplace se puede resolver por separación de variables.

1 ∂ 2 (rΨ) ∂ 2Ψ
 
1 ∂ ∂Ψ 1 U1 (r)
+ 2 sin θ + 2 2 =0 → Ψ∼ U2 (θ)U3 (φ)
r ∂r2 r sin θ ∂θ ∂θ r sin θ ∂φ2 r
con

15
 d2 U1 l(l + 1)
− U1 = 0 → U1 = A` r`+1 + B` r−`
dr2 r2
m2
   
1 d dU2
sin θ + `(` + 1) − U2 = 0 → U2 = C`m P`m (cos θ) + D`m Qm
` (cos θ)
sin θ dθ dθ sin2 θ

d2 U 3
+ m2 U3 = 0 → U3 = e±imφ
dφ2
donde P`m (x) y Qm
` (x) son las funciones de Legendre asociadas.

5.3.1. Funciones de Legendre

La ecuación de Legendre regular, usando x = cos θ, es

d 2 P`
 
d 2 dP` dP`
(1 − x ) + `(` + 1)P` = (1 − x2 ) 2 − 2x + `(` + 1)P` = 0 .
dx dx dx dx

Si hacemos una expansión en series

∞
X
f (x) = an xn+k
n=0

obtenemos
∞
X ∞
X
(k + n)(k + n − 1)an xk+n−2 + [`(` + 1) − 2(k + n) − (k + n)(k + n − 1)]an xk+n = 0
n=0 n=0

por lo que la ecuación indicial es

k(k − 1) = 0 ,
que tiene solución k = 0 y k = 1. Con k = 0, a0 = 1 y a1 = 0 tenemos

` + n + 1)(` − n)
an+2 = − an
(n + 1)(n + 2)
la cual no converge para x = ±1, por lo tanto exigimos que ` sea un entero par y truncamos la serie
generando el polinomio p` (x). Con k = 1, a0 = 1 y a1 = 0 tenemos

` + n + 2)(` − n − 1)
an+2 = − an
(n + 2)(n + 3)
la cual no converge para x = ±1, por lo tanto exigimos que ` sea un entero impar y truncamos la serie
generando el polinomio q` (x). Estas dos soluciones producen P` (x) para valores pares e impares. Para el
caso de p` (x) con ` impar y q` (x) con ` construimos la segunda solución de Legendre Q` (x). En este caso
no podemos truncar la serie, con lo cual diverge para x = ±1.

16
La función de Legendre P` (x) se puede escribir como

1 d` 2
P` (x) = (x − 1)` ,
2` `! dx`
la cual satisface la relación de ortogonalidad
Z 1
2
P`¯(x)P` (x) dx = δ¯ .
−1 2` + 1 `,`
Los primeros polinomios de Legendre son

P0 = 1
P1 = x
1
P2 = (3x2 − 1)
2
Además tenemos la relación

dP`+1 (x) dP`−1 (x)

(2` + 1)P` (x) = −
dx dx

5.3.2. Funciones de Legendre asociadas

La ecuación asociada de Legendre

d 2 P` m2
 
2 dP`
(1 − x ) 2 − 2x + `(` + 1) − P` = 0.
dx dx 1 − x2

Es fácil darse cuenta que

dm 1 2 m/2 d
m+`
P`m (x) = (1 − x2 )m/2 m
P ` (x) = `
(1 − x ) m+`
(x2 − 1)` ,
dx 2 `! dx
satisface la ecuación asociada de Legendre, y por lo tanto corresponden a las soluciones regulares. Estas
soluciones satisfacen la relación de ortogonalidad
Z 1
2 `+m
P`¯m (x)P`m (x) dx = δ¯ .
−1 2` + 1 ` − m `,`
Es importante notar que es un polinomio, con una expansión finita en x.
Para construir una solución completa, debemos tomar la sumatoria sobre todos los m y `. La restricción en
m y ` dependen de las condiciones de borde de interés. Por ejemplo, m es entero si φ ∈ [0, 2π]. Además, si
nos restringimos al rango −1 ≤ cos θ ≤ 1, entonces nos quedamos con los P`m donde `, m son enteros que
satisfacen

`≥0
−` ≤ m ≤ `

17
Para esta situación podemos definir los armónicos esféricos
s
1 2` + 1 (` − m)! m
Y`,m (θ, φ) = C`m P`m (cos θ) √ eimφ = (−1)m P (cos θ)eimφ
2π 4π (` + m)! `
que satisfacen

1 ∂ 2Ψ
   
1 1 1 ∂ ∂Ψ `(` + 1)
2
(Lθ + Lφ )Y`,m = 2 sin θ + 2 2
Y`,m = − Y`,m
r r sin θ ∂θ ∂θ sin θ ∂φ r2
con las condiciones de ortogonalidad

Z 2π Z π
dφ sin θ dθ Y`,¯∗m̄ (θ, φ)Y`,m (θ, φ) = δ`,`
¯ δm̄,m
0 0
∞ m=`
X X
∗
Y`,m (θ0 , φ0 )Y`,m (θ, φ) = δ(φ − φ0 )δ(cos θ − cos θ0 )
`=0 m=−`

Los primeros armónicos esféricos son

r
1 15
Y0,0 =√ Y2,2 = sin2 θe2iφ
4π 32π
r r
3 15
Y1,1 = − sin θeiφ Y2,1 = − sin θ cos θeiφ
8π 8π
r r
3 5
Y1,0 = cos θ Y2,0 = (3 cos2 θ − 1)
4π 16π
Notemos que esto tiene sentido solo si −` ≤ m ≤ ` si ` es entero, ya que P` (x) es un polinomio de orden `.
Es muy útil derivar

(` − m)! m
P`−m (x) = (−1)m P (x)
(` + m)! `
por lo tanto

18
 Y`,−m (θ, φ) = C`,m P`−m (cos θ)e−imφ
s
2` + 1 (` + m)! −m
= (−1)m P (cos θ)e−imφ
4π (` − m)! `
s  
2` + 1 (` + m)! m (` − m)! m
m
 ∗
= (−1) (−1) P` (cos θ) eimφ
4π (` − m)! (` + m)!
" s #∗
2` + 1 (` + m)!
= (−1)m (−1)m P m (cos θ)eimφ
4π (` − m)! `
∗
= (−1)m Y`,m (θ, φ)
o lo que es equivalente

C`,m P`m (cos θ) = (−1)m C`,−m P`−m (cos θ)

La solución general al potencial se puede escribir

∞ m=`
X  
X
` B`,m
Ψ= A`,m r + `+1 Y`,m (θ, φ)
`=0 m=−`
r

En el caso especial de la simetrı́a azimutal tenemos sólo el caso m = 0, lo cual significa que la solución
general es
∞  
X B`,m
`
Ψ= A`,m r + `+1 P` (cos θ)
`=0
r
ya que
r
2` + 1
Y`,0 (θ, φ) = P` (cos θ).
4π

Problema: Encuentre el potencial al interior de un casquete esférico donde V (θ) = Vo para 0 ≤ θ ≤ α y

cero para α ≤ θ ≤ π.

Por la simetrı́a sabemos que el potencial se puede escribir como

∞
X  r `
Ψ= A` P` (cos θ)
`=0
a
Por lo tanto en el casquete tenemos

19
 ∞
X
V (θ) = A` P` (cos θ)
`=0

Utilizando la ortogonalizacion, multiplicamos por Pm (cos θ) e integramos 0 ≤ θ ≤ π.

2m + 1 R π
Am = 0
V (θ)Pm (θ) sin θdθ
2
2m + 1 R α
= Vo 0
Pm (θ) sin θdθ
2
2m + 1 R 1
= Vo P (x)dx
cos α m
2
Vo
= [Pm−1 (cos α) − Pm+1 (cos α)]
2
y para el caso m = 0 tenemos
Vo
A0 = (1 − cos α)
2
Por lo tanto la solución es
" ∞
#
Vo X  r `
Ψ= (1 − cos α) + (P`−1 (cos α) − P`+1 (cos α)) P` (cos θ)
2 `=1
a
De aquı́ podemos calcular el campo eléctrico en el origen, y obtener
 
Vo 3 2

E= (cos α − 1) cos θr̂ − sin θθ̂
2 2
y claramente tiende a cero cuando α → π.

En general no es difı́cil de demostrar que

∞ `
1 X r<
= P (cos γ)
`+1 `
→ cos γ = cos θ cos θ0 + cos(φ − φ0 ) sin θ sin θ0
| x − x0 | `=0
r>
∞ ` `
1 X 1 r< X
∗
0
= 4π `+1
Y`,m (θ0 , φ0 )Y`,m (θ, φ).
|x−x | `=0
2` + 1 r> m=−`

20
 +V
b
Problema: Encuentre el potencial entre dos
esferas de radio a < r < b que están a poten-
cial Ψ(r = a) = −V y Ψ(r = b) = V
a
−V

5.4. Bases cilı́ndricas

En bases cilı́ndricas los operadores diferenciales para un campo escalar Ψ son

1 ∂ 2Ψ ∂ 2Ψ
 
2 1 ∂ ∂Ψ
∇Ψ = ρ + 2 2 +
ρ ∂ρ ∂ρ ρ ∂φ ∂z 2

∂Ψ 1 ∂Ψ ∂Ψ
∇Ψ = ρ̂ + φ̂ + ẑ
∂ρ ρ ∂φ ∂z
Para un campo vectorial A la divergencia es
1 ∂ 1 ∂Aφ ∂Az
∇·A= (ρAρ ) + + ,
ρ ∂ρ ρ ∂φ ∂z
el rotor es

ρ̂ ρφ̂ ẑ
1 ∂ ∂ ∂
∇×A=
ρ ∂ρ ∂φ ∂z
Aρ ρAφ Az
y el Laplaciano es
1 2 ∂Aφ
∇2 A|ρ = ∇2 Aρ − 2
Aρ − 2
ρ ρ ∂φ
1 2 ∂Aρ
∇2 A|φ = ∇2 Aφ − 2 Aφ + 2
ρ ρ ∂φ
∇2 A|z = ∇2 Az

La ecuación de Laplace se puede resolver por separación de variables:

1 ∂ 2Ψ ∂ 2Ψ
 
1 ∂ ∂Ψ
ρ + 2 2 + = 0 → Ψ ∼ U1 (ρ)U2 (φ)U3 (z)
ρ ∂ρ ∂ρ ρ ∂φ ∂z 2

21
la solución por separación de variable es

d2 U1 1 dU1 m2
 
2
+ + k − 2 U1 = 0 → U1 = Ak,m Jm (kρ) + Bk,m Nm (kρ)
dρ2 ρ dρ ρ
d2 U2
+ m2 U2 = 0 → U2 = Cm cos mφ + Dm sin mφ
dφ2
d2 U3
− k 2 U3 = 0 → U3 = Dk e+kz + Ek e−kz
dz 2
Las funciones Jm (x) y Nm (x) son conocidas como las funciones de Bessel y Neumann respectivamente. Para
el caso m = 0 podemos tratar de construir una expansión en series
∞
X
J0 (x) = an xk+n
n=0
que genera

∞
X
k k+1
a0 [k(k − 1) + k]x + a1 [(k + 1)k + (k + 1)]x + [an (n + k)(n + k − 1) + (n + k) + an−2 ]xk+n = 0.
n=2

La ecuación indicial es r(r − 1) + r = 0 por lo que tenemos una solución desgenerada r = 0. Usando r = 0,
a0 = 0 y a1 = 0, con la ecuación de recursión
1
an+2 = − an ,
(n + 2 + k)2
podemos construir la expansión en serie
" ∞
#
X (−1)n x2n
J0 (x) = a0 1 +
n=1
22n (n!)2
la cual converge para todo x. De la misma forma podemos construir las funciones Jm (x). Las funciones de
Neumann divergen para x → 0 y no son consideradas para problemas en que involucren el origen. Es útil
en problemas de borde resolver el número infinito de raı́ces xm,n de

Jm (xm,n ) = 0 n = 1, 2, 3, . . .
con lo cual se puede demostrar que

a
ρ  a2
Z  ρ 
ρ dρJm xm,r Jm xm,s = [Jm+1 (xm,r )]2 δr,s
0 a a 2
y el sistema de raı́ces xm,n forman un sistema completo en el intervalo [0, a]. De hecho esto restringe los
valores de k para sistemas discretos (pero infinitos)
xm,n
km,n = n = 1, 2, 3, . . .
a
22
para cada m. Si el rango de ρ se hace infinito, entonces requerimos de una representación continua, similar
a las series infinitas de Fourier, donde la suma discreta en k se convierte en una integral

Z ∞
dk e±kz (Cm (k)Jm (kρ) + Dm (k)Nm (kρ))(Am (k) sin mφ + Bm (k) cos mφ)
0
Z ∞
1
x dx Jm (kx)Jm (k 0 x) = δ(k 0 − k)
0 k
donde el ± depende de la restricción de z.

+V
Problema: Encuentre el potencial dentro de
un cilindro de largo L, donde la tapa superior
esta a un potencial V, la tapa inferior esta a
un potencial -V y el resto a potencial cero.

−V
En muchas situaciones es útil tomar una conjunción mas apropiada para ciertas condiciones de borde en el
infinito en termino de funciones de Hankel

(1) (2)
Hm (x) = Jm (x) + iNm (x) Hm (x) = Jm (x) − iNm (x)
que son útiles para tener condiciones de salida (radiación) en el infinito. En otras situaciones estamos
interesados en que la dependencia de z sea oscilatoria, lo que implica que tenemos que hacer la transformación

k → ik
y por lo tanto definimos las funciones de Bessel modificadas
π m+1 (1)
Im (kρ) = i−m Jm (ikρ) Km (kρ) =
i Hm (ikρ)
2
La elección de la normalización es para que estas funciones sean reales para ρ real.

La expansiones asintóticas

23
 r
2  mπ π 
lı́mx→∞ Jm (x) → cos x − −
πx 2 4
r
2  mπ π 
lı́mx→∞ Nm (x) → sin x − −
πx 2 4
m
x
lı́mx→0 Jm (x) → m
2 m!
 m
(n − 1)! 2
lı́mx→0 Nm (x) →
π x
lı́mx→0 J0 (x) → 1
2
lı́mx→0 N0 (x) → (ln x + γ − ln 2)
π
son útiles cuando tenemos que forzar ciertas condiciones de borde en el infinito.

Veremos mas adelante la importancia del Wronskiano, que en ecuaciones diferenciales determinan si dos
funciones son linealmente independientes

y1 y2
W (y1 , y2 ) =
y10 y20
Para el caso de ecuaciones diferenciales del tipo

y 00 + p1 (x)y 0 + p2 (x)y(x) = 0
Dado que y1 e y2 satisfacen esta ecuación, se puede probar que

W 0 = −p1 W
y por lo tanto
 Z x 
W (y1 , y2 ) = C exp − p1 (t)dt

Para el caso de

y (n) + p1 (x)y (n−1) + ... + pn (x)y = 0

tenemos que
 Z x 
W (y1 , y2 , ..., yn ) = C exp − p1 (t)dt

Las expresiones asintóticas se puede usar para evaluar la constante c, tal que

24
 0 0 2
Jm Nm − Jm Nm =
πx

0 0 1
Im Km − Im Km = −
x

Problema: Encuentre el potencial dentro de un cilindro

5.5. Proyectores en general

En general tenemos la relación de completitud
Z
w(x)Ψ∗n (x)Ψm (x)dx = Cm δn,m .

Podemos expandir cualquier función en este set completo

X
f (x) = an Ψn (x)
n

y tenemos que encontrar an . Por lo cual definimos el proyector

Z
w(x)Ψ∗m (x)

y lo aplicamos a ambos lados de la expresióon de arriba. Esto es

Z XZ
∗
w(x)Ψm (x)f (x)dx = w(x)Ψ∗m Ψn (x)dx,
n

por lo cual
Z
am C m = w(x)Ψ∗m (x)f (x)dx.

Si es necesario, esta expresión se puede evaluar en forma numérica.

25
6. La Ecuación de Poisson y la función de Green
Estamos interesados en resolver la ecuación de Poisson

∇2 Ψ(x) = −4πρ(x) x ∈ Ω
sujeta a la condición de borde

Ψ(x)|δΩ
La función de green en un dominio Ω corresponde al potencial producido por una carga puntual en la
posición y,

ρ(y) = δ (3) (x − y) x, y ∈ Ω
La función de Green G(x, y) satisface

∇2 G(x, y) = −4πδ (3) (x − y) x, y ∈ Ω (1)

con lo cual podemos notar que

G(x, y) = G(y, x)
esta propiedad de simetrı́a será relevante en futuros estudios.

6.1. Función de Green en el espacio infinito

Ahora construiremos la función de Green para el espacio infinito.

6.1.1. En 3D
El Laplaciano es

1 ∂2 ∂ 2G
 
2 1 ∂ ∂G 1
∇ G= 2
(rG) + 2 sin θ + 2 2
r ∂r r sin θ ∂θ ∂θ r sin θ ∂φ2
y para el problema con simetrı́a esférica G(r) tenemos que

1 d2
(rG) = 0 r 6= 0
r dr2
y por lo tanto
1
G=
r
Ya hemos demostrado que efectivamente tenemos

∇2 G = −4πδ (3) (x)

Para el caso arbitrario tenemos

26
 1
G(x, y) =
|x − y|
es la función de green para el espacio infinito, ya que

∇2 G(x, y) = −4πδ (3) (x − y)

6.1.2. En 2D
El Laplaciano es

1 ∂ 2G
 
2 1 ∂ ∂G
∇ G= ρ + 2 2
ρ ∂ρ ∂ρ ρ ∂φ
y para el problema con simetrı́a cilı́ndrica G(ρ) tenemos que
 
1 d dG
ρ =0 ρ 6= 0
ρ dρ dρ
y por lo tanto
a
G = a ln ρ = ln(x2 + y 2 )
2

δ (2) (x) = δ(x1 )δ(x2 )

por lo tanto integrando en la superficie x1 = 0 ± 1 con x2 = 0 ± 2 y i → 0
Z + Z +  2
d G d2 G

dx1 dx2 + 2 = −4π
− − dx21 dx2
podemos encontrar que

2aπ = −4π
Entonces efectivamente obtenemos que

G = −2 ln ρ
Para el caso arbitrario tenemos

G(x, y) = −2 ln |x − y|
es la función de green para el espacio infinito, ya que

∇2 G(x, y) = −4πδ (2) (x − y)

Se entiendo que las funciones deltas y las ∇2 son en 2D.

27
6.1.3. En 1D
El Laplaciano es

2 d2 G
∇ G=
dx2
Tenemos que

d2 G
=0 x 6= 0
dx2
y por lo tanto

ax + b x<0
G=
cx + d x>0
Es fácil darse cuenta que G(x) deberı́a ser continua en x = 0, con lo cual b = d = 0. Podemos elegir
esta constante como cero, ya que una constante no cambia el campo eléctrico que el potencial genera. Su
derivada tiene que ser discontinua, ya que si integramos

d2 G
= −4πδ(x)
dx2
de x = 0 ±  con  → 0, tenemos que

a − c = −4π
Vemos que la solución es

G(x) = −2π|x|
Es fácil de encontrar que efectivamente tenemos

∇2 G = −4πδ(x)
Para el caso arbitrario tenemos

G(x, y) = −2π|x − y|
es la función de green para el espacio infinito, ya que

∇2 G(x, y) = −4πδ(x − y)
Se entiendo que las funciones deltas y las ∇2 son en 1D. Aca podemos observar que la funcion G es
continua, mientras que la derivada es discontinua, lo cual tiene sentido si pensmaos que vamos integrando
desde x → −∞ hacia x → ∞.

28
 f’(x) f(x)

x x

Figura 4: (a) Derivada discontinua, (b) Funcion continua

6.2. Función de Green para un volumen arbitrario

De ahora en adelante nos concentraremos en la función de Green para regiones Ω arbitrarias en 3D y
dejaremos los problemas en otras dimensiones para la tarea. Dado que
1
∇2 = −4πδ (3) (x − y)
|x − y|
podemos notar que la función de Green para un volumen arbitrario se puede escribir como
1
G(x, y) = + F (x, y) → ∇2 F (x, y) = 0
|x − y|
siempre y cuando ∇2 F (x, y) = 0. Esto significa que F (x, y) puede ser cualquier función que satisfaga la
ecuación de Laplace. Por lo tanto para tener una solución única dentro de x, yΩ, vamos a tener que forzar
condiciones de borde sobre G, de la misma forma que forzamos condiciones de borde sobre la ecuación de
Laplace para tener soluciones únicas.

Es mas o menos intuitivo, dado que el problema a resolver es lineal, que una vez que tenemos el potencial
para la carga puntual, es factible obtener el potencial por superposición de una distribución continua de
cargas. Veremos que primero tenemos que tener cuidado con las condiciones de borde tanto para Ψ como
para G. Comenzaremos con el teorema de Green. Estamos interesados en resolver el problema

∇2 Ψ(x) = −4πρ(x) Ψ|∂Ω

en xΩ con alguna condición de borde para Ψ|∂Ω . Dado que

∇ · (Ψ∇G) = Ψ∇2 G + ∇Ψ · ∇G
∇ · (G∇Ψ) = G∇2 Ψ + ∇Ψ · ∇G

podemos obtener

Ψ∇2 G − G∇2 Ψ = ∇ · (Ψ∇G) − ∇ · (G∇Ψ)

29
Luego al integrar en todo el espacio de interes
Z I  
2 2

3 ∂G ∂Ψ
Ψ∇y G − G∇y Ψ d y = Ψ −G dSy
Ω δΩ ∂ny ∂ny
donde la normal apunta hacia afuera del volumen Ω de interés. Ahora elegimos G como la función de Green
que satisface

∇2y G(x, y) = −4πδ (3) (x − y)

con condiciones de borde que especificaremos mas abajo. Obtenemos
Z I  
3 1 ∂ ∂
Ψ(x) = ρ(y) G(x, y) d y + G(x, y) Ψ(y) − Ψ(y) G(x, y) dSy
Ω 4π δΩ ∂ny ∂ny
donde los dos términos de superficie en la ecuación representan la carga superficial y la capa de un dipolo
en la superficie. El potencial y el campo afuera del volumen es cero debido a las discontinuidades. Las
condiciones de borde de Dirichlet se usan cuando Ψ está definido en el borde, y las condiciones borde de
Neumann cuando ∂Ψ/∂n está definido en el borde.

Todavı́a no hemos forzado las condiciones de borde sobre G (es una libertad dada por F ), las cuales
tendremos que forzar para obtener una solución única a la ecuación de Poisson. Asumamos dos soluciones

U = Ψ2 − Ψ1
tal que

∇2 Ψ1 (x) = −4πρ(x) Ψ1 |δΩ

y

∇2 Ψ2 (x) = −4πρ(x) Ψ2 |δΩ .

Ambas satisfacen la ecuación de Poisson con la misma densidad de carta y las mismas condiciones de borde.
Usando las igualdades dadas más arriba tenemos
Z I  
2 2

3 ∂U
U ∇ U + | ∇U | d x = U dS
Ω δΩ ∂n
por lo tanto existe una sola solución, esto es U = 0, si requerimos que
I  
∂U
U dS = 0
δΩ ∂n
lo que implica poner condiciones de borde sobre la superficie ∂Ω. Por lo tanto las dos condiciones naturales
son

especificar el potencial en el borde, Ψ|∂Ω

30
 ∂Ψ
especificar la carga en el borde,
∂n ∂Ω

Vemos inmediatamente que especificar los dos tipos de condiciones seria redundante, y generarı́a soluciones
no reales. Esto significa que dada un tipo de condición de borde para Ψ, tenemos que hacer cero una de las
integrales de superficie de arriba, y para eso utilizamos la flexibilidad en la condiciono de borde de G que
aun nos queda.

1. En el caso de condiciones de borde tipo Dirichlet, donde el potencial Ψ(x)|x∈δΩ esta definido en el
borde, forzamos G(x, y) =|y∈δΩ 0 para y en el borde tal que
Z I
13 ∂
Ψ(x) = ρ(y) G(x, y) d x − Ψ G(x, y) dSy
Ω 4π δΩ ∂ny

Esto es equivalente a encontrar la carga en la superficie.

2. En el caso de condiciones borde de Neumann, donde la densidad de carga ∂Ψ(x)/∂nx |x∈δΩ esta definida
en el borde, forzamos ∂G(x, y)/∂ny |y∈δΩ = −4π/S para y en el borde tal que
Z I
1 3 ∂Ψ
Ψ(x) =< Ψ >S + ρ(y) G(x, y) d y + G(x, y) dSy
Ω 4π δΩ ∂ny

La razón para esta última condición de borde es para satisfacer que la integral sobre el volumen Ω de
la ecuacion de green
∇2 G(x, y) = −4πδ (x) (x − y)

require
I
∂G(x, y)
dSy = −4π 6= 0
δΩ ∂ny

La solución general para Ψ se construye con el uso de funciones de Green, las cuales construiremos a partir
de

Método de las imágenes.

Expansión ortonormal de la función de Green.

Expansión de los valores propios del problema relacionado.

Solución numérica (ver apéndice)

31
Problema: Cual es el teorema de Green equivalente en 2D

Problema: Cual es el teorema de Green equivalente en 1D

6.3. Método de las imágenes

En particular, si situamos cualquier set de cargas fuera de Ω, entonces al superposición de los potenciales
satisfacen la condición ya nombrada

∇2 G(x, y) = −4πδ (3) (x − y) x, y ∈ Ω (2)

dentro de Ω, esto quiere decir que la contribución de las cargas exteriores a la densidad de carga interior es
cero. Si logramos con este set de cargas exteriores, aparte de la carga interior, satisfacer las condiciones de
borde, entonces tendremos una solución dentro de Ω. Pronto veremos que esta solución es la única solución
al problema de borde. Por lo tanto, en orden a resolver la ecuación de Poisson, dentro Ω, tenemos la libertad
de usar cualquier distribución de carga exterior al borde para satisfacer las condiciones de borde requeridas.
Este es el método llamado de las imágenes.

Problema: Una carga frente a un plano conductor (a tierra).

Es fácil hallar la carga imagen. Podemos entonces calcular la carga superficial y la fuerza en la carga real
q. Es interesante notar que la fuerza está también relacionada con
Z
F = 2πσ 2 dS

lo cual es una declaración general.

Problema: tiene ese signo?

Problema: Una carga dentro de una esfera conductora (a tierra).

El borde se fija en Ψ = 0. Ponemos otra carga exterior, en diferentes q2 y d2 , pero en la misma linea. Esta
carga q2 está en la realidad inducida sobre la superficie.
q1 q2
Ψ = Ψin + Ψout = +
|x − d1 k| |x − d2 k|
En Mathematica podemos escribir este problema como:

32
 q q1
Ψ[x− , y− , z− ] = p −p
(x − d)2 + y 2 + z 2 (x − d1)2 + y 2 + z 2

Podemos buscar la solucion

s = Simplif y[Solve[{Ψ[R, 0, 0] == 0, Ψ[−R, 0, 0] == 0}, {q1, d1}], d < R&&d > 0]

que, dado los parámetros y las condiciones de borde, obtenemos

qR R2
{q1 → , d1 → }
d d

El potencial es

Ψ[x , y , z ] = Ψ[x, y, z]/.s

lo que da

q qR
p − s
2 2
(−d + x) + y + z 2 2
R2
d − + x + y2 + z2
d

Podemos graficar los contornos como se muestran en Fig. 5, donde también se muestra las lineas de campo.

CountourPlot[Ψ[x, 0, z]/.{R → 1, d → 0,8, q → 1}, {x, −1, 1}, {z, −1, 1}, CountourShading→False,
AxesLabel→{“x”,“z”}, PlotPoints→50]

Podemos calcular la carga superficial inducida en la superficie de la esfera

1 ∂Ψ
σ= (3)
4π ∂r r=a
y la fuerza que siente la carga real.

Problem: El origen es un punto estable, o inestable?

Las lineas de campo se definen como las lineas que su tangente es paralela al campo electrico, esto es

33
 FieldyLines FieldyLines
1
1

0.75

0.5 0.5

0.25

x x
-1 -0.5 0.5 1 -1 -0.75 -0.5 -0.25 0.25 0.5 0.75 1

-0.25

-0.5 -0.5

-0.75

-1 -1

Figura 5: (a) Equipotencial, (b) Campo eléctrico dentro de una esfera conductora a tierra.

dr E
=
ds |E|
las cuales se pueden ver en la Fig 5, y claramente son ortogonales a la superficie, como deberia ser. Resul-
tados similares pueden ser obtenidos para una carga puntual fuera de una esfera conductora.

Problema: Esferas conductoras con hemisferios a diferente potencial. Hemos construido la función de Green
por el método de las imágenes para el potencial fuera de las esferas conductoras con Ψ = 0 en el borde de
la esfera. Necesitamos integrar con las condiciones de borde de Dirichlet. Note que esto es equivalente a la
densidad de carga en la superficie.

Problema: Una carga q fuera de una esfera aislada con carga Q, pero conductora. Una esfera a tierra
(equipotencial) tiene una carga q2 . Sumando Q − q2 , distribuido uniformemente sobre la superficie, alcan-
zamos una carga total Q. Esto significa, que es el equivalente a una carga puntual de magnitud Q − q2 en
el origen. El potencial y la fuerza pueden ser construidos por superposición
q1 q2 Q − q2
Ψ = Ψin + Ψout = + +
| x − d1 k | | x − d2 k | |x|
Lo mismo se puede hacer con una esfera conductora en la tierra en un potencial V especifico.

Problema: Una esfera conductora (a tierra) en un campo eléctrico uniforme asintótico. Es construido en
una forma similar, con 2 cargas de distinto signo a una distancia R >> a, y sus dos imágenes dentro de la
esfera. El limite tiene sentido cuando 2q/R2 → Eo .

Problema: Construya por imágenes la función de Green para un plano conductor en 1D

34
6.4. Expansión de la función de Green en bases ortonormales
Hemos visto en que en general, todo lo que necesitamos hacer es resolver la ecuación de Laplace, y esta es
la solución generada para construir la función de Green. Supongamos que queremos solucionar la ecuación
de Laplace dada una condición de borde. Arriba se explica la expansión de la solución en términos de las
bases ortonormales de la ecuación de Laplace en distintos sistemas de coordenadas:

Rectangulares.

Polares.

Esféricas.

Cilı́ndricas.

De las cuales podemos solucionar problemas de valor de borde en diferentes geometrı́as.

Para encontrar la función de Green en una geometrı́a dada, expandamos la función de Green en términos
de los modos normales de la ecuación de Laplace para esa geometrı́a, y resolvemos directamente para los
coeficientes desde

∇2 G(x, y) = −4πδ (3) (x − y) x, y ∈ Ω

con la condición de borde Dirichlet

G(x, y)|∂Ω = 0

6.4.1. G en coordenadas esféricas para espacio infinito por método directo

Podemos expandir la función de Green en términos de un set completo de soluciones de la ecuación de
Laplace en coordenadas esféricas
∞ m=`
X X
G(x, y) = A`,m (rx , ry , θy , φy ) Y`,m (θx , φx )
`=0 m=−`

donde
∞ m=`
(3) δ(rx − ry ) X X ∗
δ (x − y) = Y`,m (θy , φy ) Y`,m (θx , φx )
rx2 `=0 m=−`
∗
Dado que el problema es lineal y la función delta ya tiene el coeficiente para Y`,m (θy , φy ), por simetrı́a
podemos expandir G como
∞ m=−`
X X
∗
G(x, y) = g`,m (rx , ry )Y`,m (θy , φy )Y`,m (θx , φx )
`=0 m=−`

35
Podemos utilizar
1 ∂ 2 (rΨ) L (θ, φ)
∇2 Ψ = − Ψ,
r ∂r2 r2
donde
L (θ, φ)Y`,m (θ, φ) = −`(` + 1))Y`,m (θ, φ).
Usando la ortogonalidad de los armónicos esféricos, al introducir estas expresiones en la ecuación de Green
podemos encontrar

1 d2 (rx g`,m ) `(` + 1) 4π

2
− 2
g`,m = − 2 δ(rx − ry )
rx drx rx rx
de donde obtenemos la solución
"
−(`+1)
Arx` + Brx rx < ry
g`,m (rx , ry ) = 0 ` 0 −(`+1)
A rx + B rx rx > ry
Ahora necesitamos involucrar las condiciones de borde. Como estamos hablando del espacio infinito, tenemos
 
A(r0 )rx` rx < ry `
r<
g`,m (rx , ry ) =  1 = C
B(r0 ) `+1 rx > rx0

`+1
r r>
x

donde en la última expresión hemos hecho cumplir la continuidad en rx = ry , y la simetrı́a G(x, y) = G(y, x).
El valor de C es obtenido integrando desde

rx`+1
rx = ry − ε r< = rx r> = ry rx g`,m = C
ry`+1

hasta

ry`
rx = ry + ε r< = ry r> = rx rx g`,m =C `.
rx
Por lo tanto

rx =ry +ε
1 d2 (rx g` ) `(` + 1)
Z  
4π
rx drx × 2
− 2
g`,m = − 2 δ(rx − ry )
rx =ry −ε rx drx rx rx
rx =ry +ε
d 4π 4π
(rx g`,m ) =− → C=
drx rx =ry −ε ry 2` + 1

Por lo tanto la función de Green es

∞ m=−`
X 4π r`
1 X
< ∗
G(x, y) = 0
= `+1
Y`,m (θy , φy ) Y`,m (θx , φx )
|x−x | `=0 m=−`
2` + 1 r>
donde

36
 r> = max[rx , ry ] r< = min[rx , ry ]
En el caso de simetrı́a azimutal tenemos

∞
1 1 X
` −(`+1)


= p 2 = A` rx + B` r x P` (cos γ)
|x−y | rx + ry2 − 2rx ry cos γ `=0
∞ `
1 X r<
= `+1
P` (cos γ) → cos γ = cos θx cos θx + cos(φx − φy ) sin θx sin θy
|x−y | l=0
r>

La última expresión puede también ser encontrada evaluando x y y en el eje z (donde γ = 0), y expandiendo
1/|x − y|. En esta notación r> es el más largo entre x y y, y r< el opuesto. Es interesante notar que
somos capaces de encontrar toda la solución sólo evaluando la expresión anterior donde cos γ = 1 y luego
multiplicando cada término por P` (cos γ).

6.4.2. G en coordenadas esféricas entre dos esferas conductoras por método directo
Podemos expandir la función de Green en términos de los modos normales de la ecuación de Laplace en
coordenadas esféricas con radios de las esferas a < b. Nuevamente escribimos
∞ m=−`
X X
G(x, y) = A`,m (rx , ry , θy , φy ) Y`,m (θx , φx )
`=0 m=−`

donde
∞ m=`
(3) δ(rx − ry ) X X ∗
δ (x − y) = Y`,m (θy , φy ) Y`,m (θx , φx )
r2 `=0 m=−`
∗
Notemos ya que el problema es lineal y la función delta ya tiene el coeficiente Y`,m (θy , φy ). Por lo tanto
podemos expandir G como
∞ m=`
X X
G(x, y) = g`,m (r, r0 )Y`,m
∗
(θ0 , φ0 )Y`,m (θ, φ)
`=0 m=−`

donde obtenemos

1 d2 (rx g` ) `(` + 1) 4π
2
− 2
g`,m = − 2 δ(rx − ry )
rx drx rx rx
La solución es
"
−(`+1)
A`,m rx` + B`,m rx rx < ry
g`,m (rx , ry ) = ` −(`+1)
C`,m rx + D`,m rx rx > ry

37
donde hemos usado la ortogononalidad de los armónicos esféricos. Ahora necesitamos involucrar las condi-
ciones de borde. Tomamos el espacio vacı́o entre dos esferas conductoras de radios a < b y recordando que
rY está en el medio, tenemos

a2`+1
   
`
 Ā`,m (ry ) rx − r`+1 rx < ry 
a2`+1 `
  
` 1 r>
 = C̄`,m r< − r`+1 − 2`+1
g`,m (rx , ry ) = 
 x 
rx`
 
1 < r `+1
> b
B̄`,m (ry ) `+1 − 2`+1 rx > ry
rx b
donde en la última expresión hemos hecho cumplir la continuidad en rx = ry , y la simetrı́a G(x, y) = G(y, x).
El valor de C es obtenido integrando desde rx = ry − ε hasta rx = ry + ε la ecuación de arriba para g

rx =ry +ε
1 d2 (rx g`,m ) `(` + 1)
Z  
4π
rx drx × − g`,m = − 2 δ(rx − ry )
rx =ry −ε rx drx2 rx2 rx
rx =ry +ε
d(rx g`,m ) 4π 4π 1
=− → C=  
drx rx =ry −ε ry 2` + 1  a 2`+1
1−
b

y la función de Green es
∞ m=` `
a2`+1
  
X X 4π 1 ` 1 r> ∗
G(x, y) =   a 2`+1  r< − `+1 `+1
− 2`+1 Y`,m (θy , φy )Y`,m (θx , φy )
2` + 1 r< r> b
`=0 m=−` 1−
b

6.4.3. G en coordenadas esféricas fuera de una esfera conductora por método indirecto
Supongamos que sabemos la solución de G en una geometrı́a simple. Podemos usar esta solución

∇2 Ĝ(x, y) = −4πδ (3) (x − y)

para construir G para geometrı́as más complicadas

G(x, y) = Ĝ(x, y) + F (x, y)

con

∇2 F (x, y) = 0 F |δΩ = 0
Por ejemplo, tomemos la función de Green para el espacio infinito

∞ X ` `
1 X 4π r< ∗
Ĝ(x, y) = = `+1
Y`,m (θy , φy )Y`,m (θx , φx )
|x−y | `=0 m=−`
2` + 1 r>

luego el potencial G puede ser escrito como

38
 ∞ X
`  `

X 4π r< B` ∗
G = Ĝ + F = `+1
+ A` rx` + `+1 Y`,m (θy , φy )Y`,m (θx , φx )
`=0 m=−`
2` + 1 r> rx

Donde F puede ser escrita en la forma más general, y hemos incluido la linealidad del problema en los
coeficientes A y B. Notemos que las condiciones de salto en rx = ry esta incluida en la expresión de Ĝ.
Por lo tanto sólo nos tenemos que preocupar de las condiciones de borde en rx = a (note que a < ry ), y la
condición en rx → ∞, la cual es A = 0. La solución luego

∞ m=` `
a2`+1
 
X X 4π r< ∗
G(x, y) = `+1
− Y`,m (θy , φy )Y`,m (θx , φx ),
`=0 m=−`
2` + 1 r> ry`+1 rx`+1

que es el mismo resultado para b = 0.

Problema: Encuentre el potencial produ-

cido por un anillo de carga de densidad λ λ
(concéntrico con la esfera) fuera de una esfe-
ra conductora a potencial Vo cos θ sin φ.

Problema: Encuentre el potencial produci-

do por un anillo de carga de densidad λ de
radio R < a dentro de un cilindro de radio a L
y largo L. El anillo esta puesto perpendicular λ
al eje de simetrı́a del cilindro en L/2.

39
 Problema: Encuentre el potencial externo
producido por una esfera de radio a que tiene
2
a
un potencial Ψ = e−θ sobre su superficie.

Problema: Encuentre el potencial producido por un anillo de carga entre dos esferas, asumiendo G(r, θ) y
G(r, θ, φ).

40
6.4.4. Expansión en bases de un problema relacionado
Otro método para obtener las funciones de Green es expandiendo en un sistema completo de los valores
propios de un problema relacionado. El problema relacionado es

∇2 Ψ + (f (x) + λ)Ψ = 0

el cual es satisfecho para ciertos valores {λn , Ψn }

Z
2
∇ Ψn + (f (x) + λn )Ψn = 0 → Ψn Ψm d3 x = δnm

Otros valores dan las soluciones no convergentes o las singulares. La función de Green puede, entonces, ser
encontrada por

∇2 G(x, y) + (f (x) + λ)G(x, y) = −4πδ (3) (x − y)

P P Ψ∗n (y)Ψn (x)

G(x, y) = n an (y)Ψn (x) = 4π
n
λn − λ
Notemos que la ecuación para las funciones propias pueden en principio ser diferentes desde la ecuación
satisfecha por la función de Green.

Problema: Funciones propias en el espacio infinito.

Asumamos la ecuación de ondas

1
(∇2 + k 2 )Ψk (x) = 0 →Ψk (x) = eik·x
(2π)3/2
Z ik·(x−y)
1 1 e
λ=0 → = 2 2
d3 k
|x−y | 2π k

41
7. Magnetostática
Biot, Savart (∼ 1820) y Ampere (∼ 1820) propusieron que la fuerza inducida por un cable que lleva una
corriente I1 con un elemento de longitud d`1 en un cable que también lleva una corriente I2 con un elemento
de longitud d`2 es

I1 I2 d`2 × [d`1 × (r2 − r1 )] I1

F2,1 ∼ 2 3
= d`1 × B2 (r1 )
c |r1 − r2 | c

Jean baptiste biot (1774-1862)

André-Marie Ampère (1775-1836) Sacado de la Wikipedia,
Feliz Savart (1791-1841)
Sacado de la Wikipedia, The Dibner collection at the http://www.sil.si.edu/digitalcollections/hst/scientific- Sacado de la Wikipedia, Moonik:Buste de Félix Sa-

Smithsonian Institution (USA) identity/cf/display− results.cfm?alpha− sort=b

vart à l’Institut de France situé dans le 6e arrondis-

sement de Paris en France.

Esta dependencia fue descubierta por cuidadosos experimentos. En un sentido esto define la “extraña
cantidad” denominada corriente. En el sistema CGS la constante de proporcionalidad es k = 1. Las corrientes
son medidas en statamperes (10−9 /3 Amperes).

Las corrientes fueron una entidad misteriosa hasta que fueron descubiertas la conexión entre las cargas y
las corrientes. Si definimos

X
ρ(x, t) = qi δ (3) (x − xi (t))
i
X
J(x, t) = qi vi (t)δ (3) (x − xi (t))
i

claramente tenemos la ecuación de continuidad macroscópica

∂ρ
+∇·J=0
∂t
Para casos magnetostáticos, la dependencia del tiempo no está incluida.

42
Si la corriente no está concentrada en un cable conductor infinitesimal, definimos la densidad de corriente
como la corriente por el área representativa

J = ne < v >x → J d3 x = I d`
Veremos más adelante que esta función continua es un promedio espacial de la velocidad de la carga < v >x
en alguna escala particular intermedia, si está existe. Esto asume que hay tanto una escala espacial (y
temporal) donde tal promedio es significativo, i.e., donde las desviaciones son pequeñas en todo tiempo,
mirar Fig. 4 y análisis más adelante.

8. Campo magnético y potencial magnético

Un campo magnético puede ser definido en una forma similar como en el caso de la electrostática, pero se
requiere que

d`y × (x − y)
I I
I I
F= d`y × B(y) → B(x) =
c c | x − y |3
y pasando a una formulación macroscópica

J(y) × (x − y) 3
Z Z
1 3 1
F= J(x) × B(x) d x → B(x) = dy
c c | x − y |3
Esta definición implica la bien conocida “regla de la mano derecha”. Debemos siempre ser cuidadosos sobre
esta definición, ya que incluye una “fuerza propia”, similar a una “energia propia”, efecto que se produce
cuando x = y. Por lo tanto, solo tiene sentido la fuerza de un “grupo de cargas” sobre otro “grupo de
cargas”.

Hemos asumido la idea de superposición, lo cual está bien siempre y cuando no vayamos al régimen de
mecánica cuántica, óptica no lineal, etc.

Podemos fácilmente definir el potencial vectorial

 Z 
1 J(y) 3
B(x) = ∇ × dy
c |x−y |
con lo cual tenemos la ecuación de Maxwell

∇·B =0
Acá hemos usado
1 x−y 1
∇x =− 3
= −∇y
|x − y| |x − y| |x − y|
y también
 
J 1 1
∇x × = ijk Jk ∂j → −J(y) × ∇x .
|x − y| i |x − y| |x − y|

43
Como habı́amos visto en el caso electrostático

∇ · E = 4πρ
∇×E=0

fue útil tener el rotor y la divergencia del campo. Entonces, qué es ∇ × B? Tomemos
 
1 R J(y) 3
∇ × B(x) = ∇ × ∇ × dy
c |x−y |
       
1R 1 3 1R 2 1 3
= ∇∇ · J(y) dy − J(y)∇ dy
c |x−y | c |x−y |
 
4π 1 R ∇ · J(y) 3
= J(x) + ∇ dy
c c |x−y |
R ρ(y)d3 y
 
4π 1∂
= J(x) − ∇
c c ∂t |x−y |

y por lo tanto

4π 1 ∂E
∇×B= J(x) +
c c ∂t
Obtenemos el resultado integrando por partes los términos de la forma
Z Z Z
∇ · (f A)dV = f ∇ · JdV + J · ∇f dV.

Para el caso magnetostático el último término es cero, este último término fue propuesto por Maxwell
completando las llamadas ecuaciones de Maxwell. Por lo tanto lo vamos a mantener en las ecuaciones.

Podemos definir el potencial vectorial como

Z
1 J(y)
A(x) = d3 y + ∇Ψ → B=∇×A
c |x−y |
La elección de Ψ depende del gauge que elijamos. Es útil utilizar el gauge ∇ · A = 0, ya que a partir de

∇ × (∇ × A) = −∇2 A − ∇(∇ · A)
obtenemos
4π
∇2 A = − J
c
lo cual es lo mismo que la ecuación de Poisson, pero en forma vectorial. En el caso de la electrostática, en
general, no conocemos las corrientes inducidas en los bordes.

44
Supongamos que tenemos un campo vectorial que

∇ · Ā 6= 0
entonces podemos elegir otro campo vectorial

A = Ā + ∇Ψ
que da

∇ · A = ∇ · Ā + ∇2 Ψ
Por lo tanto es factible encontrar un Ψ tal que

∇2 Ψ = −∇ · Ā
y por lo tanto podemos asumir

∇·A=0

9. Método de solución
En el caso de la electrostática el borde contribuyó al potencial vı́a la generación de cargas superficiales.
Estamos interesados en los problemas de valor de borde donde el potencial fue definido en la superficie. En
el caso presente

∇·B = 0
4π
∇×B = J
c
la construcción de la función de Green es un poco mas complicada en relación con los conceptos equivalentes a
los conductores de la electrostática. Lo que requiere entender el comportamiento del medio bajo la presencia
del campo magnético, lo cual genera corrientes microscópicas, lo que depende del tiempo y es en general un
problema tiempo dependiente. ?cuál es el campo magnético dentro de un conductor?

Tomemos estas ecuaciones cerca de una interfase y tratemos de deducir como cambia el campo al cruzar la
interfase. Tomemos una caja infinitesimal de tapa ∆S, y lado dh → 0 (ver Fig. 3a). Usando el teorema de
Gauss para la integral de volumen obtenemos
H
Ω
∇ · B dx3 = 0
H
δΩ
B · n̂ ∆S = 0
(B2 − B1 ) · n̂ ∆S ≈ 0
con lo cual obtenemos una condición de borde para los campos

45
 E2 E2
n
dl
∆S dh
dh

E1 n E1
Figura 6: Condición de borde en la superficie con dh → 0.

(B2 − B1 ) · n̂ = 0

Para el caso de campo tangencial, tomemos un rectángulo perpendicular, como se muestra en la Fig. 3b.
Integrando sobre esta superficie que atraviesa la interfase tenemos

(∇ × B) · n̂ dS = 4π J~ · n̂ dS
H H
H H
B · d` = 04π J · n̂ S.
R
(E2 − E1 ) × n̂ ∆` ≈ 4π J⊥ dh ∆`
en el limite dh → 0, con
Z
K = lı́m J⊥ dh
dh→0

Por lo tanto, la carga de la superficie genera una discontinuidad en el campo magnético tangencial debido
a una corriente superficial. Con lo cual obtenemos la segunda condición de borde

4π
(B2 − B1 ) × n̂ = K
c
donde K es la corriente superficial que se induce por la discontinuidad en el campo magnético tangencial.

En algunas situaciones es posible construir una corriente de imagen para situaciones de condiciones de borde
simples. Estudiaremos tales casos en el próximo capı́tulo. Por ahora, nos concentraremos en la situación del
espacio infinito.

Problema: El campo magnético producido por un loop de corriente de radio R.

46
 I

Figura 7: Loop de corriente

Este problema puede ser solucionado por integración directa,

d` × (x − y)
I
I
B=
c |x − y|3
con

y = R[cos φ, sin φ, 0]
y
dy
d` = dφ = [−R sin φ, R cos φ, 0]dφ
dφ
Pero ahora estudiaremos usando la función de Green en el espacio infinito. La densidad de corriente en
coordenadas esféricas es

δ(y − a)
J=I δ(θy − π/2)[− sin φy , cos φy , 0]
a
con lo que podemos integrar usando Y`,m (θ, φ) = C`,m P`m (cos θ)eimφ ,

∞ Z `
4π X Y`,m (θ, φ) r< ∗
A(r, θ, φ) = dy 3 J(y) `+1 Y`,m (θy , φy )
c `,m 2` + 1 r>
∞ ` Z 2π  −iφy
`
− eiφy eiφy + e−iφy

4πIa X X r< Y`,m (θ, φ) m −imφy e
= `+1 2` + 1
C`,m P` (0) dφy e , ,0
c `=0 m=−` r> 0 2i 2
∞ ` `
8π 2 Ia X X r<
 
1 2 m m imφ δm,−1 − δm,1 δm,1 + δm,−1
= `+1 2` + 1
|C`,m | P` (cos θ)P` (0)e , ,0
c `=0 m=−` r> 2i 2

Dado que
∗
Y`,−m (θ, φ) = (−1)m Y`,m (θ, φ)
o lo que es equivalente

C`,m P`m (cos θ) = (−1)m C`,−m P`−m (cos θ)

Para el caso m = ±1, tenemos finalmente

47
 ∞ `
4πIa X 1 r<
A= `+1
|C`,1 |2 P`1 (0)P`1 (cos θ)φ̂
c `=0 2` + 1 r>
donde ahora r> = M ax[r, a] y r< = M in[r, a]. Dado que

(` + 1)!

(−1)(`−1)/2     ` impar
1

`
n−1 n+1
P` (0) =   2 ! !
2 2
0 ` par
De esta expresión podemos calcular el campo magnético lejos del loop r >> a como

" ∞
#
πIa2 (−1)n (2n + 1)! a2n 1

2 cos θ X
Br = +2 P
2n+2 2n+1
(cos θ)
c r3 n=1
2n n! r
2
" ∞
#
sin θ X (−1)n (2n + 1)! a2n 1

πIa
Bθ = − 3 + n (n + 1)!
P
2n+3 2n+1
(cos θ)
c r n=1
2 r

El término en los paréntesis redondos es el momento del dipolo magnético para el loop plano (área por la
corriente). Este resultado será útil cuando discutamos el medio magnético en el próximo capitulo.

48
10. Fuerza de Lorentz
Ahora podemos calcular la fuerza que siente un carga q en la presencia de un campo magnético y un campo
eléctrico

ρ(x) = qδ (3) (x − x(t)) J(x) = qv(t)δ (3) (x − x(t))

utilizando las definiciones de arriba, tenemos
Z Z
dv 3 1
m = ρ(y)E(y)d y + J(y) × B(y)d3 y
dt c
que nos da
dv v
m = qE(x) + q × B(x)
dt c
Esta es la fuerza de Lorentz, aunque algunos sugieren que fue propuesta por Maxwell antes.

James Clerk Maxwell (1831-1879)

Hendrik Antoon Lorentz (1853- Sacado de la Wikipedia, Frontpiece in James Maxwell,

1928) The Scientific Papers of James Clerk Maxwell. Ed: W.

Sacado de la Wikipedia, http://th.physik.uni- D. Niven. New York: Dover, 1890

frankfurt.de/ jr/physpictheo.html

49
11. Solución numérica
Estamos interesados en resolver numéricamente la ecuación de Laplace y la de Poisson

∇ · E = −∇2 Ψ = 4πρ
En general muchos de estos problemas quedan explı́citos como problemas con condiciones de borde, no
como un problema inicial que se puede integrar con métodos estándar.

11.1. Problemas en 1D
Para el caso de una dimensión tenemos que resolver

d2 Ψ
= −4πρ(x)
dx2
que en principio se puede encontrar como
Z x Z z
0
Ψ(x) = Ψ(0) + Ψ (0)x − 4π dz ρ(y)dy
0 0
Podemos notar que para integrar esta ecuación necesitamos ser capaces de hacer la integral y además tener
dos condiciones de borde en x = 0.

Muchas veces la integral no se puede resolver en forma analı́tica y en otras no tenemos las condiciones en
el mismo punto, sino mas bien tenemos condiciones de borde Ψ(0) y Ψ(L).

11.1.1. “Shooting Method”

Para esto debemos formular un “shooting method”.

Problema: Supongamos que tenemos que resolver

d2 Ψ
= −4πρ(x)
dx2
que en principio se puede encontrar como
Z x Z z
0
Ψ(x) = Ψ(0) + Ψ (0)x − 4π dz ρ(y)dy
0 0
Podemos notar que para integrar esta ecuación necesitamos ser capaces de hacer la integral y además tener
las condiciones de borde en x = 0.

Muchas veces la integral no se puede resolver en forma analı́tica y en otras no tenemos las condiciones en el
mismo punto, sino mas bien tenemos condiciones de borde Ψ(0) y Ψ(L). Para estos casos, podemos utilizar
un “método de disparo” para encontrar la solución en forma numérica.

50
Por ejemplo, si queremos resolver el problema de

ρ(x) = sin(cos[x]) Ψ(0) = 1 Ψ(1) = 4

Para eso definimos una función en MATHEMATICA

f un[xo ] := (p[x]/.N DSolve[{p00 [x] == −4P iρ[x], p[0] == 1, p0 [0] == xo}, p[x], {x, 0, 1}])/.x− > 1

con lo cual podemos encontrar el valor de Ψ(1) dado un valor inicial xo = Ψ0 (0). Ahora tenemos que
encontrar el cero de la función

g(xo) = f un[xo] − Ψ(1)

Para encontrar el valor de Ψ0 (0) que de Ψ(1) = 4. Esto se puede hacer con el método de la secante dados
dos valores iniciales s0 = 1, s1 = 3,
sn − sn−1
sn+1 = sn − g(sn )
g(sn ) − g(sn−1 )
el cual converge rápidamente a s → 7,2, y la solución se muestra en la Fig. 8.
5

0.2 0.4 0.6 0.8 1

Figura 8: Solución al problema en una dimensión

En MATHEMATICA podemos hacer uso de un método mas eficiente que esta dado por FindRoot.

11.1.2. Elementos finitos

También es posible construir directamente una aproximación numérica del funcional
Z
S[f ] = L dx

discretizando la función f (x) en termino de funciones localizadas

X
f (x) = fi hi (x)
i

donde hi (x) es el elemento finito alrededor del punto de la grilla (xi ). Por ejemplo,

51
 
|x − xi |
1− |x − xi | ≤ ∆x

hi (x) = ∆x
0 |x − xi | > ∆x
con xi = i∆x. Pero hay otras posibilidades como una cuadrática centrada en xi o una Gaussiana centrada
en xi . Podemos calcular la integral de arriba en función de fi . Luego de optimizar esta expresión sobre los
valores desconocidos de fi , obtenemos la solución esperada. En MATHEMATICA podemos utilizar NSolve
si el problema es polinomial, o FindMinimum si el problema es no-lineal en fi . Para discretizaciones mas
avanzadas, es posible que tengamos que usar un método mas global como algoritmos genéticos o CSA, ya
que pueden existir múltiples mı́nimos.

Ejemplo: Supongamos que queremos resolver el mismo problema de arriba

d2 Ψ
= −4πρ(x)
dx2
con

ρ(x) = sin(cos[x]) Ψ(0) = 1 Ψ(1) = 4

En este caso sabemos que este problema se puede escribir como un funcional
1
L = Ψ02 − 4πρΨ
2
Por lo tanto re-escribmos la acción como
N Z 1 N Z 1
1X 0 0
X
S[{Ψi }] = Ψi Ψj dxhi (x)hj (x) − 4π Ψi dxh0i (x)ρ(x)
2 i,j=0 0 i=0 0

Con la definición
Z 1
Ai,j = dxh0i (x)h0j (x)
0
y
Z 1
bi = 4π dxh0i (x)ρ(x)
0
donde i = 0, . . . , N y j = 0, . . . , N . Entonces se puede escribir como una matriz actuando sobre el vector
Ψ̄ = (Ψ0 , . . . , ΨN ),
1
S[Ψ̄] = Ψ̄T AΨ̄ − Ψ̄T b
2
Notemos que las condiciones de borde tienen que ser satisfechas por
X X
Ψ(0) = Ψi hi (0) Ψ(1) = Ψi hi (1)
i i

52
que se usan para despejar dos Ψi . Para nuestros elementos finitos triangulares tenemos

Ψ0 = 1 ΨN = 4
Con estos elementos finitos podemos calcular inmediatamente
 
1 −1 0 . . . 0
 −1 2 −1 . . . 0 
1  ... ... ... ... ...

A= 
∆x  
 0 0 −1 2 −1 
0 0 −1 1
por lo tanto
N N
1 X X
S= Ψn An,m Ψm − bn Ψn
2 n,m=0 n=0
Si ahora buscamos el óptimo
dS
=0 i = 1, .., N − 1
dΨi
obtenemos
N
X
Ai,n Ψn − bi = 0 i = 1, .., N − 1
n=0
porque la matriz A es simétrica. Notemos que las condiciones de borde están incluidas aquı́ ya que la
sumatoria va de n = 0 hasta m = N . Separemos los términos correspondientes a condiciones de borde
N
X −1
Ai,n Ψn = bi − Ai,0 Ψ0 − Ai,N ΨN i = 1, .., N − 1
n=1
Podemos re-escribir entonces esto como

ĀΨ = b̄
y la solución es entonces

Ψ̄ = Ā−1 b̄
Problemas mas complicados también se pueden enfrentar.

11.2. Problema en 2-D y 3-D

Para el caso de una dimensión tenemos que resolver

d2 Ψ d2 Ψ
+ 2 = −4πρ(x, y) Ψ|Ω
dx2 dy
en el borde.

53
 Y

x
0.2 0.4 0.6 0.8 1

Figura 9: Solución al problema en una dimensión con elementos finitos

11.2.1. Solución por relajación

Es posible probar que el teorema del valor medio para el potencial, en el cual el potencial en el centro de
una esfera, sin ninguna carga, es igual al potencial promedio en la superficie de la esfera. Si asumimos una
discretización esto puede ser usado por iteración para obtener el valor convergente para el potencial en una
región sin ninguna carga, y donde las condiciones lı́mite para el potencial es conocida

(n+1) 1 h (n) (n) (n) (n) (n) (n)

i 4π
Ψi,j,k = Ψi+1,j,k + Ψi−1,j,k + Ψi,j+1,k + Ψi,j−1,k + Ψi,j,k+1 + Ψi,j,k−1 + ∆x2 ∆y 2 ∆z 2 ρi,j,k
6 6
asumiendo una condición de borde conocida, la solución converge a la actual solución de la ecuación de
Poisson. En otros sistemas de coordenadas es necesario discretizar apropiadamente. La extensión a 1D o
2D es trivial.

Problema: Descargas fractales

Asumamos que definimos una descarga fractal como un sistema de puntos adyacentes que están en un
potencial fijo Ψ = 1, y el infinito (en este caso, un circulo de radio R) en Ψ = 0. Dada esta configuración
de potencial, podemos usar el método de relajación para calcular el potencial el cualquier lado. En dos
dimensiones el método de relajación es

(n+1) 1 h (n) (n) (n) (n)

i
Ψi,j = Ψi+1,j + Ψi−1,j + Ψi,j+1 + Ψi,j−1
4
Por lo tanto conocemos el campo eléctrico en cualquier lado, y especialmente en los puntos adyacentes a la
descarga. La descarga se desarrollará agregando en los puntos adyacentes a la descarga. Un punto tiene la
probabilidad de ser agregado a la descarga proporcional a un poder η en el campo eléctrico local
η
Pi,j ∼ Ei,j
Una vez que el punto es agregado, recalculamos el potencial y repetimos la operación, Usamos el código de
potencial escrito en C++ para calcular el fractal siguiente. El resultado está dado para η = 1,0
El potencial es mostrado como los contornos de color. Obviamente si η = 0 tenemos una estructura en dos
dimensiones, en cambio si η va hacia infinito tenemos una lı́nea de una dimensión.

54
 Figura 10: Descarga Fractal

11.2.2. Solución directa

También es posible construir una solución numérica resolviendo la ecuación de Laplace en forma directa
desde su discretización. En forma implı́cita podemos escribir en la aproximación a primer orden de las
derivadas
Ψi+1,j + Ψi−1,j − 2Ψi,j Ψi,j+1 + Ψi,j−1 − 2Ψi,j
+ = −4πρi,j
∆x2 ∆y 2
Tenemos que resolver este set implı́cito de relaciones, pero donde están mezclados condiciones de borde
con valores que desconocemos. La idea se separar los valores que conocemos de los valores desconocidos en
forma explicita

AΨ = B → Ψ = A−1 B
Ası́ se puede resolver el problema. De hecho la subrutina “Solve” de MATHEMATICA puede ser muy
útil en estas situaciones. Notemos que otras discretizaciones son también posibles. Lo interesante, es que
en principio una vez que resolvemos un problema, todos los otros problemas (con condiciones de borde
similares) están resueltos.

Problema: Lente electrostatico. Pretendemos resolver el problema

Ψ(x < 1, y = ±0,5) = Ψ(x > 4, y = ±0,5) = −1 Ψ(2 < x < 3, y = ±0,5) = 0
en [0, 5] × [−0,5, 0,5] con una interpolación lineal entre los electrodos en el borde. Utilizaremos la función
de MATHEMATICA NDSolve para encontrar una solución numérica.

55
 ny = 20; nx = 50; pl = -1; pr = -1;
stp = {pl, pl, 0, 0, pr, pr};
stx = {0, nx/5, 2*nx/5, 3*nx/5, 4*nx/5, nx + 1};
f= Interpolation[Transpose[stx, stp], InterpolationOrder → 1];

eqna = Table[-4*a[i,j]+a[i+1,j]+a[i-1,j]+a[i,j+1]+a[i,j-1]==0,{i,2,nx},{j,2,ny}];
eqn1 = Table[a[i, 1] == f[i], {i, 1, nx + 1}];
eqn2 = Table[a[i, ny + 1] == f[i], {i, 1, nx + 1}];
eqn3 = Table[a[1, j] == pl, {j, 2, ny}];
eqn4 = Table[a[nx + 1, j] == pr, {j, 2, ny}];

eq = Flatten[{eqna, eqn1, eqn2, eqn3, eqn4}];

var = Flatten[Table[a[i, j], {i, 1, nx + 1}, {j, 1, ny + 1}]];

sol = Flatten[NSolve[eq, var]];

Dada la solución, ahora podemos reconstruir el potencial como

∆x = 0,1; xmax = nx∆x; ymax = ny∆x/2;

Ψ = Interpolation[F latten[T able[{i∆x, j∆x − ymax, a[i + 1, j + 1]/.sol}, {i, 0, nx}, {j, 0, ny}], 1]];

g1 = ContourP lot[Ψ[x, y], {x, 0, xmax}, {y, −ymax, ymax}];

el cual se muestra en la Fig. 11. Ahora podemos calcular el campo eléctrico

tmax = 8;
ef x[x− , y− ] = D[Ψ[x, y], x];
ef y[x− , y− ] = D[Ψ[x, y], y];
ex[x− , y− ] = If [x >= 0&&x <= xmax&&y >= −ymax&&y <= ymax, ef x[x, y], 0];
ey[x− , y− ] = If [x >= 0&&x <= xmax&&y >= −ymax&&y <= ymax, ef y[x, y], 0];

y programar el calculo de la trayectoria de electrones

tra[xo− , yo− , eo− , θ− ]:= (

solu =NDSolve[{x”[t]==ex[x[t],y[t]],y”[t]==ey[x[t],y[t]],
x[0]==xo,y[0]==yo,x’[0]==eo*Cos[θ],y’[0]==eo*Sin[θ]}, {x[t], y[t]}, {t, 0, tmax}];
x1[t− ] = x[t] /. Flatten[solu];
y1[t− ] = y[t] /. Flatten[solu];
Return[ParametricPlot[{x1[t], y1[t]}, {t, 0, tmax},PlotStyle → {RGBColor[1, 0, 0]}]];
)

s = Table[tra[0.1, yv, 0.01, 0], {yv, -0.5, 0.5, 0.1}];

56
Finalmente, podemos grafical todo junto como se muestra en la Fig. 11, que muestra que el lente elec-
trostático es convergente.

gra = Show[g1, s, DisplayFunction → $DisplayFunction, PlotRange − > {{0, xmax}, {-ymax, ymax}},
Frame − > True, FrameLabel − > {’‘x”,”y”,””,””}];

0.75

0.5

0.25

0
y

-0.25

-0.5

-0.75

1 2 3 4 5
x

Figura 11: Solución al problema en dos dimensiones con las trayectorias de electrones superpuestos. La
trayectoria de los electrones en un lente electrostático definen la distancia focal del lente. Esto es un lente
convergente.

Las trayectorias en rojo se calcularon numéricamente utilizando

d2 r q
2
= − ∇Ψ
dt m
Notemos que podemos normalizar
Ψ x
Ψ̄ = y= τ = vo t
Vo L
con lo cual todo el problema queda determinado en termino de un parámetro adimensional. Es factible
entonces graficar la distancia focal y la magnificación del lente electrostático como función de este parámetro.

11.2.3. Solución por elementos finitos

En realidad el método de elementos finitos tiene que ver con un procedimiento mas general para resolver
problemas de este tipo. Notemos que la ecuación de Laplace viene de un principio de optimización
Z
(∇Ψ)2 − 4πρΨ dx3 δL = 0 → ∇2 Ψ = −4πρ


L[Ψ] = →

Por lo tanto uno puede asumir una forma de interpolación

X
Ψ(x, y) = Ψi,j h(x, y, xi , yj )
i,j

donde h(x, y, xi , yi ) es el elemento finito alrededor del punto de la grilla (xi , yj ). Usando una aproximación
lineal

57
 h(x, y, xi , yi ) = h(x, xi )h(y, yi )
donde

|x − xi |
1− |x − xi | ≤ ∆x

h(x, xi ) = ∆x
0 |x − xi | > ∆x
con xi = i∆x, se puede calcular la integral de arriba. Luego optimizamos esta expresión sobre los valores
desconocidos de Ψi,j . Aquı́ lo importante es ordenar apropiadamente los ı́ndices para poder escribir este
problema como

ĀΨ = b̄
y la solución es entonces

Ψ̄ = Ā−1 b̄
donde b contiene información sobre las condiciones de borde, y A es una matriz conocida.

Si tenemos condiciones de borde donde otro sistema de coordenadas es razonable, tenemos que evaluar el
funcional en forma apropiada.

Cuando las condiciones de borde no son regulares o rectangulares, es conveniente usar una grilla no rectan-
gular. Para estos casos es muy común utilizar una grilla triangular (la mayorı́a de las veces se utiliza una
grilla variable con en la Fig. 12

Figura 12: Grilla triangular variable

Para estos casos, conviene utilizar una interpolación lineal dentro de cada triangulo

Ψi (x, y) ≈ ai,0 + ai,1 x + ai,2 y

donde los coeficientes van variando en cada triangulo de area Ai y dependen de la posición de los tres
puntos que definen cada triangulo. Es posible hacer una transformación de una interpolación dentro de

58
los triángulos a una descripción en termino de elementos finitos sobre cada punto de la grilla, pero estos
claramente dependen de la posición de los otros puntos del triangulo. El caso de una grilla regular es mas
fácil. Notemos que el Lagrangiano en dos dimensiones es
Z X
(∇Ψ)2 dx2 = Ai (a2i,1 + a2i,2 )

donde la suma es sobre los triángulos. Notemos que puede ser mas fácil organizar computacionalmente los
triángulos en termino de un método de refinamiento de la grilla, en lo que se llama un método adaptivo.

59
Appendix
A. Convención de las sumatorias de Einstein
Desde ahora, usaremos la convención de las sumatorias de Einstein. Esto significa que dos ı́ndices en una
expresión implican una doble sumatoria explı́cita. Definamos

(A × B)k = εi,j,k Ai Bj
(
P (i, j, k) i 6= j 6= k
εi,j,k =
0 otro caso

con P (i, j, k) = (−1)n , con n el número de permutaciones de elementos adyacentes que necesitamos operar
sobre (i, j, k) para producir (1, 2, 3), el cual vale 1. Por ejemplo,

P (2, 3, 1) = P (2, 1, 3) = P (1, 2, 3) = 1

P (3, 2, 1) = −P (3, 1, 2) = P (1, 3, 2) = −P (1, 2, 3) = −1

Esto es equivalente a ver si la rotación cı́clica de (i, j, k) da (1, 2, 3) que tiene valer +1 o (2, 1, 3) que tiene
valor -1. Por ejemplo,

P (2, 3, 1) → P (1, 2, 3) → P (2, 3, 1) = 1

P (3, 2, 1) → P (2, 1, 3) → P (3, 2, 1) = −1

Con estas definiciones no es difı́cil darse cuenta que

εi,j,k εl,m,k = δi,l δj,m − δi,m δj,l

con δi,j = 0 si i = j, y cero en otro caso. Esta relación fundamental ya que permite calcular productos como

(A × B)2 = εi,j,k Ai Bj εl,m,k Al Bm

= (δi,l δj,m − δi,m δj,l )Ai Bj Al Bm
= A2 B 2 − (A · B)2
A × (A × B) = A2 B − A(A · B)

60
B. Funciones Delta
Una función delta puede ser construida de muchas maneras

Como un limite de una Gaussiana con σ → 0.

x2
1 −
δ(x) = lı́m √ e 2σ 2
σ→0 2πσ
Como un limite de δi,j en un espacio de dimensión N → ∞.

La definición que nos interesa es

δ(x − a) = 0 x 6= a
Z
δ(x − a)f (x) dx = f (a) a ∈ I
I

Es en realidad una distribución. Para definir distribuciones apropiadamente es necesario definir medidas
(por ejemplo Lebesgue) Notemos que la integral de Riemman no tiene mucho sentido acá (excepto como
limité de algo mas continuo).

Entonces es claro que

Z
0
f (a) = − δ 0 (x − a)f (x) dx
I
y haciendo un cambio de variable tenemos
X 1
δ(g(x)) = δ(x − xi )
i
| g 0 (x i) |
Notemos que aquı́ asumimos que las funciones son bastante suaves con ceros simples.

Otro resultado importante es

1
∇2 = −4πδ (3) (x − y)
|x−y |
Primero es fácil darse cuenta que en coordenadas esféricas tenemos

1 ∂2 ∂ 2Ψ
 
2 1 ∂ ∂Ψ 1
∇ Ψ= (rΨ) + sin θ +
r ∂r2 r2 sin θ ∂θ ∂θ r2 sin2 θ ∂φ2
con lo cual

61
  
12
∇ =0 r 6= 0
r
pero por el teorema de Gauss tenemos que
Z Z
2 3
∇ Ψd y= ∇Ψ · ndS
Ω ∂Ω
y usando
∂Ψ 1 ∂Ψ 1 ∂Ψ
∇Ψ = r̂ + θ̂ + φ̂
∂r r ∂θ r sin θ ∂φ
con lo cual
 
1 1
∇ = − 2 r̂
r r
y por lo tanto si integramos una esfera de radio R, obtenemos

  (
x∈
Z
1 0 /Ω
∇2 d3 y =
Ω |x−y | −4π x∈Ω

para cualquier R. Esto implica la expresión de arriba desde la definición de la función delta.

Es también posible probar esto usando un procedimiento de limites con un potencial de la forma
Z
ρ(y)
Ψa (x) = p dy
Ω |x − y|2 + a2
Esto tiene que ver con el concepto de distribución, la cual es un funcional lineal tal que
 
lı́m S (ψn ) = S lı́m ψn
n→∞ n→∞

para una secuencia convergente ψn . También se puede utilizar el concepto de medida (por ejemplo Lebesgue),
tal que
Z −∞
f (x)δ [dx] = f (0)
∞
o en forma mas “suave” como
Z −∞ Z −∞
f (x)δ [dx] = f (x)dH(x)
∞ ∞
donde

62
 
1 x≥0
H(x) =
0 x < 00
Ası́, podemos definir la derivada de la distribución
Z −∞ Z −∞ Z −∞
0 0
δ[Ψ] = − Ψ (x)H(x)dx = Ψ(x)H (x)dx = Ψ(x)dH(x) = Ψ(0)
∞ ∞ ∞

En d dimensiones, tenemos

δ (d) [cx] = |c|−n δ (d) (x)

y en el caso de la composición
Z Z
(d) d f (x)
f (x)δ [g(x)]d x = dS(x)
Ω ∂Ω |∇g|

63