Realizando una transformada de Fourier, se encuentra:
2ϵμνρσpρ˜Aσ(p)=pμAν(p)−pνAμ(p) .
Puedes ver cada lado de la ecuación como un operador que mapea un vector a un tensor antisimétrico, actuando sobre el vector ˜A en la izquierda y en el vector A a la derecha. Explícitamente, al definir
$\tilde{M}_{\mu\nu}^{~~\sigma}(p) = 2\epsilon_{\mu\nu\rho}^{~\sigma}p^\rho~,\\ M_{\mu\nu}^{~~~~\lambda}(p) = p_\mu \delta_ \nu^\lambda -p_\nu \delta_\mu^\lambda~,$
su ecuación se convierte en
$ \tilde{M}_{\mu\nu}^{~~\sigma}(p) \tilde{A}_\sigma(p) = M_{\mu\nu}^{~~\lambda}(p) A_\lambda(p)~.$
Ahora su problema se convierte en el de invertir uno de estos dos operadores, digamos ˜M para la definición (están relacionados entre sí por la contracción con el tensor épsilon, por lo que no hay gran diferencia en cuál de ellos se elige). Desgraciadamente no son invertibles, porque se puede comprobar fácilmente que el vector pμ se mapea al tensor antisimétrico evanescente. Esto no es tan sorprendente, sólo es equivalente a la afirmación de que las intensidades de campo son invariantes bajo tranformaciones gauge, pero aún así significa que estrictamente hablando no se pueden invertir estos operadores y encontrar ˜A en términos de A (o viceversa). A lo sumo se puede encontrar una solución hasta una transformación gauge es decir, se puede encontrar ˜Aμ hasta un desplazamiento de la forma pμC(p) donde C(p) es una función arbitraria de p .
Una vez que acepte que esto es lo máximo que puede hacer, podrá comprobar fácilmente que el operador
$\tilde{N}_\omega^{\mu\nu}(p)\,= \,b \,\frac{1}{p^2}\epsilon_\omega^{\mu\nu\tau}p_\tau$
-- con una elección adecuada del coeficiente numérico b que no trataré de arreglar pero que se puede arreglar con un cálculo directo -- cumple la tarea de casi invertir ˜M(p) es decir, existe un C(p) tal que
$\tilde{N}_\omega^{\mu\nu}(p)\tilde{M}_{\mu\nu}^{~~\sigma}(p) \tilde{A}_\sigma(p) = \tilde{A}_\omega(p) + p_\omega C(p)~.$
(Para fijar la constante numérica b hay que utilizar la identidad que relaciona la contracción de dos tensores épsilon con la suma de un montón de productos de deltas de Kronecher).
Así que la mejor respuesta posible a tu pregunta es -en el espacio del momento-
$\tilde{A}_\omega(p) = \tilde{N}_\omega^{\mu\nu}(p) M_{\mu\nu}^{~~\lambda}(p) A_\lambda(p) - p_\omega C(p)~.$
Nótese que esta relación NO es polinómica en el momento, porque el operador ˜N contiene un 1/p2 y, por lo tanto, no es local. Si se expresa en el espacio de posición contiene algo así como la inversa de un Laplaciano (o D'Alembertiano si está en firma lorentziana). No existe ninguna relación local entre un campo gauge y su dual EM, básicamente porque necesitas invertir una derivada para obtener uno del otro.