Integral de energía potencial electrostática

Estoy tratando de calcular la energía total de un sistema simple de dos cargas a través de la integral de energía electrostática de un sistema dado en el libro de Griffiths:

tu = ϵ 0 2 V mi 2 d V .

Donde el volumen se integra en todo el espacio, por lo que el término límite que no se muestra aquí se reduce a cero. Creo que esto debería dar la misma respuesta que la fórmula estándar dada para cargas puntuales:

tu = 1 4 π ε 0 q 1 q 2 R .

Pero tengo problemas para evaluar la propia integral. Coloqué q 1 en el origen de los ejes de coordenadas y q 2 sobre el z -eje una distancia R lejos de la primera carga, y amplió el mi 2 término:

mi = mi 1 + mi 2
entonces
mi 2 = mi 1 2 + 2 mi 1 mi 2 + mi 2 2 .

Encontré que la integral de los términos propios diverge cuando se evalúa y, después de leer Griffiths, decidí descartar los términos de energía propia y solo retener la energía debido al término de intercambio.

Alquiler r = X 2 + y 2 + z 2 y r = X 2 + y 2 + ( z R ) 2 , encontré que la integral del término de interacción es:

mi 1 = 1 4 π ε 0 q 1 r 3 r y mi 2 1 4 π ε 0 q 2 r 3 r

tu = ϵ 0 V mi 1 mi 2   d V = q 1 q 2 dieciséis π 2 ε 0 V X 2 + y 2 + z 2 z R ( X 2 + y 2 + z 2 ) 3 2   ( X 2 + y 2 + ( z R ) 2 ) 3 2   d V .

Convirtiendo a coordenadas esféricas, con r = X 2 + y 2 + z 2 , θ el ángulo desde el eje z y φ el ángulo azimutal, donde he evaluado la integral azimutal:

tu = q 1 q 2 8 π ε 0 0 0 2 π r R porque ( θ ) ( r 2 2 R r porque ( θ ) + R 2 ) 3 2 pecado ( θ )   d θ   d r .

Me topé con una pared de ladrillos al tratar de evaluar la integral: normalmente usaría una sustitución en el caso de una integral simple, pero no estoy seguro de cómo hacerlo para una integral doble cuando las variables están todas mezcladas. ¿Estoy en el camino correcto?

No estoy seguro de que esta integral converja, dado que las otras dos divergen, ¿esta fórmula se aplica a cargas puntuales o solo a distribuciones de carga continuas?

Simplemente puede elegir un cuadro como origen (0,0,0) y tomar otras coordenadas como X , y , z o r , θ , ϕ . Entonces la integral se vuelve más simple. En caso de carga puntual, hice algunos argumentos en la respuesta a continuación.

Respuestas (2)

De los comentarios de la sección 2.4.4 de Griffith sobre energía electrostática, puede obtener su respuesta. Si considera las cargas puntuales, entonces, en realidad, esta integral está relacionada con la energía propia, que es infinita por lo general, para hacerla finita, a menudo introducimos un radio de corte d . (En física de partículas, a menudo usamos terminología simple y renormalizada, la renormalización es un proceso que hace infinito a finito) La integral relevante se describe bien en Griner's Electrodynamics y Jackson's ch1. Estos dos libros de texto contienen tanto el cálculo como su interpretación física también.

La fórmula de Poynting para la energía electrostática en volumen V

mi = V 1 2 ϵ 0 mi 2 d V

se puede derivar de la ley de Coulomb solo para los casos en que el campo que actúa sobre las partículas está definido en todas partes. Sin embargo, la partícula puntual tiene una densidad de carga infinita en el punto en el que está presente y el campo no está definido en ese punto. Entonces la derivación falla.

Para dos partículas puntuales en reposo, el trabajo necesario para llevar estas partículas a sus posiciones r 1 , r 2 se sabe que es

W = 1 4 π ϵ 0 q 1 q 2 | r 1 r 2 |
Si desea expresar esta energía en términos de campos EM únicamente, esto se puede escribir como

w h o yo mi   s pag a C mi ϵ 0 mi 1 ( X ) mi 2 ( X ) d 3 X
dónde mi 1 ( X ) = ϕ 1 ( X ) es el campo debido a la primera partícula y mi 2 ( X ) = ϕ 2 es el campo debido a la segunda partícula.

Prueba:

El potencial ϕ 1 es

ϕ 1 ( X ) = 1 4 π ϵ 0 q 1 | X r 1 |
y el potencial ϕ 2 ( X ) es
ϕ 2 ( X ) = 1 4 π ϵ 0 q 2 | X r 2 | .
La integral se convierte
w h o yo mi   s pag a C mi ϵ 0 mi 1 ( X ) mi 2 ( X ) d 3 X = w h o yo mi   s pag a C mi ϵ 0 ϕ 1 ( X ) ϕ 2 ( X ) d 3 X =
= w h o yo mi   s pag a C mi ϵ 0 ( ϕ 1 ϕ 2 ) d 3 X w h o yo mi   s pag a C mi ϵ 0 ϕ 1 Δ ϕ 2 d 3 X .
Para el campo electrostático, la primera integral es cero (esto se puede demostrar usando el teorema de Gauss). Para el segundo potencial, la ecuación de Poisson

Δ ϕ 2 = q 2 ϵ 0 d ( X r 2 )
se cumple por lo que llegamos a la integral

w h o yo mi   s pag a C mi 1 4 π ϵ 0 q 1 | X r 1 | q 2 ϵ 0 d ( X r 2 ) d 3 X
que tiene el valor

1 4 π ϵ 0 q 1 q 2 | r 2 r 1 | ,

que es lo mismo que W arriba.

Esta fórmula para energía EM tiene una versión general para campos dependientes del tiempo.

mi mi metro = ϵ 0 mi 1 mi 2 + 1 m 0 B 1 B 2 d 3 X

En caso de que haya más partículas involucradas, se pueden derivar fórmulas similares, con la suma de cada par de partículas.

RC Stabler, A Possible Modification of Classical Electrodynamics, Physics Letters, 8, 3, (1964), p. 185-187. http://dx.doi.org/10.1016/S0031-9163(64)91989-4

JA Wheeler, RP Feynman, Electrodinámica clásica en términos de interacción directa entre partículas, Rev. Mod. Phys., 21, 3, (1949), pág. 425-433. http://dx.doi.org/10.1103/RevModPhys.21.425

J. Frenkel, Zur Elektrodynamik punktförmiger Elektronen, Zeits. F. Phys., 32, (1925), pág. 518-534. http://dx.doi.org/10.1007/BF01331692

¿Este método es solo tu = ϵ o 2 mi neto 2 d 3 X ϵ o 2 mi 1 2 d 3 X ϵ o 2 mi 2 2 d 3 X , es decir, restando las singularidades?
No, esos términos son infinitos y no se pueden restar de forma matemáticamente válida. La fórmula que escribí anteriormente se puede derivar de una manera sencilla y matemáticamente válida del teorema del trabajo y la energía, que a su vez se puede derivar de las ecuaciones de Maxwell, la fórmula de la fuerza de Lorentz y la suposición de que las partículas actúan sobre otras partículas pero nunca sobre sí mismas.
¿Existe una fórmula para el producto escalar de 2 gradientes? ϕ 1 ϕ 2 = ( ϕ 1 ϕ 2 ) ϕ 1 Δ ϕ 2 ? ¿El lado izquierdo es un escalar mientras que el lado derecho es una matriz menos una función escalar? Probablemente me estoy perdiendo algo. Definitivamente veo cómo mi 1 mi 2 d V es igual a la conocida W calculando la integral. Para W con más de una partícula, puedo ver cómo la integral mi a mi b d V sigue siendo igual a W (nuevamente "calculándolo"). Gracias por esta linda prueba entre los 2
@ DWade64, sí, lo hay, pero tiene razón, la forma en que fue escrito no tenía sentido. La fórmula real es ϕ 1 ϕ 2 = ( ϕ 1 ϕ 2 ) ϕ 1 Δ ϕ 2 En palabras, en realidad hay una divergencia en lugar de un gradiente en el primer término. Gracias por reportar el error".
Ah, debería haber podido resolver eso, especialmente con el comentario sobre el Teorema de Gauss. Gracias por la actualización