sign in sign up
<<
Home , Skyrmion

arXiv:1912.05085v1 [quant-ph] 11 Dec 2019 hteatyi h eainbtenBryconnec- Berry between relation the are, is questions exactly our what So space adi- approximation[11–13]. from real derived abatic as is explained generally which is connection field Berry this emer- field, its prop- of gauge existence gent important the is An transforma- skyrmion magnetic gauge of 10]. erty local magnetization[9, a magneti- for as local tion treated of be rotation can of theoret- zation kind in Essentially, displace- this physics, a 6–8]. ical not location[4, its but di- of magnetization the ment of local rotation of the from rection come methods these 5]. All processing[4, magnetic and storage of mag- information application to potential skyrmions of the techni- for important position individual support provide cal given delete which a materials[4], and netic at create skyrmion to magnetic are Dzyaloshinskii-Moriya methods there many chiral Experimentally, usually the which interactions(DMI)[1–3]. from materials, magnetic originate in figurations ∗ † orsodn uhr [email protected] Author: Corresponding [email protected] antcsymosaetplgclsi con- topological are skyrmions Magnetic fte( the of h xc eainbtenW-agptnil n er con Berry and general potentials Wu-Yang between magnetic relation that means exact pote this The Wu-Yang connection, the Berry the that to verify portional we adiabati transformation, the gauge taking local by Wu- Moreover, and structure. skyrmion algebraic magnetic same that find we vector, magnetization ag rnfrain ytkn the taking By transformation. gauge h uYn oetaso general a of potentials Wu-Yang the Using .INTRODUCTION I. N SU su − ( antcSymosaeQaiMgei Monopoles Quasi-Magnetic are Skyrmions Magnetic ( antcMtraso h O,LnhuUiest,P.R.China, University, Lanzhou MOE, the of Materials Magnetic N N )W-agpotentials. Wu-Yang 1) nttt fTertclPyis&KyLbrtr o for Laboratory Key & Physics Theoretical of Institute atnsbler oa ae aaerzto fdensity of parametrization bases local subalgebra Cartan ) ae ..teBrycneto o h ueo ie tt i state mixed or pure the for connection Berry the i.e. case, ) iRn Ren, Ji-Rong nToDmninlMgei Materials Magnetic Two-Dimensional in ∗ h Wang, Zhi N SU lvlqatmsse r opeeyepesdby expressed completely are system quantum -level 2 atnsbler oa ai salclnormalized local a as basis local subalgebra Cartan (2) † n hte hr xssagnrlrlto for relation general a exists there whether and for approxima- tion adiabatic and transformation gauge tion, case. unu ytmwt n tt xrse sden- as expressed state operator general any sity the with for system matrices quantum rotation are tions opsto hoy1] unadZagfud( found de- Zhang and potential Duan theory[15], gauge composition Duan-Ge magnetic using of monopoles[14], tensor electromagnetic gauge-invariant agptnil[6,o h extended the of potentials[16], Yang 1) etaso hsqatmsse r h rjcinof projection the are su system quantum this po- Wu-Yang of the tentials prove can We system[20]. general a quantum of matrix density the for We parametriza- method bases tion local also. subalgebra state Cartan a mixed develop the includes operator density sigetninfor extension esting for general connection a Berry and transformations potentials Wu-Yang tween gauge local U by com- expressed are potentials pletely Wu-Yang These theory[17–19]. e u n a Wang Hao and Qu, Fei ( ( x U N ) ewl n htthe that find will We 1 antcmnpl ag ed,ie Wu- i.e. fields, gauge monopole magnetic (1) a connection flat ) ∈ ntr vlto fmgeiainas magnetization of evolution unitary c SU krini us-antcmonopole. quasi-magnetic is skyrmion SU eto sdsusdi ealfrthe for detail in discussed is nection ta fmgei krin spro- is skyrmions magnetic of ntial N agmgei oooehv the have monopole magnetic Yang ( 2 antzto nmgei skyrmions, magnetic in magnetization (2) N lvlqatmsse,ti sa inter- an is this system, quantum -level .W ilsuyteeatrlto be- relation exact the study will We ). ρ nlgu ote’ Hooft ’t the to Analogous . operator SU h egtdaverage weighted the s s and ism 2 antzto eas the because magnetization (2) ig 1 730000 ∂ µ SU ρ UU epoethat prove we , ( N † ag transforma- gauge ) n( on SU N ( N SU ) − ( )Cartan 1) N gauge ) N N SU SU -level -level N ( N (2) − ) 2 subalgebra local bases ni, which theoretically relate II. THE WU-YANG POTENTIALS OF N-LEVEL to (N − 1) topological quasi- similar to mag- QUANTUM SYSTEM netic skyrmions in SU(2), this provides us a new clue for creating topological quasi-particles. In this section, we analyze the emergent gauge fields for a general N-level quantum system by using su(N) Cartan subalgebra local bases parametrization. In general, despite U(x) adiabatically evolution For a general N-level quantum system, its quan- or not, the (N − 1) Wu-Yang potentials always exist tum states are described by density matrix ρ. The and relate to general geometric phase of pure or mixed density matrix ρ must satisfy three properties ρ† = ρ, state[21]. However, when the adiabatic approxima- Trρ = 1 and ρ ≥ 0 , where ρ ≥ 0 is called the positive tion is satisfied, by taking the adiabatic unitary evo- semidefinite condition of ρ, and it means all eigenval- lution as a local gauge transformation U(x), we will ues of ρ are nonnegative. We know any (N × N) Her- find that N Berry connections of the N-level quantum mitian matrix can be expanded in su(N) system can be expressed by U(x). The weighted av- bases T and the unit matrix I . Considering the nor- erage for these Berry connections with respect to the a N malization property Trρ = 1, the su(N) Lie algebra eigenvalues of density operator is equal to the expec- vector representation of density matrix ρ is[24, 25] tation value of flat connection with respect to the final 1 state of the adiabatic unitary evolution. An exact re- ρ = (I + 2N(N − 1)vaT ), (1) N N a lation between Wu-Yang potentials and the weighted p 2 average of Berry connections is obtained in the con- where the Lie algebra bases Ta(a = 1, 2, · · · , N − 1) dition of adiabatic approximation. We will find in are generators of SU(N) in fundamen- the case of SU(2) magnetic skyrmions, the Berry con- tal representation and the Lie algebra vector is v ≡ a nection is proportional to Wu-Yang potential exactly, v Ta ∈ su(N). We use the convention of summation and the and Wu-Yang magnetic over repeated indices. monopole have the same algebra structure. Therefore, The vector space for su(N) Lie algebra is com- the magnetic skyrmions are quasi-magnetic monopoles posed of vector bases Ta, when we use the conventions in magnetic materials, which is related to the abun- of normalization relations dant topological structure in Skyrme theory[22, 23]. 1 Tr(T T )= δ (2) a b 2 ab

This paper is organized as follows. In section II, and the commutation and anticommutation relations we study the Wu-Yang potentials of a general N-level 1 [Ta, Tb]= ifabcTc, {Ta, Tb} = δabI + dabcTc, (3) quantum system, which is based on the parametriza- N tion of density operator by su(N) Cartan subalgebra the constant in eq.(1) ensure that (v, v) = 1 for pure local bases. In section III, we discuss the of states and (v, v) < 1 for mixed states, where the defi- Wu-Yang for SU(2) case, it nat- nition of the inner product of Lie algebra vectors is urally relates to a magnetic skyrmion charge. In sec- (Ta, Tb) ≡ 2T r(TaTb). (4) tion IV, we take the su(2) Cartan subalgebra local basis as the local magnetization of magnetic material, The restricted region of v from the positive so the magnetic skyrmion is analyzed by gauge trans- semidefinite condition of ρ was analytically obtained formations. In section V, we study the exact relation in ref.[24, 25]. More specifically, since the eigenvalues between SU(N) Wu-Yang potentials and Berry con- of ρ, we denote it as an(n = 1, 2, · · · , N), are invari- nection in adiabatic approximation. ant under unitary transformations, the necessary and 3 sufficient condition for density matrix ρ, as described where the components ui(i =1, 2, · · · , N − 1) are well in ref.[20, 26], is determined by an as † ρ = Uρ U , (5) i d N 1 ui = ( ak − iai+1). (11) n N − 1 where ρd ≡ a ρn. The ρn(n =1, 2, · · · , N) is the pure r i(i + 1) k=1 X state density matrix So the density matrixp of the N-level quantum system is parameterized by using su(N) Cartan subalgebra ρn = Diag(0, · · · , 0, 1, 0, · · · , 0) (6) local bases as given in eq.(10) and eq.(11) i.e. u ≡ n i with the single 1 in the n-th diagonal. The a satisfies u ni. n n n a = 1 and 0 ≤ a ≤ 1. In the following, we Each local basis ni is a general su(N) Lie algebra n studyP the case of non-degenerate for the non-zero a , vector in parameter space, its covariant derivative is for convenience, we assume a1 = a2 = · · · = am−1 = defined as 0(m ∈ Z, 1 ≤ m ≤ N) and Dµni = ∂µni − ig[Aµ,ni], (12)

m m+1 N 0 < a < a < · · · < a ≤ 1. (7) a where g is coupling constant and Aµ ≡ AµTa is SU(N) gauge potential, it is also a Lie algebra vector. We can The U in eq.(5) is any special unitary transformation project A to the directions of local bases of Cartan U ∈ SU(N). µ subalgebra and other directions perpendicular to these For a given element U(x) ∈ SU(N), the local local bases[17–19] bases of su(N) Lie algebra are defined as

† Aµ =(Aµ,ni)ni + [[Aµ,ni],ni]. (13) ua(x)= U(x) (x), (8)

On the other hand, using the definition of ni, it can where we use x = {xµ}(µ = 1, 2, · · · , D) to repre- be easily proved that sent the coordinate of the N-level quantum system in † any D-dimensional parameter space, which can be the ∂µni = [∂µUU ,ni], (14) real space, the momentum space and other parameter which means that spaces, and we assume that the density operator is defined smoothly and single-valued in these parame- Dµni = ∂µni − ig[Kµ,ni]=0, (15) ter spaces. 1 † where Kµ ≡ ig ∂µUU and we can verify that it is just Moreover, the Cartan subalgebra of su(N) Lie the flat connection of SU(N). algebra is composed of diagonal generators Hi(i = Analogous to the ’t Hooft SU(2) gauge-invariant 1, 2, · · · , N − 1) which satisfy [Hi,Hj] = 0. In this electromagnetic tensor[14] of magnetic monopoles, vector space of su(N) Cartan subalgebra, the local Duan and Zhang use the gauge potential decompo- bases are defined as sition theory to define the extended SU(N) gauge-

† invariant electromagnetic tensor as[17–19] ni(x)= U(x)HiU (x). (9) i 1 fµν =(Fµν ,ni) − (ni, [Dµnk, Dν nk]), (16) It can be proved that ni still satisfy the property of ig

Cartan subalgebra [ni,nj ]=0. where Fµν = ∂µAν −∂ν Aµ −ig[Aµ,Aν ] is SU(N) gauge Using Cartan subalgebra local bases, the density field strength. Using eq.(13) one can prove(Appendix matrix of eq.(5) is expressed as(Appendix A) B) 1 † 1 i i i i ρ = Uρ U = (I + 2N(N − 1)u n ), (10) fµν = ∂µAν − ∂ν Aµ − (ni, [∂µnk, ∂ν nk]), (17) d N N i ig p 4

i i where Aµ ≡ (Aµ,ni). When we take Kµ as Aµ, the fµν the relation between Wu-Yang potentials and Berry in eq.(16) equals 0 because Dµni = 0and Fµν (K)=0. connection. So according to eq.(17), we obtain

i i i 1 III. THE TOPOLOGICAL CHARGE OF Kµν ≡ ∂µaν − ∂ν aµ = (ni, [∂µnk, ∂ν nk]), (18) ig 2-LEVEL QUANTUM SYSTEM where 1 ai ≡ (K ,n )= (∂ UU †,n ) , (i =1, 2, · · · , N−1) In this section, we analyze the topological charge µ µ i ig µ i (19) of a general 2-level quantum system. For SU(2) case, τ = 1 σ (a = 1, 2, 3) are generators of SU(2) Lie is named as i-th Wu-Yang potential which is used to a 2 a describe magnetic monopole[16]. group, where σa is Pauli matrix. Its Cartan subal- Now, thanks to the Cartan subalgebra local bases gebra index i only equals to 3, so the SU(2) Wu-Yang 3 1 † potential from eq.(19) becomes a = (∂µUU ,n ), parametrization of density matrix ρ, the physical µ iqe 3 meanings of these Wu-Yang potentials for a general where U ∈ SU(2) and g is substituted with an emer- N-level quantum system are clear. At first, it is a gent charge qe, which will be explained later. The remarkable feature that the Wu-Yang potentials in Cartan subalgebra indices i, k in eq.(18) also should eq.(19) are completely expressed by local gauge trans- be 3, so the Wu-Yang curvature tensor becomes formations U(x). For N-level quantum system, the 3 3 3 1 Kµν = ∂µaν − ∂ν aµ = (n3, [∂µn3, ∂ν n3]). (20) unitary evolution of density matrix in parameter space iqe is completely described by U(x) ∈ SU(N), so these Because n3 is a su(2) Lie algebra vector, so we a a Wu-Yang potentials must be produced in any unitary can expand it as n3 = n3τa, where n3 (a =1, 2, 3) are evolution path. Moreover, the parallel transport char- its components. Using the commutation relation of acter for ni defines a parallel transport of density ma- and the definition of the inner product trix state for this unitary evolution, however there are of Lie algebra,

(N − 1) Wu-Yang potentials seen as emergent gauge 3 1 a b c Kµν = (n3τa, [∂µn3τb, ∂ν n3τc]) fields that always exist in this parallel transport state. iqe 1 On the other hand, ai is the projection of SU(N) = na∂ nb ∂ nc(τ , [τ , τ ]) µ iq 3 µ 3 ν 3 a b c 1 † e flat connection Kµ ≡ ∂µUU on the i-th Cartan ig 1 a b c i = n3∂µn3∂ν n3iǫbcd(τa, τd) subalgebra local basis ni and Kµν is Wu-Yang curva- iqe ture tensor corresponding to i-th Wu-Yang potential 1 a b c = n3∂µn3∂ν n32iǫbcdT r(τaτd) i i iqe aµ. Kµν is a U(1) magnetic monopole electromagnetic 1 a b c 1 field tensor, so for N-level pure or mixed quantum sys- = n3∂µn3∂ν n32iǫbcd( δad) iqe 2 tem, its emergent gauge fields are magnetic monopole 1 a b c = εabcn3∂µn3∂ν n3. (21) electromagnetic fields, but note that we put no con- qe straints on the dimension of parameter space now. So the magnetic monopole charge G of SU(2) Wu- The definition of Wu-Yang potentials and Wu-Yang Yang potential in R3 is[15] curvature tensors provide us a method for the calcula- 1 a b c µ ν tion of emergent electromagnetic fields, it is indepen- G = εabcn3∂µn3∂ν n3dx ∧ dx , (µ, ν =1, 2). S2 2qe dent on state and gauge condition but only depend Z (22) on local gauge transformation U(x). In the follow- which has a electromagnetic field distribution corre- ing, we mainly use this method to analyze the mag- sponding to its positive or negative charge. When netic monopole property of magnetic skyrmions and we consider S2 as a single-point compactification of a 5 two-dimensional disk D2 with a uniform field on the expressed as D ∼ 2 boundary as Σ ≡ 2 ∪ {∞} = S , r a r n3( ,t)= n3( ,t)τa. (24) 1 G = ε na∂ nb ∂ ncdxµ ∧ dxν , (µ, ν =1, 2). From this point of view, we assume 2q abc 3 µ 3 ν 3 ZΣ e (23) ~m(r,t) ≡ n3(r,t), (25) Then G corresponds to a quantization magnetic flux that is to say, we describe the local unit magnetiza- in Σ. tion as a su(2) Cartan subalgebra local basis. After Furthermore, we can see that the Wu-Yang mag- replacing the local magnetization ~m with the su(2) netic monopole must have an n3 field distribution. Cartan subalgebra local basis n3, the SU(2) Wu-Yang † Considering n3 comes from Uτ3U , n3 can be viewed potential we discussed in previous section will has a as spin operator in magnetic material and its direc- real physical effect that we can detect it in magnetic tion is determined by local gauge transformation U. materials, and we can use it to explain the theoretical So far, physicists still have not found any real mag- origin of the emergent electromagnetic field and Berry netic monopole, but eq.(23) provides us one clue to phase phenomenal produced in magnetic materials. explore quasi-magnetic monopole in magnetic system. According to eq.(25) and the replacing ma with In fact, this configuration of eq.(23) appears in many a 3 n3, the Wu-Yang curvature tensor Kµν of (2 + 1)- different materials[27]. In the following section, we dimensional magnetic material can be expressed as will discuss the physical realization of quasi-magnetic 3 1 a b c Kµν = εabcm ∂µm ∂ν m , (µ, ν =1, 2). (26) monopole in (2 + 1)-dimensional magnetic material, qe in which the magnetization configuration can be pre- According to eq.(23), it is obvious that the Wu-Yang cisely manipulated in laboratory. magnetic monopole charge can be expressed as

1 µ ν G = ~m · (∂µ ~m × ∂ν ~m)dx ∧ dx , (µ, ν =1, 2). Σ 2qe Z (27) IV. MAGNETIC SKYRMIONS IN Many theoretical and experimental studies have TWO-DIMANSIONAL MATERIALS AS QUASI-MAGNETIC MONOPOLES proved[12, 27] that skyrmions in (2 + 1)-dimensional magnetic materials have the following structure

1 µ ν In this section, using the Cartan subalgebra lo- S = ~m · (∂µ ~m × ∂ν ~m)dx ∧ dx , (µ, ν =1, 2), 8π Σ cal basis method introduced in the previous sec- Z (28) tion, we discuss the SU(2) Wu-Yang magnetic where ~m is unit magnetization and the value of single monopole charge of magnetic material. In a 2- skyrmion charge S is +1 or −1. Now by comparing dimensional base manifold M of magnetic mate- eq.(27) and eq.(28), we find ~ rial, a smooth magnetization configuration M = q S = e G. (29) M~ (r,t) is a cross section of associative vector bun- 4π dle E(M, su(2),SU(2), P ) and its unit magnetization that is to say, a skyrmion number S always correspond

vector ~m(r,t) = M~ (r,t)/|M~ (r,t)| correspond to a to a Wu-Yang magnetic monopole charge G, and qe is a unit su(2) Lie algebra vector ~m(r,t)= m (r,t)τa with the emergent electric charge of magnetic skyrmions a r a r a m ( ,t)m ( ,t) = 1(a = 1, 2, 3). We know that which satisfies the non-Abelian analog of Dirac quan- suP(2) Cartan subalgebra is τ3, and its correspond- tization condition[28] i.e. eq.(29). The magnetiza- † ing local basis is n3(r,t) = U(r,t)τ3U (r,t), where tion configuration itself can be regarded as a quasi-

n3(r,t) remains a su(2) Lie algebra vector and can be magnetic monopole. 6

The magnetization configuration of magnetic explicitly skyrmions originates from the competition between ~m(r,t)= U(r,t)~m (r)U †(r,t). (31) exchange interaction and DM interaction in a external 0 magnetic field[3, 29]. As a topological protect quasi- where ~l is the rotation axis with many advantages for technological spin- tronic applications[4, 5], now magnetic skyrmion is ~ℓ =(− sin φ, cos φ, 0), (32) the active research subject of theoretical and exper- and θ(r,t), φ(r,t) are polar angle and azimuth of imental. The experiment researchings like magnetic ~m(r,t). The spin coherent states corresponding to skyrmion topological Hall effect[30] et. verified that ~m(r,t) are skyrmion has emergent electromagnetic field[12]. The- oretically, the emergent electromagnetic field of mag- cos θ |~m(r,t)i = 2 , netic skyrmion originates from a Berry phase, it is ob- ↑ iφ θ e sin 2 ! tained by an adiabatic process. However the Wu-Yang −e−iφ sin θ |~m(r,t)i = 2 . (33) potential can be expressed as a general gauge transfor- ↓ θ cos 2 ! mation, in spite of adiabatic or not, so we think that r r the SU(2) Wu-Yang potential is a theoretical origin So the gauge transformation from |~m0( )i to |~m( ,t)i of spin Berry phase, we will further explain it in the is expressed as following section. |~m(r,t)is = U(r,t)|~m0(r)is , s =↑ or ↓, (34)

1 0 V. BERRY CONNECTION IN GAUGE where |~m0(r)i↑ = and |~m0(r)i↓ = . It is 0 1 TRANSFORMATION VIEWPOINT ! ! well known that the matrix form of U is

0 3 1 2 In the preceding section, we find the relation be- 0 a U + iU iU + U U = U I + iU σa = , (35) tween SU(2) Wu-Yang monopole charge and skyrmion iU 1 − U 2 U 0 − iU 3 ! number. In this section, we connect our Cartan sub- with algebra local bases method to spin Berry connec- θ θ tion, and then we discuss the general relation between U 0 = cos , U 1 = sin sin φ, 2 2 SU(N) Wu-Yang potentials and Berry connection of θ U 2 = − sin cos φ, U 3 =0. (36) pure or mixed state. 2 In ferromagnetic materials, the spin Berry con- Moreover, we identify this gauge transformation nection is defined as the overlap of the neighbor spin eq.(30) with a unitary time-evolution operator which

coherent states[10]. Now we express the spin Berry evolve the initial state |~m0(r)i to |~m(r,t)i with a connection of a (2+1)-dimensional magnetization con- proper Hamiltonian h(t), and it is expressed as[31] figuration ~m(r,t) using local SU(2) gauge transforma- 1 t ′ ′ i R0 dt h(t ) tion, and discuss the relation between Wu-Yang pon- U = T(e ), (37) tential and spin Berry connection. At first, we assume where T is the time-ordering operator. The adiabatic a uniform magnetization configuration ~m (r) = τ , 0 3 approximation condition in time-evolution process de- so the local gauge transformation related ~m (r) to 0 fine a process without jumping from a given eigenvalue ~m(r,t) is to another eigenvalue, that is θ θ θ U(r,t) = exp(−i ~σ · ~ℓ) = cos I − i sin ~σ · ℓ,~ (30) r r′ ′ ′ 2 2 2 sh~m( ,t)|~m( ,t )is′ =0 , for all s 6= s . (38) 7

Under this approximation conditon, the spin Berry calculated as connection[10, 32] is calculated as

3 aµ =(Kµ,n3)

1 † † = (∂µUU , Uτ3U ) iqe r r′ ′ Aµ↑ = −i↑h~m( ,t)|~m( ,t )i↑ 2 † † = Tr(∂µUU Uτ3U ) † ′ iqe = −i↑h~m0(r)|U (r,t)∂µU(r,t)|~m0(r )i↑ 2 † 0 3 0 3 = Tr(∂µUτ3U ) = −i(∂µ(U + iU )(U − iU ) iqe 1 2 1 2 0 3 1 2 +∂µ(iU − U )(−iU − U )) 1 ∂µ(U + iU ) ∂µ(iU + U ) = Tr 2 θ iq 1 2 0 3 = sin ∂ φ e " ∂µ(iU − U ) ∂µ(U − iU ) ! 2 µ 1 U 0 − iU 3 −iU 1 − U 2 = (1 − cos θ)∂µφ. (39) 2 iU 1 − U 2 −U 0 − iU 3 !#

1 0 3 0 3 = [∂µ(U + iU )(U − iU ) iqe 1 2 1 2 +∂µ(iU + U )(iU − U ) 1 2 1 2 and +∂µ(iU − U )(−iU − U ) 0 3 0 3 +∂µ(U − iU )(−U − iU )]

2 2 θ = sin ∂µφ qe 2 1 ′ ′ = (1 − cos θ)∂ φ. (42) Aµ↓ = −i↓h~m(r,t)|~m(r ,t )i↓ µ qe † ′ = −i↓h~m0(r)|U (r,t)∂µU(r,t)|~m0(r )i↓ = i(∂ (iU 1 + U 2)(iU 1 − U 2) µ So obviously 0 3 0 3 +∂µ(U − iU )(−U − iU ))

2 θ = − sin ∂µφ 3 2 2 3 2 aµ = Aµ↑ = − Aµ↓ and Aµ = qeaµτ3. (43) 1 qe qe = (cos θ − 1)∂ φ. (40) 2 µ So the spin Berry connection of a smooth magneti- zation configuration is proportional to Wu-Yang po- 2 tential with a factor ± qe . Particularly, to a static So we define a SU(2) Berry connection matrix as skyrmion ~m(r), a real space spin Berry connection serves as emergent gauge field which originates from SU(2) Wu-Yang potential exactly.

Aµ ≡ Diag(Aµ↑, Aµ↓). (41) Now we study the SU(N) Berry connection and Wu-Yang poteneials in adiabatic unitary evolution. Similar to the case of SU(2) magnetization, the uni- form field in parameter space that we need is[20, 21]

On the other hand, for the same gauge transfor- n mation U(r,t), the SU(2) Wu-Yang pontential can be u0(x) ≡ ρd = a ρn. (44) 8

There are N orthogonal eigenstates for u0(x) to the final density matrix state in adiabatic unitary evolution. 1 0 Comparing eq.(52) with eq.(19) of SU(N) Wu-  0   1  |u0(x)i1 = . , |u0(x)i2 = . , · · · , Yang potentials, using eq.(10), the exact relation be-  .   .      tween Berry connection of weighted average for pure      0   0      or mixed state and Wu-Yang potentials is obtained     0 0 g . Aµ = (Kµ, ρ)  .   0  2 |u0(x)iN−1 = , |u0(x)iN = . . (45) g i  1   .  = (Kµ, (IN + 2N(N − 1)u ni))     2N      0   1  g 2(N − 1) p     = ui(K ,n )     2 N µ i For an explicit U(x,t) ∈ SU(N), the unitary evolution r g 2(N − 1) path is = uiai , (53) 2 N µ † r u(x,t)= U(x,t)u0(x)U (x,t), (46) where the ui is expressed as an by eq.(11). In the and the orthogonal eigenstates for u(x,t) are case of SU(2) magnetic skyrmion, N = 2 and the unit magnetization vector correspond to pure state i.e. |u(x,t)in = U(x,t)|u0(x)in. (47) a1 = 1, a2 = 0 or a1 = 0, a2 = 1, so using eq.(11) we Under the adiabatic approximation conditon eq.(38), obtain u1 = a1 =1 or u1 = a2 = −1, that is to say the Berry connections for SU(N) are calculated as g g A = a3 or A = − a3 , (54) µ 2 µ µ 2 µ ′ ′ Aµn = −inhm(x,t)|m(x ,t )in that is exactly the eq.(43). † ′ = −inhm0(x)|U (x,t)∂µU(x,t)|m0(x )in To summarize the discussion above, the methods i † = − (U (x,t)∂µU(x,t), ρn). (48) of creating or deleting magnetic skyrmion in magnetic 2 material[4] are equivalent to creating or deleting quasi- Now the SU(N) Berry connection matrix can be de- magnetic monopole by taking SU(2) gauge transfor- fined as mation for the magnetization vector. For SU(N) case,

Aµ = Diag(Aµ1, Aµ2, · · · , AµN ). (49) we take the gauge transformation on ρd, then the (N−1) Wu-Yang potentials are theoretically related to The weighted average of Aµ for pure or mixed (N − 1) topological quasi-particles for pure or mixed n state is an average of Aµn with probability a [21], state, these topological quasi-particles are similar to

i † magnetic skyrmion in SU(2), so the result of SU(N) Aµ = − (U (x,t)∂µU(x,t), ρd). (50) 2 provide us a new clue for creating topological quasi- Reversing eq.(10) we obtain particles.

† ρd = U ρU, (51) VI. CONCLUSION then i A = − (∂ U(x,t)U †(x,t), ρ) Since magnetic skyrmions were observed in fer- µ 2 µ g romagnetic materials[1, 2], various novel properties of = (Kµ, ρ). (52) 2 magnetic skyrmion have been detected which mainly So the explicit meaning of the weighted average of owing to its topological characteristics. Usually, the

Aµ is the expectation value of flat connection respect emergent electromagnetic field of magnetic skyrmions 9 is considered to be a useful tool for understanding where −i is the (i +1)-th diagonal element. The ρd in these properties[12]. To obtain this emergent electro- eq.(5) is magnetic field from the magnetization configuration n 1 2 N ρd ≡ a ρn = Diag(a , a , · · · , a ), (56) of skyrmion itself, using the method of gauge trans- formation, we take the local magnetization in mag- with ρn = Diag(0, 0, · · · , 1, 0, · · · , 0), where 1 is the netic material as a su(2) Cartan subalgebra local ba- n-th diagonal element. A recursive relation between sis. Then we verified that the SU(2) Wu-Yang poten- ρn(n =1, 2, · · · , N) and Hi(i =1, 2, · · · , N − 1) is tial equals to the adiabatical Berry connection multi- 2(n − 1) 2(n + 1) ρ − ρ = H − H (57). plied by a factor ± 2 , so the same algebraic structure n+1 n n n−1 n n qe r r between magnetic skyrmion and Wu-Yang monopole After summation, the ρn is expressed as suggest that the magnetic skyrmion is quasi-magnetic 2n 2(n − 2) 2(n − 1) ρ =− H +( − )H monopole. n n − 1 n−1 n − 1 n − 2 n−2 r r r Moreover, by extending the discussion from + · · · + (1 − 2)H1 + ρ1. (58) SU(2) case to SU(N), we find that the magneti- We calculate the ρ1 as zation is naturally replaced by density operator ρ. N− 1 1 2 And then, by applying the su(N) Cartan subalgebra ρ = I + H . (59) 1 N N k(k + 1) k k=1 s parametrization method to this density operator, we X obtain (N − 1) Wu-Yang potentials of a general N- Then level quantum system. After considering the adiabatic n ρd = a ρn approximation, we obtain N Berry connections for N− 1 2i 2(i + 1) the adiabatic unitary evolution, the weighted average = [( − )(aN + aN−1 + · · · i +1 i i=1 r r for these Berry connections is exactly the expectation X i+2 2(i + 1) i+1 value of flat connection with respect to the final state +a ) − a ]Hi + ρ1 r i of adiabatic unitary evolution. Finally, we find that N− 1 2 the (N − 1) Wu-Yang potentials can be used to calcu- = [(− )(1 − (a1 + a2 + · · · + ai+1)) i(i + 1) i=1 s late the weighted average of Berry connections, it is a X 2 new representation different from ref.[21]. Analogous −(i + 1) ai+1]H + ρ i(i + 1) i 1 to SU(2) magnetic skyrmions, the result of SU(N) s N− case provide us a new clue for creating topological 1 2 = [a1 + a2 + · · · + ai i(i + 1) quasi-particles. i=1 s X i+1 −ia − 1]Hi + ρ1 N− 1 1 2 = I + [a1 + a2 + · · · N N i(i + 1) VII. APPENDIX A i=1 s i Xi+1 +a − ia ]Hi. (60)

In our conventions of normalization relation So eq.(2) and commutate and anticommutate relations † 1 i ρ = UρdU = (IN + 2N(N − 1)u ni), (61) eq.(3) for su(N) Lie algebra, the generators Hi(i = N p 1, 2, · · · , N − 1) of su(N) Cartan subalgebra are with

i 1 i N 1 k i H = Diag(1, 1, · · · , 1, −i, 0, · · · , 0),(55) u = ( a − ia +1). (62) i N − 1 2i(i + 1) r i(i + 1) k=1 X p p 10

VIII. APPENDIX B where we used (ni, [nk,L]) = 0 for any L ∈ su(N) in the second equal sign, and we used eq.(63) in the fifth We begin from the projection formula eq.(13) equal sign. According to

Aµ =(Aµ,nk) + [[Aµ,nk],nk]. (63)

The differential of Aµ is Dµnk = ∂µnk − ig[Aµ,nk], (69)

k k ∂µAν = ∂µAν nk + Aν ∂µnk + [∂µ[Aν ,nk],nk]

+[[Aν ,nk], ∂µnk]. (64) we obtain and the projection of ∂µAν to ni is

((∂µAν − ∂µAν ),ni) ([Dµnk, Dν nk],ni) (70) k k = ∂µAν (nk,ni)+ Aν (∂µnk,ni) = ([(∂µnk − ig[Aµ,nk]), (∂ν nk − ig[Aν ,nk])],ni) +([∂µ[Aν ,nk],nk],ni) + ([[Aν ,nk], ∂µnk],ni) 2 = ([∂µnk, ∂ν nk],ni)+(ig) ([Aµ,nk], [Aν ,nk],ni) k k −∂ν Aµ(nk,ni) − Aµ(∂ν nk,ni) −ig([∂µnk, [Aν ,nk]],ni) − ig([[Aµ,nk], ∂ν nk],ni). −([∂ν [Aµ,nk],nk],ni) − ([[Aµ,nk], ∂ν nk],ni)

= ([[Aν ,nk], ∂µnk],ni) − ([[Aµ,nk], ∂ν nk],ni) i i +(∂µAν − ∂ν Aµ), (65) After transposition, we obtain where we used (ni, [nk,L]) = 0 for any L ∈ su(N) in the second equal sign.

Any Lie algebra vectors A,B,C satisfy Jacobbi ([[Aν ,nk], ∂µnk],ni) − ([[Aµ,nk], ∂ν nk],ni) 1 1 identity = ([D n , D n ],n ) − ([∂ n , ∂ n ],n ) ig µ k ν k i ig µ k ν k i

[A, [B,C]] + [B, [C,A]] + [C, [A, B]]=0. (66) −ig([Aµ,nk], [Aν ,nk],ni)) 1 1 = ([D n , D n ],n ) − ([∂ n , ∂ n ],n ) We take A = [Aµ,nk], B = Aν ,C = nk, so the Jacobbi ig µ k ν k i ig µ k ν k i identity is +ig([Aµ,Aν ],ni). (71)

[[Aµ,nk], [Aν ,nk]]

= −[Aν , [nk, [Aµ,nk]]] − [nk, [[Aµ,nk],Aν ]]. (67) where we used eq.(68) in second equal sign. After Then the projection of [[Aµ,nk], [Aν ,nk]] to ni is transposition, and using eq.(65) we obtain

([[Aµ,nk], [Aν ,nk]],ni)

= −([Aν , [nk, [Aµ,nk]]],ni) − ([nk, [[Aµ,nk],Aν ]],ni) i 1 fµν ≡ (Fµν ,ni) − ([Dµnk, Dν nk],ni) = −([Aν , [nk, [Aµ,nk]]],ni) ig 1 = ([Aν , [[Aµ,nk],nk]],ni) =(∂ Ai − ∂ Ai ) − ([∂ n , ∂ n ],n ), (72) µ ν ν µ ig µ k ν k i = −([[[Aµ,nk],nk],Aν ],ni) k = ([(Aµ − Aµnk),Aν],ni)

= −([Aµ,Aν ],ni), (68) where Fµν = ∂µAν − ∂ν Aµ − ig[Aµ,Aν]. 11

[1] S. M¨uhlbauer, B. Binz, F. Jonietz, C. Pfleiderer, A. [17] Y. S. Duan and P. M. Zhang, Mod. Phys. Lett. A 17, Rosch, A. Neubauer, R. Georgii and P. B¨oni, Science 2283 (2002). 323, 915 (2009). [18] X. Liu, Y. S. Duan, and P. M. Zhang, Commun. Theor. [2] X. Z. Yu, Y. Onose, N. Kanazawa, J. H. Park, J. H. Han, Phys. 44, 371 (2005). Y. Matsui, N. Nagaosa and Y. Tokura, 465, 901 [19] Y. S. Duan and J. R. Ren, Topological gauge field theory (2010). (LanZhou University Press, 2019). [3] U. K. R¨oβler, A. N. Bogdanov and C. Pfleiderer, Nature [20] L. J. Boya and K. Dixit. Phys. Rev. A 78, 042108. 442, 797 (2006). [21] S. Chaturvedi, E. Ercolessi, G. Marmo, G. Morandi, N. [4] A. Fert, N. Reyren, and V. Cros, Nat. Rev. Mat. 2, 17031 Mukunda and R. Simon, Eur. Phys. J. C. 35, 413 (2004). (2017). [22] Y. M. Cho. Phys. Rev. Lett. 87, 252001 (2001). [5] A. Fert, V. Cros and J. Sampaio, Nat. Nanotech. 8, 152 [23] Y. M. Cho, Kyoungtae Kimm, J. H. Yoon and P. M. (2013). Zhang, Eur. Phys. J. C. 77, 88 (2017). [6] J. R. Ren, H. Wang, Z. Wang and F. Qu, Sci. China [24] G. Kimura, Phys. Lett. A 314, 339 (2003). Phys. Mech. Astron. 62, 950021 (2019). [25] M. S. Byrd and N. Khaneja, Phys. Rev. A 68, 062322 [7] X. Zhang, M. Ezawa and Y. Zhou, Sci. Rep. 5, 9400 (2003). (2015). [26] M. S. Byrd, L. J. Boya, M. Mims and E. C. G. Sudarshan, [8] N. Romming, A. Kubetzka, C. Hanneken, K. von J. Phys.: Conf. Ser. 87, 012006 (2007). Bergmann and R. Wiesendanger, Phys. Rev. Lett. 114, [27] J. H. Han, Skyrmions in condensed (Springer, 177203 (2015). 2017). [9] P. Bruno, V. K. Dugaev and M. Taillefumier, Phys. Rev. [28] Y. M. Shnir, Magnetic monopoles (Springer, 2006). Lett. 93, 096806 (2004). [29] S. Rohart and A. Thiaville, Phys. Rev. B 88, 184422. [10] G. Tatara, Physica E 106, 208 (2019). [30] A. Neubauer, C. Pfleiderer, B. Binz, A. Rosch, R. Ritz, P. [11] T. Schulz, R. Ritz, A. Bauer, M. Halder, M. Wagner, G. Niklowitz, and P. B¨oni, Phys. Rev. Lett. 102, 186602. C. Franz, C. Pfleiderer, K. Everschor, M. Garst and A. [31] A. Bohm, A. Mostafazadeh, H. Koizumi, Q. Niu and J. Rosch, Nat. Phys. 8, 301 (2012). Zwanziger, The geometric phase in quantum systems: [12] N. Nagaosa and Y. Tokura, Nat. Nanotech. 8, 899 (2013). foundations, mathematical concepts, and applications [13] K. Everschor-Sitte and K. Sitte, J. Appl. Phys 115, in molecular and condensed matter physics (Springer, 172602 (2014). 2003). [14] G. ’t Hooft, Nucl. Phys. B 79, 276 (1974). [32] X. G. Wen, Quantum field theory of many-body systems: [15] Y. S. Duan and M. L. Ge, Sci. Sinica. 11, 1072 (1979). from the origin of sound to an origin of light and [16] T. T. Wu, C. N. Yang, Phys. Rev. D 12, 3845 (1975); (Oxford University Press, 2004). Phys. Rev. D 12, 3843 (1975); Phys. Rev. D 14, 437 (1975).