Nonuniform currents and spins of relativistic electron vortices in a magnetic field.
Koen van Kruining
Armen G. Hayrapetyan
J\"org B. G\"otte
We present a relativistic description of electron vortex beams in a homogeneous magnetic field. Including spin from the beginning reveals that spin-polarized electron vortex beams have a complicated azimuthal current structure, containing small rings of counterrotating current between rings of stronger corotating current. Contrary to many other problems in relativistic quantum mechanics, there exists a set of vortex beams with exactly zero spin-orbit mixing in the highly relativistic and nonparaxial regime. The well defined phase structure of these beams is analogous to simpler scalar vortex beams, owing to the protection by the Zeeman effect. For states that do show spin-orbit mixing, the spin polarization across the beam is nonuniform rendering the spin and orbital degrees of freedom inherently inseparable.
LandaubeamsarXiv2.tex
\documentclass[twocolumn,superscriptaddress,prl]{revtex4-1}
\usepackage{amsmath}
\usepackage{natbib}
\usepackage{graphicx}
\usepackage{cancel}
\begin{document}
\title{Nonuniform currents and spins of relativistic electron vortices in a magnetic field.}
\author{Koen van Kruining}
\email{koen@pks.mpg.de}
\author{Armen G. Hayrapetyan}
\affiliation{Max-Planck-Institut f\"ur Physik komplexer Systeme, 01187 Dresden, Germany}
\author{J\"org B. G\"otte}
\affiliation{Nanjing University, Nanjing 210093, China}
\affiliation{University of Glasgow, Glasgow G12 8QQ, United Kingdom}
\begin{abstract}
We present a relativistic description of electron vortex beams in a homogeneous magnetic field. Including spin from the beginning reveals that spin-polarized electron vortex beams have a complicated azimuthal current structure, containing small rings of counterrotating current between rings of stronger corotating current. Contrary to many other problems in relativistic quantum mechanics, there exists a set of vortex beams with exactly zero spin-orbit mixing in the highly relativistic and nonparaxial regime. The well defined phase structure of these beams is analogous to simpler scalar vortex beams, owing to the protection by the Zeeman effect. For states that do show spin-orbit mixing, the spin polarization across the beam is nonuniform rendering the spin and orbital degrees of freedom inherently inseparable.
\end{abstract}
\maketitle
\textit{Introduction}---
The concept of light beams carrying orbital angular momentum along the propagation axis has been widely utilized in modern optics \cite{Allen+:PiO39:1999, MolinaTerriza+:NatPhys3:2007,FrankeArnold+:LPR2:2008}. Based on analogies of the governing wave equations, vortex beams have also been predicted and generated for electrons \cite{BliokhBliokhSNori07,*BliokhDennisNori11, * BliokhNori12b,UchidaTonomura:Nature464:2010,Verbeeck+:Nature467:2010,McMorran+:Science331:2011,VerbeeckSchatschneiderLS-PLS-TT11,vanBoxemVerbeeckPartoens13,GGMFKB14,*GGMFKB15,*GrilloKarimiGFDennisBoyd14,Laroqueteal16,TTMousleyBabikerY17} and neutrons \cite{Clark+:Nature525:2015}, as well as proposed for atoms \cite{HayrapetyanMatulaSurzhykovFritzsche13,LEBabikerAlD14}. This promises the ability to probe and manipulate matter on smaller length scales, but also opens up the possibility to consider the interaction of vortex beams with external fields \cite{BliokhSchattschneiderVNori12,GreenshieldsStampsFranke-Arnold12, *GreenshieldsFranke-ArnoldStamps15,Karlovets12,HayrapetyanMAielloSFritzsche14, GSchattshneiderBliokhNVerbeeck13, *SchattschneiderSS-PLofflerS-TBliokhNori14}, other vortex beams \cite{Bialynicki-Birula04, *Bialynicki-BirulaChmura05,*Bialynicki-BirulaRadozycki06, Ivanov11, *IvanovSerbo11, *Ivanov12, *IvanofSSurzhykovFritzsche16} and atoms \cite{Matulaetal14, *Serboetal15, *ZaytsevSerboShabajev17}. 

In the simplest description these vortex beams are scalar and obey the paraxial Schr\"odinger equation. Going beyond the paraxial approximation reveals a linking between the spin and orbital degrees of freedom arising whenever the beam is tightly confined, complicating the vortex structure \cite{Barnett17, Bialynicki-Birula17}. And whereas light beams as solutions of Maxwell's equation are naturally relativistic, for particles it is important to distinguish between the nonrelativisitic regime based on Schr\"odinger's equation and the relativistic regime covered by the Dirac equation. 

Whether or not a nonrelativistic description suffices depends not only on the energy of the electron beam involved, but also on the importance the spin of the particle in the interaction in question, as spin is naturally included in the Dirac equation \cite{Dirac28, LL4}. For electrons traveling through a magnetic field it is of particular importance to take the spin into account, because it interacts strongly with the field.

%Electron beams carrying orbital angular momentum have been studied both non-relativistically \cite{BliokhBliokhSNori07, BliokhSchattschneiderVNori12} and relativistically \cite{BliokhDennisNori11}, also for electrons interacting with light \cite{HayrapetyanMAielloSFritzsche14}. Whether a non-relativistic description suffices depends not only on the energy of the electron beam involved, but also on whether the spin polarization of the electrons plays an important role, since spin is naturally included in the Dirac equation, but has to be included by hand into the Schr\"odinger equation. For electrons traveling through a magnetic field it is of paricular importance to take the spin into account, because it interacts with the field. 

We analytically solve the Dirac equation for an electron in a homogeneous magnetic field, a problem first considered by Landau \cite{LL3,LL4}. The interaction with the magnetic field confines the beam and gives rise to a set of discrete energy levels (Landau levels) \cite{LL3, BliokhSchattschneiderVNori12}. On top of that the Zeeman effect shifts the energy of the positive and negative spin states relative to each other. The quantized Landau and Zeeman contributions to the energy determine which states  undergo spin-orbit mixing with each other and completely forbid spin-orbit mixing for some of them. The inclusion of spin also leads to a (for some states large) redistribution of the azimuthal current within the beam, revealing a pattern of concentric rings of clockwise and counterclockwise rotating current. Our results and conclusions are not only applicable to electrons propagating in beams, but also for electrons confined in Penning traps. 

Throughout this letter we set $c=\hbar=1$, use the standard representation for the Dirac matrices, slashes to denote contraction with Dirac matrices, the positive z-axis as quantization axis for angular momentum and the metric signature diag($+$ $-$ $-$ $-$). 

\textit{Electron beams in a magnetic field and their spin-orbit structure}---
A magnetic field can be incorporated in the Dirac equation using the gauge covariant momentum operator $P_\mu=p_\mu-eA_\mu=i\partial_\mu-eA_\mu$, with $A_\mu$ the vector potential and $e$ the electron charge. Choosing the magnetic field in the positive z-direction we take the vector potential $A_\mu=\frac 12 B(0,-y,x,0)$, with $B$ the magnitude of the magnetic field. Using cylindrical coordinates and first solving the `squared' Dirac equation $(\cancel P+m)(\cancel P-m)\Psi=0$, we assume a solution of the form $\Psi=e^{i(kz-\mathcal E t\pm l\phi)}\psi(r)u$, with $\mathcal E$ the total energy, $u$ a bispinor and $l$ positive. We rescale the radial coordinate $r$ as $\tilde r=\sqrt{|e|B/2}r$.  At a field strength of one Tesla $\tilde r=1$ corresponds to 36 nanometer. The rescaled equation for the spin and radial parts becomes
\begin{multline}
\frac{B|e|}2\left(\frac 1{\tilde r}\partial_{\tilde r}\tilde r\partial_{\tilde r}-\frac{l^2}{\tilde r^2}\mp 2 l-\tilde r^2-2\Sigma_z\right)\psi(\tilde r)u=\\
-(\mathcal E_L^2+\mathcal E_Z^2)\psi(\tilde r)u,\nonumber
\end{multline}
with $\Sigma_i=\mbox{diag}(\hat \sigma_i,\hat \sigma_i)$, and $\hat \sigma_i$ the Pauli matrices. The interaction energy of the electrons spin magnetic moment is $\mathcal E^2_Z=2\sigma_z B|e|$(=Zeeman energy, $\sigma_z=\pm\frac 12$). $\mathcal E_L$ Is the sum of the electrons orbital kinetic energy and the interaction energy of the orbital magnetic moment (Landau energy). The radial differential equation has the well-known solution \cite{LL3, BliokhSchattschneiderVNori12}
$$
\psi(\tilde r)=\tilde r^le^{-\frac{\tilde r^2}2}L_p^l(\tilde r^2),\quad \mathcal E_L^2=B|e|(2p+l(1\pm 1)+1),
$$
with $L_p^l$ an associated Laguerre polynomial. Here the $\pm$-sign is the sign of the orbital angular momentum. For negative orbital angular momentum $\mathcal E_L^2$ is independent of $l$ because the kinetic and magnetic contributions cancel (FIG.~\ref{niveauschema}). These solutions are nondiffracting Laguerre-Gauss beams, with $p$ the radial quantum number indicating how many  rings surround the central spot or ring.  The solutions of the squared Dirac equation describe superpositions of positive and negative energy states. Applying $\cancel P+m$ to the wave functions projects out the positive energy part  (the full calculation is in the supplementary material). The physical solutions are
\begin{widetext}\begin{eqnarray}
&\Psi=&e^{i(kz-\mathcal E t+l\phi)-\frac{\tilde r^2}2}\left(
\tilde r^l L_p^l(\tilde r^2)\left[\begin{array}{c}m+\mathcal E\\0\\k\\0\end{array}\right]+
\sqrt 2i e^{i\phi}\tilde r^{l+1}L^{l+1}_p(\tilde r^2)\left[\begin{array}{c}0\\0\\0\\ \sqrt{B|e|}\end{array}\right]
\right),\qquad \mbox{spin}>0\mbox{, OAM}\ge 0, \nonumber\\
&\Psi=&e^{i(kz-\mathcal E t+l\phi)-\frac{\tilde r^2}2}\left(
\tilde r^l L_p^l(\tilde r^2)\left[\begin{array}{c}0\\m+\mathcal E\\0\\-k\end{array}\right]-
\sqrt 2(p+l)i e^{-i\phi}\tilde r^{l-1}L^{l-1}_p(\tilde r^2)\left[\begin{array}{c}0\\0\\ \sqrt{B|e|}\\0\end{array}\right]
\right),\qquad \mbox{spin}<0\mbox{, OAM}> 0,\nonumber\\
&\Psi=&e^{i(kz-\mathcal E t-l\phi)-\frac{\tilde r^2}2}\left(
\tilde r^l L_p^l(\tilde r^2)\left[\begin{array}{c}m+\mathcal E\\0\\k\\0\end{array}\right]-
\sqrt 2(p+1)i e^{i\phi}\tilde r^{l-1}L^{l-1}_{p+1}(\tilde r^2)\left[\begin{array}{c}0\\0\\0\\ \sqrt{B|e|}\end{array}\right]
\right),\qquad \mbox{spin}>0\mbox{, OAM}< 0,\nonumber\\
&\Psi=&e^{i(kz-\mathcal E t-l\phi)-\frac{\tilde r^2}2}\left(
\tilde r^l L_p^l(\tilde r^2)\left[\begin{array}{c}0\\m+\mathcal E\\0\\-k\end{array}\right]+
\sqrt 2i e^{-i\phi}\tilde r^{l+1}L^{l+1}_{p-1}(\tilde r^2)\left[\begin{array}{c}0\\0\\ \sqrt{B|e|}\\0\end{array}\right]
\right),\qquad \mbox{spin}<0\mbox{, OAM}\le 0,\nonumber
\end{eqnarray}\end{widetext}
for each of the four combinations of positive and negative spin and orbital angular momentum. Whenever we derive an expression which is different for these four solutions, we put the corresponding expressions in the same order. The second term in the brackets is the spin-orbit mixing term, which appears because orbital angular momentum is not a good quantum number \cite{Dirac28}. 

Of particular interest is the last expression (negative spin and orbital angular momentum). Rewriting $L_{p-1}^{l+1}=-{L'}_p^l$ \cite{RHB}, one sees that the spin-orbit term is zero for $p=0$. The lack of spin-orbit mixing for these states stems from all states having a well defined angular momentum and squared energy. The Zeeman effect shifts the squared energy upwards by $\mathcal E_Z^2=B|e|$ for the states with positive spin and downwards by the same amount for the states with negative spin. The Landau quantization generates a squared energy ladder with level spacing $\Delta\mathcal E_L^2=2B|e|$, twice the Zeeman shift. So the positive spin states are shifted upward one level compared to the negative spin states (FIG.~\ref{niveauschema}) and for the lowest lying states with negative spin there is no positive spin state with equal squared energy they can spin-orbit mix with.
\begin{figure}
\includegraphics[width=\columnwidth]{niveauschemaplat.png}
\caption{The energy levels for a fixed value of $k$ sorted by their total angular momentum. The states with positive spin (red) have one quantum of squared energy more than the states with negative spin (blue). Thus the ground states are not degenerate with any opposite spin states and cannot spin-orbit mix as indicated by the arrows. \label{niveauschema}}
\end{figure}
Without the spin-orbit mixing term, the wave function factorizes into a product state of a constant bispinor and a scalar function. Typically, both for light and electrons, such a simple separation in a spin part and a spatial part is not possible, making these negative angular momentum $p=0$ states quite special. This clean separation of spin and orbital angular momentum also makes the ground states perfectly spin polarized, a condition which otherwise has only been achieved with a more complicated combination of magnetic and electric fields \cite{Karimietal12, Nsofinietal16} high loss of beam intensity \cite{SchattschneiderGrilloAubry17} or extremely high laser intensities \cite{DellwegMuller17}. That they are (for a given $k$) the lowest energy states suggests that there should be a way to selectively populate these `scalar like' unperturbed nonparaxial vortex states. 

\textit{Detailed analysis of the current structure}---
The detailed charge flow within the beam can be computed using the four current $j_\mu=\Psi^\dagger\gamma_0\gamma_\mu\Psi$. Integrating its zeroth component over the entire transverse plane gives a  useful normalization factor. Using $\int_0^\infty x^l L_p^l(x)^2e^{-x}dx=\frac{(l+p)!}{p!}$, the integrated probability density is evaluated to be resp.
\begin{eqnarray}
&\textstyle{\int} j_0=& \pi{\textstyle\frac{(l+p)!}{p!}}\left(m^2+\mathcal E^2+2m\mathcal E+k^2+2B|e|(l+p+1)\right),\nonumber\\
&\textstyle{\int} j_0=& \pi{\textstyle\frac{(l+p)!}{p!}}\left(m^2+\mathcal E^2+2m\mathcal E+k^2+2B|e|(l+p)\right),\nonumber\\
&\textstyle{\int} j_0=& \pi{\textstyle\frac{(l+p)!}{p!}}\left(m^2+\mathcal E^2+2m\mathcal E+k^2+2B|e|(p+1)\right),\nonumber\\
&\textstyle{\int} j_0=& \pi{\textstyle\frac{(l+p)!}{p!}}\left(m^2+\mathcal E^2+2m\mathcal E+k^2+2B|e|p\right).\nonumber
\end{eqnarray}
The last term in the brackets is in each case $\mathcal E_L^2+\mathcal E_Z^2$. Using $\mathcal E=\sqrt{m^2+k^2+\mathcal E_L^2+\mathcal E_Z^2}$, the integrated probability density can in each case be written as ${\int }j_0=2\pi \mathcal E(\mathcal E+m)\frac{(l+p)!}{p!}$. The total current in the z-direction through the transverse plane is
$$
\textstyle{\int }j_z={ \int} j_0\displaystyle \frac k{\mathcal E},
$$
so the electrons have the same speed as particles with mass $\sqrt{m^2+\mathcal E_L^2+\mathcal E_Z^2}$. For the transverse current components one can transform the Dirac matrices into 
\begin{eqnarray}
&\gamma_r=&\cos\phi\gamma_x+\sin\phi\gamma_y,\nonumber\\
&\gamma_\phi=&-\sin\phi\gamma_x+\cos\phi\gamma_y.\nonumber
\end{eqnarray}
The radial component is always zero and the azimuthal component is
\begin{eqnarray}
&j_\phi=&2\sqrt 2\tilde r^{2l+1}e^{-\tilde r^2}L_p^l(\tilde r^2)L_p^{l+1}(\tilde r^2)(\mathcal E+m)\sqrt{B|e|},\nonumber\\
&j_\phi=&2\sqrt 2(p+l)\tilde r^{2l-1}e^{-\tilde r^2}L_p^l(\tilde r^2)L_p^{l-1}(\tilde r^2)(\mathcal E+m)\sqrt{B|e|},\nonumber\\
&j_\phi=-&2\sqrt 2(p+1)\tilde r^{2l-1}e^{-\tilde r^2}L_p^l(\tilde r^2)L_{p+1}^{l-1}(\tilde r^2)(\mathcal E+m)\sqrt{B|e|},\nonumber\\
&j_\phi=-&2\sqrt 2\tilde r^{2l+1}e^{-\tilde r^2}L_p^l(\tilde r^2)L_{p-1}^{l+1}(\tilde r^2)(\mathcal E+m)\sqrt{B|e|},\nonumber
\end{eqnarray}
\begin{figure}[h!]
\includegraphics[width=\columnwidth]{fig2arb.pdf}
\includegraphics[width=\columnwidth]{fig2bcrb.pdf}
\caption{$j_\phi$ For positive spin (red), negative spin (blue) and a spin 0 beam for comparison (dashed) for $l=2,p=3$ (a), $l=-2,p=3$ (b) and $l=-2,p=0$ (c). The spin part of the current gives rise to a series of dips where the current flows in the opposite direction, which are absent when spin is neglected. For negative $l$, the azimuthal current is negative near the center but positive on the outside due to the interaction with the magnetic field. The most striking difference from a spin 0 vortex beam occurs for negative $l$ and $p=0$ where negative spin is a Landau-Zeeman ground state and the azimuthal current is exactly zero everywhere.\label{azc}}
\end{figure}
where we have used the surface element $dz\times d\tilde r$. Rescaling $d\tilde r$ back to $dr$ gives a current proportional to $B|e|$. These expressions are quite different from the azimuthal currents for scalar vortex beams in a magnetic field \cite{BliokhSchattschneiderVNori12}, because the spin contribution is included in them as well \cite{Ohanian85}. As can be seen in FIG.~\ref{azc} the inclusion of the spin current reveals complicated patterns of flows and counterflows, which are absent if spin is neglected. These keep their shape even for magnetic field strengths at which there is no appreciable spin-orbit induced change in the beam profile (FIG.~\ref{fig3})
\begin{figure}
\begin{minipage}{0.5\columnwidth}
\includegraphics[width=\columnwidth]{fig3arb.pdf}
\end{minipage}\begin{minipage}{0.5\columnwidth}
\includegraphics[width=\columnwidth]{fig3brb.pdf}
\end{minipage}
\caption{The regions of negative(=clockwise) azimuthal electron current marked for the $p=3$, $l=2$ state (a) and $p=3$, $l=-2$ (b) for negative (left side, in blue) and positive spin (right side, in red) superposed on the beam profiles for a magnetic field of 1 Tesla. The negative currents occur on the inner side of the dark fringes for positive spin and on the outer side for negative spin. Visible rearrangement of the electron density due to spin-orbit mixing only appears around 1 gigatesla.} \label{fig3}
\end{figure}

\textit{Nonuniform spin}---
As a consequence of spin-orbit mixing, the spin polarization of an electron becomes nonuniform, similar to the nonuniform spin appearing in structured light \cite{BliokhNori12, *BliokhBekshaevNori14, *BekshaevBliokhNori15,NBAB15,Antognozzietal16}, which is used for direction sensitive optical switching \cite{KienHakuta08, *KienRauschenbeutel14a, PhysRevLett.110.213604, *Petersen67,*MSASRauschenbeutel14, *PhysRevX.5.041036, PhysRevLett.113.237203,*PhysRevA.91.042116,*PhysRevA.93.062104, *PhysRevA.93.063830, Kopenhagen15}. Its existence can be inferred decomposing the probability current in a spin and an orbital part \cite{Gordon1928, Ohanian85} and comparing the z-components of the orbital part and the total current, finding that $\mathrm{Re}(\Psi^*P_z/m\Psi)\neq\Psi^*\gamma_0\gamma_z\Psi$. The difference has to be made up for by a spin current $\nabla\times \vec{ \mathcal S}$ caused by a spin component perpendicular to $\hat z$ \footnote{Both Gordon and Ohanian only consider the free case, but because the probability current is gauge invariant, one can see that one must use the gauge covariant momentum operator if fields are present.}. Using 
$\Sigma_r=\cos\phi\Sigma_x+\sin\phi\Sigma_y\mbox{, }
\Sigma_\phi=-\sin\phi\Sigma_x+\cos\phi\Sigma_y,$ 
it can be shown that the radial spin is zero and the azimuthal spin is $\pm\frac 12 k/(\mathcal E+m)j_\phi$, where the sign is given by the sign of the total spin in the z-direction. The ground states' spin polarization is uniform because their spin-orbit mixing is zero. This in contrast to structured light, where the nonuniformity inevitably appears in any finite width beam. 

The difference between a uniformly and a nonuniformly spin polarized state is that for a uniformly polarized state one can always choose a direction along which a spin measurement will certainly give the outcome spin up whereas this is impossible for a nonuniformly polarized state, because spin and spatial degrees of freedom are entangled. For our electron beams this entanglement can be shown by taking their density matrices and tracing out everything except the spin. The remaining mixed spin state is for positive spin
$$
\rho_s=\frac{((m+\mathcal E)^2+k^2)|\uparrow\rangle\langle\uparrow|+(\mathcal E_L^2+\mathcal E_Z^2)|\downarrow\rangle\langle\downarrow|}{2\mathcal E(\mathcal E+m)}
$$
and the same with the spins interchanged for negative spin, showing that one cannot separate the spin and orbital degrees of freedom.

\textit{Gauge covariant angular momentum operator}---
With our choice of gauge, the exact solutions of the Dirac equation are eigenfunctions of the canonical angular momentum operator ($\hat J_z=-i\partial_\phi+\frac 12\Sigma_z$) with eigenvalues resp. $l+\frac 12,\;l-\frac 12,\;-l+\frac 12$ and $-l-\frac 12$. The canonical momentum is not gauge covariant but can be made so by the usual minimal substitution, yielding: $\mathcal {\hat J}_z=-i\partial_\phi-erA_\phi+\frac 12\Sigma_z$. This operator does not have any stationary solution of the Dirac equation or the `squared' Dirac equation as its eigenstate, as can be verified by applying it to any (linear combination of degenerate) basis state. Its expectation value can be computed by adding $\displaystyle\frac{\int\Psi^\dagger |e|rA_\phi\Psi}{\int \Psi^\dagger\Psi}=\displaystyle\frac{\int\Psi^\dagger\tilde r^2\Psi}{\int \Psi^\dagger\Psi}$ to the canonical angular momentum, the result is (suppl. mat.)
\begin{eqnarray}
&\mathcal J_z=&2p+2l+\frac 32+\frac{\mathcal E^2_L+\mathcal E^2_Z}{2\mathcal E(\mathcal E+m)},\nonumber\\
&\mathcal J_z=&2p+2l+\frac 12-\frac{\mathcal E^2_L+\mathcal E^2_Z}{2\mathcal E(\mathcal E+m)},\nonumber \end{eqnarray}\begin{eqnarray}
&\mathcal J_z=&2p+\frac 32+\frac{\mathcal E^2_L+\mathcal E^2_Z}{2\mathcal E(\mathcal E+m)},\nonumber\\
&\mathcal J_z=&2p+\frac 12-\frac{\mathcal E^2_L+\mathcal E^2_Z}{2\mathcal E(\mathcal E+m)}.\nonumber
\end{eqnarray}
If one would neglect the spin-orbit term, one would always get a half-integer expectation value for the gauge covariant angular momentum, although the states, even without spin-orbit term are \emph{not} eigenstates of $\mathcal {\hat J}_z$. This fortuitous coincidence has been overlooked in the literature until now, to the best of our knowledge. The reason that the expectation value of $\mathcal{\hat J}_z$ is not a half integer number is that the orbital contribution changes by two quanta when $l$ or $p$ is changed by one whereas the spin contribution changes by the usual one quantum upon spin flip. Therefore the main term and the spin orbit term have different expectation values for $\mathcal{\hat J}_z$ and one takes the probability weighted average of the both terms. $\mathcal{\hat L}_z+2\mathcal{\hat S}_z$ Does have half-integer expectation values. This last quantity determines the z-component of the magnetic moment, $M_z$, of the electron as can be verified by computing (details in suppl. mat.)
\begin{multline}
M_z=\int \frac e2 rj_\phi=-\frac{\textstyle{\int}j_0}{\mathcal E}\frac{\mathcal E_L^2+\mathcal E_Z^2}{2B}=\\
\frac{e}{2\mathcal E}{\textstyle\int}j_0(2p+l(1\pm1)+1+2\mathcal S)=\frac e{2\mathcal E}\textstyle{\int}j_0(\mathcal L+2\mathcal S)\nonumber
\end{multline}
showing that the gauge covariant operators are the ones determining the magnetic moment.

Apart from not having any stationary eigenfunctions, the gauge covariant angular momentum operators also do not generate a Lie group. These two properties can be proven more generally. Taking the commutator of two of these operators gives (suppl. mat.)
$$
[\mathcal{\hat J}_j,\mathcal{\hat J}_k]=-i\epsilon_{jkl}(\mathcal{\hat J}_l-x_l\mathbf{x\cdot B}),
$$
showing that they violate the closure axiom for Lie algebras if there is any magnetic field present. For the existence of stationary solutions we change notation and write the components of the gauge covariant momenta and `boost' operators as an antisymmetric tensor $\mathcal{\hat J}_{\mu\nu}=x_{[\mu} P_{\nu]}+\frac i2\sigma_{\mu\nu}$. The brackets on the indices indicate antisymmetrization, $T_{[\mu\nu]}=T_{\mu\nu}-T_{\nu\mu}$. With this notation, $\mathcal{\hat J}_{12}=\mathcal{\hat J}_z$. Now the existence of physical states that are eigenstates of $\mathcal{\hat J}_{\mu\nu}$ is only possible if the commutator $[\cancel P-m,\mathcal{\hat J}_{\mu\nu}]$ vanishes. This commutator is (suppl. mat.)
$$
[\cancel P-m,\mathcal{\hat J}_{\mu\nu}]=iex_{[\mu}F_{\nu]\lambda}\gamma^{\lambda},
$$
which vanishes only for an extremely restricted class of possible fields. Taking $\mathcal{\hat J}_{12}$ and writing out the field components explicitly, we have
$$
[\cancel P-m,\mathcal{\hat J}_{12}]=ie\left((xE_y-yE_x)\gamma_0-B_z \mathbf x\cdot\vec\gamma+\gamma_3 \mathbf x\cdot \mathbf B\right).
$$
So the only possible field that would allow for physical eigenstates of $\mathcal{\hat J}_{12}$, is a constant electric field in the z-direction.

\textit{Conclusion}---
We have shown that in a homogeneous magnetic field there exist electron vortex beams without spin-orbit mixing and thus with a very `clean' vortex core. For these beams, spin-orbit mixing remains absent even for strong magnetic fields and relativistic speeds. Including the effect of spin reveals an internal rearrangement of the azimuthal current which is quite substantial if the orbital angular momentum and magnetic field point in opposite directions. For electron vortex beams the current scales linearly with the beam intensity and the spin rearrangement of the azimuthal current can be magnified by using a strong enough electron beam. If an electron vortex beam is wide enough, a suitable test particle can probe these current rearrangements similar to how a small dielectric particle can probe the local Poynting vector of a light beam \cite{HeFrieseHeckenbergRubinsztein-Dunlop95, PhysRevLett.88.053601,NBAB15,Antognozzietal16}. 

\textit{Acknowledgements}--- KvK thanks Valentin Walther for the useful discussion helping him clarify the nature of the nonuniform electron spin. This work was supported by the Engineering and Physical Sciences Research Council of the United Kingdom with grants Nos. EP/I012451/1 and EP/M01326X/1 and the the National Key Research and Development Program of China under contract number 2017YFA0303700.
\bibliography{/home/koen/Documents/Publicaties/Landaubundels/Landaubeams.bib}
%%%% supplementary material %%%%

\onecolumngrid
\pagebreak
\begin{center}
\Large Supplementary material
\end{center}
\subsection*{Obtaining the exact solutions of the Dirac equation in a magnetic field}
The solutions of te squared Dirac equation in a constant magnetic field (symmetric gauge) are, using the rescaled coordinate $\tilde r=\sqrt{\frac{|e|B}2}r$ and taking $l$ positive
$$\Psi=e^{i(kz-\mathcal E t\pm l\phi)}\tilde r^{l}e^{-\frac{\tilde r^2}2}L_p^{l}(\tilde r^2)\left(\left[\begin{array}{c}1\\0\\0\\0\end{array}\right]\vee\left[\begin{array}{c}0\\1\\0\\0\end{array}\right]\right).$$
The exact solutions pf the first order Dirac equation can be obtained by applying $\cancel P+m$ to the solutions of the squared Dirac equation. Using
$$
(\cancel P+m)\left[\begin{array}{c}1\\0\\0\\0\end{array}\right]=\left[\begin{array}{c}
i\partial_t+m\\0\\-i\partial_z\\ -i\partial_x+\partial_y+\frac{|e|B}2(i x-y)
\end{array}\right],\qquad (\cancel P+m)\left[\begin{array}{c}0\\1\\0\\0\end{array}\right]=\left[\begin{array}{c}
0\\ i\partial_t+m\\-i\partial_x-\partial_y-\frac{|e|B}2(ix+y)\\-i\partial_z
\end{array}\right],
$$
the de derivatives with respect to $t$ and $z$ are easy to compute and give resp. $\mathcal E$ and $k$. For the transverse derivatives, one can use the rescaled coordinates $\tilde x=\sqrt{\frac{|e|B}2}x,\; \tilde y=\sqrt{\frac{|e|B}2}y$ to rewrite them as
$$
(\cancel P+m)\left[\begin{array}{c}1\\0\\0\\0\end{array}\right]=\left[\begin{array}{c}
m+\mathcal E\\0\\k\\ \sqrt{\frac{|e|B}2}(-i\partial_{\tilde x}+\partial_{\tilde y}+(i \tilde x-\tilde y))
\end{array}\right],\qquad (\cancel P+m)\left[\begin{array}{c}0\\1\\0\\0\end{array}\right]=\left[\begin{array}{c}
0\\m+\mathcal E\\\sqrt{\frac{|e|B}2}(-i\partial_{\tilde x}-\partial_{\tilde y}-(i\tilde x+\tilde y))\\k
\end{array}\right].
$$
The components $\sqrt{\frac{|e|B}2}(-i\partial_{\tilde x}+\partial_{\tilde y}+(i \tilde x-\tilde y))$ and $\sqrt{\frac{|e|B}2}(-i\partial_{\tilde x}-\partial_{\tilde y}-(i\tilde x+\tilde y))$ give rise to the spin-orbit mixing terms, whose explicit computation is rather lengthy. The following three identities will be of use
\begin{eqnarray}
(\partial_{\tilde x}\pm i\partial_{\tilde y})\tilde r^n&=&n\tilde r^{n-1}e^{\pm i\phi},\nonumber\\
(\partial_{\tilde x}\pm\partial_{\tilde y})\tilde r^{|n|}e^{\pm i|n|\phi}&=&0,\nonumber\\
(\partial_{\tilde x}\pm\partial_{\tilde y})\tilde r^{|n|}e^{\mp i|n|\phi}&=&2|n|\tilde r^{|n|-1}e^{\pm i(|n|-1)\phi}.\nonumber
\end{eqnarray}
The form of the spin-orbit term depends on the signs of the spin and orbital angular momentum. For spin and orbital angular momentum positive, one has
\begin{multline}
 \sqrt{\frac{|e|B}2}(-i\partial_{\tilde x}+\partial_{\tilde y}+(i \tilde x-\tilde y))e^{i(kz-\mathcal E t+ l\phi)}\tilde r^{l}e^{-\frac{\tilde r^2}2}L_p^{l}(\tilde r^2)=\\
i \sqrt{2|e|B}\left(e^{i(kz-\mathcal E t+ (l+1)\phi)}\tilde r^{l+1}e^{-\frac{\tilde r^2}2}\left(L_p^l(\tilde r^2)-{L'}_p^{l}(\tilde r^2)\right)\right).\nonumber
\end{multline}
Using the recurrence relations for Laguerre polynomials ${L'}_p^l(\tilde r^2=-L_{p-1}^{l+1}(\tilde r^2)$ (prime denotes differentiation with respect to $\tilde r^2$) and $L_p^l(\tilde r^2)=L_p^{l+1}(\tilde r^2)-L_{p-1}^{l+1}(\tilde r^2)$, the Laguerre polynomials in the brackets become simply $L_p^{l+1}(\tilde r^2)$, thus the overall solution of the first order Dirac equation in this case becomes
$$
\Psi=e^{i(kz-\mathcal E t+l\phi)-\frac{\tilde r^2}2}\left(
\tilde r^l L_p^l(\tilde r^2)\left[\begin{array}{c}m+\mathcal E\\0\\k\\0\end{array}\right]+
\sqrt 2i e^{i\phi}\tilde r^{l+1}L^{l+1}_p(\tilde r^2)\left[\begin{array}{c}0\\0\\0\\ \sqrt{B|e|}\end{array}\right]
\right). \nonumber\\
$$


Now for positive orbital angular momentum and negative spin, the spin-orbit term is
\begin{multline}
 \sqrt{\frac{|e|B}2}(-i\partial_{\tilde x}-\partial_{\tilde y}-(i \tilde x+\tilde y))e^{i(kz-\mathcal E t+ l\phi)}\tilde r^{l}e^{-\frac{\tilde r^2}2}L_p^{l}(\tilde r^2)=\\
i \sqrt{2|e|B}\left(e^{i(kz-\mathcal E t+ (l-1)\phi)}\tilde r^{l-1}e^{-\frac{\tilde r^2}2}\left(-lL_p^l(\tilde r^2)-\tilde r^2{L'}_p^{l}(\tilde r^2)\right)\right).\nonumber
\end{multline}
Using the recurrence relation
$$
\tilde r^2 {L'}_p^{l}(\tilde r^2)=pL_p^{l}(\tilde r^2)-(p+l)L_{p-1}^{l}(\tilde r^2)
$$
to rewrite the derivatives of the Laguerre polynomial, one gets
$$
-\left(lL_p^l(\tilde r^2)+pL_p^l(\tilde r^2)-(p+l)L_{p-1}^l(\tilde r^2)\right).
$$
With the relation $L_p^{l-1}(\tilde r^2)=L_p^l(\tilde r^2)-L_{p-1}^l(\tilde r^2)$ this expression simplifies to $-(p+l)L_{p}^{l-1}(\tilde r^2)$ and the solution of the first order Dirac equation for positive orbital angular momentum and negative spin becomes
$$
\Psi=e^{i(kz-\mathcal E t+l\phi)-\frac{\tilde r^2}2}\left(
\tilde r^l L_p^l(\tilde r^2)\left[\begin{array}{c}0\\m+\mathcal E\\0\\-k\end{array}\right]-
\sqrt 2(p+l)i e^{-i\phi}\tilde r^{l-1}L^{l-1}_p(\tilde r^2)\left[\begin{array}{c}0\\0\\ \sqrt{B|e|}\\0\end{array}\right]
\right).
$$
For negative orbital angular momentum and positive spin, one has
\begin{multline}
 \sqrt{\frac{|e|B}2}(-i\partial_{\tilde x}+\partial_{\tilde y}+(i \tilde x-\tilde y))e^{i(kz-\mathcal E t-l\phi)}\tilde r^{l}e^{-\frac{\tilde r^2}2}L_p^{l}(\tilde r^2)=\\
i \sqrt{2|e|B}\left(e^{i(kz-\mathcal E t+ (l+1)\phi)}\tilde r^{l-1}e^{-\frac{\tilde r^2}2}\left(\tilde r^2 L_p^l(\tilde r^2)-lL_p^l(\tilde r^2)-\tilde r^2{L'}_p^{l}(\tilde r^2)\right)\right).\nonumber
\end{multline}
Using $L_p^l(\tilde r^2)=-{L'}_{p-1}^{l+1}(\tilde r^2)$, where the prime denotes differentiation with respect to $\tilde r^2$, one can rewrite the first Laguerre polynomial:
$$
-\left(r^2 {L'}^{l-1}_{p+1}(\tilde r^2)+lL_p^l(\tilde r^2)+\tilde r^2{L'}_p^{l}(\tilde r^2)\right).
$$
Using again $\tilde r^2 {L'}_p^{l}(\tilde r^2)=pL_p^{l}(\tilde r^2)-(p+l)L_{p-1}^{l}(\tilde r^2)$ this expression becomes
$$
-\left((p+1)L_{p+1}^{l-1}(\tilde r^2)-(p+l)L_p^{l-1}(\tilde r^2)+lL_p^l(\tilde r^2)+pL_p^l(\tilde r^2)-(p+l)L_{p-1}^l(\tilde r^2)\right).
$$
Because of $L_p^{l-1}(\tilde r^2)=L_p^l(\tilde r^2)-L_{p-1}^l(\tilde r^2)$, everything but the first term cancels and the solution for the first order Dirac equation for negative orbital angular momentum and positive spin becomes
$$
\Psi=e^{i(kz-\mathcal E t-l\phi)-\frac{\tilde r^2}2}\left(
\tilde r^l L_p^l(\tilde r^2)\left[\begin{array}{c}m+\mathcal E\\0\\k\\0\end{array}\right]-
\sqrt 2(p+1)i e^{i\phi}\tilde r^{l-1}L^{l-1}_{p+1}(\tilde r^2)\left[\begin{array}{c}0\\0\\0\\ \sqrt{B|e|}\end{array}\right]
\right)
$$
The case of negative orbital angular momentum and spin is simple, one has
$$
\sqrt{\frac{|e|B}2}(-i\partial_{\tilde x}-\partial_{\tilde y}-(i \tilde x+\tilde y))e^{i(kz-\mathcal E t+ l\phi)}\tilde r^{l}e^{-\frac{\tilde r^2}2}L_p^{l}(\tilde r^2)=
-i \sqrt{2|e|B}\left(e^{i(kz-\mathcal E t+ (l-1)\phi)}\tilde r^{l-1}e^{-\frac{\tilde r^2}2}\tilde r^2{L'}_p^{l}(\tilde r^2)\right),\nonumber
$$
and using again $L_p^l(\tilde r^2)=-{L'}_{p-1}^{l+1}(\tilde r^2)$, the overall solution of the first order Dirac equation becomes
$$
\Psi=e^{i(kz-\mathcal E t-l\phi)-\frac{\tilde r^2}2}\left(
\tilde r^l L_p^l(\tilde r^2)\left[\begin{array}{c}0\\m+\mathcal E\\0\\-k\end{array}\right]+
\sqrt 2i e^{-i\phi}\tilde r^{l+1}L^{l+1}_{p-1}(\tilde r^2)\left[\begin{array}{c}0\\0\\ \sqrt{B|e|}\\0\end{array}\right]
\right)
$$
\subsection*{Explicit forms of the radial and azimuthal gamma and spin matrices}
\begin{eqnarray}
\gamma^r=\cos \phi \gamma^x+\sin \phi \gamma^y=\left[\begin{array}{c c c c}
0 & 0 & 0& e^{-i\phi}\\ 0& 0 & e^{i\phi}& 0\\ 0 & -e^{-i\phi} & 0 & 0\\ -e^{i\phi}& 0 & 0 & 0
\end{array}\right],\nonumber\\
\gamma^\phi=-\sin\phi \gamma^x+\cos\phi\gamma^y=\left[\begin{array}{c c c c}
0 & 0 & 0& -ie^{-i\phi}\\ 0& 0 & ie^{i\phi}& 0\\ 0 & ie^{-i\phi} & 0 & 0\\ -ie^{i\phi}& 0 & 0 & 0
\end{array}\right],\nonumber\\
\Sigma_r=\cos\phi\Sigma_x+\sin\phi\Sigma_y=\left[\begin{array}{c c c c}
0&e^{-i\phi} &0&0\\ e^{i\phi} &0 & 0&0\\ 0&0&0&e^{-i\phi}\\ 0&0&e^{i\phi}&0
\end{array}\right],\nonumber\\
\Sigma_\phi=-\sin\phi\Sigma_x+\cos\phi\Sigma_y=\left[\begin{array}{c c c c}
0& -i e^{-i\phi}&0&0\\ ie^{i\phi}& 0 & 0&0\\ 0&0&0& -ie^{-i\phi}\\ 0&0&ie^{i\phi}&0
\end{array}\right].\nonumber
\end{eqnarray}
\subsection*{Explicit evaluation of $\displaystyle\frac{\int\Psi^\dagger\tilde r^2\Psi}{\int \Psi^\dagger\Psi}$.}
The quantity $\displaystyle\frac{\int\Psi^\dagger\tilde r^2\Psi}{\int \Psi^\dagger\Psi}$ appears in the computation of the gauge covariant angular momentum of the electron vortex states. For the four different combinations of positive and negative orbital angular momentum and spin (order the same as in the main text)  $\int\Psi^\dagger\tilde r^2\Psi$ can be shown to be resp.
\begin{eqnarray}
\int\Psi^\dagger\tilde r^2\Psi=2\pi\int_0^\infty\left((\mathcal E+m)^2+k^2\right)\tilde r^{2l+2}{L^l_p}^2(\tilde r^2) e^{-\tilde r^2}+2 B|e| \tilde r^{2l+4}{L_p^{l+1}}^2(\tilde r^2) e^{-\tilde r^2} d\tilde r,\nonumber\\
\int\Psi^\dagger\tilde r^2\Psi=2\pi\int_0^\infty\left((\mathcal E+m)^2+k^2\right)\tilde r^{2l+2}{L^l_p}^2(\tilde r^2) e^{-\tilde r^2}+2B|e|(p+l)^2\tilde r^{2l}{L_p^{l-1}}^2(\tilde r^2)e^{-\tilde r^2}d\tilde r,\nonumber\\
\int\Psi^\dagger\tilde r^2\Psi=2\pi\int_0^\infty\left((\mathcal E+m)^2+k^2\right)\tilde r^{2l+2}{L^l_p}^2(\tilde r^2) e^{-\tilde r^2}+2B|e|(p+1)^2\tilde r^{2l}{L_{p+1}^{l-1}}^2(\tilde r^2)e^{-\tilde r^2}d\tilde r\nonumber\\
\int\Psi^\dagger\tilde r^2\Psi=2\pi\int_0^\infty\left((\mathcal E+m)^2+k^2\right)\tilde r^{2l+2}{L^l_p}^2(\tilde r^2) e^{-\tilde r^2}+2B|e|\tilde r^{2l+4}{L_{p-1}^{l+1}}^2(\tilde r^2)e^{-\tilde r^2}d\tilde r.\nonumber
\end{eqnarray}
Substituting $x=\tilde r^2$ turns these integrals to
\begin{eqnarray}
\pi\int_0^\infty\left((\mathcal E+m)^2+k^2\right)x^{l+1}{L^l_p}^2(x) e^{-x}+2 B|e| x^{l+2}{L_p^{l+1}}^2(x) e^{-x} dx,\nonumber\\
\pi\int_0^\infty\left((\mathcal E+m)^2+k^2\right)x^{l+1}{L^l_p}^2(x) e^{-x}+2B|e|(p+l)^2x^l{L_p^{l-1}}^2(x)e^{-x}dx,\nonumber\\
\pi\int_0^\infty\left((\mathcal E+m)^2+k^2\right)x^{l+1}{L^l_p}^2(x) e^{-x}+2B|e|(p+1)^2x^l{L_{p+1}^{l-1}}^2(x)e^{-x}dx\nonumber\\
\pi\int_0^\infty\left((\mathcal E+m)^2+k^2\right)x^{l+1}{L^l_p}^2(x) e^{-x}+2B|e|x^{l+2}{L_{p-1}^{l+1}}^2(x)e^{-x}dx.\nonumber
\end{eqnarray}
For Laguerre polynomials, we have the orthogonality relation $\int_0^\infty L_p^l(x)L_{p'}^l(x)e^{-x}dx=\frac{(l+p)!}{p!}\delta_{pp'}$. Using this relation and $L_p^{l-1}(x)=L_p^l(x)-L_{p-1}^l(x)$, we obtain the following integral identity
$$
\int_0^\infty x^{l+1}{L_p^l}^2(x)e^{-x}dx=\int_0^\infty x^{l+1}\left(L_p^{l+1}(x)-L_{p-1}^{l+1}(x)\right)^2e^{-x}dx=\frac{(l+p+1)!}{p!}+\frac{l+p!}{(p-1)!}=\frac{l+p!}{p!}(2p+l+1),
$$
which can be used to evaluate all the integals and obtain
\begin{eqnarray}
\int\Psi^\dagger\tilde r^2\Psi=\pi\frac{(l+p)!}{p!}\left(\left((\mathcal E+m)^2+k^2\right)(2p+l+1)+2B|e|(p+l+1)(2p+l+2)\right),\nonumber\\
\int\Psi^\dagger\tilde r^2\Psi=\pi\frac{(l+p)!}{p!}\left(\left((\mathcal E+m)^2+k^2\right)(2p+l+1)+2B|e|(p+l)(2p+l)\right),\nonumber\\
\int\Psi^\dagger\tilde r^2\Psi=\pi\frac{(l+p)!}{p!}\left(\left((\mathcal E+m)^2+k^2\right)(2p+l+1)+2B|e|(p+1)(2p+l+2)\right),\nonumber\\
\int\Psi^\dagger\tilde r^2\Psi=\pi\frac{(l+p)!}{p!}\left(\left((\mathcal E+m)^2+k^2\right)(2p+l+1)+2B|e|p(2p+l)\right).\nonumber
\end{eqnarray}
By noting that $\mathcal E_L^2+\mathcal E_Z^2$ is $2B|e|$ times resp. $p+l+1$, $p+l$, $p+1$ and $p$ and rearranging, one gets
\begin{eqnarray}
\int\Psi^\dagger\tilde r^2\Psi=\pi\frac{(l+p)!}{p!}\left(2\mathcal E(\mathcal E+m)(2p+l+1)+\mathcal E_L^2+\mathcal E_Z^2\right),\nonumber\\
\int\Psi^\dagger\tilde r^2\Psi=\pi\frac{(l+p)!}{p!}\left(2\mathcal E(\mathcal E+m)(2p+l+1)-\mathcal E_L^2-\mathcal E_Z^2\right),\nonumber\\
\int\Psi^\dagger\tilde r^2\Psi=\pi\frac{(l+p)!}{p!}\left(2\mathcal E(\mathcal E+m)(2p+l+1)+\mathcal E_L^2+\mathcal E_Z^2\right),\nonumber\\
\int\Psi^\dagger\tilde r^2\Psi=\pi\frac{(l+p)!}{p!}\left(2\mathcal E(\mathcal E+m)(2p+l+1)-\mathcal E_L^2-\mathcal E_Z^2\right).\nonumber
\end{eqnarray}
Dividing by $\int\Psi^\dagger\Psi={\textstyle \int j_0}=2\pi\mathcal E(\mathcal E+m)\frac{(p+l)!}{p!}$ and adding the canonical angular momentum, resp. $l+\frac 12,\;l-\frac 12,\;-l+\frac 12$ and $-l-\frac 12$, one gets for the gauge covariant angular momentum
\begin{eqnarray}
&\mathcal J_z=&2p+2l+\frac 32+\frac{\mathcal E^2_L+\mathcal E^2_Z}{2\mathcal E(\mathcal E+m)},\nonumber\\
&\mathcal J_z=&2p+2l+\frac 12-\frac{\mathcal E^2_L+\mathcal E^2_Z}{2\mathcal E(\mathcal E+m)},\nonumber\\
&\mathcal J_z=&2p+\frac 32+\frac{\mathcal E^2_L+\mathcal E^2_Z}{2\mathcal E(\mathcal E+m)},\nonumber\\
&\mathcal J_z=&2p+\frac 12-\frac{\mathcal E^2_L+\mathcal E^2_Z}{2\mathcal E(\mathcal E+m)}.\nonumber
\end{eqnarray}
\section*{Computation of the magnetic moment}
Using the explicit form of the azimuthal Dirac matrix, it is easy to see that only the crossterms between the main and spin-orbit parts contribute to the azimuthal current and these can be computed to be
\begin{eqnarray}
&j_\phi=&2\sqrt 2\tilde r^{2l+1}e^{-\tilde r^2}L_p^l(\tilde r^2)L_p^{l+1}(\tilde r^2)(\mathcal E+m)\sqrt{B|e|},\nonumber\\
&j_\phi=&2\sqrt 2(p+l)\tilde r^{2l-1}e^{-\tilde r^2}L_p^l(\tilde r^2)L_p^{l-1}(\tilde r^2)(\mathcal E+m)\sqrt{B|e|},\nonumber\\
&j_\phi=&-2\sqrt 2(p+1)\tilde r^{2l-1}e^{-\tilde r^2}L_p^l(\tilde r^2)L_{p+1}^{l-1}(\tilde r^2)(\mathcal E+m)\sqrt{B|e|},\nonumber\\
&j_\phi=&-2\sqrt 2\tilde r^{2l+1}e^{-\tilde r^2}L_p^l(\tilde r^2)L_{p-1}^{l+1}(\tilde r^2)(\mathcal E+m)\sqrt{B|e|}.\nonumber
\end{eqnarray}
Now $M_z=\int\frac e2 rj_\phi \tilde r d\phi d\tilde r=\sqrt{\frac 2{B|e|}}\int\frac e2 \tilde rj_\phi r d\phi d\tilde r$. Substituting the explicit currents into this integral, using $x=\tilde r^2$ and performing the angular integration yields
\begin{eqnarray}
M_z=-2\pi(\mathcal E+m)|e|\int_0^\infty x^{l+1}L_p^l(x)L_p^{l+1}(x) e^{-x}dx,\nonumber\\
M_z=-2\pi(\mathcal E+m)|e|(p+l)\int_0^\infty x^lL_p^l(x)L_p^{l-1}(x) e^{-x}dx,\nonumber\\
M_z=2\pi(\mathcal E+m)|e|(p+1)\int_0^\infty x^lL_p^l(x)L_{p+1}^{l-1}(x) e^{-x}dx,\nonumber\\
M_z=2\pi(\mathcal E+m)|e|\int_0^\infty x^{l+1}L_p^l(x)L_{p-1}^{l+1}(x) e^{-x}dx.\nonumber
\end{eqnarray}
Now one can again use $L_p^l(x)=L_p^{l+1}(x)-L_{p-1}^{l+1}(x)$ and the orthogonality relation of associated Laguerre polynomials to evaluate these integrals
\begin{eqnarray}
M_z=-2\pi(\mathcal E+m)|e|\frac{(l+p+1)!}{p!}=-\frac{\textstyle{\int}j_0}{\mathcal E}\frac{\mathcal E_L^2+\mathcal E_Z^2}{2B},\nonumber\\
M_z=-2\pi(\mathcal E+m)|e|(p+l)\frac{(p+l)!}{p!}=-\frac{\textstyle{\int}j_0}{\mathcal E}\frac{\mathcal E_L^2+\mathcal E_Z^2}{2B},\nonumber\\
M_z=-2\pi(\mathcal E+m)|e|(p+1)\frac{(l+p)!}{p!}=-\frac{\textstyle{\int}j_0}{\mathcal E}\frac{\mathcal E_L^2+\mathcal E_Z^2}{2B},\nonumber\\
M_z=-2\pi(\mathcal E+m)|e|\frac{(l+p)!}{(p-1)!}=-\frac{\textstyle{\int}j_0}{\mathcal E}\frac{\mathcal E_L^2+\mathcal E_Z^2}{2B}.\nonumber
\end{eqnarray}
\section*{Commutator identities for the gauge covariant angular momentum operators}
For this section we write the angular omentum operators in antisymmetric tensor form. The gauge covariant angular momentum can be split in a spin and an orbital part like
$$
\mathcal J_{\mu\nu}=\mathcal L_{\mu\nu}+\frac i2\sigma_{\mu\nu},\qquad\mbox{with } \sigma_{\mu\nu }=\frac 12[\gamma_\mu,\gamma_\nu]\mbox{ and }\mathcal L_{\mu\nu}=x_{[\mu}P_{\nu]}\equiv x_\mu P_\nu-x_\nu P_\mu.
$$
Because $\mathcal L_{\mu\nu}$ contains no Dirac matrices, one obviously has $[\mathcal L_{\mu\nu},\sigma_{\rho\sigma}]=0$, so $[\mathcal J_{\mu\nu},\mathcal J_{\rho\sigma}]=[\mathcal L_{\mu\nu},\mathcal L_{\rho\sigma}]-\frac 14[\sigma_{\mu\nu},\sigma_{\rho\sigma}]$. Writing out the commutator for the $\sigma$-tensor gives
$$
[\sigma_{\mu\nu},\sigma_{\rho\sigma}]=\frac 14\left((\gamma_\mu\gamma_\nu-\gamma_\nu\gamma_\mu)(\gamma_\rho\gamma_\sigma-\gamma_\sigma\gamma_\rho)-(\gamma_\rho\gamma_\sigma-\gamma_\sigma\gamma_\rho)(\gamma_\mu\gamma_\nu-\gamma_\nu\gamma_\mu)\right).
$$
One can check that if all four indices are different, this commutator is zero. If two indices are the same one can eliminate the identical Dirac matrices and obtain after some algebra
$$
[\sigma_{\mu\nu},\sigma_{\rho\sigma}]=2(-\eta_{\mu\rho}\sigma_{\nu\sigma}+\eta_{\mu\sigma}\sigma_{\nu\rho}+\eta_{\nu\rho}\sigma_{\mu\sigma}-\eta_{\nu\sigma}\sigma_{\mu\rho} ).
$$
For the orbital part, we need the commutation relations of the gauge covariant momentum $P_\mu=i\partial_\mu-eA_\mu$. It is easy to check that
$$
[ix_{[\mu}\partial_{\nu]},ix_{[\rho}\partial_{\sigma]}]=-\eta_{\mu\rho}x_{[\nu}\partial_{\sigma]}+\eta_{\mu\sigma}x_{[\nu}\partial_{\rho]}
+\eta_{\nu\rho}x_{[\mu}\partial_{\sigma]}-\eta_{\nu\sigma}x_{[\mu}\partial_{\rho]}.
$$
To get the commutators for the gauge coveriant orbital angular momenta, we need to add $[ix_{[\mu}\partial_{\nu]},-ex_{[\rho}A_{\sigma]}]+[-ex_{[\mu}A_{\nu]},ix_{[\rho}\partial_{\sigma]}]=[ix_{[\mu}\partial_{\nu]},-ex_{[\rho}A_{\sigma]}]-[ix_{[\rho}\partial_{\sigma]},-ex_{[\mu}A_{\nu]}]$ (the vector potentias commute with each other). Using that both terms are the same up to the index swap $\mu\leftrightarrow\rho$, $\nu\leftrightarrow\sigma$ and using $\partial_\mu A_\rho-(\mu\leftrightarrow\rho)=F_{\mu\rho}$ these terms can be evaluated to be
\begin{multline}
[ix_{[\mu}\partial_{\nu]},-ex_{[\rho}A_{\sigma]}]+[-ex_{[\mu}A_{\nu]},ix_{[\rho}\partial_{\sigma]}]=-ie(\eta_{\mu\rho}x_{[\nu}A_{\sigma]}-\eta_{\nu\rho}x_{[\mu}A_{\sigma]}-\eta_{\mu\sigma} x_{[\nu}A_{\rho]}+\eta_{\nu\sigma}x_{[\mu}A_{\rho]})\\
-ie(x_\mu x_\rho F_{\nu\sigma}-x_\mu x_\sigma F_{\nu\rho}-x_\nu x_\rho F_{\mu\sigma}+x_\nu x_\sigma F_{\mu\rho}).\nonumber
\end{multline}
Using $-x_{[\mu}\partial_{\rho]}-ie x_{[\mu}A_{\rho]}=ix_{[\mu}P_{\rho]}$, and putting things together gives
\begin{eqnarray}
-i[\mathcal L_{\mu\nu},\mathcal L_{\rho\sigma}]=\eta_{\mu\rho}\mathcal L_{\nu\sigma}-\eta _{\mu\sigma}\mathcal L_{\nu\rho}-\eta_{\nu\rho}\mathcal L_{\mu\sigma}+\eta_{\nu\sigma}\mathcal L_{\mu\rho}+e(x_\mu x_\rho F_{\nu\sigma}-x_\mu x_\sigma F_{\nu\rho}-x_\nu x_\rho F_{\mu\sigma}+x_\nu x_\sigma F_{\mu\rho}),\nonumber\\
-i[\mathcal J_{\mu\nu},\mathcal J_{\rho\sigma}]=\eta_{\mu\rho}\mathcal J_{\nu\sigma}-\eta _{\mu\sigma}\mathcal J_{\nu\rho}-\eta_{\nu\rho}\mathcal J_{\mu\sigma}+\eta_{\nu\sigma}\mathcal J_{\mu\rho}+e(x_\mu x_\rho F_{\nu\sigma}-x_\mu x_\sigma F_{\nu\rho}-x_\nu x_\rho F_{\mu\sigma}+x_\nu x_\sigma F_{\mu\rho}).\nonumber
\end{eqnarray}
Then using $\mathcal J_x=\mathcal J_{23}$, $\mathcal J_y=\mathcal J_{31}$ and $\mathcal J_z=\mathcal J_{12}$, one gets 
$$
[\mathcal J_j,\mathcal J_k]=-i\epsilon_{jkl}(\mathcal J_l-x_l\mathbf{x\cdot B}),
$$
For the commutator $[\cancel P-m,J_{\mu\nu}]$, one can first note that $m$ commutes with any operator. Again using $\mathcal J_{\mu\nu}=\mathcal L_{\mu\nu}+\frac i2\sigma_{\mu\nu}$, one can compute the commutators of the spin and orbital parts seperately using $[P_\mu,P_\nu]=-ie F_{\mu\nu}$:
\begin{eqnarray}
[\cancel P,\sigma_{\mu\nu}]=\frac 12P^\lambda\left[ \gamma_\lambda,[\gamma_\mu,\gamma_\nu]\right]=2P^\lambda\eta_{\lambda[\mu}\gamma_{\nu]}=2P_{[\mu}\gamma_{\nu]}=-2 \gamma_{[\mu}P_{\nu]},\nonumber\\
\left[\gamma^\lambda P_\lambda,x_{[\mu}P_{\nu]}\right]=i\gamma^{\lambda}\eta_{\lambda[\mu}P_{\nu]}+iex_{[\mu}F_{\nu]\lambda}\gamma^{\lambda}=i\gamma_{[\mu}P_{\nu]}+iex_{[\mu}F_{\nu]\lambda}\gamma^{\lambda},\nonumber\\
\left[\gamma^\lambda P_\lambda,x_{[\mu}P_{\nu]}+\frac i2\sigma_{\mu\nu}\right]=[\cancel P,J_{\mu\nu}]=iex_{[\mu}F_{\nu]\lambda}\gamma^{\lambda}.\nonumber
\end{eqnarray}
\end{document}
LandaubeamsarXiv2.bbl
%merlin.mbs apsrev4-1.bst 2010-07-25 4.21a (PWD, AO, DPC) hacked
%Control: key (0)
%Control: author (8) initials jnrlst
%Control: editor formatted (1) identically to author
%Control: production of article title (-1) disabled
%Control: page (0) single
%Control: year (1) truncated
%Control: production of eprint (0) enabled
\begin{thebibliography}{67}%
\makeatletter
\providecommand \@ifxundefined [1]{%
 \@ifx{#1\undefined}
}%
\providecommand \@ifnum [1]{%
 \ifnum #1\expandafter \@firstoftwo
 \else \expandafter \@secondoftwo
 \fi
}%
\providecommand \@ifx [1]{%
 \ifx #1\expandafter \@firstoftwo
 \else \expandafter \@secondoftwo
 \fi
}%
\providecommand \natexlab [1]{#1}%
\providecommand \enquote  [1]{``#1''}%
\providecommand \bibnamefont  [1]{#1}%
\providecommand \bibfnamefont [1]{#1}%
\providecommand \citenamefont [1]{#1}%
\providecommand \href@noop [0]{\@secondoftwo}%
\providecommand \href [0]{\begingroup \@sanitize@url \@href}%
\providecommand \@href[1]{\@@startlink{#1}\@@href}%
\providecommand \@@href[1]{\endgroup#1\@@endlink}%
\providecommand \@sanitize@url [0]{\catcode `\\12\catcode `\$12\catcode
  `\&12\catcode `\#12\catcode `\^12\catcode `\_12\catcode `\%12\relax}%
\providecommand \@@startlink[1]{}%
\providecommand \@@endlink[0]{}%
\providecommand \url  [0]{\begingroup\@sanitize@url \@url }%
\providecommand \@url [1]{\endgroup\@href {#1}{\urlprefix }}%
\providecommand \urlprefix  [0]{URL }%
\providecommand \Eprint [0]{\href }%
\providecommand \doibase [0]{http://dx.doi.org/}%
\providecommand \selectlanguage [0]{\@gobble}%
\providecommand \bibinfo  [0]{\@secondoftwo}%
\providecommand \bibfield  [0]{\@secondoftwo}%
\providecommand \translation [1]{[#1]}%
\providecommand \BibitemOpen [0]{}%
\providecommand \bibitemStop [0]{}%
\providecommand \bibitemNoStop [0]{.\EOS\space}%
\providecommand \EOS [0]{\spacefactor3000\relax}%
\providecommand \BibitemShut  [1]{\csname bibitem#1\endcsname}%
\let\auto@bib@innerbib\@empty
%</preamble>
\bibitem [{\citenamefont {Allen}\ \emph {et~al.}(1999)\citenamefont {Allen},
  \citenamefont {Padgett},\ and\ \citenamefont {Babiker}}]{Allen+:PiO39:1999}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {L.}~\bibnamefont
  {Allen}}, \bibinfo {author} {\bibfnamefont {M.~J.}\ \bibnamefont {Padgett}},
  \ and\ \bibinfo {author} {\bibfnamefont {M.}~\bibnamefont {Babiker}},\
  }\href@noop {} {\bibfield  {journal} {\bibinfo  {journal} {Prog. Opt.}\
  }\textbf {\bibinfo {volume} {XXXIX}},\ \bibinfo {pages} {291} (\bibinfo
  {year} {1999})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Molina-Terriza}\ \emph {et~al.}(2007)\citenamefont
  {Molina-Terriza}, \citenamefont {Torres},\ and\ \citenamefont
  {Torner}}]{MolinaTerriza+:NatPhys3:2007}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {G.}~\bibnamefont
  {Molina-Terriza}}, \bibinfo {author} {\bibfnamefont {J.~P.}\ \bibnamefont
  {Torres}}, \ and\ \bibinfo {author} {\bibfnamefont {L.}~\bibnamefont
  {Torner}},\ }\href {\doibase 10.1038/nphys607} {\bibfield  {journal}
  {\bibinfo  {journal} {Nat. Phys.}\ }\textbf {\bibinfo {volume} {3}},\
  \bibinfo {pages} {305} (\bibinfo {year} {2007})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Franke-Arnold}\ \emph {et~al.}(2008)\citenamefont
  {Franke-Arnold}, \citenamefont {Allen},\ and\ \citenamefont
  {Padgett}}]{FrankeArnold+:LPR2:2008}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {S.}~\bibnamefont
  {Franke-Arnold}}, \bibinfo {author} {\bibfnamefont {L.}~\bibnamefont
  {Allen}}, \ and\ \bibinfo {author} {\bibfnamefont {M.~J.}\ \bibnamefont
  {Padgett}},\ }\href {\doibase 10.1002/lpor.200810007} {\bibfield  {journal}
  {\bibinfo  {journal} {Laser \& Phot. Rev.}\ }\textbf {\bibinfo {volume}
  {2}},\ \bibinfo {pages} {299} (\bibinfo {year} {2008})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Bliokh}\ \emph {et~al.}(2007)\citenamefont {Bliokh},
  \citenamefont {Bliokh}, \citenamefont {Savel'ev},\ and\ \citenamefont
  {Nori}}]{BliokhBliokhSNori07}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {K.~Y.}\ \bibnamefont
  {Bliokh}}, \bibinfo {author} {\bibfnamefont {Y.~P.}\ \bibnamefont {Bliokh}},
  \bibinfo {author} {\bibfnamefont {S.}~\bibnamefont {Savel'ev}}, \ and\
  \bibinfo {author} {\bibfnamefont {F.}~\bibnamefont {Nori}},\ }\href {\doibase
  10.1103/PhysRevLett.99.190404} {\bibfield  {journal} {\bibinfo  {journal}
  {Phys. Rev. Lett.}\ }\textbf {\bibinfo {volume} {99}},\ \bibinfo {pages}
  {190404} (\bibinfo {year} {2007})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Bliokh}\ \emph {et~al.}(2011)\citenamefont {Bliokh},
  \citenamefont {Dennis},\ and\ \citenamefont {Nori}}]{BliokhDennisNori11}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {K.~Y.}\ \bibnamefont
  {Bliokh}}, \bibinfo {author} {\bibfnamefont {M.~R.}\ \bibnamefont {Dennis}},
  \ and\ \bibinfo {author} {\bibfnamefont {F.}~\bibnamefont {Nori}},\ }\href
  {\doibase 10.1103/PhysRevLett.107.174802} {\bibfield  {journal} {\bibinfo
  {journal} {Phys. Rev Lett}\ }\textbf {\bibinfo {volume} {107}},\ \bibinfo
  {pages} {174802} (\bibinfo {year} {2011})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Bliokh}\ and\ \citenamefont
  {Nori}(2012{\natexlab{a}})}]{BliokhNori12b}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {K.~Y.}\ \bibnamefont
  {Bliokh}}\ and\ \bibinfo {author} {\bibfnamefont {F.}~\bibnamefont {Nori}},\
  }\href {\doibase 10.1103/PhysRevA.86.033824} {\bibfield  {journal} {\bibinfo
  {journal} {Phys. Rev. A}\ }\textbf {\bibinfo {volume} {86}},\ \bibinfo
  {pages} {033824} (\bibinfo {year} {2012}{\natexlab{a}})}\BibitemShut
  {NoStop}%
\bibitem [{\citenamefont {Uchida}\ and\ \citenamefont
  {Tonomura}(2010)}]{UchidaTonomura:Nature464:2010}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {M.}~\bibnamefont
  {Uchida}}\ and\ \bibinfo {author} {\bibfnamefont {A.}~\bibnamefont
  {Tonomura}},\ }\href {\doibase 10.1038/nature08904} {\bibfield  {journal}
  {\bibinfo  {journal} {Nature}\ }\textbf {\bibinfo {volume} {464}},\ \bibinfo
  {pages} {737} (\bibinfo {year} {2010})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Verbeeck}\ \emph {et~al.}(2010)\citenamefont
  {Verbeeck}, \citenamefont {Tian},\ and\ \citenamefont
  {Schattschneider}}]{Verbeeck+:Nature467:2010}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {J.}~\bibnamefont
  {Verbeeck}}, \bibinfo {author} {\bibfnamefont {H.}~\bibnamefont {Tian}}, \
  and\ \bibinfo {author} {\bibfnamefont {P.}~\bibnamefont {Schattschneider}},\
  }\href {\doibase 10.1038/nature09366} {\bibfield  {journal} {\bibinfo
  {journal} {Nature}\ }\textbf {\bibinfo {volume} {467}},\ \bibinfo {pages}
  {301} (\bibinfo {year} {2010})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {McMorran}\ \emph {et~al.}(2011)\citenamefont
  {McMorran}, \citenamefont {Agrawal}, \citenamefont {Anderson}, \citenamefont
  {Herzing}, \citenamefont {Lezec}, \citenamefont {McClelland},\ and\
  \citenamefont {Unguris}}]{McMorran+:Science331:2011}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {B.~J.}\ \bibnamefont
  {McMorran}}, \bibinfo {author} {\bibfnamefont {A.}~\bibnamefont {Agrawal}},
  \bibinfo {author} {\bibfnamefont {I.~M.}\ \bibnamefont {Anderson}}, \bibinfo
  {author} {\bibfnamefont {A.~A.}\ \bibnamefont {Herzing}}, \bibinfo {author}
  {\bibfnamefont {H.~J.}\ \bibnamefont {Lezec}}, \bibinfo {author}
  {\bibfnamefont {J.~J.}\ \bibnamefont {McClelland}}, \ and\ \bibinfo {author}
  {\bibfnamefont {J.}~\bibnamefont {Unguris}},\ }\href {\doibase
  10.1126/science.1198804} {\bibfield  {journal} {\bibinfo  {journal}
  {Science}\ }\textbf {\bibinfo {volume} {331}},\ \bibinfo {pages} {192}
  (\bibinfo {year} {2011})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Verbeeck}\ \emph {et~al.}(2011)\citenamefont
  {Verbeeck}, \citenamefont {Schattschneider}, \citenamefont {Lazar},
  \citenamefont {St\"oger-Pollach}, \citenamefont {L\"offler}, \citenamefont
  {Steiger-Thirsfeld},\ and\ \citenamefont
  {Tendeloo}}]{VerbeeckSchatschneiderLS-PLS-TT11}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {J.}~\bibnamefont
  {Verbeeck}}, \bibinfo {author} {\bibfnamefont {P.}~\bibnamefont
  {Schattschneider}}, \bibinfo {author} {\bibfnamefont {S.}~\bibnamefont
  {Lazar}}, \bibinfo {author} {\bibfnamefont {M.}~\bibnamefont
  {St\"oger-Pollach}}, \bibinfo {author} {\bibfnamefont {S.}~\bibnamefont
  {L\"offler}}, \bibinfo {author} {\bibfnamefont {A.}~\bibnamefont
  {Steiger-Thirsfeld}}, \ and\ \bibinfo {author} {\bibfnamefont {G.~V.}\
  \bibnamefont {Tendeloo}},\ }\href {\doibase 10.1063/1.3662012} {\bibfield
  {journal} {\bibinfo  {journal} {Appl. Phys. Lett.}\ }\textbf {\bibinfo
  {volume} {99}},\ \bibinfo {pages} {203109} (\bibinfo {year}
  {2011})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {van Boxem}\ \emph {et~al.}(2013)\citenamefont {van
  Boxem}, \citenamefont {Verbeeck},\ and\ \citenamefont
  {Partoens}}]{vanBoxemVerbeeckPartoens13}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {R.}~\bibnamefont {van
  Boxem}}, \bibinfo {author} {\bibfnamefont {J.}~\bibnamefont {Verbeeck}}, \
  and\ \bibinfo {author} {\bibfnamefont {B.}~\bibnamefont {Partoens}},\ }\href
  {http://stacks.iop.org/0295-5075/102/i=4/a=40010} {\bibfield  {journal}
  {\bibinfo  {journal} {EPL (Europhysics Letters)}\ }\textbf {\bibinfo {volume}
  {102}},\ \bibinfo {pages} {40010} (\bibinfo {year} {2013})}\BibitemShut
  {NoStop}%
\bibitem [{\citenamefont {Grillo}\ \emph
  {et~al.}(2014{\natexlab{a}})\citenamefont {Grillo}, \citenamefont {Gazzadi},
  \citenamefont {Karimi}, \citenamefont {Mafakheri}, \citenamefont {Boyd},\
  and\ \citenamefont {Frabboni}}]{GGMFKB14}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {V.}~\bibnamefont
  {Grillo}}, \bibinfo {author} {\bibfnamefont {G.~C.}\ \bibnamefont {Gazzadi}},
  \bibinfo {author} {\bibfnamefont {E.}~\bibnamefont {Karimi}}, \bibinfo
  {author} {\bibfnamefont {E.}~\bibnamefont {Mafakheri}}, \bibinfo {author}
  {\bibfnamefont {R.~W.}\ \bibnamefont {Boyd}}, \ and\ \bibinfo {author}
  {\bibfnamefont {S.}~\bibnamefont {Frabboni}},\ }\href {\doibase
  10.1063/1.4863564} {\bibfield  {journal} {\bibinfo  {journal} {Appl. Phys.
  Lett.}\ }\textbf {\bibinfo {volume} {104}},\ \bibinfo {pages} {043109}
  (\bibinfo {year} {2014}{\natexlab{a}})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Grillo}\ \emph {et~al.}(2015)\citenamefont {Grillo},
  \citenamefont {Gazzadi}, \citenamefont {Mafakheri}, \citenamefont {Frabboni},
  \citenamefont {Karimi},\ and\ \citenamefont {Boyd}}]{GGMFKB15}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {V.}~\bibnamefont
  {Grillo}}, \bibinfo {author} {\bibfnamefont {G.~C.}\ \bibnamefont {Gazzadi}},
  \bibinfo {author} {\bibfnamefont {E.}~\bibnamefont {Mafakheri}}, \bibinfo
  {author} {\bibfnamefont {S.}~\bibnamefont {Frabboni}}, \bibinfo {author}
  {\bibfnamefont {E.}~\bibnamefont {Karimi}}, \ and\ \bibinfo {author}
  {\bibfnamefont {R.~W.}\ \bibnamefont {Boyd}},\ }\href {\doibase
  10.1103/PhysRevLett.114.034801} {\bibfield  {journal} {\bibinfo  {journal}
  {Phys Rev. Lett.}\ }\textbf {\bibinfo {volume} {114}},\ \bibinfo {pages}
  {034801} (\bibinfo {year} {2015})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Grillo}\ \emph
  {et~al.}(2014{\natexlab{b}})\citenamefont {Grillo}, \citenamefont {Karimi},
  \citenamefont {Gazzadi}, \citenamefont {Frabboni}, \citenamefont {Dennis},\
  and\ \citenamefont {Boyd}}]{GrilloKarimiGFDennisBoyd14}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {V.}~\bibnamefont
  {Grillo}}, \bibinfo {author} {\bibfnamefont {E.}~\bibnamefont {Karimi}},
  \bibinfo {author} {\bibfnamefont {G.~C.}\ \bibnamefont {Gazzadi}}, \bibinfo
  {author} {\bibfnamefont {S.}~\bibnamefont {Frabboni}}, \bibinfo {author}
  {\bibfnamefont {M.~R.}\ \bibnamefont {Dennis}}, \ and\ \bibinfo {author}
  {\bibfnamefont {R.~W.}\ \bibnamefont {Boyd}},\ }\href {\doibase
  10.1103/PhysRevX.4.011013} {\bibfield  {journal} {\bibinfo  {journal} {Phys.
  Rev. X}\ }\textbf {\bibinfo {volume} {4}},\ \bibinfo {pages} {011013}
  (\bibinfo {year} {2014}{\natexlab{b}})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Larocque}\ \emph {et~al.}(2016)\citenamefont
  {Larocque}, \citenamefont {Bouchard}, \citenamefont {Grillo}, \citenamefont
  {Sit}, \citenamefont {Frabboni}, \citenamefont {Dunin-Borkowski},
  \citenamefont {Padgett}, \citenamefont {Boyd},\ and\ \citenamefont
  {Karimi}}]{Laroqueteal16}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {H.}~\bibnamefont
  {Larocque}}, \bibinfo {author} {\bibfnamefont {F.}~\bibnamefont {Bouchard}},
  \bibinfo {author} {\bibfnamefont {V.}~\bibnamefont {Grillo}}, \bibinfo
  {author} {\bibfnamefont {A.}~\bibnamefont {Sit}}, \bibinfo {author}
  {\bibfnamefont {S.}~\bibnamefont {Frabboni}}, \bibinfo {author}
  {\bibfnamefont {R.~E.}\ \bibnamefont {Dunin-Borkowski}}, \bibinfo {author}
  {\bibfnamefont {M.~J.}\ \bibnamefont {Padgett}}, \bibinfo {author}
  {\bibfnamefont {R.~W.}\ \bibnamefont {Boyd}}, \ and\ \bibinfo {author}
  {\bibfnamefont {E.}~\bibnamefont {Karimi}},\ }\href {\doibase
  10.1103/PhysRevLett.117.154801} {\bibfield  {journal} {\bibinfo  {journal}
  {Phys. Rev. Lett.}\ }\textbf {\bibinfo {volume} {117}},\ \bibinfo {pages}
  {154801} (\bibinfo {year} {2016})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Thirunavukkarasu}\ \emph {et~al.}(2017)\citenamefont
  {Thirunavukkarasu}, \citenamefont {Mousley}, \citenamefont {Babiker},\ and\
  \citenamefont {Yuan}}]{TTMousleyBabikerY17}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {G.}~\bibnamefont
  {Thirunavukkarasu}}, \bibinfo {author} {\bibfnamefont {M.}~\bibnamefont
  {Mousley}}, \bibinfo {author} {\bibfnamefont {M.}~\bibnamefont {Babiker}}, \
  and\ \bibinfo {author} {\bibfnamefont {J.}~\bibnamefont {Yuan}},\ }\href
  {\doibase 10.1098/rsta.2015.0438} {\bibfield  {journal} {\bibinfo  {journal}
  {Phil. Trans. R. Soc. Lond. A}\ }\textbf {\bibinfo {volume} {375}} (\bibinfo
  {year} {2017}),\ 10.1098/rsta.2015.0438}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Clark}\ \emph {et~al.}(2015)\citenamefont {Clark},
  \citenamefont {Barankov}, \citenamefont {Huber}, \citenamefont {Arif},
  \citenamefont {Cory},\ and\ \citenamefont {Pushin}}]{Clark+:Nature525:2015}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {C.~W.}\ \bibnamefont
  {Clark}}, \bibinfo {author} {\bibfnamefont {R.}~\bibnamefont {Barankov}},
  \bibinfo {author} {\bibfnamefont {M.~G.}\ \bibnamefont {Huber}}, \bibinfo
  {author} {\bibfnamefont {M.}~\bibnamefont {Arif}}, \bibinfo {author}
  {\bibfnamefont {D.~G.}\ \bibnamefont {Cory}}, \ and\ \bibinfo {author}
  {\bibfnamefont {D.~A.}\ \bibnamefont {Pushin}},\ }\href {\doibase
  10.1038/nature15265} {\bibfield  {journal} {\bibinfo  {journal} {Nature}\
  }\textbf {\bibinfo {volume} {525}},\ \bibinfo {pages} {504} (\bibinfo {year}
  {2015})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {{Hayrapetyan, Armen G.}}\ \emph
  {et~al.}(2013)\citenamefont {{Hayrapetyan, Armen G.}}, \citenamefont
  {{Matula, Oliver}}, \citenamefont {{Surzhykov, Andrey}},\ and\ \citenamefont
  {{Fritzsche, Stephan}}}]{HayrapetyanMatulaSurzhykovFritzsche13}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibnamefont {{Hayrapetyan, Armen
  G.}}}, \bibinfo {author} {\bibnamefont {{Matula, Oliver}}}, \bibinfo {author}
  {\bibnamefont {{Surzhykov, Andrey}}}, \ and\ \bibinfo {author} {\bibnamefont
  {{Fritzsche, Stephan}}},\ }\href {\doibase 10.1140/epjd/e2013-30191-x}
  {\bibfield  {journal} {\bibinfo  {journal} {Eur. Phys. J. D}\ }\textbf
  {\bibinfo {volume} {67}},\ \bibinfo {pages} {167} (\bibinfo {year}
  {2013})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Lembessis}\ \emph {et~al.}(2014)\citenamefont
  {Lembessis}, \citenamefont {Ellinas}, \citenamefont {Babiker},\ and\
  \citenamefont {Al-Dossary}}]{LEBabikerAlD14}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {V.~E.}\ \bibnamefont
  {Lembessis}}, \bibinfo {author} {\bibfnamefont {D.}~\bibnamefont {Ellinas}},
  \bibinfo {author} {\bibfnamefont {M.}~\bibnamefont {Babiker}}, \ and\
  \bibinfo {author} {\bibfnamefont {O.}~\bibnamefont {Al-Dossary}},\ }\href
  {\doibase 10.1103/PhysRevA.89.053616} {\bibfield  {journal} {\bibinfo
  {journal} {Phys. Rev. A}\ }\textbf {\bibinfo {volume} {89}},\ \bibinfo
  {pages} {053616} (\bibinfo {year} {2014})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Bliokh}\ \emph {et~al.}(2012)\citenamefont {Bliokh},
  \citenamefont {Schattschneider}, \citenamefont {Verbeeck},\ and\
  \citenamefont {Nori}}]{BliokhSchattschneiderVNori12}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {K.~Y.}\ \bibnamefont
  {Bliokh}}, \bibinfo {author} {\bibfnamefont {P.}~\bibnamefont
  {Schattschneider}}, \bibinfo {author} {\bibfnamefont {J.}~\bibnamefont
  {Verbeeck}}, \ and\ \bibinfo {author} {\bibfnamefont {F.}~\bibnamefont
  {Nori}},\ }\href {\doibase 10.1103/PhysRevX.2.041011} {\bibfield  {journal}
  {\bibinfo  {journal} {Phys. Rev. X}\ }\textbf {\bibinfo {volume} {2}},\
  \bibinfo {pages} {041011} (\bibinfo {year} {2012})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Greenshields}\ \emph {et~al.}(2012)\citenamefont
  {Greenshields}, \citenamefont {Stamps},\ and\ \citenamefont
  {Franke-Arnold}}]{GreenshieldsStampsFranke-Arnold12}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {C.}~\bibnamefont
  {Greenshields}}, \bibinfo {author} {\bibfnamefont {R.~L.}\ \bibnamefont
  {Stamps}}, \ and\ \bibinfo {author} {\bibfnamefont {S.}~\bibnamefont
  {Franke-Arnold}},\ }\href {http://stacks.iop.org/1367-2630/14/i=10/a=103040}
  {\bibfield  {journal} {\bibinfo  {journal} {New J. Phys.}\ }\textbf {\bibinfo
  {volume} {14}},\ \bibinfo {pages} {103040} (\bibinfo {year}
  {2012})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Greenshields}\ \emph {et~al.}(2015)\citenamefont
  {Greenshields}, \citenamefont {Franke-Arnold},\ and\ \citenamefont
  {Stamps}}]{GreenshieldsFranke-ArnoldStamps15}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {C.~R.}\ \bibnamefont
  {Greenshields}}, \bibinfo {author} {\bibfnamefont {S.}~\bibnamefont
  {Franke-Arnold}}, \ and\ \bibinfo {author} {\bibfnamefont {R.~L.}\
  \bibnamefont {Stamps}},\ }\href
  {http://stacks.iop.org/1367-2630/17/i=9/a=093015} {\bibfield  {journal}
  {\bibinfo  {journal} {New J. Phys.}\ }\textbf {\bibinfo {volume} {17}},\
  \bibinfo {pages} {093015} (\bibinfo {year} {2015})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Karlovets}(2012)}]{Karlovets12}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {D.~V.}\ \bibnamefont
  {Karlovets}},\ }\href {\doibase 10.1103/PhysRevA.86.062102} {\bibfield
  {journal} {\bibinfo  {journal} {Phys. Rev. A}\ }\textbf {\bibinfo {volume}
  {86}},\ \bibinfo {pages} {062102} (\bibinfo {year} {2012})}\BibitemShut
  {NoStop}%
\bibitem [{\citenamefont {Hayrapetyan}\ \emph {et~al.}(2014)\citenamefont
  {Hayrapetyan}, \citenamefont {Matula}, \citenamefont {Aiello}, \citenamefont
  {Surzhikov},\ and\ \citenamefont
  {Fritzsche}}]{HayrapetyanMAielloSFritzsche14}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {A.~G.}\ \bibnamefont
  {Hayrapetyan}}, \bibinfo {author} {\bibfnamefont {O.}~\bibnamefont {Matula}},
  \bibinfo {author} {\bibfnamefont {A.}~\bibnamefont {Aiello}}, \bibinfo
  {author} {\bibfnamefont {A.}~\bibnamefont {Surzhikov}}, \ and\ \bibinfo
  {author} {\bibfnamefont {S.}~\bibnamefont {Fritzsche}},\ }\href {\doibase
  10.1103/PhysRevLett.112.134801} {\bibfield  {journal} {\bibinfo  {journal}
  {Phys. Rev. Lett}\ }\textbf {\bibinfo {volume} {112}},\ \bibinfo {pages}
  {134801} (\bibinfo {year} {2014})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Guzzinati}\ \emph {et~al.}(2013)\citenamefont
  {Guzzinati}, \citenamefont {Schattschneider}, \citenamefont {Bliokh},
  \citenamefont {Nori},\ and\ \citenamefont
  {Verbeeck}}]{GSchattshneiderBliokhNVerbeeck13}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {G.}~\bibnamefont
  {Guzzinati}}, \bibinfo {author} {\bibfnamefont {P.}~\bibnamefont
  {Schattschneider}}, \bibinfo {author} {\bibfnamefont {K.~Y.}\ \bibnamefont
  {Bliokh}}, \bibinfo {author} {\bibfnamefont {F.}~\bibnamefont {Nori}}, \ and\
  \bibinfo {author} {\bibfnamefont {J.}~\bibnamefont {Verbeeck}},\ }\href
  {\doibase 10.1103/PhysRevLett.110.093601} {\bibfield  {journal} {\bibinfo
  {journal} {Phys. Rev. Lett.}\ }\textbf {\bibinfo {volume} {110}},\ \bibinfo
  {pages} {093601} (\bibinfo {year} {2013})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Schattschneider}\ \emph {et~al.}(2014)\citenamefont
  {Schattschneider}, \citenamefont {Schachinger}, \citenamefont
  {St\"oger-Pollach}, \citenamefont {L\"offler}, \citenamefont
  {Steiger-Thirsfeld}, \citenamefont {Bliokh},\ and\ \citenamefont
  {Nori}}]{SchattschneiderSS-PLofflerS-TBliokhNori14}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {P.}~\bibnamefont
  {Schattschneider}}, \bibinfo {author} {\bibfnamefont {T.}~\bibnamefont
  {Schachinger}}, \bibinfo {author} {\bibfnamefont {M.}~\bibnamefont
  {St\"oger-Pollach}}, \bibinfo {author} {\bibfnamefont {S.}~\bibnamefont
  {L\"offler}}, \bibinfo {author} {\bibfnamefont {A.}~\bibnamefont
  {Steiger-Thirsfeld}}, \bibinfo {author} {\bibfnamefont {K.~Y.}\ \bibnamefont
  {Bliokh}}, \ and\ \bibinfo {author} {\bibfnamefont {F.}~\bibnamefont
  {Nori}},\ }\href {\doibase 10.1038/ncomms5586} {\bibfield  {journal}
  {\bibinfo  {journal} {Nat. Comm.}\ }\textbf {\bibinfo {volume} {5}},\
  \bibinfo {pages} {4586} (\bibinfo {year} {2014})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Bialynicki-Birula}(2004)}]{Bialynicki-Birula04}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {I.}~\bibnamefont
  {Bialynicki-Birula}},\ }\href {\doibase 10.1103/PhysRevLett.93.020402}
  {\bibfield  {journal} {\bibinfo  {journal} {Phys. Rev. Lett.}\ }\textbf
  {\bibinfo {volume} {93}},\ \bibinfo {pages} {020402} (\bibinfo {year}
  {2004})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Bialynicki-Birula}\ \emph {et~al.}(2005)\citenamefont
  {Bialynicki-Birula}, \citenamefont {Bialynicka-Birula},\ and\ \citenamefont
  {Chmura}}]{Bialynicki-BirulaChmura05}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibnamefont {Bialynicki-Birula}},
  \bibinfo {author} {\bibfnamefont {Z.}~\bibnamefont {Bialynicka-Birula}}, \
  and\ \bibinfo {author} {\bibfnamefont {B.}~\bibnamefont {Chmura}},\
  }\href@noop {} {\bibfield  {journal} {\bibinfo  {journal} {Laser physics}\
  }\textbf {\bibinfo {volume} {15}},\ \bibinfo {pages} {1371} (\bibinfo {year}
  {2005})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Bialynicki-Birula}\ and\ \citenamefont
  {Radozycki}(2006)}]{Bialynicki-BirulaRadozycki06}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {I.}~\bibnamefont
  {Bialynicki-Birula}}\ and\ \bibinfo {author} {\bibfnamefont {T.}~\bibnamefont
  {Radozycki}},\ }\href {\doibase 10.1103/PhysRevA.73.052114} {\bibfield
  {journal} {\bibinfo  {journal} {Phys. Rev. A}\ }\textbf {\bibinfo {volume}
  {73}},\ \bibinfo {pages} {052114} (\bibinfo {year} {2006})}\BibitemShut
  {NoStop}%
\bibitem [{\citenamefont {Ivanov}(2011)}]{Ivanov11}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {I.~P.}\ \bibnamefont
  {Ivanov}},\ }\href {\doibase 10.1103/PhysRevD.83.093001} {\bibfield
  {journal} {\bibinfo  {journal} {Phys. Rev. D}\ }\textbf {\bibinfo {volume}
  {83}},\ \bibinfo {pages} {093001} (\bibinfo {year} {2011})}\BibitemShut
  {NoStop}%
\bibitem [{\citenamefont {Ivanov}\ and\ \citenamefont
  {Serbo}(2011)}]{IvanovSerbo11}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {I.~P.}\ \bibnamefont
  {Ivanov}}\ and\ \bibinfo {author} {\bibfnamefont {V.~G.}\ \bibnamefont
  {Serbo}},\ }\href {\doibase 10.1103/PhysRevA.84.033804} {\bibfield  {journal}
  {\bibinfo  {journal} {Phys. Rev. A}\ }\textbf {\bibinfo {volume} {84}},\
  \bibinfo {pages} {033804} (\bibinfo {year} {2011})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Ivanov}(2012)}]{Ivanov12}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {I.~P.}\ \bibnamefont
  {Ivanov}},\ }\href {\doibase 10.1103/PhysRevD.85.076001} {\bibfield
  {journal} {\bibinfo  {journal} {Phys. Rev. D}\ }\textbf {\bibinfo {volume}
  {85}},\ \bibinfo {pages} {076001} (\bibinfo {year} {2012})}\BibitemShut
  {NoStop}%
\bibitem [{\citenamefont {Ivanov}\ \emph {et~al.}(2016)\citenamefont {Ivanov},
  \citenamefont {Seipt}, \citenamefont {Surzhykov},\ and\ \citenamefont
  {Fritzsche}}]{IvanofSSurzhykovFritzsche16}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {I.~P.}\ \bibnamefont
  {Ivanov}}, \bibinfo {author} {\bibfnamefont {D.}~\bibnamefont {Seipt}},
  \bibinfo {author} {\bibfnamefont {A.}~\bibnamefont {Surzhykov}}, \ and\
  \bibinfo {author} {\bibfnamefont {S.}~\bibnamefont {Fritzsche}},\ }\href
  {\doibase 10.1103/PhysRevD.94.076001} {\bibfield  {journal} {\bibinfo
  {journal} {Phys. Rev. D}\ }\textbf {\bibinfo {volume} {94}},\ \bibinfo
  {pages} {076001} (\bibinfo {year} {2016})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Matula}\ \emph {et~al.}(2014)\citenamefont {Matula},
  \citenamefont {Hayrapetyan}, \citenamefont {Serbo}, \citenamefont
  {Surzhykov},\ and\ \citenamefont {Fritzsche}}]{Matulaetal14}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {O.}~\bibnamefont
  {Matula}}, \bibinfo {author} {\bibfnamefont {A.~G.}\ \bibnamefont
  {Hayrapetyan}}, \bibinfo {author} {\bibfnamefont {V.~G.}\ \bibnamefont
  {Serbo}}, \bibinfo {author} {\bibfnamefont {A.}~\bibnamefont {Surzhykov}}, \
  and\ \bibinfo {author} {\bibfnamefont {S.}~\bibnamefont {Fritzsche}},\ }\href
  {http://stacks.iop.org/1367-2630/16/i=5/a=053024} {\bibfield  {journal}
  {\bibinfo  {journal} {New Journal of Physics}\ }\textbf {\bibinfo {volume}
  {16}},\ \bibinfo {pages} {053024} (\bibinfo {year} {2014})}\BibitemShut
  {NoStop}%
\bibitem [{\citenamefont {Serbo}\ \emph {et~al.}(2015)\citenamefont {Serbo},
  \citenamefont {Ivanov}, \citenamefont {Fritzsche}, \citenamefont {Seipt},\
  and\ \citenamefont {Surzhykov}}]{Serboetal15}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {V.~G.}\ \bibnamefont
  {Serbo}}, \bibinfo {author} {\bibfnamefont {I.~P.}\ \bibnamefont {Ivanov}},
  \bibinfo {author} {\bibfnamefont {S.}~\bibnamefont {Fritzsche}}, \bibinfo
  {author} {\bibfnamefont {D.}~\bibnamefont {Seipt}}, \ and\ \bibinfo {author}
  {\bibfnamefont {A.}~\bibnamefont {Surzhykov}},\ }\href {\doibase
  10.1103/PhysRevA.92.012705} {\bibfield  {journal} {\bibinfo  {journal} {Phys.
  Rev. A}\ }\textbf {\bibinfo {volume} {92}},\ \bibinfo {pages} {012705}
  (\bibinfo {year} {2015})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Zaytsev}\ \emph {et~al.}(2017)\citenamefont
  {Zaytsev}, \citenamefont {Serbo},\ and\ \citenamefont
  {Shabajev}}]{ZaytsevSerboShabajev17}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {V.~A.}\ \bibnamefont
  {Zaytsev}}, \bibinfo {author} {\bibfnamefont {V.~G.}\ \bibnamefont {Serbo}},
  \ and\ \bibinfo {author} {\bibfnamefont {V.~M.}\ \bibnamefont {Shabajev}},\
  }\href {\doibase 10.1103/PhysRevA.95.012702} {\bibfield  {journal} {\bibinfo
  {journal} {Phys. Rev. A}\ }\textbf {\bibinfo {volume} {95}},\ \bibinfo
  {pages} {012702} (\bibinfo {year} {2017})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Barnett}(2017)}]{Barnett17}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {S.~M.}\ \bibnamefont
  {Barnett}},\ }\href {\doibase 10.1103/PhysRevLett.118.114802} {\bibfield
  {journal} {\bibinfo  {journal} {Phys. Rev. Lett.}\ }\textbf {\bibinfo
  {volume} {118}},\ \bibinfo {pages} {114802} (\bibinfo {year}
  {2017})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Bialynicki-Birula}\ and\ \citenamefont
  {Bialynicka-Birula}(2017)}]{Bialynicki-Birula17}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {I.}~\bibnamefont
  {Bialynicki-Birula}}\ and\ \bibinfo {author} {\bibfnamefont {Z.}~\bibnamefont
  {Bialynicka-Birula}},\ }\href {\doibase 10.1103/PhysRevLett.118.114801}
  {\bibfield  {journal} {\bibinfo  {journal} {Phys. Rev. Lett.}\ }\textbf
  {\bibinfo {volume} {118}},\ \bibinfo {pages} {114801} (\bibinfo {year}
  {2017})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Dirac}(1928)}]{Dirac28}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {P.~A.~M.}\
  \bibnamefont {Dirac}},\ }\href {\doibase 10.1098/rspa.1928.0023} {\bibfield
  {journal} {\bibinfo  {journal} {Proc. Roy. Soc. Lond. A}\ }\textbf {\bibinfo
  {volume} {117}},\ \bibinfo {pages} {610} (\bibinfo {year}
  {1928})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Lifshitz}\ \emph {et~al.}(1991)\citenamefont
  {Lifshitz}, \citenamefont {Pitajevski},\ and\ \citenamefont
  {Berertstetstki}}]{LL4}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {J.~M.}\ \bibnamefont
  {Lifshitz}}, \bibinfo {author} {\bibfnamefont {L.~P.}\ \bibnamefont
  {Pitajevski}}, \ and\ \bibinfo {author} {\bibfnamefont {W.~B.}\ \bibnamefont
  {Berertstetstki}},\ }\href@noop {} {\emph {\bibinfo {title}
  {Quantenelektrodynamik}}},\ \bibinfo {edition} {7th}\ ed.,\ \bibinfo {series}
  {Lehrbuch der Theoretischen Physik}, Vol.~\bibinfo {volume} {4}\ (\bibinfo
  {publisher} {Harri Deutsch},\ \bibinfo {year} {1991})\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Landau}\ and\ \citenamefont {Lifshitz}(1991)}]{LL3}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {L.~D.}\ \bibnamefont
  {Landau}}\ and\ \bibinfo {author} {\bibfnamefont {J.~M.}\ \bibnamefont
  {Lifshitz}},\ }\href@noop {} {\emph {\bibinfo {title} {Quantum mechanics
  (non-relativistic theory)}}},\ \bibinfo {edition} {3rd}\ ed.,\ \bibinfo
  {series} {Course of theoretical physics}, Vol.~\bibinfo {volume} {3}\
  (\bibinfo  {publisher} {Butterworth-heineman},\ \bibinfo {year}
  {1991})\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Riley}\ \emph {et~al.}(2006)\citenamefont {Riley},
  \citenamefont {Hobson},\ and\ \citenamefont {Bence}}]{RHB}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {K.~F.}\ \bibnamefont
  {Riley}}, \bibinfo {author} {\bibfnamefont {M.~P.}\ \bibnamefont {Hobson}}, \
  and\ \bibinfo {author} {\bibfnamefont {S.~J.}\ \bibnamefont {Bence}},\
  }\href@noop {} {\emph {\bibinfo {title} {Mathematical methods for physics and
  engineering}}},\ \bibinfo {edition} {3rd}\ ed.\ (\bibinfo  {publisher}
  {Cambridge university press},\ \bibinfo {year} {2006})\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Karimi}\ \emph {et~al.}(2012)\citenamefont {Karimi},
  \citenamefont {Marrucci}, \citenamefont {Grillo},\ and\ \citenamefont
  {Santamato}}]{Karimietal12}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {E.}~\bibnamefont
  {Karimi}}, \bibinfo {author} {\bibfnamefont {L.}~\bibnamefont {Marrucci}},
  \bibinfo {author} {\bibfnamefont {V.}~\bibnamefont {Grillo}}, \ and\ \bibinfo
  {author} {\bibfnamefont {E.}~\bibnamefont {Santamato}},\ }\href {\doibase
  10.1103/PhysRevLett.108.044801} {\bibfield  {journal} {\bibinfo  {journal}
  {Phys. Rev. Lett.}\ }\textbf {\bibinfo {volume} {108}},\ \bibinfo {pages}
  {044801} (\bibinfo {year} {2012})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Nsofini}\ \emph {et~al.}(2016)\citenamefont
  {Nsofini}, \citenamefont {Sarenac}, \citenamefont {Wood}, \citenamefont
  {Cory}, \citenamefont {Arif}, \citenamefont {Clark}, \citenamefont {Huber},\
  and\ \citenamefont {Pushin}}]{Nsofinietal16}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {J.}~\bibnamefont
  {Nsofini}}, \bibinfo {author} {\bibfnamefont {D.}~\bibnamefont {Sarenac}},
  \bibinfo {author} {\bibfnamefont {C.~J.}\ \bibnamefont {Wood}}, \bibinfo
  {author} {\bibfnamefont {D.~G.}\ \bibnamefont {Cory}}, \bibinfo {author}
  {\bibfnamefont {M.}~\bibnamefont {Arif}}, \bibinfo {author} {\bibfnamefont
  {C.~W.}\ \bibnamefont {Clark}}, \bibinfo {author} {\bibfnamefont {M.~G.}\
  \bibnamefont {Huber}}, \ and\ \bibinfo {author} {\bibfnamefont {D.~A.}\
  \bibnamefont {Pushin}},\ }\href {\doibase 10.1103/PhysRevA.94.013605}
  {\bibfield  {journal} {\bibinfo  {journal} {Phys. Rev. A}\ }\textbf {\bibinfo
  {volume} {94}},\ \bibinfo {pages} {013605} (\bibinfo {year}
  {2016})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Schattschneider}\ \emph {et~al.}(2017)\citenamefont
  {Schattschneider}, \citenamefont {Grillo},\ and\ \citenamefont
  {Aubry}}]{SchattschneiderGrilloAubry17}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {P.}~\bibnamefont
  {Schattschneider}}, \bibinfo {author} {\bibfnamefont {V.}~\bibnamefont
  {Grillo}}, \ and\ \bibinfo {author} {\bibfnamefont {D.}~\bibnamefont
  {Aubry}},\ }\href {\doibase https://doi.org/10.1016/j.ultramic.2016.11.029}
  {\bibfield  {journal} {\bibinfo  {journal} {Ultramicroscopy}\ }\textbf
  {\bibinfo {volume} {176}},\ \bibinfo {pages} {188 } (\bibinfo {year}
  {2017})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Dellweg}\ and\ \citenamefont
  {M\"uller}(2017)}]{DellwegMuller17}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {M.~M.}\ \bibnamefont
  {Dellweg}}\ and\ \bibinfo {author} {\bibfnamefont {C.}~\bibnamefont
  {M\"uller}},\ }\href {\doibase 10.1103/PhysRevLett.118.070403} {\bibfield
  {journal} {\bibinfo  {journal} {Phys. Rev. Lett.}\ }\textbf {\bibinfo
  {volume} {118}},\ \bibinfo {pages} {070403} (\bibinfo {year}
  {2017})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Ohanian}(1986)}]{Ohanian85}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {H.~C.}\ \bibnamefont
  {Ohanian}},\ }\href {\doibase 10.1119/1.14580} {\bibfield  {journal}
  {\bibinfo  {journal} {Am. J. Phys.}\ }\textbf {\bibinfo {volume} {54}},\
  \bibinfo {pages} {500} (\bibinfo {year} {1986})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Bliokh}\ and\ \citenamefont
  {Nori}(2012{\natexlab{b}})}]{BliokhNori12}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {K.~Y.}\ \bibnamefont
  {Bliokh}}\ and\ \bibinfo {author} {\bibfnamefont {F.}~\bibnamefont {Nori}},\
  }\href {\doibase 10.1103/PhysRevA.85.061801} {\bibfield  {journal} {\bibinfo
  {journal} {Phys. Rev. A}\ }\textbf {\bibinfo {volume} {85}},\ \bibinfo
  {pages} {061801} (\bibinfo {year} {2012}{\natexlab{b}})}\BibitemShut
  {NoStop}%
\bibitem [{\citenamefont {Bliokh}\ \emph {et~al.}(2014)\citenamefont {Bliokh},
  \citenamefont {Bekshaev},\ and\ \citenamefont {Nori}}]{BliokhBekshaevNori14}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {K.~Y.}\ \bibnamefont
  {Bliokh}}, \bibinfo {author} {\bibfnamefont {A.~Y.}\ \bibnamefont
  {Bekshaev}}, \ and\ \bibinfo {author} {\bibfnamefont {F.}~\bibnamefont
  {Nori}},\ }\href {\doibase 10.1038/ncomms4300} {\bibfield  {journal}
  {\bibinfo  {journal} {Nat. Comm.}\ }\textbf {\bibinfo {volume} {5}},\
  \bibinfo {pages} {3300} (\bibinfo {year} {2014})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Bekshaev}\ \emph {et~al.}(2015)\citenamefont
  {Bekshaev}, \citenamefont {Bliokh},\ and\ \citenamefont
  {Nori}}]{BekshaevBliokhNori15}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {A.~Y.}\ \bibnamefont
  {Bekshaev}}, \bibinfo {author} {\bibfnamefont {K.~Y.}\ \bibnamefont
  {Bliokh}}, \ and\ \bibinfo {author} {\bibfnamefont {F.}~\bibnamefont
  {Nori}},\ }\href {\doibase 10.1103/PhysRevX.5.011039} {\bibfield  {journal}
  {\bibinfo  {journal} {Phys. Rev. X}\ }\textbf {\bibinfo {volume} {5}},\
  \bibinfo {pages} {011039} (\bibinfo {year} {2015})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Neugebauer}\ \emph {et~al.}(2015)\citenamefont
  {Neugebauer}, \citenamefont {Bauer}, \citenamefont {Aiello},\ and\
  \citenamefont {Banzer}}]{NBAB15}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {M.}~\bibnamefont
  {Neugebauer}}, \bibinfo {author} {\bibfnamefont {T.}~\bibnamefont {Bauer}},
  \bibinfo {author} {\bibfnamefont {A.}~\bibnamefont {Aiello}}, \ and\ \bibinfo
  {author} {\bibfnamefont {P.}~\bibnamefont {Banzer}},\ }\href {\doibase
  10.1103/PhysRevLett.114.063901} {\bibfield  {journal} {\bibinfo  {journal}
  {Phys. Rev. Lett.}\ }\textbf {\bibinfo {volume} {114}},\ \bibinfo {pages}
  {063901} (\bibinfo {year} {2015})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Antognozzi}\ \emph {et~al.}(2016)\citenamefont
  {Antognozzi}, \citenamefont {Berminham}, \citenamefont {Harniman},
  \citenamefont {Simpson}, \citenamefont {Senior}, \citenamefont {Hayward},
  \citenamefont {Hoerber}, \citenamefont {Dennis}, \citenamefont {Bekshaev},
  \citenamefont {Bliokh},\ and\ \citenamefont {Nori}}]{Antognozzietal16}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {M.}~\bibnamefont
  {Antognozzi}}, \bibinfo {author} {\bibfnamefont {C.~R.}\ \bibnamefont
  {Berminham}}, \bibinfo {author} {\bibfnamefont {R.~l.}\ \bibnamefont
  {Harniman}}, \bibinfo {author} {\bibfnamefont {S.}~\bibnamefont {Simpson}},
  \bibinfo {author} {\bibfnamefont {J.}~\bibnamefont {Senior}}, \bibinfo
  {author} {\bibfnamefont {R.}~\bibnamefont {Hayward}}, \bibinfo {author}
  {\bibfnamefont {H.}~\bibnamefont {Hoerber}}, \bibinfo {author} {\bibfnamefont
  {M.~R.}\ \bibnamefont {Dennis}}, \bibinfo {author} {\bibfnamefont
  {Y.}~\bibnamefont {Bekshaev}, \bibfnamefont {Alexandr}}, \bibinfo {author}
  {\bibfnamefont {K.~Y.}\ \bibnamefont {Bliokh}}, \ and\ \bibinfo {author}
  {\bibfnamefont {F.}~\bibnamefont {Nori}},\ }\href {\doibase
  10.1038/NPHYS3732} {\bibfield  {journal} {\bibinfo  {journal} {Nat. Phys.}\
  }\textbf {\bibinfo {volume} {12}},\ \bibinfo {pages} {731} (\bibinfo {year}
  {2016})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Kien}\ and\ \citenamefont
  {Hakuta}(2008)}]{KienHakuta08}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {F.~L.}\ \bibnamefont
  {Kien}}\ and\ \bibinfo {author} {\bibfnamefont {K.}~\bibnamefont {Hakuta}},\
  }\href {\doibase 10.1103/PhysRevA.77.013801} {\bibfield  {journal} {\bibinfo
  {journal} {Phys. Rev. A}\ }\textbf {\bibinfo {volume} {77}},\ \bibinfo
  {pages} {013801} (\bibinfo {year} {2008})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Kien}\ and\ \citenamefont
  {Rauschenbeutel}(2014)}]{KienRauschenbeutel14a}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {F.~L.}\ \bibnamefont
  {Kien}}\ and\ \bibinfo {author} {\bibfnamefont {A.}~\bibnamefont
  {Rauschenbeutel}},\ }\href {\doibase 10.1103/PhysRevA.90.023805} {\bibfield
  {journal} {\bibinfo  {journal} {Phys. Rev. A}\ }\textbf {\bibinfo {volume}
  {90}},\ \bibinfo {pages} {023805} (\bibinfo {year} {2014})}\BibitemShut
  {NoStop}%
\bibitem [{\citenamefont {Junge}\ \emph {et~al.}(2013)\citenamefont {Junge},
  \citenamefont {O'Shea}, \citenamefont {Volz},\ and\ \citenamefont
  {Rauschenbeutel}}]{PhysRevLett.110.213604}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {C.}~\bibnamefont
  {Junge}}, \bibinfo {author} {\bibfnamefont {D.}~\bibnamefont {O'Shea}},
  \bibinfo {author} {\bibfnamefont {J.}~\bibnamefont {Volz}}, \ and\ \bibinfo
  {author} {\bibfnamefont {A.}~\bibnamefont {Rauschenbeutel}},\ }\href
  {\doibase 10.1103/PhysRevLett.110.213604} {\bibfield  {journal} {\bibinfo
  {journal} {Phys. Rev. Lett.}\ }\textbf {\bibinfo {volume} {110}},\ \bibinfo
  {pages} {213604} (\bibinfo {year} {2013})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Petersen}\ \emph {et~al.}(2014)\citenamefont
  {Petersen}, \citenamefont {Volz},\ and\ \citenamefont
  {Rauschenbeutel}}]{Petersen67}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {J.}~\bibnamefont
  {Petersen}}, \bibinfo {author} {\bibfnamefont {J.}~\bibnamefont {Volz}}, \
  and\ \bibinfo {author} {\bibfnamefont {A.}~\bibnamefont {Rauschenbeutel}},\
  }\href {\doibase 10.1126/science.1257671} {\bibfield  {journal} {\bibinfo
  {journal} {Science}\ }\textbf {\bibinfo {volume} {346}},\ \bibinfo {pages}
  {67} (\bibinfo {year} {2014})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Mitsch}\ \emph {et~al.}(2014)\citenamefont {Mitsch},
  \citenamefont {Sayrin}, \citenamefont {Albrecht}, \citenamefont
  {Schneewei\ss{}},\ and\ \citenamefont
  {Rauschenbeutel}}]{MSASRauschenbeutel14}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {R.}~\bibnamefont
  {Mitsch}}, \bibinfo {author} {\bibfnamefont {C.}~\bibnamefont {Sayrin}},
  \bibinfo {author} {\bibfnamefont {B.}~\bibnamefont {Albrecht}}, \bibinfo
  {author} {\bibfnamefont {P.}~\bibnamefont {Schneewei\ss{}}}, \ and\ \bibinfo
  {author} {\bibfnamefont {A.}~\bibnamefont {Rauschenbeutel}},\ }\href
  {\doibase 10.1038/ncomms6713} {\bibfield  {journal} {\bibinfo  {journal}
  {Nat. Comm.}\ }\textbf {\bibinfo {volume} {5}},\ \bibinfo {pages} {5713}
  (\bibinfo {year} {2014})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Sayrin}\ \emph {et~al.}(2015)\citenamefont {Sayrin},
  \citenamefont {Junge}, \citenamefont {Mitsch}, \citenamefont {Albrecht},
  \citenamefont {O'Shea}, \citenamefont {Schneeweiss}, \citenamefont {Volz},\
  and\ \citenamefont {Rauschenbeutel}}]{PhysRevX.5.041036}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {C.}~\bibnamefont
  {Sayrin}}, \bibinfo {author} {\bibfnamefont {C.}~\bibnamefont {Junge}},
  \bibinfo {author} {\bibfnamefont {R.}~\bibnamefont {Mitsch}}, \bibinfo
  {author} {\bibfnamefont {B.}~\bibnamefont {Albrecht}}, \bibinfo {author}
  {\bibfnamefont {D.}~\bibnamefont {O'Shea}}, \bibinfo {author} {\bibfnamefont
  {P.}~\bibnamefont {Schneeweiss}}, \bibinfo {author} {\bibfnamefont
  {J.}~\bibnamefont {Volz}}, \ and\ \bibinfo {author} {\bibfnamefont
  {A.}~\bibnamefont {Rauschenbeutel}},\ }\href {\doibase
  10.1103/PhysRevX.5.041036} {\bibfield  {journal} {\bibinfo  {journal} {Phys.
  Rev. X}\ }\textbf {\bibinfo {volume} {5}},\ \bibinfo {pages} {041036}
  (\bibinfo {year} {2015})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Ramos}\ \emph {et~al.}(2014)\citenamefont {Ramos},
  \citenamefont {Pichler}, \citenamefont {Daley},\ and\ \citenamefont
  {Zoller}}]{PhysRevLett.113.237203}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {T.}~\bibnamefont
  {Ramos}}, \bibinfo {author} {\bibfnamefont {H.}~\bibnamefont {Pichler}},
  \bibinfo {author} {\bibfnamefont {A.~J.}\ \bibnamefont {Daley}}, \ and\
  \bibinfo {author} {\bibfnamefont {P.}~\bibnamefont {Zoller}},\ }\href
  {\doibase 10.1103/PhysRevLett.113.237203} {\bibfield  {journal} {\bibinfo
  {journal} {Phys. Rev. Lett.}\ }\textbf {\bibinfo {volume} {113}},\ \bibinfo
  {pages} {237203} (\bibinfo {year} {2014})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Pichler}\ \emph {et~al.}(2015)\citenamefont
  {Pichler}, \citenamefont {Ramos}, \citenamefont {Daley},\ and\ \citenamefont
  {Zoller}}]{PhysRevA.91.042116}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {H.}~\bibnamefont
  {Pichler}}, \bibinfo {author} {\bibfnamefont {T.}~\bibnamefont {Ramos}},
  \bibinfo {author} {\bibfnamefont {A.~J.}\ \bibnamefont {Daley}}, \ and\
  \bibinfo {author} {\bibfnamefont {P.}~\bibnamefont {Zoller}},\ }\href
  {\doibase 10.1103/PhysRevA.91.042116} {\bibfield  {journal} {\bibinfo
  {journal} {Phys. Rev. A}\ }\textbf {\bibinfo {volume} {91}},\ \bibinfo
  {pages} {042116} (\bibinfo {year} {2015})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Ramos}\ \emph {et~al.}(2016)\citenamefont {Ramos},
  \citenamefont {Vermersch}, \citenamefont {Hauke}, \citenamefont {Pichler},\
  and\ \citenamefont {Zoller}}]{PhysRevA.93.062104}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {T.}~\bibnamefont
  {Ramos}}, \bibinfo {author} {\bibfnamefont {B.}~\bibnamefont {Vermersch}},
  \bibinfo {author} {\bibfnamefont {P.}~\bibnamefont {Hauke}}, \bibinfo
  {author} {\bibfnamefont {H.}~\bibnamefont {Pichler}}, \ and\ \bibinfo
  {author} {\bibfnamefont {P.}~\bibnamefont {Zoller}},\ }\href {\doibase
  10.1103/PhysRevA.93.062104} {\bibfield  {journal} {\bibinfo  {journal} {Phys.
  Rev. A}\ }\textbf {\bibinfo {volume} {93}},\ \bibinfo {pages} {062104}
  (\bibinfo {year} {2016})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Vermersch}\ \emph {et~al.}(2016)\citenamefont
  {Vermersch}, \citenamefont {Ramos}, \citenamefont {Hauke},\ and\
  \citenamefont {Zoller}}]{PhysRevA.93.063830}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {B.}~\bibnamefont
  {Vermersch}}, \bibinfo {author} {\bibfnamefont {T.}~\bibnamefont {Ramos}},
  \bibinfo {author} {\bibfnamefont {P.}~\bibnamefont {Hauke}}, \ and\ \bibinfo
  {author} {\bibfnamefont {P.}~\bibnamefont {Zoller}},\ }\href {\doibase
  10.1103/PhysRevA.93.063830} {\bibfield  {journal} {\bibinfo  {journal} {Phys.
  Rev. A}\ }\textbf {\bibinfo {volume} {93}},\ \bibinfo {pages} {063830}
  (\bibinfo {year} {2016})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {S\"ollner}\ \emph {et~al.}(2015)\citenamefont
  {S\"ollner}, \citenamefont {Mahmoodian}, \citenamefont {Hansen},
  \citenamefont {Lidolo}, \citenamefont {Javadi}, \citenamefont {Kirsanske},
  \citenamefont {Pregnolato}, \citenamefont {El-Ella}, \citenamefont {Lee},
  \citenamefont {Song}, \citenamefont {Stobbe},\ and\ \citenamefont
  {Lohdahl}}]{Kopenhagen15}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {I.}~\bibnamefont
  {S\"ollner}}, \bibinfo {author} {\bibfnamefont {S.}~\bibnamefont
  {Mahmoodian}}, \bibinfo {author} {\bibfnamefont {S.~L.}\ \bibnamefont
  {Hansen}}, \bibinfo {author} {\bibfnamefont {L.}~\bibnamefont {Lidolo}},
  \bibinfo {author} {\bibfnamefont {A.}~\bibnamefont {Javadi}}, \bibinfo
  {author} {\bibfnamefont {G.}~\bibnamefont {Kirsanske}}, \bibinfo {author}
  {\bibfnamefont {T.}~\bibnamefont {Pregnolato}}, \bibinfo {author}
  {\bibfnamefont {H.}~\bibnamefont {El-Ella}}, \bibinfo {author} {\bibfnamefont
  {E.~H.}\ \bibnamefont {Lee}}, \bibinfo {author} {\bibfnamefont {J.~D.}\
  \bibnamefont {Song}}, \bibinfo {author} {\bibfnamefont {S.}~\bibnamefont
  {Stobbe}}, \ and\ \bibinfo {author} {\bibfnamefont {P.}~\bibnamefont
  {Lohdahl}},\ }\href {\doibase 10.1038/NNANO.2015.159} {\bibfield  {journal}
  {\bibinfo  {journal} {Nat. Nanotech.}\ }\textbf {\bibinfo {volume} {10}},\
  \bibinfo {pages} {775} (\bibinfo {year} {2015})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {Gordon}(1928)}]{Gordon1928}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {W.}~\bibnamefont
  {Gordon}},\ }\href {\doibase 10.1007/BF01327881} {\bibfield  {journal}
  {\bibinfo  {journal} {Zeitschrift f\"ur Physik}\ }\textbf {\bibinfo {volume}
  {50}},\ \bibinfo {pages} {630} (\bibinfo {year} {1928})}\BibitemShut
  {NoStop}%
\bibitem [{Note1()}]{Note1}%
  \BibitemOpen
  \bibinfo {note} {Both Gordon and Ohanian only consider the free case, but
  because the probability current is gauge invariant, one can see that one must
  use the gauge covariant momentum operator if fields are present.}\BibitemShut
  {Stop}%
\bibitem [{\citenamefont {He}\ \emph {et~al.}(1995)\citenamefont {He},
  \citenamefont {Friese}, \citenamefont {Heckenberg},\ and\ \citenamefont
  {Rubinsztein-Dunlop}}]{HeFrieseHeckenbergRubinsztein-Dunlop95}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {H.}~\bibnamefont
  {He}}, \bibinfo {author} {\bibfnamefont {M.~E.~J.}\ \bibnamefont {Friese}},
  \bibinfo {author} {\bibfnamefont {N.~R.}\ \bibnamefont {Heckenberg}}, \ and\
  \bibinfo {author} {\bibfnamefont {H.}~\bibnamefont {Rubinsztein-Dunlop}},\
  }\href {\doibase 10.1103/PhysRevLett.75.826} {\bibfield  {journal} {\bibinfo
  {journal} {Phys. Rev. Lett}\ }\textbf {\bibinfo {volume} {75}},\ \bibinfo
  {pages} {826} (\bibinfo {year} {1995})}\BibitemShut {NoStop}%
\bibitem [{\citenamefont {O'Neil}\ \emph {et~al.}(2002)\citenamefont {O'Neil},
  \citenamefont {MacVicar}, \citenamefont {Allen},\ and\ \citenamefont
  {Padgett}}]{PhysRevLett.88.053601}%
  \BibitemOpen
  \bibfield  {author} {\bibinfo {author} {\bibfnamefont {A.~T.}\ \bibnamefont
  {O'Neil}}, \bibinfo {author} {\bibfnamefont {I.}~\bibnamefont {MacVicar}},
  \bibinfo {author} {\bibfnamefont {L.}~\bibnamefont {Allen}}, \ and\ \bibinfo
  {author} {\bibfnamefont {M.~J.}\ \bibnamefont {Padgett}},\ }\href {\doibase
  10.1103/PhysRevLett.88.053601} {\bibfield  {journal} {\bibinfo  {journal}
  {Phys. Rev. Lett.}\ }\textbf {\bibinfo {volume} {88}},\ \bibinfo {pages}
  {053601} (\bibinfo {year} {2002})}\BibitemShut {NoStop}%
\end{thebibliography}%
