Let's consider free Maxwell theory:
$$ S_{Maxwell} = \int d^dx \; -\frac{1}{4e^2}F_{\mu\nu}F^{\mu\nu} $$
In such theory one can define monopole operator using path integral via correlation functions with other operators:
$$ \langle\mathcal{M}(x) \mathcal{O}_1(x_1)\dots \mathcal{O}_n(x_n)\rangle = \int D A_{\mu}\;\mathcal{O}_1(x_1)\dots \mathcal{O}_n(x_n) \;e^{-S_{Maxwell}} $$
Here path integration goes through configurations, where in point x included magnetic flux across surface, that surrounds point x of insertion of monopole.
I wanna to understand, how calculate correlation functions between monopole operators (in arbitrary d!).
For example:
$$ \langle\mathcal{M}(x)\rangle = ??? $$ $$ \langle\mathcal{M}(x)\mathcal{M}(y)\rangle = ??? $$
How to calculate such correlators? I think that simplest example in Maxwell theory in 3d (monopole in such theory rather instantons), due to duality with compact scalar.
I will be very appreciated for any suggestions/answers!
Edit(3D case):
Following section 8.1.1 David Tong: Lectures on Gauge Theory:
$$ \langle\mathcal{M}(x)\rangle = \int D\sigma \;e^{i\sigma(x)}\; e^{-\int d^3x\; \frac{e^2}{8\pi^2}\partial_\mu \sigma \partial^\mu \sigma} $$ $$ \langle\mathcal{M}(x)\mathcal{M}(y)\rangle = \int D\sigma \;e^{i\sigma(x)}e^{i\sigma(y)} \;e^{-\int d^3x\; \frac{e^2}{8\pi^2}\partial_\mu \sigma \partial^\mu \sigma} $$
Taking this Gaussian integrals using: $$ \int D\sigma \;e^{i\int d^3x\;J\sigma}\; e^{-\int d^3x\; \frac{e^2}{8\pi^2}\partial_\mu \sigma \partial^\mu \sigma} = e^{\frac{2\pi^2}{e^2}\int d^3x d^3y\; J(x)\Box^{-1}(x-y)J(y)} $$
I obtained:
$$ \langle\mathcal{M}(x)\rangle = e^{\frac{2\pi^2}{e^2}\Box^{-1}(0)} $$ $$ \langle\mathcal{M}(x)\mathcal{M}(y)\rangle =e^{\frac{4\pi^2}{e^2}\Box^{-1}(0)} e^{\frac{4\pi^2}{e^2} \Box^{-1}(x-y)} $$
$$ \Box^{-1}(x-y) = \frac{1}{|x-y|} $$
Maxwell theory in 3d is scale invariant.. Why such correlators doesn't have such property?
Is my calculations true? How to generalize this to higher dimmensions?