Artificial magnetic field induced by an evanescent wave

Cold atomic gases are perfect laboratories for realization of quantum simulators. In order to simulate solid state systems in the presence of magnetic fields special effort has to be made because atoms are charge neutral. There are different methods for realization of artificial magnetic fields, that is the creation of specific conditions so that the motion of neutral particles mimics the dynamics of charged particles in an effective magnetic field. Here, we consider adiabatic motion of atoms in the presence of an evanescent wave. Theoretical description of the adiabatic motion involves artificial vector and scalar potentials related to the Berry phases. Due to the large gradient of the evanescent field amplitude, the potentials can be strong enough to induce measurable effects in cold atomic gases. We show that the resulting artificial magnetic field is able to induce vortices in a Bose-Einstein condensate trapped close to a surface of a prism where the evanescent wave is created. We also analyze motion of an atomic cloud released from a magneto-optical trap that falls down on the surface of the prism. The artificial magnetic field is able to reflect falling atoms that can be observed experimentally.

Cold atomic gases are perfect laboratories for realization of quantum simulators. In order to simulate solid state systems in the presence of magnetic fields special effort has to be made because atoms are charge neutral. There are different methods for realization of artificial magnetic fields, that is the creation of specific conditions so that the motion of neutral particles mimics the dynamics of charged particles in an effective magnetic field. Here, we consider adiabatic motion of atoms in the presence of an evanescent wave. Theoretical description of the adiabatic motion involves artificial vector and scalar potentials related to the Berry phases. Due to the large gradient of the evanescent field amplitude, the potentials can be strong enough to induce measurable effects in cold atomic gases. We show that the resulting artificial magnetic field is able to induce vortices in a Bose-Einstein condensate trapped close to a surface of a prism where the evanescent wave is created. We also analyze motion of an atomic cloud released from a magneto-optical trap that falls down on the surface of the prism. The artificial magnetic field is able to reflect falling atoms that can be observed experimentally.
C old atomic gases are flexible laboratories with a great potential to investigate a variety of problems from many fields of physics 1,2 . Trapping potentials for atoms and mutual atom interactions can be controlled and engineered nearly at will. Mixtures of different atomic species of the Fermi and Bose statistics can be prepared and investigated experimentally.
Atoms are charge neutral and seem not suitable to simulate orbital magnetism. However, there are also methods that allow for a generation of artificial gauge fields, i.e. the creation of specific conditions such that the motion of neutral particles mimics the dynamics of charged particles in an effective magnetic field, see refs 3-5 and references therein. There are also proposals to simulate the spin-orbit coupling with the help of non-Abelian artificial gauge fields [5][6][7][8][9][10][11][12][13][14] . Ultra-cold atoms in synthetic gauge fields open possibilities for deep understanding of fundamental problems like high-temperature superconductivity or strongly interacting counterparts of topological insulators 15 .
In the present article we focus on a method of the synthetic magnetic field creation which involves adiabatic motion of atoms in a laser radiation whose intensity changes in space 4,16 . We assume that atoms are placed close to a dielectric surface and experience an evanescent field that penetrates the region of the atoms location. Large gradient of the generalized Rabi frequency is responsible for the presence of geometrical vector and scalar potentials that are related to the Berry phases 17 . These potentials generate an artificial magnetic field for the charge neutral particles which is able to induce vortices in an atomic cloud cooled down to the quantum degeneracy.
Mirrors are basic elements of optical devices. Similarly atomic mirrors that can reflect beams of atoms are of great experimental interest. They can base on the interaction between magnetic moments of atoms with magnetic fields 18 or take advantage of a spatially dependent light-shift of atomic energy levels. An evanescent field is able to create a strong optical dipole potential for atoms and is a convenient tool to build atomic mirrors [19][20][21][22][23] . Evanescent fields can be created by shining laser radiation on a dielectric surface but they are also connected with surface plasmons that are collective oscillations of free eletrons at a metal surface 24 . Plasmonically tailored dipole potentials and nanofibers [25][26][27][28] offer promising tools for trapping, guiding and manipulating atoms. In the present publication we investigate also under what conditions the artificial magnetic fields created by evanescent waves can be responsible for reflection of atoms falling down on a surface of a prism at temperatures higher than the critical value for the quantum degeneracy.

Results
Adiabatic motion of atoms. We consider a two-level atom interacting with an external laser field. Assume that the atomic energy level difference is v 0 and the atom at rest is located at r. The Hamiltonian of the internal degrees of freedom of the atom interacting with an electromagnetic field is time periodic due to the periodicity of the electromagnetic wave, i.e. H in (t1 2p/ v)5H in (t). According to the Floquet theorem [29][30][31] (an analogue of the Bloch theorem in solid-state physics), the operator (the socalled Floquet Hamiltonian) H in~Hin {i hL t possesses time periodic eigenstates. The corresponding eigenvalues are defined modulo v and are called quasi-energies. Just as in the solid-state physics, one can reduce considerations to a single Floquet zone (the equivalent of the Brillouin zone). Eigenstates jx 1 (r)ae and jx 2 (r)ae of the Floquet Hamiltonian for the two-level atom problem, i.e. the socalled dressed states in the atomic optics context, can be found analytically if the rotating wave approximation (RWA) is applied. The quasi-energy splitting equals is the generalized Rabi frequency. D 5 v 0 2v is the detuning from the resonance frequency and k(r)5d?E(r)/ is the Rabi frequency where d and E(r) stand for the atomic dipole moment and the electric field vector, respectively. In the present paper we consider small detuning and neglect spontaneous emission of an atom. Therefore, the presented results are relevant to, e.g., the long-lived clock transition in ytterbium atoms 32 .
If we assume that the atom is initially prepared in, e.g., the jx 1 (r)ae dressed state and the atom moves sufficiently slowly, the dressed state is adiabatically followed by the atom. Then, non-trivial vector A and scalar W potentials can appear in the Hamiltonian that describes the motion of the center of mass of the atom and the system can mimic dynamics of a charged particle in the presence of a magnetic field 4,5,16 . The vector potential A is the consequence of the Berry phase 17 that emerges due to the adiabatic approximation, i.e.
Thus, the vector potential A can be expressed in the form and the corresponding magnetic field, B 5 = 3 A, depends on the gradient of the phase of the external electromagnetic wave and the gradient of the generalized Rabi frequency V(r). The scalar potential W reads The artificial gauge potentials are called 'geometrical' because they depend only on the spatial variation of the dressed state as one can see in (2) and (3).
In the following we consider a total internal reflection of light in a prism that creates an evanescent wave. An evanescent wave seems to have all important properties necessary to generate artificial gauge potentials for adiabatically moving atoms, i.e. it has a gradient of the phase and a large gradient of the amplitude.
Bose-Einstein condensate in an artificial magnetic field. Let us consider a prism made of dielectric material with the refractive index n . 1 and an electromagnetic plane wave that propagates in the prism. The wave strikes the boundary between the dielectric medium and the vacuum at the angle of incidence h greater than the critical angle h 0 5 arcsin(1/n) for the total internal reflection. The evanescent wave, that appears in the vacuum, propagates along the boundary (x direction) and decays exponentially with increasing distance from the boundary (z direction), see Fig. 1, where E 0 describes the amplitude and the direction of the electric field vector, w(x) 5 xk 0 n sin h is the running phase, is the penetration depth and k 0 5 2p/l is the wavenumber. We have chosen the TE polarization but similar analysis can be performed for the TM polarization. The transmission The dressed states of a two-level atom in the presence of such an evanescent field, obtained within the RWA, read x with energies V(z)/2, 2 V(z)/2 respectively and W(z)5 arctg(jk(x, z)j/D). We assume that slowly moving atoms follow one of the dressed states, e.g. jx 1 (x, z)ae. It is possible because the energy split in the dressed atom picture is V which leads to the separation of the dynamics of each dressed state and allows for the adiabatic elimination of one of them. The condition for the applicability of the adiabatic approximation can be obtained by rewriting the internal state iand solving the Schrödinger equation as a power series in velocity 33 . Then the adiabatic motion requires c 2 j j=1 that gives the range of velocities v=V= x 2 j+x 1 h i j j . The vector potential associated with the adiabatic motion is then The calculation of the curl of the vector potential A allows one to obtain the artificial magnetic field vector, which has a nonzero component in the y direction only, where B 0~ hk 2 0 2 and a(z) 5 jt TE (h)j 2 e 22z/d , and we have introduced the parameter The magnetic field B(z) can be shaped by changing h. The angle of incidence determines both the maximal value of the magnetic field and a range Dz over which B(z) is significant. The strongest field is created if the incident angle of the plane wave h is much greater than the critical angle for the total internal reflection h 0 . Then, however, the magnetic field is present on a small region Dz < d < 1/k 0 . Integrating the artificial magnetic field, Eq. (8), we obtain Thus, the maximal magnetic field is of the order of hk 2 0 . When h approaches h 0 , the penetration depth d increases but the artificial magnetic field becomes weaker B / 1/d.
In order to decide which values of h are suitable for an experiment one should analyze how many vortices in ultra-cold atomic gases the artificial magnetic field is able to create. The vortex density can be expressed by r v 5 B/(2p ). If B keeps significant value in a square of area (Dz) 2 , the number of vortices in this square is (Dz) 2 r v . The field B(z) depends on z coordinate only, thus, the space where the magnetic field is significant forms a layer of width Dz. Therefore, it is more instructive to estimate the number of vortex rows N rows~D z ffiffiffiffi ffi r v p which, for h close to the critical angle h 0 , can be approximated by For n 5 1.4 and l 5 578 nm we obtain N rows 5 1 and Dz < 2.3 mm if h 2 h 0 < 8 ? 10 24 rad while if h 2 h 0 < 10 25 rad, N rows 5 3 and Dz < 20.8 mm. The adjustment of the incidence angle with accuracy of 10 24 rad has been demonstrated in a laboratory 22 but much better accuracy can be attainable in evanescent wave experiments.
In Fig. 2 we show the plots of the magnetic fields B(z) that correspond to these two examples. By changing h one can control how many vortex rows are realized experimentally and then investigate how the spatial arrangement of vortices changes with an increase of N rows . Atoms are charge neutral but to get a sense of the order of magnitude of the artificial magnetic field created by means of the evanescent wave let us assume that atoms possess the elementary charge e. Then, the black curve in Fig. 2 corresponds to magnetic field B/e < 0.3 mT that is present on the region Dz < 10 mm.
An increase of the parameter s, Eq. (9), causes essentially a shift of B(z) towards greater values of z, see Fig. 3. Thus, the location of the region where the artificial magnetic field is present can be suitably chosen by a change of the parameter s. It allows one to trap atomic clouds sufficiently far away from the surface of the prism and consequently eliminate the influence of the van der Waals interaction between the atom and the dielectric wall 19 .
In order to trap atoms close to the surface of the prism an external magnetic trap or additional laser beams have to be applied 22 . The geometric potential W is too weak to overcome the gravitational attraction. In Fig. 4 we show for parameters corresponding to those used in Fig. 2. For example for ytterbium atoms, the maximal force created by this potential is 0.17 mg where g is the gravitational acceleration. Also the optical dipole potential created by the considered evanescent waves can be too weak to keep an atomic cloud above the surface of a prism. Indeed, the artificial magnetic fields suitable for experiments require large penetration depth d. However, by increasing d we decrease the optical dipole force because =V becomes smaller. If the incidence angle of a laser beam is significantly greater than the critical angle h 0 , it is possible to approximate the beam by a plane wave. The situation changes for the incident angles close to h 0 and such a situation is investigated in the present publication. Therefore  we have to consider the realistic laser radiation in the analysis of the artificial magnetic fields induced by evanescent waves.
We consider a laser beam incident with an angle h in on a boundary between a dielectric medium and the vacuum. The beam is represented by a Gaussian superposition of plane waves. The resulting electric field in the vacuum is a sum of two contributions The first contribution describes the evanescent field, i.e. it corresponds to the superposition of plane waves incident with angles h . h 0 , where the last exponential term describes the profile of the beam. l is the distance of the waist of the incident beam from the surface, Dh 5 2/(nk 0 w) describes the Gaussian distribution of the incident angles, where w is the waist of the beam, and w y is the radius of the transverse distribution 22 . The second contribution in Eq. (13), i.e. E 2 (r, t), describes the propagation of waves that strike the surface with h , h 0 and it is given by the similar formula as Eq. (14) but the range of the integration is between 0 and h 0 and izk 0 ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 1{n 2 sin 2 h p substitutes for e 2z/d . Let us assume that the laser beam has experimentally realistic parameters: l 5 578 nm, l 5 680 mm, w 5 440 mm and w y 5 440 mm 22 . We consider the incident angles for which d h in ð Þ=w y and consequently the Rabi frequency changes much more slowly with the change of y than z. Thus, the dominant component of the resulting artificial magnetic field vector is the y-component, i.e. B r ð Þ<{B r ð Þŷ similarly as in the previous section. A new feature for the present configuration is that the shape of B(r) depends both on the incident angle h in and on the parameter s.
In Fig. 5 we show B(0, 0, z) versus z for s 5 5 and for the same incident angles as in the case of a single plane wave, cf. Fig. 2. For h in 2 h 0 5 8 ? 10 24 rad the artificial magnetic field is nearly the same as in Fig. 2 and consequently the plane wave approximation of a laser beam is valid. For h in 2 h 0 5 10 25 rad there is a noticeable difference between the cases of the single plane wave and the realistic laser beam. That is, the maximal field is greater but the range of Dz where the magnetic field is significant is smaller.
So far we have analyzed the creation of the artificial magnetic field by means of an evanescent wave and estimated number of vortices that can be formed in ultra-cold atomic gases in the presence of such fields. In order to confirm the predictions we switch to the numerical simulations within the mean field approximation. We consider a Bose-Einstein condensate trapped in a harmonic potential in the presence of the vector potential A in the two-dimensional (2D) approximation. For the parameters we are going to use, the geometrical scalar potential (3) and the optical dipole potential are very weak and therefore are neglected. The Gross-Pitaevskii equation (GPE) in the units of the harmonic trapping potential reads where m is the chemical potential of the system and g stands for the atomic interaction strength (we assume AEyjyae 5 1). In the numerical simulation we discretize 2D space. A naive approximation of, e.g., h x y(x, z) < [y(x 1 dx, z)2y(x2dx, z)]/(2dx) leads to a discrete version of the GPE which is not gauge invariant. This problem can be overcome by adopting the Schwinger line integral used in the lattice gauge theories 34 . That is, in order to make a discrete version of the energy functional gauge invariant, terms like y*(x, z)y(x 1 dx, z) have to be exchanged by gauge invariant terms where the Schwinger line integral U(x, z; x 1 dx, z) 5 exp(iA x dx). It corresponds to the following substitution in the GPE and similarly for (h z 1 iA z ) 2 y. The resulting discrete GPE is gauge invariant and recovers Eq. (15) in the limit where dx and dz go to zero. The ground states of the GPE found numerically for the artificial magnetic field created by means of the evanescent wave with the   Gaussian profile are presented in Fig. 6. Top and bottom panel in Fig. 6 corresponds to the black and the red curve in Fig. 5, respectively. The interaction coefficient g has been chosen so that the Thomas-Fermi radius of a cloud of ytterbium atoms, in the harmonic potential with the frequency v trap 5 2p 3 16 Hz, is R TF 5 10 mm (top panel) and R TF 5 15 mm (bottom panel). Many different combinations of the total particle number and the trap frequency along the third direction lead to the same coefficient g in the 2D approximation. Therefore, it is meaningful to provide g, or R TF equivalently, only and do not consider a particular choice of the particle number and the trap frequency in the y direction. We analyze the behaviour of the system in the Thomas-Fermi regime. Thus, the healing length of the system, and consequently the size of the vortex cores, is much smaller than R TF and also smaller than the region of space Dz where the artificial magnetic field is significant. The numbers of vortex rows that can be estimated according to where B is half of the maximum value of the magnetic field are 1 and 2, respectively, which is smaller than the actual numbers visible in Fig. 6. It means that our estimate is even too pessimistic. The reason is that in the estimation we take the minimal value of the magnetic field in the considered region but in fact atoms feel stronger field.
Cold atoms in an artificial magnetic field. We would like to analyze now if the presence of artificial magnetic fields induced by evanescent waves can be observed in experiments with cold atomic gases at temperatures higher than the critical value for the quantum degeneracy. We will deal with velocities of atoms of the order of a meter per second which are much greater than in the case of ultracold gases. To fulfill the adiabaticity condition it is necessary to apply a sufficiently large detuning D. This in turn implies a large Rabi frequency because in order to have a significant artificial magnetic field the parameter s in Eq. (9) must be at least of the order of unity.
Such conditions can be difficult to fulfill for ytterbium atoms. Therefore, in the present section we consider a different arrangement for creation of artificial magnetic fields. Let us consider two degenerate internal states j1ae and j2ae that belong to the ground-state manifold of, e.g., 87 Rb atoms and two laser beams, characterized by the Rabi frequencies k 1 (x, z) and k 2 (x) and the wave vectors k 1 < k 2 ; k 0 , which can couple these states to an excited state. We assume large detuning of the laser radiations from the resonant transition and c~k 1 j j 2 z k 2 j j 2 À Á D 2 =1. The first beam strikes a surface of a prism with an incident angle h greater but very close to the critical angle h 0 for the total internal reflection and produces an evanescent wave, see Fig. 1. The other beam propagates in the vacuum along the prism surface, i.e. with the wave vector k 2~{ k 0x . The laser beams can induce Raman transitions between the two ground states. We assume that atoms follow adiabatically the dressed state where f~{k Ã and it is significant if the parameters is of order of unity or greater. We consider a cloud of atoms prepared in a magneto-optical trap 1 mm above the horizontally oriented prism surface (cf. Fig. 1), at temperature T 5 10 mK -in Fig. 1, the gravitational force is oriented along z axis and points downwards. We assume that after the release of atoms from the trap they are accelerated to the average velocity v~{1x m=s. With such initial conditions the cloud of atoms expands in the gravitational field and falls down on the surface of the prism. Artificial magnetic field is created by an evanescent wave which penetrates space above the prism. The dressed state (19) uncouples from an atomic excited state and there is no light shift of this level. The geometrical scalar potential W is too weak to reflect the falling atoms. Thus, an atom does not hit the dielectric surface only if the artificial magnetic field is able to bend its trajectory due to the artificial Lorentz force. Figure 7 shows density of atoms 35 ms after the release from the trap obtained in classical trajectory simulations. If the artificial magnetic field was not present, such a time would be sufficient for all atoms to fall down and strike the dielectric surface. Due to the dispersion of initial atomic velocities (i.e. Dv < 0.03 m/s for T 5 10 mK), many atoms hit the prism surface and are lost. However, some fraction of atoms (i.e. those which enter the region of the magnetic field with a sufficiently small angle with respect to the region boundary) are able to bounce repeatedly on the prism surface due to the presence of the synthetic magnetic field. Figure 7 shows the atomic density at the moment of time corresponding to the turning point of atoms. At t 5 35 ms atoms reach the highest positions, return, bounce off the dielectric surface and again reach the highest position at t < 70 ms and so on.
In the example illustrated in Fig. 7 we have assumed that after the atomic trap is turned off, atoms are accelerated to the average velocity v~{1x m=s. It is necessary in order to ensure that atoms will enter the region of the artificial magnetic field with small angles with respect to the region boundary. Then, the Lorentz force is able to bend atomic trajectories and reflect atoms from the prism surface. Such an additional acceleration can be omitted in an experiment if the surface of the prism is not oriented horizontally but forms angle of about 0.1 rad with the gravitational force vector. Then, the initial acceleration is not necessary because the gravitational field accelerates atoms to the suitable velocities.

Discussion
We have considered atoms that move slowly in the presence of an evanescent wave. The theoretical description of the adiabatic atomic motion involves the geometrical Berry phases that can be represented by vector and scalar potentials experienced by atoms. Such artificial gauge potentials are the stronger the greater gradients of the phase and amplitude of an external electromagnetic field 4 . An evanescent wave is a good candidate for the realization of the gauge fields due to an exponential decay of its amplitude.
We analyze two laser beam configurations that lead to synthetic magnetic fields. The first configuration assumes creation of an evanescent wave by means of laser radiation which is nearly resonant with an electronic transition. This method can be applied to atoms with a large radiative lifetime. If a Bose-Einstein condensate is placed in such an evanescent wave, the synthetic magnetic field can induce vortices in the condensate. In order to create a large number of vortices the angle of incidence has to be very close to the critical angle for the total internal reflection. Then, realistic profile of a laser beam has to be taken into account. We show that an evanescent wave with the Gaussian profile with experimentally attainable parameters is able to create an unidirectional magnetic field.
We also analyze artificial magnetic fields induced by two laser beams. The first laser radiation creates an evanescent wave while the other beam propagates in the vacuum and together with the first one can induce Raman transitions. Atoms prepared in a dressed state which uncouples from an excited state experience a synthetic magnetic field and are not affected by a radiative decay. We illustrate this method with a system of a cold atomic gas at temperature higher than the critical value for the quantum degeneracy. The presence of the artificial magnetic field influences motion of an atomic cloud that is released from a magneto-optical trap and falls down on a dielectric surface. In such a way one can realize a new type of the atomic mirror that is based on the artificial Lorentz force and this effect can be measured experimentally.
After the present article appeared in the arXiv basis, an article on a similar topic was submitted to arXiv and recently published 36 . Rb atoms, at T 5 10 mK, are prepared in a magneto-optical trap (v trap 5 2p 3 100 Hz) at z 5 1 mm above the surface of the prism and at x 5 0. It is assumed that after the trap is turned off, atoms are accelerated to the average velocity v~{1x m=s. Next, the cloud of atoms moves in the gravitational field and falls down on the dielectric surface. Atoms feel an artificial magnetic field induced by two red-detuned (D 5 10 GHz, l 5 795.5 nm) laser beams with the ratio of the Rabi frequencys~5. The first beam propagates in the prism and strikes its surface at the incident angle h 5 h 0 1 3 ? 10 24 rad and creates an evanescent wave. The other one propagates in the vacuum along the x direction. We assume that atoms are prepared in the dressed state (19) which uncouples from an atomic excited state. The only scalar potentials experienced by atoms correspond to the geometrical potential W and the gravitational one. The former creates too weak force (smaller than 0.2 mg) to overcome the gravitational field and is not able to reflect atoms. Note that in the absence of the artificial magnetic field all atoms are lost at the prism surface after evolution time shorter than 35 ms. If the artificial field is on, 5% of atoms are reflected from the prism surface.