arxiv: v1 [cond-mat.mes-hall] 22 Jan 2019

Transcription

1 Spin Vortices of a two-dimensional Electron in inhomogeneous Magnetic Fields Amin Naseri, Jesko Sirker Department of Physics and Astronomy, University of Manitoba, Winnipeg, Canada R3T N (Dated: January 3, 019) arxiv: v1 [cond-matmes-hall] Jan 019 We show that an inhomogeneous magnetic field can turn the spin texture of a two-dimensional single electron into a spin vortex An analytical relation between the phases of the asymptotic eigenfunction and the components of the magnetic field is derived which reveals that the net topological charge q of the vortex can take on any q Z The ground-state spin textures with q = 1, are illustrated numerically Pursuing topological phenomena has become one of the main themes in condensed matter physics 1 In this context topology refers to the existence of discrete quantum numbers (topological charges) which are related to the global topological properties of the system and are thus insensitive to small local perturbations The existence of a Berry phase is fundamental for the topological quantization of many physical quantities In the approach introduced by Michael Berry, 3 the geometrical phase of a wave function acquired under cyclic adiabatic evolution is studied which unravels the topology hidden in the parameter space of the system 4 We demonstrate below that for an electron confined to two dimensions the Berry connection is directly related to the winding number of vortices in the spin texture of the single electron in the presence of an inhomogeneous magnetic field How do the spin fields of a two-dimensional (D) electron arrange themselves if the spin is not a good quantum number? This question is the main motivation of our study In reference 5 it has been reported that the ground-state spin texture of an electron in a spin-orbit coupled quantum dot can develop a vortex with a topological charge q = 1 From a different perspective, the spin and the orbital angular momentum of electron vortex beams 6,7 propagating in a 3D space subject to electric and magnetic fields have been studied theoretically and experimentally 8,9 In a setting more relevant to our study, electromagnetic fields with topological patterns are used to manipulate the spin of electron beams 10 In the following, the spin texture of an electron strongly confined in a D space is explored in the presence of an inhomogeneous magnetic field B = (B 1, B, B 3 ) 11 where B can represent an external magnetic field or the local magnetic moments of the material which hosts the electron 1 Our analytical study shows that the spin texture of a D electron can be a vortex and its topological charge can be any integer q Z depending on the inplane magnetic field components B 1, In the main text, we focus on spatially dependent inhomogeneous fields The possible spin textures due to a homogeneous magnetic field and spin-orbit couplings are addressed analytically in the Appendix The analytical results are restricted to magnetic fields which are a continuous function of position r Ground-state spin textures of specific models with continuous inhomogeneous magnetic fields and winding numbers q = 1, are illustrated numerically Other numerical examples are also given concerning discontinuous magnetic fields not covered by our analytical results with q = 1, Although our analysis is aiming for the spin texture of an electron, it is applicable to any spin-1/ particle A generic Hamiltonian H which describes a nonrelativistic electron in a magnetic field is H = H 0 + δ B Σ, (1) where = δb 0 is the Zeeman energy and B 0 is a constant with the dimension of the magnetic field Σ = (σ 1, σ, σ 3 ) is the vector of Pauli matrices, and H 0 commutes with all Pauli matrices, [H 0, σ j ] = 0 for j = 1,, 3 If all energy scales of the model are negligible in comparison with the Zeeman energy, the spin dynamics can be described by the semi-classical Landau-Lifshitz-Gilbert equations 13 Here we study the spin texture of an electron within the quantum realm, in which eigenfunctions of H govern the physics of the electron The spin fields are defined by σ j (r) = Ψ (r)σ j Ψ(r), () and Ψ(r) is an eigenfunction of H in the coordinate representation We use a cylindrical coordinate system, r = (ρ, θ, z), in the following It is worth noting that a proper angular momentum in a 3D space is the sum of the orbital and the spin angular momentum However, the orbital angular momentum does not contribute to the in-plane spin texture of a D electron Investigating non-trivial topologies of in-plane spin textures σ(r) = (σ 1 (r), σ (r)), (3) is the purpose of our study σ(r) is an observable quantity and a unitary transformation of H leaves its topological features invariant 14 We examine two conditions under which σ 1, (r) is a constant, and therefore, has no topological pattern First, if the Hamiltonian is translationally invariant [H, p] = 0, where p = i is the canonical momentum, then Ψ can be chosen to be a simultaneous eigenfunction of H and p which leads to pσ j (r) = 0 since (pψ) = pψ (ie Ψ and Ψ have opposite momenta) Second, if there is a good spin quantum number, then by a unitary transformation σ 3 can be chosen to commute with the Hamiltonian [H, σ 3 ] = 0

2 Therefore, the in-plane spin fields are identically zero σ 1, (r) = 0 We conclude that in order to have non-trivial in-plane spin textures, the system cannot be translationally invariant and [H, Σ ˆn] 0 must be satisfied for any fixed 3D unit vector ˆn In the two-dimensional electron dot studied in the following, translational symmetry is broken by a scalar confining potential V 0 The Hamiltonian is then given by H 0 = Π m + V 0 (4) where Π = p ea is the mechanical momentum, and A is a vector potential The mass and charge of the electron are m and e, respectively We also consider examples where A and B not only depend on the in-plane coordinates (ρ, θ) but are also functions of z In those cases, V 0 will contain a potential mωzz / which strongly confines the electron in the x y plane and allows to neglect the degree of freedom perpendicular to the plane That is to say, ω z is the largest energy scale in the system The in-plane components of V 0 are taken to be harmonic potentials in the numerical studies presented in the following Its symmetries and its analytical form, however, do not play any role in our general results An eigenfunction of H, which is a spinor, can be formally written as ( ) e is 1 ψ Ψ(r) = 1, (5) e is ψ where S 1, and ψ 1, are real functions The spin and charge densities become For the 3D vector field σ 1 (r) = ψ 1 ψ cos (S S 1 ), σ (r) = ψ 1 ψ sin (S S 1 ), σ 3 (r) = ψ 1 ψ, Ψ(r) = ψ 1 + ψ (6) Σ(r) = (σ 1 (r), σ (r), σ 3 (r)), (7) it follows immediately that Σ(r) = Ψ(r) which implies that the spin is conserved globally dv Σ(r) = dv Ψ(r) = 1 (8) The normalized two-component in-plane vector field σ(r) = (cos (S S 1 ), sin (S S 1 )), (9) represents a map from R to a 1-sphere if ψ 1, 0 The single-valuedness of the wave function requires (S S 1 ) ϑ=π ϑ=0 = πq, q Z, (10) where ϑ is an angle which parametrizes a circle (ie a closed loop) in the x y plane, with an arbitrary radius and an arbitrary point as the center 14 Therefore, σ(r) locally describes a vortex with charge q for q 0 Given an eigenfunction Ψ of H, two questions are of concern for our study: First, is it possible to tailor H such that q 0 and is fixed for any arbitrary loop around a certain domain in the x y plane? Second, can vortices with any given q Z be realized? We derive an analytical relation between the phase difference (S S 1 ) in the asymptotic regime and the in-plane components of B to demonstrate that the answer to both questions is positive If the eigenfunctions of H are known, any information about the spin fields are extractable readily However, here we will describe the topology of σ(r) without a priori knowledge of the exact form of Ψ, except its asymptotic behavior Our starting point is the time-dependent Schrödinger equation HΨ = i t Ψ Continuity equations for spin fields 15,16 can be derived and by approximating them, see the following, one can find a relation between the spin fields and in-plane components of B Multiplying Ψ σ j from the left to the time-dependent Schrödinger equation and taking its imaginary part, we have J j + t σ j (r) = δ [B Σ(r)] ˆx j, (11) where ˆx j are unit vectors in a 3D Cartesian coordinate system The spin-current densities J j are J j = 1 [ Ψ σ j pψ ( pψ ) σ j Ψ ] e m m Aσ j(r) (1) The spin-continuity equations are in fact just a decomposition of the ordinary continuity equation into the Pauli matrix basis As the system is stationary, the time derivative of the densities is zero in all the equations We spend some lines on the spin-continuity equations If there is a good spin quantum number eg [H, σ 3 ] = 0, all the terms in the continuity equations for σ 1, (r) vanish identically while the continuity equation for the third component of the spin density reduces to the ordinary continuity equation J 0 + t Ψ = 0 where the current density is J 0 = 1 [ Ψ pψ ( pψ ) Ψ ] e m m A Ψ (13) If the spin is not a good quantum number, then the components of the spin torque play the role of a source/sink in the continuity equations of the spin fields Since the spin is conserved globally according to Eq (8), the net contribution of the spin-torque components must vanish, which can be verified as follows i dv B Σ(r) = Ψ [Σ, H] Ψ = Ψ (ΣH HΣ) Ψ = Ψ Σ Ψ (E E) = 0, (14) where E is an eigenenergy H Ψ = E Ψ This relation implies dv J j = 0, by taking an integral of Eq (11) in the static case

3 3 In order to achieve a local relation between the inplane spin densities and the components of the magnetic field B1, analytically, we argue that the transverse current densities Jj are negligible in the continuity equations in the regime ρ Lx,y where Lx,y are the confinement lengths of the dot The confinement of the electron implies that no net current can flow out of the dot, and hence, the transverse current densities Jj have to vanish in the distances far enough from the confinement area That is to say the Zeeman energy outweighs the kinetic energy in the regime ρ Lx,y and the kinetic term can be neglected in this limit Dismissing the kinetic term in H, Eq (11) for the third component j = 3 gives a) b) >0 c) d) 0 σ1 (r) B1, σ (r) B (15) which indicates that σ(r) orients itself along the inplane magnetic field at large distances17 The relation in Eq (15) is not valid if the magnetic field is a discontinuous function of r as σ1, (r) must be continuous functions of the position This equation also does not distinguish between σ(r) and its time-reversed partner σ(r) although it can reveal the net topological charge of σ(r) Taking the gradient from both side of Eq (15) and rearranging terms, we find (S S1 ) = B1 B B B1 B1 + B (16) Taking a loop integral, we can thus connect the winding number q defined in Eq (10) with the in-plane components of the magnetic field I 1 (B1 B B B1 ) q= dl (17) π B1 + B Here dl is an element of a closed loop in the plane with ρ Lx,y encircling the center of the confinement Eq (15) counts the net winding number in the spin texture of the dot although vortices can exist in the region ρ Lx,y which have a different winding number (eg if Ψ is not the ground state of H) and Eq (15) cannot address them Furthermore, all the eigenstates of H have the same net topological charge according to Eq (15), if Ψ meets the constraints mentioned above Instead of keeping the magnetic field fixed and moving around the center of the confinement to define the winding number, we can also keep the position r fixed at which we consider the eigenfunction Ψ(r) Ignoring again the kinetic term for ρ Lx,y fixed, a cyclic adiabatic variation of the angle θ for the magnetic field B = (B1 (ρ, θ), B (ρ, θ), 0) leads to a Berry connection Aθ = ihψ θ Ψi = 1 B1 θ B B θ B1 B1 + B (18) This shows that the winding number q of the vortex is directly related to the Berry connection under cyclic adiabatic variation of the field FIG 1: (color online) The ground-state in-plane spin texture σ(r) and the charge density Ψ of a single electron in the presence of inhomogeneous magnetic fields obtained by an exact diagonalization with confinement ωz ωx,y a) A vortex with q = generated by B The parameters are /(~ω) = 14, ωc /ω 0, bl = 07, and ωx = ωy b) A vortex with q = 1 generated by B 1 with /(~ω) = 1, ωc /ω = 00, bl = 005, and ωx = ωy c) A vortex with q = 1 generated by B+ 1 The parameters are /(~ω) = 1, ωc /ω 0, bl = 01, and ωy /ωx = 07 d) A vortex with q = generated by B+ The parameters are /(~ω) = 14, ωc /ω 0, bl = 07, and ωx = ωy We exemplify in the following magnetic fields which induce vortices with q Z where q can be calculated by Eq (17) If B is an external magnetic field, it is chosen to be divergence-free B = 0 as we are not interested to analyze the spin textures due to magnetic monopoles Besides, B must be curl-free B = 0 according to the Maxwell equations since the model is taken to be timeindependent The magnetic field n n B n = B0 b ρ (cos (nθ) sin ϕ, sin (nθ) sin ϕ, cos ϕ)(19) which becomes an in-plane field for ϕ = π/ fulfills the constraints and generates a vortex with a charge q = n where n is a positive integer n Z+, and b has the dimension of an inverse length, and satisfies blx,y 1 It is straightforward to write down a vector potential with only non-zero R z-component An = (0, 0, An ), for instance, An = dyb n x, which gives B = A Alternatively, one can again conn n sider the geometric phase the ground state wave function Ψ = ( exp(inθ) sin(ϕ/), cos(ϕ/)) acquires under a cyclic adiabatic parameter change in the magnetic field The Berry connections in this case read Aθ = ihψ θ Ψi =

4 4 a) b) > 0 and B 0 ( ν, 0, 0) y l, B II = B 0 (µ, 0, 0) l < y < l, B 0 ( ν, 0, 0) y l, (3) FIG : (color online) The ground-state in-plane spin texture σ(r) and the charge density Ψ of a single electron in the fields a) B I which generates a vortex with q = 1, and b) B II inducing a vortex with q = The parameters are ω x = ω y, ω c/ω = 00, and /( ω) = 1 in both cases n(cos ϕ 1)/ and A ϕ = 0 The Berry curvature then becomes Ω = ϕ A θ θ A ϕ = n sin(ϕ)/ If we now integrate over the whole surface of the unit sphere we obtain πn so the Chern number is indeed n The magnetic field B + n = B 0 b n ρ n (cos (nθ), sin (nθ), zf(x, y)), (0) has an in-plane vortex with a charge q = n with n Z + f(x, y) is chosen in a way to give a divergence-free field, where x, y, z are Cartesian coordinates in a 3D space However, B + n can have a non-zero curl B + n 0, which generates a time-dependent electric field To remedy this problem, we consider the Hamiltonian in the limit ω c /ω 0 while /( ω) is kept finite, where ω c = B 0 /(em) is the cyclotron frequency In this limit, the induced time-dependent electric field and the vector potential vanish, and leave a stationary system We illustrate the ground-state spin textures of some models numerically via a standard exact diagonalization The scalar potential V 0 is taken to be the potential of a 3D harmonic oscillator V 0 = 1 m ( ω xx + ω yy + ω zz ) (1) ω x,y,z are confinement frequencies which satisfy ω z ω x,y and ω = ω x ω y The unit of length is set equal to l = /(mω), and ω is chosen as the unit of the energy In Fig 1, the ground-state spin textures of an electron in the presence of the magnetic fields B ± n with n = 1, and the corresponding charge densities are presented The topological charges of the spin textures are in agreement with the analytical results, Eq (17) Next, we present the ground-state spin textures of an electron experiencing a discontinuous magnetic field, a case which is not captured by Eq (15) As examples, we consider { B 0 (η, 0, 0) y 01l, B I = () B 0 ( η, 0, 0) y < 01l, 0 where the parameters η = 0001, µ = 01 and ν = 0 are used for illustration These steps in the magnetic fields can be thought of as being due to magnetic domains in the background of the plane where the electron is confined 11,18 In addition, we apply a homogenous magnetic field perpendicular to the plane B 0 = B 0 ẑ, generated by the vector potential A 0 = (B 0 /)( y, x, 0) The ground-state spin textures in the presence of B I = B 0 +B I or B II = B 0 +B II are depicted in Fig, which shows B I induces a vortex with q = 1, and B II gives rise to a vortex with q = Although Eq (15) fails to work because of the discontinuities in the magnetic fields, it gives a hint that σ(r) needs to turn to align itself, for instance, with the magnetic field at y ± in the case of B I In these set-ups, the energy densities confined in the domains need to be of the same order of magnitude in order to favor the presence of the vortices In conclusion, our study shows that the spin texture of a single electron, which is strongly confined in a D space and is subject to an inhomogeneous magnetic field can be a spin vortex with any q Z A local relation between the spin fields and the components of the magnetic field is found in the asymptotic regime which reveals the net topological charge of the spin texture if the magnetic field does not vanish at ρ L x,y and z = 0 Additionally, numerical studies were presented to illustrate the spin vortices The quantization of the topological charge relies on the single valuedness of the eigenfunctions and does not require any symmetry of the model, and hence, has merely a quantum mechanical origin Intuitively speaking, the vector of the spin fields orients itself along the magnetic field locally to minimize the Zeeman energy and gives the electron a topological character In order to favor the presence of the vortices, the Zeeman coupling needs to be comparable with the energy unit of the system (bl) q ω Our investigation does not cover the role of σ 3 (r) and is focused on the in-plane spin texture We have also restricted ourselves to a stationary system and to the regime k B T ω where thermal fluctuations can be neglected Here k B is the Boltzmann constant and T is the temperature Exploring and classifying the spin textures induced by an arbitrary inhomogeneous magnetic field can be an interesting direction to generalize the results of our study and we leave it as an open future problem AN acknowledges fruitful discussions with Robert Stamps, Igor Proskurin and Wenchen Luo He also acknowledges support by the Canada Research Chairs through Tapash Chakraborty AN and JS acknowledge support by the Natural Sciences and Engineering Research Council (NSERC) through the Discovery Grants program

5 5 Appendix: Spin-orbit Couplings The spin texture of a confined electron in a D space subject to Rashba or Dresselhaus spin-orbit coupling (SOC) hosts a spin vortex with q = 1 5 We use Eq (17) to count the winding number of the vortices The following Hamiltonian describes a D electron in the presence of SOCs and a homogeneous magnetic field H SOC = H 0 + g 1 ([p y ea y ] σ x [p x ea x ] σ y )(A1) + g ([p y ea y ] σ y [p x ea x ] σ x ) + σ 3, where H 0 is given in Eq (4), A = (A x, A y, 0) = (B/)( y, x, 0) is the vector potential and the strength of Rashba (Dresselhaus) term is determined by g 1 (g ) Due to the broken translational symmetry, the electron has confinement lengths L x,y Therefore, the density of momentum decays to zero for distances larger than the confinement lengths ρ L x,y and can be neglected in comparison to the other terms in the Hamiltonian Neglecting the derivative terms in H SOC in the asymptotic regime, we have H SOC = H 0 eg 1 (A y σ x A x σ y ) (A) eg (A y σ y A x σ x ) + σ 3, where H 0 also does not include derivative terms The effective magnetic fields at ρ L x,y associated with the SOCs are then B R = eg 1 (A y, A x, 0), (A3) B D = eg ( A x, A y, 0) Putting these effective fields into Eq (17), the winding number for the case of Rashba SOC (ie g 1 0 and g = 0) is q R = +1, and for the case of Dresselhaus SOC (ie g 1 = 0 and g 0) is q D = 1 Since the spin densities are observables, the results are independent of the gauge of the vector potential 1 Davide Castelvecchi, Nature 547, 7-74 (017) D J Thouless, Topological Quantum Numbers in Nonrelativistic Physics (World Scientific, River Edge, NJ, 1998) 3 M V Berry, Proc R Soc London A 39, 45 (1984) 4 D Xiao, M-C Chang, and Q Niu, Rev Mod Phys 8, 1959 (010) 5 Wenchen Luo, Amin Naseri, Jesko Sirker, Tapash Chakraborty, arxiv: (018) 6 K Y Bliokh, I P Ivanov, G Guzzinati, L Clark, R Van Boxem, A Béché, R Juchtmans, M A Alonso, P Schattschneider, F Nori, JVerbeeck, Phys Rep 690 (017) 7 S M Lloyd, M Babiker, G Thirunavukkarasu, and J Yuan, Rev Mod Phys 89, (017) 8 K Y Bliokh, Y P Bliokh, S Savel ev, and F Nori, Phys Rev Lett 99, (007) 9 M Uchida and A Tonomura, Nature 464, 737 (010) 10 E Karimi, L Marrucci, V Grillo, and E Santamato, Phys Rev Lett 108, (01) 11 A Nogaret, J Phys Condens Matter, 5301 (010) 1 A C Hewson, The Kondo Problem to Heavy Fermions (Cambridge University Press, Cambridge, UK, 1993) 13 Jung Hoon Han, Skyrmions in Condensed Matter, (Springer International Publishing, 017) 14 M Nakahara, Geometry, Topology and Physics (IOP, Bristol, 1990) 15 Junren Shi, Ping Zhang, Di Xiao, and Qian Niu, Phys Rev Lett 96, (006) 16 Z An, F Q Liu, Y Lin and C Liu, Sci Rep, 388 (01) 17 For a 3D electron, one can derive B Σ(r) 0 if B j are in the same order of magnitude 18 F M Peeters and A Matulis, Phys Rev B 48, (1993)