La introducción del cuadripotencial nos permite escribir los campos eléctrico y magnético en términos tensoriales, lo que será muy útil para derivar la forma de las ecuaciones de Maxwell manifiestamente covariantes.
Por definición de los potenciales:
\begin{eqnarray} \nonumber \vec B & = & \vec{\nabla}\times\vec A\\ \nonumber \vec E & = & -\vec{\nabla}\phi-\frac{1}{c}\frac{\partial\vec A}{\partial t} \end{eqnarray}Si lo escribimos en términos de las componentes explícitas del cuadripotencial obtenemos:
$$ \begin{array}{rcl} B^i & = & \sum_{j,k} \epsilon_{ijk} \partial_j A^k\\ E^i & = & -\partial_i A^0-\partial_0 A^i \end{array} $$donde $\epsilon_{ijk}$ es el símbolo de Levi-Civita:
$$ \epsilon_{ijk}= \left\{ \begin{array}{ll} 0 & \mathrm{Si\;hay\;un\;índice\;repetido}\\ +1 & \mathrm{Si}\;ijk\;\mathrm{es\; permutación\;par}\\ -1 & \mathrm{Si}\;ijk\;\mathrm{es\; permutación\;impar} \end{array} \right. $$Una interesante simetría se produce si escribimos las componentes explícitas:
$$ \begin{array}{rcl} B^1 & = & \partial_2 A^3-\partial_3 A^2\\ E^1 & = & -\partial_1 A^0-\partial_0 A^1 \end{array} $$Si usamos las componentes contravariantes de la derivada (que implican simplemente un cambio en el signo de las componentes espaciales de las derivadas) podemos escribir:
$$ \begin{array}{rcl} B^1 & = & -\partial^2 A^3+\partial^2 A^2\\ E^1 & = & \partial^1 A^0-\partial^0 A^1 \end{array} $$que se puede escribir como:
$$ \begin{array}{rcl} B^1 & = & -(\partial^2 A^3-\partial^2 A^2)\\ E^1 & = & -(\partial^0 A^1-\partial^1 A^0) \end{array} $$Estas relaciones dejan adivinar que la combinación de cantidades $\partial^mu A^\nu$ parece jugar un papel unificador en la nueva formulación de la teoría electromagnética, de la misma manera que los campos eléctrico y magnético lo jugaron en el tiempo de Maxwell. Esta nueva cantidad es tensorial y tiene segundo rango (es la combinación simultánea del tensor derivada y el tensor potencial).
De aquí, si definimos el tensor contravariante de segundo rango:
$$ F^{\mu\nu}\equiv\partial^\mu A^\nu - \partial^\nu A^\mu $$Vemos que el campo eléctrico y el campo magnético se pueden escribir en términos de este tensor como:
$$ \begin{array}{rcl} B^i & = & - \epsilon_{ijk} F^{jk}\\ E^i & = & - F^{0i} \end{array} $$donde es importante anotar que en la primera expresión no se usa la convención de suma de Einstein.
Otra manera de verlo podría ser escribir explícitamente las componentes del tensor en coordenadas cartesianas:
$$ F^{\mu\nu}: \left( \begin{array}{cccc} 0 & -E_x & -E_y & -E_z\\ E_x & 0 & -B_z & +B_y \\ E_y & +B_z & 0 & -B_x \\ E_z & -B_y & +B_x & 0 \\ \end{array} \right) $$Llamamos a $F^{\mu\nu}$ el tensor electromagnético o tensor de Faraday.
Proposición: Propiedades del tensor de Faraday. El tensor de Faraday tiene las siguientes propiedades:
- Antisimetría: $F^{\mu\nu}=-F^{\nu\mu}$
- Traza nula: $\sum_\mu F^{\mu\mu}=0$
- Propiedad cíclica de la derivada contravariante: $$\partial^{[\alpha}F^{\mu\nu]}=F^{[\mu\nu,\alpha]}\equiv\partial^\alpha F^{\mu\nu}+\partial^\mu F^{\nu\alpha}+\partial^\nu F^{\alpha\mu}=0$$ donde en lo sucesivo $[\alpha\beta\gamma]=\alpha\beta\gamma,\beta\gamma\alpha,\gamma\alpha\beta$ es la permutación cíclica de los índices.
- Determinante invariante: $\det F^{\mu\nu}=(\vec E\cdot \vec B)^2$
- Escalar: $F^{\mu\nu}F_{\mu\nu}=2(B^2-E^2)$.
Las ecuaciones de Maxwell inhomogenéas las habíamos escrito como:
$$ \partial_\mu \partial^\mu A^\nu = \frac{4\pi}{c} J^\nu $$Podemos sumar y restar $\partial_\mu\partial^\nu A^\mu$ para obtener:
$$ \partial_\mu(\partial^\mu A^\nu-\partial^\nu A^\mu)+ \partial_\mu \partial^\nu A^\mu= \frac{4\pi}{c} J^\nu $$Como para escribir la ecuación anterior asumimos el Gauge de Lorentz, $\partial_\mu^\mu A^\nu=0$ las ecuaciones de Maxwell inhomogeneas se pueden escribir en términos del tensor de Faraday como:
$$ {F^{\mu\nu}}_{,\mu}=\frac{4\pi}{c}J^\nu $$Qué es la forma manifiestamente covariante de las ecuaciones inhomogéneas de Maxwell.
¿Y las ecuaciones de Maxwell homogéneas?. Estas ecuaciones son:
\begin{eqnarray} \nonumber \vec{\nabla}\cdot\vec{B}&=&0\\ \nonumber \vec{\nabla}\times\vec{E}+\frac{1}{c}\frac{\partial\vec B}{\partial t}&=&0\\ \end{eqnarray}Vemos que se trata de 4 ecuaciones que se pueden escribir también en la forma:
$$ \begin{array}{rcl} \partial_1 B_x + \partial_2 B_y + \partial_3 B_z & = & 0 \\ \partial_0 B_x + \partial_2 E_z - \partial_3 E_y & = & 0 \\ \partial_0 B_y + \partial_3 E_x - \partial_1 E_z & = & 0 \\ \partial_0 B_z + \partial_1 E_y - \partial_2 E_x & = & 0 \\ \end{array} $$Reemplazando los valores del campo por las respectivas entradas del vector de Faraday estas ecuaciones se pueden escribir como:
$$ \begin{array}{rcl} \partial_1 F^{32} + \partial_2 F^{13} + \partial_3 F^{21} & = & 0 \\ \partial_0 F^{32} + \partial_2 F^{30} - \partial_3 F^{20} & = & 0 \\ \partial_0 F^{13} + \partial_3 F^{10} - \partial_1 F^{30} & = & 0 \\ \partial_0 F^{21} + \partial_1 F^{20} - \partial_2 F^{10} & = & 0 \\ \end{array} $$Si se cambian las derivadas covariantes por derivadas contravariantes y se reorganizan los índices, se obtienen expresiones que recuerdan una propiedad introducida antes:
$$ \begin{array}{rcl} \partial^1 F^{23} + \partial^2 F^{31} + \partial^3 F^{12} & = & 0 \\ \partial^0 F^{32} + \partial^3 F^{20} + \partial^2 F^{03} & = & 0 \\ \partial^0 F^{13} + \partial^1 F^{30} + \partial^3 F^{10} & = & 0 \\ \partial^0 F^{21} + \partial^2 F^{10} + \partial^1 F^{02} & = & 0 \\ \end{array} $$que no son otra cosa que la propiedad cíclica de la derivada contravariante del tensor de Faraday:
$$ \partial^{[\alpha}F^{\mu\nu]}=0 $$En síntesis, en términos explícitamente covariantes, las ecuaciones de Maxwell se pueden escribir como:
\begin{eqnarray} {F^{\mu\nu}}_{,\mu} & = & \frac{4\pi}{c}J^\nu\;\;\mathrm{(Gauss-E)\;y\;(Ampère-Maxwell)}\\ \partial^{[\alpha}F^{\mu\nu]} & = & 0\;\;\mathrm{(Gauss-B)\;y\;(Maxwell-Faraday)}\\ \end{eqnarray}Algunas de las propiedades mencionadas antes se pueden expresar y desarrollar en sistemas de álgebra computacional como SymPy. Para hacerlo primero debermos definir las componentes de los campos:
import sympy as sp
Ex,Ey,Ez,Bx,By,Bz=sp.symbols("E_x,E_y,E_z,B_x,B_y,B_z")
La representación matricial del tensor electromagnético será:
F=sp.Matrix([[0,-Ex,-Ey,-Ez],[Ex,0,-Bz,By],[Ey,Bz,0,-Bx],[Ez,-By,+Bx,0]])
F
Con esta representación podemos obtener el tensor de Faraday completamente covariante usando:
$$ F_{\mu\nu}=\eta_{\mu\alpha}\eta_{\nu\beta}F^{\alpha\beta} $$¿Cómo puede calcularse esta operación matricialmente? El producto de dos matrices $A$ y $B$ se puede escribir como:
$$ (AB)_{ij}=\sum_n A_{in}B_{nj} $$Usando la representación matricial de los tensores como:
$$ F_c = \eta F \eta^\mathrm{T} $$Proposición: Transformación de tensores de segundo rango. Si $A^{\mu\nu}$ es un tensor de segundo rango y $A$ es su representación matricial (una matriz $4\times 4$ con entradas iguales a las componentes del tensor en una base dada) en el sistema de referencia del laboratorio, la representación matricial del tensor en el sistema de referencia vijero es:
$$A'=\Lambda A \Lambda^\mathrm{T}$$
Para calcular la representación matricial del tensor de Faraday covariante, debemos introducir el tensor métrico:
eta=sp.Matrix([[1,0,0,0],[0,-1,0,0],[0,0,-1,0],[0,0,0,-1]])
Y hacer el producto matricial:
Fcov=eta*F*sp.transpose(eta)
Fcov
En términos de las componentes matriciales del tensor el producto del tensor de Faraday contravariante y covariante se puede escribir como:
$$ F^{\mu\nu}F_{\mu\nu}=-\sum_\mu \left(\sum_\nu F^{\mu\nu}F_{\nu\mu}\right) $$El producto entre paréntesis se puede considerar las diagonales de la matriz:
$$ A=FF_c $$que a su vez es la representación matricial del tensor ${A^\mu}_\nu=F^{\mu\alpha}F_{\alpha\nu}$. Así el escalar buscado es la traza, es decir la suma de las componentes en la diagonal, de la matriz A:
$$ F^{\mu\nu}F_{\mu\nu}=-\sum_\mu{A^\mu}_\mu=-{A^\mu}_\mu\equiv -\mathrm{tr}({A^\mu}_\nu) $$O en términos matriciales:
$$ FF_c=-\mathrm{tr}(A) $$Proposición: Módulo de un tensor de segundo rango. Si $A^{\mu\nu}$ es un tensor de segundo rango y $A$ es su representación matricial (una matriz $4\times 4$ con entradas iguales a las componentes del tensor en una base dada), $A_{\mu\nu}$ es el correspondiente tensor covariante con $A_c$ su representación matricial, el módulo del tensor se define como:
$$|A|\equiv A^{\mu\nu}A_{\mu\nu}$$Si definimos el tensor ${\bar{A}^\mu}_\nu=A^{\mu\alpha}A_{\alpha\nu}$, con representación matricial $\bar{A}=AA_c$, la magnitud de un tensor simétrico sera: $$|A|=\mathrm{tr}(\bar{A})$$ y si es antisimétrico: $$|A|=-\mathrm{tr}(\bar{A})$$
En SymPy:
A=-(F*Fcov).trace()
A
Que es justamente lo que habíamos calculado previamente.
¿Cómo calcular ahora el determinante del tensor de Faraday?
sp.factor(F.det())
Que es justamente la propiedad que habíamos visto antes.
Llegados a este punto podemos encontrar las leyes de transformación de los campos eléctricos y magnético.
Siendo el tensor de Faraday, un tensor, al pasar de un sistema de referencia inercial a otro a través de las transformaciones de Lorentz-Einstein $\Lambda^{\mu'}_\nu$ el tensor transforma como:
$$ F^{\mu'\nu'}={\Lambda^{\mu'}}_\mu {\Lambda^{\nu'}}_\nu F^{\mu\nu} $$Si definimos la matriz $\bar{F}$ como:
$$ \bar{F}^{\mu'\nu}=\sum_\mu {\Lambda^{\mu'}}_\mu F^{\mu\nu} $$o bien, matricialmente:
$$ \bar{F} = \Lambda F $$Entonces la transformación original es:
$$ F^{\mu'\nu'}=\sum_\nu \bar{F}^{\mu'\nu} {\Lambda^{\nu'}}_\nu $$Por la simetría de la transformación de Lorentz-Einstein:
$$ F^{\mu'\nu'}=\sum_\nu \bar{F}^{\mu'\nu} {\Lambda_{\nu}}^{\nu'} $$que en términos matriciales se puede escribir como:
$$ F'=\Lambda F \Lambda^\mathrm{T} $$Podemos intentar realizar estas operaciones usando SymPy. Para ello definamos los símbolos que usaremos:
Ahora veamos como es la transformación. Para ello introduzcamos primero la TLE en la configuración estándar:
u=sp.symbols("u_L")
g=1/sp.sqrt(1-u**2)
D=sp.Matrix([[g,-g*u,0,0],[-g*u,g,0,0],[0,0,1,0],[0,0,0,1]])
D
La transformacion nos quedará:
sp.simplify(D*F*sp.transpose(D))
De donde resulta entonces que en la transformación estándar las reglas de transformación de los campos eléctricos se pueden escribir como:
$$ \begin{array}{rcl} E'_x & = & E_x\\ E'_y & = & \gamma(E_y - u_L B_z)\\ E'_z & = & \gamma(E_z - u_L B_y)\\ \end{array} $$Y las del campo magnético:
$$ \begin{array}{rcl} B'_x & = & B_x\\ B'_y & = & \gamma(B_y + u_L E_z)\\ B'_z & = & \gamma(B_z - u_L E_y)\\ \end{array} $$Las componentes $E_y$ y $B_z$ de esta transformación coinciden con las que habíamos obtenido en la deducción de las transformaciones de Lorentz-Einstein, con lo que la teoría es completamente consistente.
Es posible mostrar que para una configuración general, la transformación de los campos sería:
$$ \begin{array}{rcl} \vec{E} & = & \vec{E}' + (\gamma-1)\vec{E}_\perp - \gamma \vec{u}_L\times\vec B\\ \vec{B} & = & \vec{B}' + (\gamma-1)\vec{B}_\perp + \gamma \vec{u}_L\times\vec E\\ \end{array} $$donde $\vec{E}_\perp$ y $\vec{B}_\perp$ son los vectores de campo eléctrico y magnético perpendiculares a la velocidad relativa del sistema de referencia.
Si bien en la teoría electromagnética clásica, la fuerza de Lorentz se postula, en realidad esta resultad de las transformaciones del tensor de Faraday que vimos en la sección anterior.
Para ello considere por ejemplo el caso de una carga que se mueve en un campo magnético estático $\vec B$. En el sistema de referencia de la carga ($\vec{u}_L=\vec{v}_L$):
$$ \vec{E}' = \gamma_v \vec{v}_L\times\vec B\\ $$y si postulamos que la fuerza es simplemente $\vec F'=q\vec E'$ la fuerza será:
$$ \vec{F}' = \gamma_v q \vec{v}_L\times\vec B\\ $$Ahora bien, sabemos de las reglas de transformación de la fuerza que la componente perpendicular transforma como:
$$ \vec F_\perp = \frac{\gamma_u^{-1}}{1+\vec{u}_l\cdot\vec{v}'_L} F'_\perp $$Dado que en el sistema de referencia de la carga $\vec v'_L=\vec o$ y que $\vec u_L=\vec v_L$ la fuerza sobre la carga en el sistema de referencia del observador será:
$$ \vec F=q \vec{v}_L\times\vec B\\ $$Si además hay un campo eléctrico en el sistema del observador (donde la partícula se mueve) la fuerza será entonces:
$$ \vec F = q \vec E + \frac{1}{c}q\vec{v}\times \vec B $$Nota 4.19. La fuerza de Lorentz y la rotación del campo de Faraday. Una manera e interpretar el resultado anterior y el origen en general de la fuerza de Lorentz podría ser la de pensar de la siguiente manera. El campo de Faraday cuando solo hay campo magnético tiene solo componentes espaciales. Cuando una carga eléctrica está en movimiento es como si el campo de Faraday se rotará. Al rotarse aparecen componentes temporales del campo de Faraday que la carga percibe como campo eléctrico y de allí que experimente una fuerza.
Escribiendo explícitamente la componente $x$ queda:
$$ \frac{\mathrm{d}p_x}{\mathrm{d}t}=\frac{q}{c} (cE_x + v_y B_z - v_z B_y) $$Multiplicando por $\gamma$ a ambos lados podemos escribir:
$$ \frac{\mathrm{d}P^1}{\mathrm{d}\tau}=\frac{q}{c} (U^0 F^{10} - U^2 F^{12} - U^3 F^{13}) $$o en términos de las componentes covariantes de la cuadrivelocidad,
$$ \frac{\mathrm{d}P^1}{\mathrm{d}\tau}=\frac{q}{c} (U^0 F^{10} + U_2 F^{12} + U_3 F^{13}) $$De donde podemos proponer la generalización para calcular la cuadrifuerza:
$$ f^\mu=\frac{q}{c}U_\nu F^{\mu\nu} $$que llevada a un conjunto de cargas se puede escribir como:
$$ K^\mu=J_\nu F^{\mu\nu} $$donde en esta última expresión $K^\mu$ representa una densidad de fuerza (fuerza por unidad de volumen).
Nota: La componente temporal de la fuerza de Lorentz. La generalización presentada en el párrafo anterior apareció simplemente como un ejercicio de inducción, pero no lo es estrictamente. Si bien las componentes espaciales de la cuadrifuerza siquen la misma regla:
$$f^i=\frac{q}{c}F^{i\nu}U_\nu$$la componente temporal no tiene porque seguirla. ¿De dónde viene la idea de que podemos generalizarla a esta componente?. Hagamos el ejercicio inverso e intentemos determinar en términos de cantidaes convencionales que significa esta componente:
$$f^0=\frac{q}{c}F^{0\nu}U_\nu=\frac{q}{c}(F^{01} U_1+F^{02} U_2+F^{03} U_3)$$Si recordamos que $E^i=-F^{i0}$ y $U_i=-U^i$ entonces:
$$f^0=\gamma\frac{q}{c}\vec E\cdot \vec v$$en unidades luz y reconociendo que la trifuerza experimentada por una carga en presencia de un campo eléctrico es $\vec F=q\vec E$ obtenemos:
$$f_L^0=\gamma\vec F\cdot \vec v_L$$que es justamente lo que sabíamos sobre la componente temporal de la fuerza. Es decir la componente temporal de la fuerza de Lorentz no es otra cosa que la componente temporal de la cuadrifuerza de un campo eléctrico.
Ahora bien, sabemos también que $f_L^0=\mathrm{d}E/\mathrm{d}t$, siendo $E$ en este contexto, la energía total de la partícula. De allí que la componente cero de la fuerza de Lorentz es simplemente la potencia de la fuerza de Lorentz, la tasa a la cual se transfiere o remueve energía de la partícula.
Uno de los aspectos más oscuros de la teoría electromagnética es la interpretación intuitiva de cada una de las cantidades relevantes de la teoría. Si bien la interpretación de los campos es relativamente directa, están asociados con las fuerzas que experimentan partículas y corrientes eléctricas en presencia de ellos, la interpretación de los potenciales es un poco más complicada.
Si eso pasa con la electrodinámica vectorial en tres dimensiones, la electrodinámica covariante en cuatro dimensiones no es menos sencilla.
La fuerza de Lorentz sin embargo nos ofrece un raro vistazo a la geometría del tensor de Faraday en 4 dimensiones. Para entenderlo mejor podemos apelar a una cantidad clásica: el tensor de esfuerzos $\sigma_{ij}$.
Definición: Tensor de esfuerzos. El tensor de esfuerzos es una cantidad clásica usada en la mecánica del continuo para calcular la tracción producida en cualquier punto dentro de un sólido por el material vecino. La tracción $\vec{T}^{(\hat{n})}$ se define como la fuerza por unidad de área que experimenta el elemento de sólido en dirección de una superficie imaginaria definida por el vector $\hat{n}$. La tracción viene dada en términos del tensor de esfuerzos $\sigma_{ij}$ por el producto punto:
$$\vec{T}^{(\hat{n})}:\left(\sum_i n_i \sigma_{i1},\sum_i n_i \sigma_{i2},\sum_i n_i \sigma_{i3}\right)\equiv \hat{n}\cdot \sigma$$
La característica más interesante y relevante para nuestra discusión es que el producto punto de un vector (dirección) y un tensor de segundo rango (tensor de esfuerzos) produce un vector (tracción).
Particularmente interesantes son las componentes de la tracción en dirección de los vectores coordenados:
$$ \begin{array}{rcl} \vec{T}^{(e_1)} & : & (\sigma_{11},\sigma_{12},\sigma_{13})\\ \vec{T}^{(e_2)} & : & (\sigma_{21},\sigma_{22},\sigma_{23})\\ \vec{T}^{(e_2)} & : & (\sigma_{31},\sigma_{32},\sigma_{33}) \end{array} $$Escrito de esta manera, nos damos cuenta que las componentes del tensor de esfuerzo en la base $\hat{e}_1, \hat{e}_2, \hat{e}_3$ no son otra cosa que las componentes de las tracciones sobre la superficie en cada una de las direcciones. Esto le da una interpretación geométrica concreta a la tracción.
Con el tensor de Faraday sucede una situación similar. La cuadrifuerza de Lorentz $f^\mu$ juega aquí el papel de la tracción y el tensor de Faraday es como el tensor de esfuerzos:
$$ f_L^\mu= q F^{\mu\nu} U_{L\nu} $$De modo que uno puede ver las componentes del tensor de Faraday como los factores por los que hace falta multiplicar la velocidad para obtener las componentes de la cuadrifuerza sobre la partícula. En otros términos, del tensor de Faraday:
$$ F^{\mu\nu}: \left( \begin{array}{cccc} 0 & -E_x & -E_y & -E_z\\ E_x & 0 & -B_z & +B_y \\ E_y & +B_z & 0 & -B_x \\ E_z & -B_y & +B_x & 0 \\ \end{array} \right) $$podemos pasar al cuadrivector de fuerza escrito de la forma:
$$ \left( \begin{array}{c} \mathrm{d}E/\mathrm{d}t\\ F_{Lx}\\ F_{Ly}\\ F_{Lz} \end{array} \right) = q \left( \begin{array}{cccc} 0+ & (-v_{Lx})(-E_x)+ & (-v_{Ly})(-E_y) +& (-v_{Lz})(-E_z)\\ E_x+ & 0+ & (-v_{Ly})(-B_z)+ & (-v_{Lz})(+B_y) \\ E_y+ & (-v_{Lx})(+B_z)+ & 0+ & (-v_{Lz})(-B_x) \\ E_z+ & (-v_{Lx})(-B_y)+ & (-v_{Ly})(+B_x)+ & 0 \\ \end{array} \right) $$donde asumimos que sumamos sobre las filas de la matriz anterior (abusando de la notación).
Es decir, las filas del tensor de Faraday nos dan una idea de cuánto vale la fuerza electromagnética en una determinada dirección.
Ahora bien. Si hacemos por ejemplo $\vec v:(v_x,0,0)$, la fuerza de Lorentz queda:
$$ \left( \begin{array}{c} \mathrm{d}E/\mathrm{d}t\\ F_{Lx}\\ F_{Ly}\\ F_{Lz} \end{array} \right) = q \left( \begin{array}{c} 0 + v_x E_x \\ E_x + 0 \\ E_y - v_x B_z \\ E_z + v_x B_y \\ \end{array} \right) $$Es decir, si se proyecta el vector de Faraday sobre una dirección de la velocidad específica el negativo de la columna de correspondiente a la velocidad nos dará las componentes de la fuerza espacial magnética experimentada por la partícula. Es como si la proyección espacial del tensor de Faraday produjera las fuerzas electromagnéticas.
Por otro lado si la partícula esta quieta $U_L^\mu:(1,0,0,0)$ entonces:
$$ \left( \begin{array}{c} \mathrm{d}E/\mathrm{d}t\\ F_{Lx}\\ F_{Ly}\\ F_{Lz} \end{array} \right) = q \left( \begin{array}{c} 0 \\ E_x \\ E_y \\ E_z \\ \end{array} \right) $$de modo que podríamos interpretar la proyección temporal del vector de Faraday como relacionada con la fuerza eléctrica experimentada por la partícula.
En síntesis la dinámica de cargas y campos electromagnéticos clásicos se describen con:
\begin{eqnarray} {F^{\mu\nu}}_{,\mu} & = & \frac{4\pi}{c}J^\nu \;\; \mathrm{(Gauss-E)\;y\;(Ampère-Maxwell)}\\ F^{[\mu\nu,\alpha]} & = & 0 \;\; \mathrm{(Gauss-B)\;y\;(Maxwell-Faraday)}\\ {J^\mu}_{,\mu} & = & 0 \;\; \mathrm{(Continuidad)}\\ K^\mu & = & F^{\mu\nu} J_\nu \;\; \mathrm{(Fuerza\;de\;Lorentz)} \end{eqnarray}donde $F^{\mu\nu}\equiv\partial^\mu A^\nu-\partial^\nu A^\mu$ y el cuadripotencial cumple ${A^\mu}_{,\mu}=0$ (gauge de Lorenz).