Magnetoelectric boundary simulated by a Chern-Simons-like model

In the literature, we can find a wide range of models where external $\delta$-like potentials coupled to fields are used to describe material boundaries. This subject has been of great interest for both experimental and theoretical reasons, once typical experimental setups are composed by material boundaries whose influence on the dinamics of fields cannot be disregarded. In this scenario, we can mention, for instance, the interactions between semi-transparent mirrors and point-like particles [1, 2, 3], investigations which concern the interaction between an atom and a $\delta$-like mirror [4, 5], the Casimir energy due to the presence of two $\delta$-like mirrors [6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17], and so on.

In the work of reference [18], in the context of the Casimir Effect, it was proposed a model composed by the Maxwell Lagrangian augmented by a term where the gauge field is coupled to $\delta$-potentials concentrated along parallel planes. The coupling exhibits a Chern-Simons-type form and describes semi-transparent mirrors. In the limit where the coupling constant between the mirrors and the field diverges, the model recovers the presence of perfectly conducting plates. In this specific limit, this model is the same as the one considered previously, in reference [19]. In reference [5] it was also considered the interaction between one of this delta-like plate (with Chern-Simons coupling) and a polarizable atom. In reference [20] this model was considered again to study the Casimir Effect.

However, there are some subjects not yet explored appropriately in the literature for such kind of Chern-Simons-type coupling with a delta-like potential. We can mention, for instance, the meaning of the modifications that the free photon propagator undergoes due to the presence of this single semi-transparent mirror, and the influence of this kind of surface in the dynamics of point-like field sources. In addition, this kind of plate exhibits magnetoelectric properties which were not explored in the literature and could be of some relevance for the analytical description of this metamaterials [21, 23, 22]. We hope that the present work could be a contribution to pave the way in this subject.

In this paper we consider some peculiarities of the model poposed in [5]. In section 2 we compute the modification undergone by the free photon propagator. We perform our analysis in a similar manner that was employed in Ref. [2], where another type of $\delta$-potential was considered. In section 3 we investigate the interaction between a stationary point-like charge and the mirror. We show that the classical image method is found as a particular case of the obtained result. In section 4 we calculate the electromagnetic field configuration induced by a stationary point-like charge in the presence of the mirror. We show that, in this setup, it emerges the field of a magnetic monopole, without the typical singularity introduced by semi-infinite Dirac strings. It is a new magnetoelectric effect and is an indication that the model could be of some relevance for the description of magnetoelectric metamaterials. In section (5) we consider the interaction between the plate and a quantum particle with spin. We show that it emerges a magnetoelectric anisotropic contribution to this energy linear in the spin. Section 6 is devoted for conclusions and final remarks.

In this paper we work in a $3+1$-dimensional Minkowski space-time with metric $\eta^{\rho\nu}=(1,-1,-1,-1)$. The Levi-Civita tensor is denoted by $\epsilon^{\rho\nu\alpha\beta}$ with $\epsilon^{0123}=1$.

In this section we consider a model composed by the Maxwell Lagrangian augmented by a $\delta$-potential concentrated along a plane and with a Chern-Simons-like coupling. The model describes a kind of partially reflective mirror. Without loss of generality, and for convenience, we will consider the mirror perpendicular to the $x^{3}$ axis and placed on the plane $x^{3}=a$. The corresponding model is given by the following Lagrangian density,

	$\displaystyle{\cal{L}}$	$\displaystyle=$	$\displaystyle-\frac{1}{4}F_{\rho\nu}F^{\rho\nu}-\frac{1}{2\xi}\left(\partial_{\rho}A^{\rho}\right)^{2}$		(1)
			$\displaystyle-\frac{1}{2}\mu v^{\rho}\epsilon_{\rho\nu\alpha\beta}A^{\nu}\partial^{\alpha}A^{\beta}\delta\left(x^{3}-a\right)-J^{\rho}A_{\rho}\ ,$

where $A^{\rho}$ is the photon field, $F^{\rho\nu}=\partial^{\rho}A^{\nu}-\partial^{\nu}A^{\rho}$ is the field strength, $\xi$ is a gauge fixing parameter, $j^{\rho}$ is an external source, $v^{\rho}=\eta^{\rho}_{\ 3}=\left(0,0,0,1\right)$ is the normal vector to the surface and $\mu>0$ is a dimensionless coupling constant which accounts for the degree of transparency of the mirror.

There is a generalization of the Chern-Simons electrodynamics in 3 + 1 dimensions, the so called Carroll-Field-Jackiw model [24], which exhibits a Lorentz symmetry breaking controlled by a single background 4-vector. The Lagrangian (1) can be seen as a kind of the Carroll-Field-Jackiw Electrodynamics where the Chern-Simons-like term is defined just along a plane perpendicular to the background vector.

The model (1) exhibits a $\delta$-type divergence on the mirror. In order to understand the role of this divergence on the field strength tensor, it is convenient to write explicitly the dynamical equations, as follows

$\displaystyle\partial_{\rho}F^{\rho\nu}+\mu\delta\left(x^{3}-a\right){\tilde{F}}^{3\nu}=J^{\nu}\ ,$

(2)

where we defined ${\tilde{F}}^{\rho\nu}=\frac{1}{2}\epsilon_{\rho\nu\alpha\beta}F^{\alpha\beta}$. For the electric and magnetic fields, we obtain

	$\displaystyle{\bf{\nabla}}\cdot{\bf{E}}$	$\displaystyle=$	$\displaystyle J^{0}-\mu\delta\left(x^{3}-a\right)\left({\bf{B}}_{\perp}\cdot{\hat{z}}\right)\ ,$		(3)
	$\displaystyle{\bf{\nabla}}\times{\bf{B}}$	$\displaystyle=$	$\displaystyle{\bf{J}}+\frac{\partial{\bf{E}}}{\partial t}-\mu\delta\left(x^{3}-a\right)\left({\bf{E}}_{\parallel}\times{\hat{z}}\right)\ ,$		(4)

where we defined the vectors perpendicular and parallel to the plate, ${\bf{B}}_{\perp}=\left(0,0,B^{3}\right)$, ${\bf{E}}_{\parallel}=\left(E^{1},E^{2},0\right)$, respectively.

At this point some comments are in order. The presence of the delta-like potential leads to additional terms in the dynamical equations, which depend on the fields and not on their derivatives. In Eq. (3) the delta-like term is proportional to the component of the magnetic field perpendicular to the mirror and can be interpreted as an extra source for the divergence of the electric field. In the same way, in Eq. (4) the delta-like term is an additional source for the curl of the magnetic field and depends on the components of the electric field parallel to the surface. For the potential considered in Ref. [2] the delta-like contributions for the dynamical field equations exhibit different behaviors in comparison with (2), (3) and (4) and depend on the derivatives of the fields.

Neglecting surface terms, one can write

$\displaystyle{\cal{L}}=\frac{1}{2}A^{\rho}{\cal{O}}_{\rho\nu}A^{\nu}-J^{\nu}A_{\nu}\ ,$

(5)

where ${\cal{O}}_{\rho\nu}$ is a differential operator.

For our purpose, it is convenient to split ${\cal{O}}_{\rho\nu}$ into two parts, one corresponding to the usual photon operator and the other one corresponding to the $\delta$-like term, as follows

$\displaystyle{\cal{O}}_{\rho\nu}={\cal{O}}^{(0)}_{\rho\nu}+\Delta{\cal{O}}_{\rho\nu}\ ,$

(6)

with

	$\displaystyle{\cal{O}}^{(0)}_{\rho\nu}$	$\displaystyle=$	$\displaystyle\eta_{\rho\nu}\Box\ ,$		(7)
	$\displaystyle\Delta{\cal{O}}_{\rho\nu}$	$\displaystyle=$	$\displaystyle-{\mu}\delta\left(x^{3}-a\right)\epsilon_{3\rho\alpha\nu}\partial_{\parallel}^{\alpha}\ ,$		(8)

where we defined the operator $\Box=\partial^{\alpha}\partial_{\alpha}$ and used a gauge where $\xi=1$. We notice that the derivative in (8) is defined only in the Minkowski coordinates parallel to the mirror, namely, $\partial_{\parallel}^{\alpha}=\left(\partial^{0},\partial^{1},\partial^{2}\right)$, because of the fixed index in the Levi-Civita tensor.

The free photon propagator satisfies the differential equation ${\cal{O}}^{(0)\rho\nu}(x)G^{(0)}_{\nu\beta}\left(x,y\right)=\eta^{\rho}_{\ \beta}\delta^{4}\left(x-y\right)$ and is given in the Feynman gauge, $(\xi=1)$, by

$\displaystyle G^{(0)}_{\rho\nu}\left(x,y\right)=-\eta_{\rho\nu}\int\frac{d^{4}p}{\left(2\pi\right)^{4}}\frac{e^{-ip\cdot\left(x-y\right)}}{p^{2}}\ .$

(9)

The propagator corresponding to the model (5) must satisfy the equation ${\cal{O}}^{\rho\nu}(x)G_{\nu\lambda}\left(x,y\right)=\eta^{\rho}_{\ \lambda}\delta^{4}\left(x-y\right)$ and can be obtained recursively [2, 1, 17, 3, 16]. It is simple to check that

	$\displaystyle G_{\rho\nu}\left(x,y\right)=G^{(0)}_{\rho\nu}\left(x,y\right)$
	$\displaystyle-\int d^{4}z\ G_{\rho\gamma}\left(x,z\right)\Delta{\cal{O}}^{\gamma\sigma}\left(z\right)G^{(0)}_{\sigma\nu}\left(z,y\right)\ ,$		(10)

For convenience let us write $G_{\rho\nu}\left(x,y\right)$ and $G^{(0)}_{\rho\nu}\left(x,y\right)$ as Fourier integrals in the parallel coordinates, as follows

	$\displaystyle G_{\rho\nu}\left(x,y\right)=\int\frac{d^{3}p_{\parallel}}{\left(2\pi\right)^{3}}\ {\cal{G}}_{\rho\nu}\left(x^{3},y^{3};p_{\parallel}\right)e^{-ip_{\parallel}\cdot\left(x_{\parallel}-y_{\parallel}\right)},$		(11)
	$\displaystyle G^{(0)}_{\rho\nu}\left(x,y\right)=\int\frac{d^{3}p_{\parallel}}{\left(2\pi\right)^{3}}\ {\cal{G}}^{(0)}_{\rho\nu}\left(x^{3},y^{3};p_{\parallel}\right)e^{-ip_{\parallel}\cdot\left(x_{\parallel}-y_{\parallel}\right)},$		(12)

where we defined $x_{\parallel}^{\rho}=\left(x^{0},x^{1},x^{2}\right)$ and $p_{\parallel}^{\rho}=\left(p^{0},p^{1},p^{2}\right)$ as being the coordinates and momenta parallel to the mirror, respectively. The functions ${\cal{G}}_{\rho\nu}\left(x^{3},y^{3};p_{\parallel}\right)$ and ${\cal{G}}^{(0)}_{\rho\nu}\left(x^{3},y^{3};p_{\parallel}\right)$ are usually known as reduced Green’s functions [2, 1, 17, 3, 16].

Using the fact that [25]

$\displaystyle\int\frac{dp^{3}}{\left(2\pi\right)}\frac{e^{ip^{3}\left(x^{3}-y^{3}\right)}}{p^{2}}=-\frac{e^{-\sigma\mid x^{3}-y^{3}\mid}}{2\sigma}\ ,$

(13)

with $p^{3}$ standing for the momentum perpendicular to the mirror and $\sigma=\sqrt{-p_{\parallel}^{2}}$, from Eq. (9) we have for the reduced photon propagator

$\displaystyle{\cal{G}}^{(0)}_{\rho\nu}\left(x^{3},y^{3};p_{\parallel}\right)$

$\displaystyle=$

$\displaystyle\eta_{\rho\nu}\frac{e^{-\sigma\mid x^{3}-y^{3}\mid}}{2\sigma}\ .$

(14)

Substituting (8) into (2), using Eqs. (11), (12) and (14), and performing some simple integrals, we obtain

		$\displaystyle{\cal{G}}_{\rho\nu}$	$\displaystyle\left(x^{3},y^{3};p_{\parallel}\right)=\eta_{\rho\nu}\frac{e^{-\sigma\mid x^{3}-y^{3}\mid}}{2\sigma}$		(15)
			$\displaystyle-i\mu\ \epsilon_{3\ \ \nu}^{\ \gamma\alpha}\ {\cal{G}}_{\rho\gamma}\left(x^{3},a;p_{\parallel}\right)p_{\parallel\alpha}\frac{e^{-\sigma\mid y^{3}-a\mid}}{2\sigma}\ .$

From the above equation, the reduced Green’s function ${\cal{G}}_{\rho\nu}\left(x^{3},y^{3};p_{\parallel}\right)$ must be obtained recursively. Evaluating Eq. (15) for $y^{3}=a$ and performing some simple manipulations, we arrive at

	$\displaystyle{\cal{G}}_{\rho\gamma}\left(x^{3},a;p_{\parallel}\right)\left(\eta_{\ \nu}^{\gamma}+\frac{i\mu}{2\sigma}\ \epsilon_{3\ \ \nu}^{\ \gamma\alpha}\ p_{\parallel\alpha}\right)=\eta_{\rho\nu}\frac{e^{-\sigma\mid x^{3}-a\mid}}{2\sigma}\ .$		(16)
			(17)

Now we multiply both sides of (16) by the inverse of the operator $\eta_{\ \nu}^{\gamma}+\frac{i\mu}{2\sigma}\ \epsilon_{3\ \ \nu}^{\ \gamma\alpha}\ p_{\parallel\alpha}$, what leads to

	$\displaystyle{\cal{G}}_{\rho\gamma}\left(x^{3},a;p_{\parallel}\right)=\frac{e^{-\sigma\mid x^{3}-a\mid}}{2\sigma}\Biggl{[}\eta_{\rho\gamma}$
	$\displaystyle-\frac{\mu^{2}}{\mu^{2}+4}\left(\eta_{\parallel\rho\gamma}-\frac{p_{\parallel\rho}p_{\parallel\gamma}}{p_{\parallel}^{2}}\right)-\frac{2i\mu}{\sigma\left(\mu^{2}+4\right)}\epsilon_{3\rho\alpha\gamma}p_{\parallel}^{\alpha}\Biggr{]}\ ,$

where we defined $\eta_{\parallel}^{\alpha\beta}=\eta^{\alpha\beta}-\eta^{\alpha 3}\eta^{\beta 3}$.

Substituting Eq. (2) in (15) and using the Eq. (11), the modified photon propagator due to the presence of the semi-transparent mirror reads

		$\displaystyle{{G}}_{\rho\nu}\left(x,y\right)$	$\displaystyle=\int\frac{d^{3}p_{\parallel}}{\left(2\pi\right)^{3}}\ e^{-ip_{\parallel}\cdot\left(x_{\parallel}-y_{\parallel}\right)}\Biggl{\{}\eta_{\rho\nu}\frac{e^{-\sigma\mid x^{3}-y^{3}\mid}}{2\sigma}$		(19)
			$\displaystyle-\frac{e^{-\sigma\left(\mid x^{3}-a\mid+\mid y^{3}-a\mid\right)}}{4+\mu^{2}}\Biggl{[}\frac{\mu^{2}}{2\sigma}\left(\eta_{\parallel\rho\nu}-\frac{p_{\parallel\rho}p_{\parallel\nu}}{p_{\parallel}^{2}}\right)$
			$\displaystyle+\frac{i\mu}{p_{\parallel}^{2}}\epsilon_{3\alpha\rho\nu}p_{\parallel}^{\alpha}\Biggr{]}\Biggr{\}}\ .$

The propagator (19) is continuous and well defined all over the space, except at coincident points, where $x^{\mu}=y^{\mu}$. The first term on the right-hand side is just the free photon propagator. The second and third terms are corrections which account for the presence of the semi-transparent mirror. As an important check, we point out that by taking the limit $\mu\rightarrow\infty$ in Eq. (19) we recover the well-known photon propagator due to the presence of a single perfect mirror [26].

In the present section we consider the interaction energy between a stationary point-like charge and the semi-transparent mirror studied in the previous section. Since we have a quadratic Lagrangian in the field variables, as discussed in references [27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38], the contribution due to the external sources to the ground state energy of the system is given by

$\displaystyle E=\frac{1}{2T}\int d^{4}x\ d^{4}y\ J^{\rho}\left(x\right){{G}}_{\rho\nu}\left(x,y\right)J^{\nu}\left(y\right)\ ,$

(20)

where $T$ is the time variable and it is implicit the limit $T\to\infty$.

With no loss of generality, we choose a charge placed at position ${\bf b}=\left(0,0,b\right)$. The corresponding external source reads

$\displaystyle J^{\rho}\left(x\right)=q\eta^{\rho 0}\delta^{3}\left({\bf x}-{\bf b}\right)\ ,$

(21)

where the parameter $q$ is the electric charge intensity.

We notice that the first term on the right-hand side of expression (19) comes from the free photon propagator (without the presence of the mirror) and, therefore, it does not contribute to the interaction energy between the charge and the mirror. Indeed, this contribution does not depend on the distance between the charge and the mirror, it is present even in the absence of the mirror and provides the charge self energy. Thus, only the last two terms of propagator contribute to the interaction energy.

Substituting Eqs. (21) and (19) in (20), discarding self interacting terms and performing the integrals in $d^{3}{\bf x}$, $d^{3}{\bf y}$, $dx^{0}$, $dp^{0}$, $dy^{0}$, we obtain

$\displaystyle E_{MC}\left(R,\mu\right)=-\frac{q^{2}}{16\pi^{2}}\left(\frac{\mu^{2}}{4+\mu^{2}}\right)\int d^{2}{\bf{p}}_{\parallel}\frac{e^{-2R\sqrt{{\bf{p}}_{\parallel}^{2}}}}{\sqrt{{\bf{p}}_{\parallel}^{2}}}\ ,$

(22)

where $R=\mid a-b\mid$ is the distance between the mirror and the charge. The sub-index $MC$ means that we have the interaction energy between the mirror and the charge.

Using polar coordinates and performing some simple integrations, the interaction energy reads

$\displaystyle E_{MC}\left(R,\mu\right)=-\frac{q^{2}}{16\pi R}\left(\frac{\mu^{2}}{4+\mu^{2}}\right)\ .$

(23)

Eq. (23) is the exact result for the interaction energy between a point-like charge and the two-dimensional semi-transparent mirror, described by the model (1). It is interesting to realize that the energy (23) still exhibits a Coulomb-like behavior. The role played by the coupling constant $\mu$ is just to attenuate the image charge ($q_{im}$) by an uniform multiplicative factor, namely, $q_{im}=q\frac{\mu^{2}}{4+\mu^{2}}$. As expected, this energy vanishes in the limit $\mu\rightarrow 0$, where we have no mirror present. In the limit $\mu\rightarrow\infty$, we have the gauge field subjected to the boundary conditions imposed by a perfect mirror and

$\displaystyle\lim_{\mu\rightarrow\infty}E_{MC}\left(R,\mu\right)=-\frac{q^{2}}{16\pi R}\ ,$

(24)

what recovers the image method. Thus, the model (1) describes the interaction obtained via image method in the limit of a perfect conductor.

The force between the point-like charge and the mirror is given by

	$\displaystyle F_{MC}\left(R,\mu\right)=-\frac{\partial}{\partial R}E_{MC}\left(R,\mu\right)=-\frac{q^{2}}{16\pi R^{2}}\frac{\mu^{2}}{4+\mu^{2}}\ ,$		(25)
			(26)

which is always negative and, therefore, has an attractive behavior. In the figure (1) we have a plot for the force (25) multiplied by $\frac{16\pi}{q^{2}}$ as a function of $R$ and $\mu$.

Here, one comment is in order. The parameter $\mu$ is dimensionless. In this case, by a dimensional analysis, we could infer the $q^{2}/R$ behavior for the energy (23) previously, just before we perform the calculations. The numerical factor $-1/(16\pi)$, including the minus sign, and the dependence on $\mu$ cannot be obtained before the calculations are performed.

The dimensionless numerical factor $\mu^{2}/(4+\mu^{2})$ lies in the range $0\leq\mu^{2}/(4+\mu^{2})\leq 1$, as we can see in figure (2), so the force (25) as well as the energy (23), attain their minimum values for a perfectlly conducting plate.

In this section we calculate the electromagnetic field configuration produced by a stationary point-like charge in the vicinity of the semi-transparent magnetoelectric mirror. We start by considering the particular solution for the electromagnetic four-potential in the presence of an external source

$\displaystyle A_{\rho}\left(x\right)=\int d^{4}y\ G_{\rho\nu}\left(x,y\right)J^{\nu}\left(y\right)\ ,$

(27)

where $G^{\rho\nu}\left(x,y\right)$ is the propagator (19) and $J_{\nu}\left(y\right)$ is given by (21).

From now on, with no loss of generality, we shall take a coordinate system where $b^{3}>0$ and $a=0$. It means that the mirror is placed along the plane $x^{3}=0$ and the particle is placed at a given point on the positive $x^{3}$ axis.

Each term on the right hand side of (19) leads to a different contribution to the 4-vector potential in (27). The first term is the one obtained from the free photon propagator (without the presence of the mirror),

	$\displaystyle\int d^{4}y\frac{d^{3}p_{\parallel}}{\left(2\pi\right)^{3}}\ e^{-ip_{\parallel}\cdot(x_{\parallel}-y_{\parallel})}\eta_{\rho\nu}\frac{e^{-\sigma\mid x^{3}-y^{3}\mid}}{2\sigma}q\eta^{\nu 0}\delta^{3}({\bf y}-{\bf b})$		(28)
			(29)
	$\displaystyle=q\eta_{\rho}^{\ 0}\int\frac{d^{2}p_{\parallel}}{\left(2\pi\right)^{2}}\ e^{i{\bf p}_{\parallel}\cdot{\bf x}_{\parallel}}\frac{e^{-\sqrt{{\bf p}_{\parallel}^{2}}\mid x^{3}-y^{3}\mid}}{2\sqrt{{\bf p}_{\parallel}^{2}}}\ .\ \ \ \$		(30)

Defining as $\varphi$ the angle between the vectors ${\bf p}_{\parallel}$ and ${\bf x}_{\parallel}$ in the ${\bf p}_{\parallel}$ space and $s=\sqrt{{\bf p}_{\parallel}^{2}}$, we can write the above integral in cylindrical coordinates as follows

	$\displaystyle\int\!\!\!$	$\displaystyle d^{4}y$	$\displaystyle\!\!\!\frac{d^{3}p_{\parallel}}{\left(2\pi\right)^{3}}\ e^{-ip_{\parallel}\cdot(x_{\parallel}-y_{\parallel})}\eta_{\rho\nu}\frac{e^{-\sigma\mid x^{3}-y^{3}\mid}}{2\sigma}q\eta^{\nu 0}\delta^{3}({\bf y}-{\bf b})$		(31)
		$\displaystyle=$	$\displaystyle\frac{q\eta_{\rho}^{\ 0}}{(2\pi)^{2}}\int_{0}^{\infty}ds\int_{0}^{2\pi}d\varphi\ \frac{e^{-s|x^{3}-b^{3}|}}{2}e^{is|{\bf x}_{\parallel}|\cos\varphi}$		(33)
		$\displaystyle=$	$\displaystyle\frac{q}{4\pi}\frac{\eta_{\rho}^{\ 0}}{\sqrt{|{\bf x}_{\parallel}|^{2}+|x^{3}-b^{3}|^{2}}}=\frac{q}{4\pi}\frac{\eta_{\rho}^{\ 0}}{|{\bf x}-{\bf b}|}\ ,$		(35)

what is the well known potential of a stationary electric charge $q$ placed at ${\bf b}$.

The second term of the propagator (19) gives the following contribution to the 4-potential (27)

	$\displaystyle\int d^{4}y\frac{d^{3}p_{\parallel}}{\left(2\pi\right)^{3}}-\frac{\mu^{2}}{4+\mu^{2}}\left(\eta_{\parallel\rho\nu}-\frac{p_{\parallel\rho}p_{\parallel\nu}}{p_{\parallel}^{2}}\right)$		(36)
			(37)
	$\displaystyle e^{-ip_{\parallel}\cdot(x_{\parallel}-y_{\parallel})}\frac{e^{-\sigma\left(\mid x^{3}-a\mid+\mid y^{3}-a\mid\right)}}{2\sigma}q\eta^{\nu 0}\delta^{3}({\bf y}-{\bf b})$		(38)
			(39)
	$\displaystyle=-\frac{q\eta_{\rho}^{\ 0}}{(2\pi)^{2}}\frac{\mu^{2}}{4+\mu^{2}}\int_{0}^{\infty}\!\!\!ds\int_{0}^{2\pi}\!\!\!d\varphi\ \frac{e^{-s(|x^{3}|+|b^{3}|)}}{2}e^{is|{\bf x}_{\parallel}|\cos\varphi}$		(40)
			(41)
	$\displaystyle=-\frac{q\eta_{\rho}^{\ 0}}{4\pi}\frac{\mu^{2}}{4+\mu^{2}}\frac{1}{|{\bf x}+sgn(x^{3}){\bf b}|}\ .\ \ \ \ \ \$		(42)

where we defined the sign function $sgn(x>0)=1$; $sgn(x<0)=-1$ and $sgn(0)=0$.

The contribution of the third term in (19) to the field (27) must be considered cautiously, as follows

	$\displaystyle-\frac{i\mu}{4+\mu^{2}}\int d^{4}y\frac{d^{3}p_{\parallel}}{\left(2\pi\right)^{3}}e^{-ip_{\parallel}\cdot\left(x_{\parallel}-y_{\parallel}\right)}\frac{e^{-\sigma\left(\mid x^{3}\mid+\mid y^{3}\mid\right)}}{2\sigma}$		(43)
			(44)
	$\displaystyle\times\epsilon_{3\alpha\rho\nu}\frac{p_{\parallel}^{\alpha}}{p_{\parallel}^{2}}q\delta^{3}({\bf y}-{\bf b})$		(45)
			(46)
	$\displaystyle=\frac{q\ \mu}{4+\mu^{2}}\epsilon_{3j\rho 0}\left(\frac{\partial|{\bf x}_{\parallel}|}{\partial x^{j}}\right)\int_{0}^{\infty}\frac{ds}{2\pi}\frac{e^{-s\left(\mid x^{3}\mid+\mid b^{3}\mid\right)}}{s}$		(47)
			(48)
	$\displaystyle\frac{\partial}{\partial|{\bf x}_{\parallel}|}\int_{0}^{2\pi}\frac{d\varphi}{2\pi}e^{is|{\bf x}_{\parallel}|\cos\varphi}$		(49)
			(50)
	$\displaystyle=\frac{q}{2\pi}\frac{\mu\ \epsilon_{0\rho j3}}{4+\mu^{2}}\frac{x_{\parallel}^{j}}{|{\bf x}_{\parallel}|^{2}}\Bigg{[}1-\frac{|x^{3}|+b^{3}}{\sqrt{{\bf x}_{\parallel}^{2}+[x^{3}+sgn(x^{3})b^{3}]^{2}}}\Bigg{]}.$		(51)

Collecting the terms (31), (36) and (43), using cylindrical coordinates: $\rho=|{\bf x}_{\parallel}|$, $\phi=arctan(\rho/x^{3})$ and performing some simple algebra, we have the potencial for the mirror-charge system

	$\displaystyle A^{0}({\bf x})$	$\displaystyle=$	$\displaystyle\frac{q}{4\pi}\Bigg{[}\frac{1}{|{\bf x}-{\bf b}|}-\frac{\mu^{2}}{4+\mu^{2}}\frac{1}{|{\bf x}+sgn(x^{3}){\bf b}|}\Bigg{]}$		(52)
	$\displaystyle{\bf A}({\bf x})$	$\displaystyle=$	$\displaystyle\frac{q}{4\pi}\frac{2\mu}{4+\mu^{2}}\Bigg{[}1-sgn(x^{3})\frac{x^{3}+sgn(x^{3})b^{3}}{|{\bf x}+sgn(x^{3}){\bf b}|}\Bigg{]}\frac{\hat{\phi}}{\rho}\ .$		(54)
					(55)

from which we obtain the fields

	$\displaystyle{\bf E}\left({\bf{x}}\right)$	$\displaystyle=$	$\displaystyle\frac{q}{4\pi}\frac{{\bf x}-{\bf b}}{|{\bf x}-{\bf b}|^{3}}-\frac{q}{4\pi}\left(\frac{\mu^{2}}{4+\mu^{2}}\right)\frac{{\bf x}+sgn(x^{3}){\bf b}}{|{\bf x}+sgn(x^{3}){\bf b}|^{3}}$		(56)
	$\displaystyle{\bf B}\left({\bf{x}}\right)$	$\displaystyle=$	$\displaystyle\frac{q}{4\pi}\left(\frac{2\mu}{4+\mu^{2}}\right)sgn(x^{3})\frac{{\bf x}+sgn(x^{3}){\bf b}}{|{\bf x}+sgn(x^{3}){\bf b}|^{3}}\ .$		(59)

Here, some comments are in order. In the side of the mirror where the charge $q$ is placed, $x^{3}>0$ and the field $A^{0}({\bf x})$ in (52) is the potential produced by the charge itself, placed at ${\bf b}$, and an image charge with intensity $-q\mu^{2}/(4+\mu^{2})$ and placed at position $-{\bf b}$. In the opposite side, $x^{3}<0$ and the field $A^{0}({\bf x})$ reduces to the one produced by a charge placed at ${\bf b}$ and with intensity $q[1-\mu^{2}/(4+\mu^{2})]$. The field $A^{0}({\bf x})$ in (52) gives the electric field in (56).

In the side of the mirror where the charge is placed, the vector potential ${\bf A}({\bf x})$ in (52) is the one related to a Dirac monopole with intensity $2q\mu/(4+\mu^{2})$, placed at $-{\bf b}$ and whose corresponding semi-infinite Dirac string is lying on the axis $-\infty<x^{3}<-b^{3}$, which is in the opposite side of the mirror. So, we do not have the typical divergence commonly associated with the string in the description of magnetic monopoles. In the opposite side of the mirror where the charge is placed, the vector potential in (52) is the same as the one produced by a Dirac monopole placed at ${\bf b}$, with intensity $-2q\mu/(4+\mu^{2})$ and related to a semi-infinite Dirac string lying along the axis $b^{3}<x^{3}<\infty$. Once again, there is no divergence for the potential imposed by the string. The field ${\bf A}({\bf x})$ in (52) gives the monopole magnetic field in (56) which does not exhibit any divergence usually imposed by Dirac strings (a delta-like magnetic monopole) [39, 40, 41].

It is important to reinforce that the flux of the magnetic field through a closed surface is always zero, for surfaces placed in any side of the mirror. The magnetic flux is also zero through surfaces which encloses a portion of the plate as well. It is due to the fact that we do not have a real magnetic monopole in any situation. All the magnetic monopoles are virtual ones, lying on the opposite side of the mirror were the field is taken.

In the limit $\mu\to\infty$, the field $A^{0}({\bf x})$ in (52) and the electric field in (56) reduce to the ones obtained for a perfect conductor, and the vector potential ${\bf A}({\bf x})$ in (52) vanishes, as well as the magnetic field in (56). So we recover the same results obtained from a perfect conductor in the limit $\mu\to\infty$.

When $\mu=0$, we recover the case with just a single point-like charge, where we have only a coulombian electric field, as it should be.

It is very interesting to notice that the modulus of the magnetic field connected to the magnetic image monopole, given by the second expression in Eq. (56), can be greater or lower than the modulus of the electric field of the image charge, given by the second term on the right hand side of the first Eq. (56) . In fact, from Eqs. (56), we have

$$|{\bf E}_{image}|=\frac{\mu}{2}|{\bf B}|\ .$$

(60)

So when $\mu=2$, the image charge and image monopole produce fields with the same intensities. When $\mu>2$ the image electric field dominates over the image magnetic field. For $\mu<2$, we have the opposite situation.

In figure (2) we have the graphic for the attenuation factor $\mu^{2}/(4+\mu^{2})$ for the image electric charge. It is a monotonic function of the coupling constant $\mu$. In figure (3) we have the attenuation factor $2\mu/(4+\mu^{2})$ for the image magnetic monopole. It is always positive, equal to zero for $\mu=0$ and $\mu=\infty$, and attains its maximum value, $1/2$, for $\mu=2$.

In this section we consider the interaction energy between a charged quantum particle with spin $1/2$ and the plate. This kind of interaction is known for standard situations, but the novelty in this case relies on the fact that we have a magnetoelectric plate. So, for a particle with spin, we shall have a cross interaction between the charge and spin of the particle due to its mirrored image in the magnetoelectric plate.

We shall take a simplified model, where the particle is described by the well known Pauli spin theory. We shall also take the particle as fixed, so its quantum state is just described by spin variables. The quantum states are denoted by $|+\rangle$ and $|-\rangle$, standing for spin up and down, respectively, with the ${\hat{z}}$ taken as the quantization axis.

Let us start by considering the interaction energy between the magnetoelectric plate and a classical charged particle with magnetic moment ${\bf m}$. In this case, the field source which accounts for the presence of such a particle is

$$J_{C,MD}^{\mu}({\bf x})=q\eta^{\mu 0}\delta^{3}({\bf{x}}-{\bf{b}})+\epsilon^{0\mu jk}m_{j}\left[\partial_{k}\delta^{3}({\bf{x}}-{\bf{b}})\right]\ .$$

(61)

The first term on the right hand side of (61) is, obviously, the source for a stationary point-like charge placed at ${\bf x}={\bf b}$. The second term is the source for a stationary magnetic dipole ${\bf m}$ placed at ${\bf x}={\bf b}$, as discussed in the appendix. From now on, with no loss of generality, we shall take ${\bf b}=(0,0,b)$.

Substituting Eq. (61) in (20) and performing some straightforward calculations, similar to the ones of the previous sections, we are taken to the interaction energy

	$\displaystyle E$	$\displaystyle=$	$\displaystyle\left[-\frac{q^{2}}{16\pi\mid b\mid}+\frac{\left[{\bf{m}}^{2}+\left({\bf{m}}\cdot{\hat{z}}\right)^{2}\right]}{64\pi\mid b\mid^{3}}\right]\left(\frac{\mu^{2}}{4+\mu^{2}}\right)$		(62)
			$\displaystyle-\frac{q\left({\bf{m}}\cdot{\hat{z}}\right)}{8\pi\mid b\mid^{2}}\left(\frac{\mu}{4+\mu^{2}}\right)\ .$		(63)

The first two terms on the right hand side of Eq. (62) are the direct interactions, that is, the charge $\times$ image-charge and the magnetic dipole $\times$ image-magnetic dipole ones, respectively. In comparison with the case of a perfect mirror, we have an atenuation factor $\mu^{2}/(4+\mu^{2})$ for these first two terms. The last term in (62) is a cross interaction due to the simultaneous presence of, both, the charge and the magnetic moment. Notice that in the limit of a perfect conducting plate, where $\mu=\infty$, the last term in (62) vanishes and the attenuation factor in the first terms equals to one.

On the right hand side of (62), the first term is proportional to the distance to the power $-1$ while the second term is proportional to the distance to the power $-3$. The last term exhibits a dependence on the distance to the power $-2$ so, for larger distances, this crossing term is more relevant in comparison with the direct magnetic dipole $\times$ image-magnetic dipole interaction.

For a quantum particle with spin $1/2$, the magnetic dipole is a quantum operator proportional to the spin operator ${\bf m}=\gamma{\bf S}$, where $\gamma$ is a constant of proportionality. In this sense, the Hamiltonian operator of the system is

	$\displaystyle H$	$\displaystyle=$	$\displaystyle\left[-\frac{q^{2}}{16\pi\mid b\mid}+\frac{\gamma^{2}\left[{\bf S}^{2}+\left(S_{z}\right)^{2}\right]}{64\pi\mid b\mid^{3}}\right]\left(\frac{\mu^{2}}{4+\mu^{2}}\right)$		(64)
			$\displaystyle-\frac{q\gamma S_{z}}{8\pi\mid b\mid^{2}}\left(\frac{\mu}{4+\mu^{2}}\right)\ .$		(66)

The constant $\gamma$ is proportional to the charge $q$ and both have the same signal.

A general spin-$1/2$ state is given by

$$|\psi\rangle=\cos\left(\theta/2\right)e^{-i\varphi/2}|+\rangle+\sin\left(\theta/2\right)e^{i\varphi/2}|-\rangle\ .$$

(67)

where $\theta$ and $\varphi$ can be interpreted as the polar and azimuthal angels, respectively, associated with the mean spin direction of the state taken the origin at the particle position.

The mean value of the Hamiltonian operator (64) in the general spin-$1/2$ state (67) is

	$\displaystyle\langle\psi|H|\psi\rangle$	$\displaystyle=$	$\displaystyle\left[-\frac{q^{2}}{16\pi\mid b\mid}+\frac{\gamma^{2}}{64\pi\mid b\mid^{3}}\right]\left(\frac{\mu^{2}}{4+\mu^{2}}\right)$		(68)
			$\displaystyle-\frac{q\gamma\cos(\theta)}{16\pi\mid b\mid^{2}}\left(\frac{\mu}{4+\mu^{2}}\right)\ .$		(70)

The first term on the right hand side of (68) leads to a Coulombian force between the particle and the plate, as discussed in section (3). There is a remarkable difference between the second term in (68) and its classical counterpart. In the case of a quantum particle with spin, the contribution to the energy coming from the second term does not depend on the mean orientation of the spin with respect to the plate. In the case of a classical particle with magnetic moment, this contribution depends on the angle between the magnetic moment of the particle and the normal vector do the plate, as one can see in (62), with the scalar product ${\bf{m}}\cdot{\hat{z}}$.

The third term on the right hand side of (68) is the most interesting one. It is due to the magnetoelectric properties of the plate and is a $\theta$-dependent contribution, the mean angle between the spin of the state (67) and the normal to the plate. From the third term in (68) we have a mean contribution to the force between the plate and particle which falls with the distance between them slower in comparison with the second term in (68). Besides, this force also depends on the $\theta$ angle.

In addition to a force, from the third term in (68), we also obtain a mean torque acting on the quantum particle with respect to the angle $\theta$. This torque is always non-positive, attains its minimum value for $\theta=\pi/2$ and vanishes when $\theta=0,\pi$.

The Maxwell electrodynamics in the vicinity of a semi-transparent mirror with magnetoelectric properties have been investigated. We have considered a mirror described by a $\delta$-potential with a Chern-Simons-like coupling. We have found exactly the modification undergone by the photon propagator due to the presence of the mirror.

We have computed the interaction energy between a stationary point-like charge and the mirror and we have shown that the role played by the coupling constant $\mu$ is just to attenuate the image charge by an uniform multiplicative factor. We have also shown that the interaction via image method is recovered in the limit case of a perfect mirror ($\mu\to\infty$).

We have shown that a stationary point-like charge produces an electric field and a magnetic field due to the presence of the mirror. The electric field is the one induced by the charge itself and the field associated to an image charge attenuated by a multiplicative factor. The magnetic field is the one associated with a magnetic monopole. The interesting magnetoelectric properties simulated by the model is an indication that it can be used to describe magnetoelectric surfaces.

We have investigated some physical phenomena which emerge from a setup composed by the plate and a stationary charged quantum particle with spin $1/2$. We have shown that, in addition a Coulombian and an isotropic dipole-dipole-like terms, the energy exhibits a magnetoelectric contribution with the features of a charge-magnetic dipole interaction. This last contribution leads also to a torque acting on the particle with spin.

We hope that the present work could pave the way to use analytical methods to describe magnetoelectric material media.

Acknowledgements:
For financial support, L.H.C. Borges and H.L. Oliveira thank to CAPES (Brazilian agency) and F.A. Barone thanks to CNPq under the grants 311514/2015-4 and 313978/2018-2.