CF 5F3 A FÍSICA NUCLEAR Y DEPARTÍCULAS

SEMESTRE 2022-1
CLASE 7 : DISPERSIÓN DE ELECTRONES MEDICIÓN
DE DISTRIBUCIÓN DE CARGA NUCLEAR. APLICACIONES
SEMESTRE 2022-1
Dispersión de Rutherford
Ahora consideraremos la sección eficaz de un electrón con energía 𝐸 que se dispersa en
un núcleo atómico con carga 𝑍𝑒. Para que el cálculo de la cinemática de reacción sea lo
suficientemente preciso, debe ser tanto relativista como mecánico cuántico. Para núcleos
pesados ​y electrones de baja energía, se puede despreciar el retroceso.
En una dispersión elástica la energía 𝐸 y el módulo de la cantidad de movimiento 𝑝 son
los mismos antes y después de la dispersión. El cálculo es similar al cálculo de la
trayectoria hiperbólica de un cometa que es desviado por el Sol al atravesar el sistema
solar.
Esto conduce a la fórmula de Rutherford: la dispersión de una partícula con carga z𝑒 y
energía cinética 𝐸𝑘𝑖𝑛 en un núcleo objetivo con carga 𝑍𝑒:
𝑑𝜎 𝑧𝑍𝑒 2 2
=
𝑑Ω 𝑅𝑢𝑡ℎ𝑒𝑟𝑓𝑜𝑟𝑑
4𝜋𝜀0 2 · 4𝐸𝑘𝑖𝑛 2 sen4 𝜃/2

En mecánica cuántica no relativista utilizando la regla de oro de Fermi se obtiene

exactamente el mismo resultado para esta sección transversal
Dispersión de una distribución de carga extendida
Considere el caso de un objeto tan pesado que el retroceso es insignificante. Entonces
podemos usar momentos tridimensionales. Si 𝑍𝑒 es pequeño, es decir, si
𝑍𝑒 ≪ 1
Se puede aplicar la aproximación de Born. Las funciones de onda 𝜓𝑖 y 𝜓𝑓 del electrón
entrante y saliente se pueden describir mediante ondas planas
1 1 ′
𝜓𝑖 = 𝑒 𝑖𝑝𝑥/ℏ 𝜓𝑓 = 𝑒 𝑖𝑝 𝑥/ℏ
𝑉 𝑉
El volumen 𝑉 debe ser grande en comparación con el centro de dispersión, y también
lo suficientemente grande como para que los estados de energía discretos en este
volumen puedan aproximarse a un continuo. Consideramos un haz de electrones con
una densidad de partículas 𝑛𝑎 (número de partículas por unidad de volumen). Se elige
un volumen de integración suficientemente grande, la condición de normalización viene
dada por
න 𝜓𝑖 2 𝑑𝑉 = 𝑛𝑎 · 𝑉 donde 𝑉 = 𝑁𝑎 /𝑛𝑎
𝑉
La velocidad de reacción 𝑊 (número de procesos físicos por unidad de tiempo) viene
dada por el producto de la sección transversal y la velocidad de la partícula del haz 𝑣𝑎
dividida por el volumen anterior.

𝜎𝑣𝑎 2𝜋 2 𝑑𝑛
=𝑊= ⟨𝜓𝑓 ℋ𝑖𝑛𝑡 𝜓𝑖 ⟩
𝑉 ℏ 𝑑𝐸𝑓

Aquí, 𝐸𝑓 es la energía total (energía cinética y masa en reposo) del estado final. Dado
que despreciamos el retroceso y que la masa en reposo es una constante, 𝑑𝐸𝑓 = 𝑑𝐸 ′
= 𝑑𝐸. La densidad 𝑛 de los posibles estados finales en el espacio de fase es

4𝜋 𝒑′ 2 𝑑 𝒑′ · 𝑉
𝑑𝑛 𝒑′ =
2𝜋ℏ 3

Por lo tanto, la sección transversal para la dispersión de un electrón en un elemento

angular sólido 𝑑Ω es
1 2𝜋 ′ 2
2 𝑉 𝒑 𝑑 𝒑′
𝑑𝜎 · 𝑣𝑎 · = ⟨𝜓𝑓 ℋ𝑖𝑛𝑡 𝜓𝑖 ⟩ 𝑑Ω
𝑉 ℏ 2𝜋ℏ 3 𝑑𝐸𝑓
 ′ 2 ′2
′
𝐸 𝑑𝜎 𝑉 𝐸 2
En buena aproximación 𝑣𝑎 ≈ 𝑐 ⇒ 𝒑 ≈ ⇒ = ⟨𝜓𝑓 ℋ𝑖𝑛𝑡 𝜓𝑖 ⟩
𝑐 𝑑Ω 2𝜋 2 ℏ𝑐 4

donde el hamiltoniano de interacción para una carga en un potencial eléctrico 𝜙 es

ℋ𝑖𝑛𝑡 = 𝑒𝜙
𝑒 𝑖𝒑′ ·𝒙 𝑖𝒑·𝒙
−
⇒ 𝜓𝑓 ℋ𝑖𝑛𝑡 𝜓𝑖 = න 𝑒 ℏ 𝜙 𝒙 𝑒 ℏ 𝑑 3 𝑥
𝑉

La transferencia de momento 𝒒 = 𝒑 − 𝒑′ . Entonces, el elemento matricial

𝑒 𝑖𝒒·𝒙
𝜓𝑓 ℋ𝑖𝑛𝑡 𝜓𝑖 = න 𝜙 𝒙 𝑒 ℏ 𝑑 3 𝑥
𝑉

El teorema de Green: para dos campos escalares 𝑣 y 𝑢 elegidos arbitrariamente, que

decrecen rápidamente a grandes distancias, la siguiente ecuación es válida para un
volumen de integración suficientemente grande:

න 𝑢𝛻 2 𝑣 − 𝑣𝛻 2 𝑢 𝑑 3 𝑥 = 0
Usando además,
𝑖𝒒·𝒙 ℏ2 2 𝑖𝒒·𝒙
𝑒 ℏ = − 2𝛻 𝑒 ℏ
𝒒

−𝑒ℏ2 2𝜙 𝒙
𝑖𝒒·𝒙
𝜓𝑓 ℋ𝑖𝑛𝑡 𝜓𝑖 = න 𝛻 𝑒 ℏ 𝑑3 𝑥
𝑉𝒒2

El potencial electrostático 𝜙 𝒙 para una distribución de carga 𝜚 𝒙 , en este caso

independiente del tiempo cumple la ecuación de Poisson

−𝜚 𝒙
𝛻2𝜙 𝒙 =
𝜀0

Se define la función 𝑓 tal que 𝜚 𝒙 = 𝑍𝑒𝑓 𝒙 , ∫ 𝑓 𝒙 𝑑 3 𝑥 = 1 condición de

normalización

𝑒ℏ2 𝑖𝒒·𝒙
3
𝑍 · 4𝜋𝛼ℏ 3𝑐 𝑖𝒒·𝒙
𝜓𝑓 ℋ𝑖𝑛𝑡 𝜓𝑖 = න 𝜚 𝒙 𝑒 ℏ 𝑑 𝑥 = න 𝑓 𝒙 𝑒 ℏ 𝑑3 𝑥
𝜀0 · 𝑉 𝒒 2 𝒒 2·𝑉
La integral
𝑖𝒒·𝒙
𝐹 𝒒 = න𝑓 𝒙 𝑒 ℏ 𝑑3 𝑥

es la transformada de Fourier de la función de carga 𝑓 𝒙 , normalizada a la carga total.

Se llama factor de forma de la distribución de carga. El factor de forma contiene toda la
información sobre la distribución espacial de la carga del objeto en estudio.

Analizaremos los factores de forma y su significado.

1) Dispersión de Rutherford, dispersión elástica: 𝐸 = 𝐸 ′ , 𝒑 = 𝒑′

𝜃
𝒒 = 2 𝒑 sen
2
 𝑑𝜎 4𝑍 2 𝛼 2 ℏ𝑐 2 𝐸 ′2 1
𝐸 = 𝒑 ·𝑐 = ∼
𝑑Ω 𝑅𝑢𝑡ℎ𝑒𝑟𝑓𝑜𝑟𝑑
𝒒𝑐 4 𝒒4

En experimentos reales de dispersión con núcleos o nucleones, se nota que las

secciones transversales, determinadas teóricamente, concuerdan con las secciones
transversales experimentales solo en el límite cuando 𝒒 → 0.

A valores más grandes de 𝒒 , las secciones transversales experimentales son

sistemáticamente más pequeñas. La razón de esto radica en la extensión espacial de
núcleos y nucleones.

La longitud de onda reducida del fotón virtual (intermediario de la interacción)

disminuye y la resolución aumenta. El efecto de la dispersión ya no se debe a la carga
total, sino solo parte de ella. Por tanto, la sección transversal disminuye.

La extensión espacial de un núcleo se describe mediante un factor de forma.

Restringiremos la discusión a los factores de forma de los sistemas esféricamente
simétricos que no tienen una orientación preferida en el espacio.
Simbolizamos este hecho escribiendo el factor de forma como

𝑑𝜎 𝑑𝜎
= · 𝐹 𝒒2 2
𝑑Ω 𝑒𝑥𝑝
𝑑Ω 𝑝𝑢𝑛𝑡𝑢𝑎𝑙

Ejemplos de factores de forma

Las primeras mediciones de factores de forma: principios de la década de 1950 en un
acelerador lineal en la Universidad de Stanford, California. Se midieron las secciones
transversales para una variedad de núcleos con electrones con energías de unos 500
MeV.
𝐹 𝒒2 = ∫ 𝑒 𝑖𝒒·𝒙/ℏ 𝑓 𝒙 𝑑 3 𝑥, න 𝑓 𝒙 𝑑3𝑥 = 1
𝑉

Si 𝑓 solo tiene dependencia radial, 𝑓 𝒙 = 𝑓 𝒙 = 𝑓 𝑟 = 𝜌 𝑟 /𝑍𝑒

∞
sen 𝒒 𝑟/ℏ 2
𝐹 𝒒2 = 4𝜋 න 𝑓 𝑟 𝑟 𝑑𝑟
0 𝒒 𝑟/ℏ
 ∞ 1 2𝜋 ∞
1 = න 𝑓 𝒙 𝑑 3 𝑥 = න න න 𝑓 𝑟 𝑟 2 𝑑𝜙𝑑 cos 𝜗 𝑑𝑟 = 4𝜋 න 𝑓 𝑟 𝑟 2 𝑑𝑟
𝑉 0 −1 0 0

La distribución radial de carga eléctrica: Transformada de Fourier del factor forma

medido experimentalmente.

1 𝑓(𝑟)
𝑓 𝑟 = න𝐹 𝒒 2 𝑒 −𝑖𝑞·𝑥ℏ 𝑑 3 𝑞
Ԧ
3 1,0
2𝜋 0,9
0,5
𝑓( 𝑟) se parametriza → 𝐹(𝒒2) 0,1
𝑟

La dispersión de un objeto con una superficie bien definida da como resultado un

patrón de difracción (sucesión de máximos y mínimos). Para una esfera homogénea
con radio 𝑅, por ejemplo, un mínimo se encuentra en

𝒒𝑅
≈ 4,5
ℏ
𝒒𝑅
≈ 4,5 𝑅 𝑝Ԧ
ℏ 𝑝Ԧ 𝑞Ԧ /2 𝜃
𝜃 ′
𝑝Ԧ − 𝑝Ԧ = 𝑞Ԧ = sen
12 𝒒 𝑝Ԧ′ 𝑝Ԧ 2
𝐶: ≈ 1,8 𝑓𝑚−1 : Mínimo 𝑝Ԧ′
ℏ

⇒ 𝑅 = 4,5ℏ/ 𝒒 ≈ 2,5 𝑓𝑚

Medición del factor de forma de 12C por dispersión de electrones.

Se muestra la sección transversal diferencial, que se midió a
una energía fija del haz de 420 𝑀𝑒𝑉 en 7 ángulos de dispersión
diferentes. La curva discontinua corresponde a la dispersión de
una onda plana en una esfera homogénea con un borde difuso
(Aproximación de Born); la curva sólida corresponde a un
análisis de fase en el método de descomposición en ondas
parciales en la que los parámetros se adaptaron a los datos de
la medición.
 1
𝐹 𝒒2 = ∫ 𝑒 𝑖𝒒·𝒙/ℏ 𝑓 𝒙 𝑑 3 𝑥, 𝑓 𝑟 = 3
න𝐹 𝒒 2 𝑒 −𝑖𝑞·𝑥ℏ 𝑑 3 𝑞
Ԧ
2𝜋

Ejercicio:
a) Demuestre que el factor de forma para una distribución de carga uniforme esta dado por
∞
sen 𝑞𝑟/ℏ 2
𝐹 𝑞Ԧ = 𝐹 𝑞Ԧ = 𝐹 𝑞 = 4𝜋 න 𝑓 𝑟 𝑟 𝑑𝑟
0 𝑞𝑟/ℏ

b) En primera aproximación, núcleos pesados de radio 𝑅, pueden describirse por una

distribución uniforme de carga 𝜌 𝑟Ԧ = 𝜌 𝑟Ԧ = 𝜌 𝑟 , normalizada en la unidad, es decir
∫ 𝜌 𝑟 𝑑3𝑟 = 1. Demuestre que el factor de forma es

3
𝐹 𝑞 = 3 sen 𝑥 − 𝑥 cos 𝑥 = 𝐹 𝑥 , donde 𝑥 = 𝑞𝑅/ℏ
𝑥

c) Encuentre, gráfica o numéricamente, los 3 primeros valores positivos de 𝑥 donde la

función 𝑓(𝑥) es cero, es decir 𝑓(𝑥𝑖 ) = 0, 𝑖 = 1,2,3.
d) En la figura se muestra la sección eficaz diferencial de difusión de electrones con
energía 𝐸 = 750 𝑀𝑒𝑉 por núcleos de 40𝐶𝑎 y 48𝐶𝑎 en función del ángulo de dispersión 𝜃.
¿Qué ángulos de difusión corresponden a las posiciones de los 3 ceros de la función,
evaluados en el item (c) y determine con estos valores el radio promedio (media
artimética) de ambos núcleos. Considere que 𝑝 ≫ 𝑚𝑐.

ƛ𝑑𝐵 =?
a) Integrando

2𝜋 𝜋 ∞
𝐹 𝑞Ԧ = න න න 𝑓 𝑟 𝑒 𝑖𝑞𝑟 cos 𝜃/ℏ 𝑟 2 sen 𝜃 𝑑𝑟𝑑𝜃𝑑𝜙 Haciendo 𝑦 = cos 𝜃
0 0 0

1 ∞ ∞ 1
𝑖𝑞𝑟𝑦/ℏ 2
ℏ 𝑖𝑞𝑟𝑦/ℏ
= −2𝜋 න න 𝑓 𝑟 𝑒 𝑟 𝑑𝑟𝑑𝑦 = 2𝜋 න 𝑓 𝑟 𝑒 𝑟 2 𝑑𝑟
−1 0 0 𝑖𝑞𝑟 −1

∞ ∞
ℏ 𝑖𝑞𝑟/ℏ −𝑖𝑞𝑟/ℏ 2
sen 𝑞𝑟/ℏ 2
= 2𝜋 න 𝑓 𝑟 𝑒 −𝑒 𝑟 𝑑𝑟 = 4𝜋 න 𝑓 𝑟 𝑟 𝑑𝑟
0 𝑖𝑞𝑟 0 𝑞𝑟/ℏ

𝑅
3 sen 𝑞𝑟/ℏ 2 3ℏ 𝑅 2 sen 𝑞𝑟/ℏ
b) 𝐹 𝑞 = 4𝜋 න 𝑟 𝑑𝑟 = න 𝑟 𝑑𝑟
4𝜋𝑅3 0 𝑞𝑟/ℏ 𝑞𝑅3 0 𝑞𝑟/ℏ

𝑅 𝑅 𝑅
3ℏ ℏ 𝑞𝑟 ℏ 𝑞𝑟 3ℏ2 𝑞𝑟 ℏ 𝑞𝑟
= −𝑟 cos − න − cos 𝑑𝑟 = 2 3 −𝑟 cos + sen
𝑞𝑅3 𝑞 ℏ 0 0 𝑞 ℏ 𝑞 𝑅 ℏ 𝑞 ℏ 0
 3ℏ2 𝑞𝑅 ℏ 𝑞𝑅 3ℏ2 𝑞𝑅 𝑞𝑅 𝑞𝑅
−𝑅 cos + sen = 3 3 sen − cos
𝑞 2 𝑅3 ℏ 𝑞 ℏ 𝑞 𝑅 ℏ ℏ ℏ

3 𝑞𝑅
= 3 sen 𝑥 − 𝑥 cos 𝑥 , 𝑥=
𝑥 ℏ

c) 𝐹 𝑥 = 0 ⇒ sen 𝑥 − 𝑥 cos 𝑥 = 0 ⇒ 𝑥 = tan 𝑥

Aplicando el método numérico de Newton-Raphson para aproximar raíces

𝑓 𝑥𝑛 sen 𝑥𝑛 − 𝑥𝑛 cos 𝑥𝑛 1
𝑥𝑛+1 = 𝑥𝑛 − ′ = 𝑥𝑛 − = 𝑥𝑛 − + cot 𝑥𝑛
𝑓 𝑥𝑛 cos 𝑥𝑛 − cos 𝑥𝑛 + 𝑥𝑛 sen 𝑥𝑛 𝑥𝑛

Otra forma es la solución gráfica

𝑥1 = 4,5; 𝑥2 = 7,7; 𝑥3 = 10,9

d) Como

𝜃 𝐸 𝜃 2𝐸 𝜃𝑖
𝑞 = 2𝑝 sen = 2 sen ⇒ 𝑥 = 𝑅𝑖 sen , donde 𝜃1 ≈ 17º, 𝜃2 ≈ 30º, 𝜃 ≈ 45º
2 𝑐 2 ℏ𝑐 2

40𝐶𝑎 ⇒ 3,73 𝑓𝑚
𝑅 promedio : 48𝐶𝑎 ⇒ 3,89 𝑓𝑚

- La sección transversal se midió en un amplio rango de 𝒒

Dentro de este rango, cambia en siete órdenes de magnitud.

- No se ven uno, sino tres mínimos en el patrón de

difracción. Este comportamiento de la sección transversal
significa que 𝐹(𝒒2) y la distribución de carga 𝜌(𝑟) se pueden
determinar con mucha precisión.

- Los mínimos de 48𝐶𝑎 se desplazan a valores de 𝒒

ligeramente más bajos que los de 40𝐶𝑎. Esto muestra que
48𝐶𝑎 es más grande.
 𝑑𝜎 2
𝑑𝜎
=𝐹 𝒒
𝑑Ω 𝑒𝑥𝑝
𝑑Ω 𝑅𝑢𝑡ℎ𝑒𝑟𝑓𝑜𝑟𝑑,𝑀𝑜𝑡𝑡

𝜌0
𝜌 𝑟 =
𝑟−𝑅1/2 /𝑎
1+𝑒
𝑟 = 𝑅1/2 − 2,2𝑎, 𝜌 = 0,9𝜌0

𝑟 = 𝑅1/2 + 2,2𝑎, 𝜌 = 0,1𝜌0

𝑑 = 4,4𝑎
𝑟2 1/2
≃ 𝑟0 · 𝐴1/3 , 𝑟0 = 0,94 𝑓𝑚
5 2
𝑅2 = 𝑟 , 𝑅 = 1,21 · 𝐴1/3 𝑓𝑚
3
𝑅 ≃ 𝑟0 · 𝐴1/3 , 𝑟0 ≃ 1,2 𝑓𝑚
Información acerca del radio nuclear se logra, analizando el factor de forma cuando
𝒒2 → 0
∞ 𝑛
𝒒 ·𝑅 1 𝑖 𝒒 𝒙 cos 𝜃
≪ 1, 𝐹 𝒒2 = න𝑓 𝒙 ෍ 𝑑 3 𝑥, 𝜗 = ∡ 𝒙, 𝒒
ℏ 𝑛! ℏ
𝑛=0

∞ 1 2𝜋 2
1 𝒒𝑟
𝐹 𝒒2 = න න න 𝑓 𝑟 1− cos 2 𝜗 + ⋯ 𝑑𝜙𝑑 cos 𝜗 𝑟 2 𝑑𝑟
0 −1 0 2 ℏ

∞ 2 ∞
1 𝒒
= 4𝜋 න 𝑓 𝑟 𝑟 2 𝑑𝑟 − 4𝜋 න 𝑓 𝑟 𝑟 4 𝑑𝑟 + ⋯
6 ℏ2
0 0 𝐹 𝒒2
𝑟2

1 𝒒2 2 𝑑𝐹 𝒒 2
𝐹 𝒒2 =1− 2
𝑟 + ⋯ ⇒ 𝑟 2 = −6ℏ2
2 อ
6ℏ 𝑑𝒒 𝑞
𝒒2 =0
En la exponencial
1 1
≈ 0,87 𝐺𝑒𝑉 → 𝑟 2 2 ≈ 0,8 𝑓𝑚
𝑎
∞
𝑟 2 = 4𝜋 න 𝑓 𝑟 𝑟 4 𝑑𝑟
0

∫ 𝜌 𝑟 𝑟 2 𝑑 3 𝑟Ԧ
𝑟2 =
∫ 𝜌 𝑟 𝑑 3 𝑟Ԧ

න 𝑓 𝑟 𝑑 3 𝑟Ԧ = 1, න𝜌 𝑟 𝑑 3 𝑟Ԧ = 𝑍𝑒
𝑉
 𝑅𝑁 ≈ 𝑟0 𝐴1/3, 𝑟0 ≈ 1,2 𝑓𝑚

3
𝑉 ∼ 𝑅𝑁 ∼ 𝐴

𝑀 𝐴
𝜌𝑀 = ≈ ≈ 1
𝑉 𝐴
- La densidad de carga 𝜌(0) en el centro del núcleo disminuye ligeramente al crecer el
número de masa. Si consideramos la presencia de neutrones al multiplicar por
𝐴/𝑍 se encuentra una densidad nuclear casi idéntica en el interior para casi todos los
núcleos.
𝜌𝑛 ≈ 0,17 𝑛𝑢𝑐𝑙𝑒𝑜𝑛𝑒𝑠/𝑓𝑚3

- Algunos núcleos tienen una ligera desviación de la forma esférica y poseen

deformaciones elipsoidales: los lantánidos (los elementos de las "tierras raras"). Su
forma exacta no puede determinarse mediante la dispersión elástica de electrones

- Los núcleos ligeros como 6,7𝐿𝑖, 8𝐵, 4𝐻𝑒, son casos especiales. En estos núcleos, no
se forma una meseta de densidad constante en el interior del núcleo y la densidad de
carga es aproximadamente gaussiana.
Teoría de Scattering, Aproximación de Born
1. Amplitud de Scattering

Procederemos a demostrar que podemos obtener la sección transversal diferencial en

el sistema referencial del centro de masa a partir de una forma asintótica de la solución
de la ecuación de Schrödinger

ℏ2 2
− 𝛻 𝜓 𝑟Ԧ + 𝑉 𝑟Ԧ 𝜓 𝑟Ԧ = 𝐸𝜓 𝑟Ԧ (1.1)
2𝜇

Observemos primero la determinación de la amplitud de scattering 𝑓(𝜃, 𝜙), que puede

ser obtenida de las soluciones de (1.1), que a su vez se puede reescribir como

2𝜇
𝛻2 + 𝑘2 𝜓 𝑟Ԧ = 2 𝑉 𝑟Ԧ 𝜓 𝑟Ԧ (1.2)
ℏ
donde 𝑘 2 = 2𝜇𝐸/ℏ2. La solución general de la ecuación (1.1) consiste de una suma de
dos componentes:

1.Una solución general de la ecuación homogénea:

2𝜇𝐸
2 2
𝛻 + 𝑘 𝜓 𝑟Ԧ = 0, con 𝑘02 = 2 (1.3)
ℏ
que resulta en la solución

𝜓ℎ 𝑟Ԧ → 𝜙𝑖 𝑟Ԧ = 𝐴𝑒 𝑖𝑘0 ·𝑟Ԧ (1.4)

que es la onda plana incidente; y

2. Una solución particular de (1.2) con el potencial de interacción.

2. Solución general de la ecuación de Schrödinger en términos de la función

de Green
La solución general de (1.2) puede expresarse en términos de la función de Green
 2𝜇
𝜓 𝑟Ԧ = 𝜙𝑖 𝑟Ԧ + Ԧ 𝑉 𝑟Ԧ ′ 𝜙 𝑟′
න𝐺 𝑟Ԧ − 𝑟′ Ԧ 𝑑3 𝑟′
Ԧ (2.1)
ℏ

Donde 𝜙𝑖 𝑟Ԧ = 𝑒 𝑖𝑘0·𝑟Ԧ y 𝐺 𝑟Ԧ − 𝑟′
Ԧ es la función de Green correspondiente al operador del
lado izquierdo de la ecuación (1.3). La función de Green se obtiene resolviendo la
ecuación para una fuente puntual

𝛻 2 + 𝑘 2 𝐺 𝑟Ԧ − 𝑟′
Ԧ = 𝛿 𝑟Ԧ − 𝑟′
Ԧ 2.2

Donde 𝐺 𝑟Ԧ − 𝑟′
Ԧ y 𝛿 𝑟Ԧ − 𝑟′
Ԧ están dadas por sus transformadas de Fourier como sigue

1 Ԧ 𝑟Ԧ′
𝐺 𝑟Ԧ − 𝑟′
Ԧ = 3
න𝑒 𝑖𝑞· 𝑟−
𝐺෨ 𝑞Ԧ 𝑑3 𝑞Ԧ 2.3
2𝜋

1 𝑖𝑞· Ԧ 𝑟Ԧ′
𝑟−
𝛿 𝑟Ԧ − 𝑟′
Ԧ = 3 න𝑒 𝑑3 𝑞Ԧ 2.4
2𝜋
La sustitución de (2.3) y (2.4) en (2.2) conduce a

1
𝑘 − 𝑞 𝐺෨ 𝑞Ԧ = 1,
2 2
𝐺෨ 𝑞Ԧ = 2 2.5
𝑘 − 𝑞2

La expresión de 𝐺 𝑟Ԧ − 𝑟′
Ԧ puede obtenerse introduciendo (2.5) en (2.3)
′
Ԧ 𝑟Ԧ ∞ 𝜋 2𝜋
1 𝑒 𝑖𝑞· 𝑟− 3
1 Ԧ 𝑟Ԧ ′
𝑞 2 𝑑𝑞
𝐺 𝑟Ԧ − 𝑟′
Ԧ = 3 න 𝑘 2 − 𝑞2 𝑑 𝑞
Ԧ= 3න
න 𝑒 𝑖𝑞 𝑟− cos 𝜃 sen 𝜃 𝑑𝜃 න 𝑑𝜙
2 2
(2.6)
2𝜋 2𝜋 0 0 0 𝑘 −𝑞

Para integrar sobre 𝜃 en (2.6) se necesita hacer el cambio de variable 𝑥 = cos 𝜃

𝜋 1
Ԧ 𝑟Ԧ′ Ԧ 𝑟Ԧ′ 𝑥
1 ′
Ԧ 𝑟Ԧ Ԧ 𝑟Ԧ ′
න 𝑒 𝑖𝑞 𝑟− cos 𝜃 sen 𝜃 𝑑𝜃 = න 𝑒 𝑖𝑞 𝑟− 𝑑𝑥 = 𝑒 𝑖𝑞 𝑟− − 𝑒 −𝑖𝑞 𝑟− (2.7)
0 −1 𝑖𝑞 𝑟Ԧ − 𝑟Ԧ ′

entonces (2.6) se expresa como

∞
1 𝑞 Ԧ 𝑟Ԧ′
𝑖𝑞 𝑟− Ԧ 𝑟Ԧ′
−𝑖𝑞 𝑟−
𝐺 𝑟Ԧ − 𝑟′
Ԧ = 2 ′
න 2 2
𝑒 − 𝑒 𝑑𝑞 2.8
4𝜋 𝑖 𝑟Ԧ − 𝑟Ԧ 0 𝑘 − 𝑞