Abstract
Reorientation of skyrmion crystal (SkX) with respect to crystallographic axes is believed to be insensitive to anisotropies of fourth order in spin-orbit coupling, for which sixth order terms are considered for explanation. Here, we show that this is wrong due to an oversimplified assumption that SkX possesses hexagonal symmetry. When the deformation of SkX is taken into account, fourth order anisotropies such as exchange anisotropy and magnetocrystalline anisotropy have pinning (in this work, the word ‘pinning’ refers to the reorientation effects of intrinsic anisotropy terms) effects on SkX. In particular, we reproduce some experiments of MnSi and Fe1−xCoxSi by considering the effect of fourth order magnetocrystalline anisotropy alone. We reproduce the rotation of SkX in Cu2OSeO3 by considering the combined effects of the exchange and magnetocrystalline anisotropies. And we use the exchange anisotropy to explain the reorientation of SkX in VOSe2O5.
1 Introduction
Helimagnets have attracted extensive interest since the first observation of magnetic skyrmions in 2009 [1]. Magnetic skyrmions in helimagnets are nontrivial spin textures, in which the spins point in all of the directions wrapping a sphere. Their topological protection [2] and facile current driven motion [3, 4] make them possible to be applied in novel spintronic and information storage devices [5, 6].
In helimagnets such as MnSi, Fe1−xCoxSi, and Cu2OSeO3, the ferromagnetic exchange interaction (for Cu2OSeO3, the exchange interaction consists of ferromagnetic and antiferromagnetic interactions, but the field-induced ground state is closer to ferromagnetic than antiferromagnetic [7, 8]) and the Dzyaloshinsky-Moriya interaction (DMI), which arises due to the broken of space inversion symmetry [9, 10], dominate the free energy when studying bulk material free from any external magnetic field. The former favors parallel spin alignment, while the latter favors the twist of the spins. They compete with each other and result in SkX at appropriate magnetic field just below the Curie temperature [1, 11–14]. In experiments, when the magnetic field is along directions with high symmetry, such as the [001], [111] and [110] directions, the wave vectors of SkX are orientated with respect to the crystallographic axes [1, 12, 15, 16]. This indicates the existence of anisotropy energy. The anisotropies of fourth order in spin-orbit coupling, such as the exchange anisotropy and fourth order magnetocrystalline anisotropy, are widely used to explain the pinning of helical phase, the transition from helical to conical phase and the appearance of tilted conical phase [1, 17–20]. However, according to the perturbation theory [1, 12, 16, 21], which treats the anisotropies perturbatively and approximates SkX by a triple-Q structure with three equivalent wave vectors forming a regular triangle, they are insensitive to the pinning of SkX. As a consequence, anisotropies with higher order are proposed. In our opinion, ignoring the deformation of SkX is oversimplified, because many experiments show that the structure of SkX is sensitive to anisotropy of the system which destroys its hexagonal symmetry [22–24].
In this work, we study the pinning effects of the exchange anisotropies and the fourth order magnetocrystalline anisotropy on deformable SkX. We apply a rescaled free-energy-density model for T point group and describe Bloch SkX by a three-order Fourier expansion with deformation-related degrees of freedom. Firstly, we study four anisotropies (three types of exchange anisotropies and a fourth order magnetocrystalline anisotropy in helimagnets with T symmetry) separately. It is found that they have different pinning effects on SkX. Then, by plotting the deformation-related parameters as functions of one exchange anisotropy, we figure out that the deformation of SkX is characterized by the change of amplitudes, lengths and azimuth angles of wave vectors. Next, we compare our results with some experiments, the fourth order magnetocrystalline anisotropy may explain the pinning of SkX in MnSi [1, 25, 26] and Fe1−xCoxSi [12, 27–29]. To reproduce the rotation of SkX in Cu2OSeO3 [16], we consider both the exchange and magnetocrystalline anisotropy, and find that at certain conditions rotation of SkX occurs when temperature or magnetic field changes. Lastly, we expand our model so that it is applicable to helimagnets hosting Néel SkX. It is found that exchange anisotropy has pinning effects on Néel SkX in helimagnets but not in or helimagnets.
2 Model
Based on the continuum spin model established by Bak and Jensen [17], we write the rescaled free-energy density [30] for helimagnets with the symmetry of T point group in the following form:
Here, is the rescaled magnetization. The first two terms in Eq. (1) represent the ferromagnetic exchange interaction and the DMI, respectively. The third term is the Zeeman energy under the rescaled magnetic field . is the Landau expansion with the rescaled temperature t, it consists of the second and the fourth order terms. The last term is the anisotropy energy. In this work, we consider only the exchange anisotropy and the fourth order magnetocrystalline anisotropy, and we express aswhere , and are the coefficients of exchange anisotropy, is the coefficient of magnetocrystalline anisotropy, and represent and , respectively.
For bulk B20 materials, the skyrmion plane rotates with respect to the applied magnetic field. To describe the configuration of SkX under magnetic field with different direction, we should choose an appropriate cartesian coordinates system O- in which the magnetic field is along the axis. Let the azimuthal and polar angles that characterize the magnetic field be θ and ψ, respectively. We rotate O- counterclockwise about the axis by angle θ, and get a new cartesian coordinates system O-. We then perform a second rotation, this time about the axis by angle ψ, and we get the final cartesian coordinates system O- [Figure 1]. In terms of orthogonal matrices, the product of the two operations can be written as . Due to the relation , we haveIn the cartesian coordinates system O-, we apply the n-order Fourier decomposition to describe the magnetization texture of SkX [31], Here, is the average magnetization over the entire SkX, and is the number of nth order waves. The ith order waves are characterized by their wave vectors and polarizations . In the presence of anisotropy energy, SkX with hexagonal symmetry will go through deformation, and the deformation-related parameters are introduced through the following equationIn reciprocal space, and are the normal strains; and reflect the shear deformation and rotation of the plane spanned by , respectively. are the undeformed wave vectors, they all can be expressed as a linear combination of and (without loss of generality, for hexagonal SkX, we set , ). As to , we decompose them along the basis vectors , , and (for the chosen of the orthogonal basis, see Ref. [32]), and we havewhere are the complex coefficients.
FIGURE 1
According to Eqs. (1) and (2), the free energy density is a functional of and ( denotes ), i.e., . Applying the following coordinate transformationthe free-energy density can be rewritten as, after averaging over a magnetic unit cell
At certain temperature t, magnetic field , rotation angles θ and ψ, exchange and magnetocrystalline anisotropies , , and , the parameters describing SkX are calculated via minimization of Eq. (9). In this work, we set the order of Fourier expansion .
Our analytical method can only deal with periodic magnetization structure. For the cases where the periodicity of skyrmions is broken, e.g., the thermal-induced disorder or the pinning from impurities, the review [33] and references therein are good to refer to.
3 Results and Discussion
We first investigate the pinning effects of anisotropies , , and on Bloch SkX, separately. The value of θ is ; thus, and correspond to the directions [001], [111] and [110], respectively. The temperature and the magnetic field are set to be and (in the O- coordinate system) so that SkX exists as a stable or metastable state. The thermodynamic parameters for MnSi [34] and Cu2OSeO3 [7] are available. Using these parameters, we have for MnSi and for Cu2OSeO3, these points are near the skyrmion stable region in the magnetic field-temperature phase diagram. The anisotropy coefficients of helimagnets are hard to get in experiments. We only find the relative exchange anisotropy for GaV4O8, which is about [35]. In this work, the values used for the anisotropy coefficients are . We think, to some extent, the values are within a realistic range.
We change the rotation-related parameter , minimize the free energy density and then plot ω as a function of , the angle between the wave vector and the axis, in Figure 2. Figures 2A,B show the effects of exchange anisotropy on SkX. For , a negative (Figure 2A) prefers a wave vector along the [100] or [010] direction; while a positive (Figure 2B) prefers a wave vector along the [110] or [10] direction. For and [110], ω reaches its minimum at and , respectively (Figures 2A,B), i.e., both negative and positive prefer a wave vector along the [110] direction for , and along the [001] direction for . Figures 1C–F show the effects of and on SkX, respectively. For , a , no matter what its sign is, pins a wave vector of SkX along the [010] direction; while a pins a wave vector of SkX along the [100] direction. For , is between and or between and , meaning that no wave vector is along any direction with high symmetry. For , both and pin a wave vector of SkX along the [001] direction. Figures 2G,H show the effects of on SkX. For and [111] , has the same the pinning effects on SkX as ; while for , is different from , it results in a wave vector along the [110] direction.
FIGURE 2
In Figure 2, we give the values of , difference between the maximum and minimum of free energy. They are much smaller than ω (about −0.1), about for and , and about - for other cases. The strength of -induced anisotropy in (001) plane is obviously smaller than that in (111) and (110) planes. Comparing the energy curves for and for , we find that they have the same and are symmetric about . The similarity between and can be inferred from their energy formula in Eq. 2, which are related by the coordinate transformation . It should be noticed that for , the periodicity of is for and , and is for and . This can be explained by symmetry analysis. The and terms in Eq. (2) have a higher symmetry than T point group, they are invariant with respect to fourfold rotations around the axes. and terms in Eq. (2) have lower symmetry and are invariant with respect ot twofold rotations around the axes, meaning the broken of the equivalence between [100] and [010].
SkX is treated as a deformable structure. To reveal how anisotropy energy deforms SkX, we take as an example and plot some deformation-related parameters as functions of in Figure 3. It can be found that for nonzero , 1) the wave amplitudes (Figure 3A), 2) the wave lengths of and are not equal to each other (Figure 3B), and 3) the angle between and deviates from (Figure 3C). We conclude that anisotropy energy breaks the hexagonal symmetry of SkX by changing the amplitudes of, the lengths of, and the angles between the wave vectors. In many of the small-angle neutron scattering (SANS) experiments (1; 27; 36), the observed Bragg spots have different intensities, this might be explained by our calculation. By energy minimization, we find that the dominant coefficients are , and , which represent the wave amplitudes of the first order waves with vectors , and . Their ratios reflect the relative intensities of the first-order Bragg spots. In the inset of the Figure 3A, two Bragg spots are brighter than the other four, because .
FIGURE 3
We now compare our results with some experiments. The SANS experiments of [12, 27–29] show that for and [110] directions, two of the six scattering spots are aligned with the [10] axis; for , two sets of six scattering spots are observed, one is aligned with one the [100] direction, the other one the [010] direction. This is compatible with the results shown in Figure 2G. Therefore, a negative may explain the pinning of SkX in . Different from , MnSi [26] is observed to have a wave vector along the [110] direction for . This may be explained by a positive (Figure 2H. We should point out that at zero magnetic field, a negative (positive) prefers the () directions for the helical state, which is indeed the case for (MnSi) [1, 28, 29, 37, 38]. In the work [26], two kinds of sixth order magnetocrystalline anisotropies and are thought to be responsible for the pinning of SkX in MnSi. However, this is contrary to other works [12, 21] which show that the second sixth order magnetocrystalline anisotropy determines the reorientation of SkX for and it pins SkX with a wave vector along the [010] or [100] direction depending on the sign of its coefficient. The SANS experiments of MnSi in Ref. [26] can not be explained by the sixth order magnetocrystalline anisotropies.
Cu2OSeO3 is another helimagnet hosting SkX, a peculiar experimental phenomenon about it is that for , SkX is reorientated with a wave vector along the [10] or [001] direction depending on the temperature and magnetic field conditions [16]. To explain this, we should consider the exchange and fourth order magnetocrystalline anisotropies at the same time. As a first step, we determine the signs of anisotropies , , and according to the fact that [100] is an easy axis for the helical state at zero field [20, 39, 40]. Then we confirm by Figure 4 that for , a “dominant” magnetocrystalline anisotropy pins SkX with a wave vector along the [001] direction, while a “dominant” exchange anisotropy pins SkX with a wave vector along the [10] direction. Here, the word “dominant” depends on the type of exchange anisotropy considered. When , the “dominant” anisotropy is for (Figures 4A,B), and is for or (Figures 4C–F). Lastly, we take the temperature and magnetic field into account and try to reproduce the rotation of SkX in Cu2OSeO3. The anisotropies we considered are and , and their values are 0.005 and −0.1, respectively, the same as that for Figure 4A. We fix the magnetic field and change the temperature from 0.5 to 0.8 (Figure 5A). It is found that the angle , for which the free energy reaches its minimum, drops suddenly from to at . Then we fix the temperature and change the magnetic field from 0.15 to 0.27 [Figure 5B]. Similar phenomenon is observed, drops from to at . Our results agree with the experiments of Cu2OSeO3 [16].
FIGURE 4
FIGURE 5
In another published work [41], we explain the electric-field-induced continuous rotation of SkX [42] by extending the present model. Unlike a previous theory [42] which explains the phenomenon by considering both the fourth and sixth order magnetocrystalline anisotropies, we find that a combination of fourth order exchange anisotropies and magnetocrystalline anisotropies dominates the phenomena. This is because the theory used in [42] obtains a positive coefficient of the fourth order magnetocrystalline anisotropy which is inconsistent with other experiments [20, 39, 40], while our model obtains a negative .
In polar magnets with symmetry, the DMI and the exchange anisotropy are different from that in Eq. (1). By applying the symmetry analysis, we derive the DMI: (in this case, the Néel SkX is stabilized, and the basis vectors in Eq. (6) are chosen to be , , and [32]), and the exchange anisotropy:Here, we have ignored the terms due to the fact that in polar magnets, SkX plane is perpendicular to the n-fold axis no matter what direction the applied magnetic field is along [43, 44]. The term with coefficient is rotationally symmetric, and it has no pinning effects on SkX.
For and point groups, is zero. As a result, the orientation of the wave vector of SkX is insensitive to the exchange anisotropy. In Ref. [44], based on a discrete model and Monte Carlo simulations, the authors attribute the pinning of Néel SkX in polar magnet GaV4Se8 to the Dzyaloshinskii-Moriya vectors. However, according to the continuum model, the DMI possesses rotational symmetry and has no pinning effects on Néel SkX. In our opinion, this contradiction is because the continuum model ignores higher order DMI terms which emerge during the process of transforming the discrete model to the continuum model. These higher order DMI terms possess lower symmetry or and might reorientate Néel SkX.
For point groups, we have . The term with coefficient will deform and thus reorientate the Néel SkX. Because it possesses symmetry, which is different from the symmetry possessed by the undeformed SkX. To study the pinning effects of on Néel SkX, we plot ω as a function of for 1) positive and 2) negative in Figure 6. It is found that a positive prefers a wave vector along the [110] or [10] direction, and a negative prefers a wave vector along the [100] or [00] direction. In experiments, very few helimagnets hosting Néel SkX have been found. VOSe2O5 [45] is one of these helimagnets, in which the Néel SkX is orientated with a wave vector along the [100] or [00] direction. In previous theories, less attention has been paid to the reorientation of Néel SkX in helimagnets. Here, a negative gives a possible explanation for the SkX-reorientation-related phenomena in VOSe2O5.
FIGURE 6
4 Conclusion
In conclusion, the exchange and fourth order magnetocrystalline anisotropies deform SkX by changing the amplitudes of, the lengths of, and the angles between wave vectors and thus show pinning effects on SkX. The results of magnetocrystalline anisotropy [exchange anisotropy] may explain some experiments of MnSi and Fe1−xCoxSi [VOSe2O5]. By considering the exchange and magnetocrystalline anisotropies at the same time, the rotation of SkX in Cu2OSeO3 is reproduced.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author contributions
XW and YH conceived the idea. XW finished the analytical deduction, and performed all the calculations. XW, YH, ZH, and BW discussed the results for revision and co-wrote the manuscript.
Funding
This work was supported by the NSFC (National Natural Science Foundation of China) through fund Nos. 11772360, 11832019, 11572355 and 51901081, the National Key Research and Development Program of China (Grant No. 2020YFA0309300), the Guangdong Basic and Applied Basic Research Foundation (Grant No. 2019A1515012016), and the Pearl River Nova Program of Guangzhou (Grant No. 201806010134).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1.
MühlbauerSBinzBJonietzFPfleidererCRoschANeubauerAet alSkyrmion Lattice in a Chiral Magnet. Science (2009) 323:915–9. 10.1126/science.1166767
2.
NagaosaNTokuraY. Topological Properties and Dynamics of Magnetic Skyrmions. Nat Nanotech (2013) 8:899–911. 10.1038/nnano.2013.243
3.
JonietzFMühlbauerSPfleidererCNeubauerAMünzerWBauerAet alSpin Transfer Torques in Mnsi at Ultralow Current Densities. Science (2010) 330:1648–51. 10.1126/science.1195709
4.
YuXZKanazawaNZhangWZNagaiTHaraTKimotoKet alSkyrmion Flow Near Room Temperature in an Ultralow Current Density. Nat Commun (2012) 3:988. 10.1038/ncomms1990
5.
ZhangXEzawaMZhouY. Magnetic Skyrmion Logic gates: Conversion, Duplication and Merging of Skyrmions. Sci Rep (2015a) 5:9400. 10.1038/srep09400
6.
ZhangXZhaoGPFangohrHLiuJPXiaWXXiaJet alSkyrmion-skyrmion and Skyrmion-Edge Repulsions in Skyrmion-Based Racetrack Memory. Sci Rep (2015b) 5:7643. 10.1038/srep07643
7.
JansonORousochatzakisITsirlinAABelesiMLeonovAARößlerUKet alThe Quantum Nature of Skyrmions and Half-Skyrmions in Cu2oseo3. Nat Commun (2014) 5:5376–11. 10.1038/ncomms6376
8.
YangJHLiZLLuXZWhangboM-HWeiS-HGongXGet alStrong Dzyaloshinskii-Moriya Interaction and Origin of Ferroelectricity inCu2OSeO3. Phys Rev Lett (2012) 109:107203. 10.1103/physrevlett.109.107203
9.
DzyaloshinskyI. A Thermodynamic Theory of "weak" Ferromagnetism of Antiferromagnetics. J Phys Chem Sol (1958) 4:241–55. 10.1016/0022-3697(58)90076-3
10.
MoriyaT. Anisotropic Superexchange Interaction and Weak Ferromagnetism. Phys Rev (1960) 120:91–8. 10.1103/physrev.120.91
11.
WilhelmHBaenitzMSchmidtMRößlerUKLeonovAABogdanovAN. Precursor Phenomena at the Magnetic Ordering of the Cubic Helimagnet Fege. Phys Rev Lett (2011) 107:127203. 10.1103/physrevlett.107.127203
12.
MünzerWNeubauerAAdamsTMühlbauerSFranzCJonietzFet alSkyrmion Lattice in the Doped Semiconductor Fe1-Xcoxsi. Phys Rev B (2010) 81:041203. 10.1103/physrevb.81.041203
13.
RuffEWidmannSLunkenheimerPTsurkanVBordácsSKézsmárkiIet alMultiferroicity and Skyrmions Carrying Electric Polarization in Gav4s8. Sci Adv (2015) 1:e1500916. 10.1126/sciadv.1500916
14.
OmraniAWhiteJPršaKŽivkovićIBergerHMagrezAet alExploration of the Helimagnetic and Skyrmion Lattice Phase Diagram in Cu2oseo3 Using Magnetoelectric Susceptibility. Phys Rev B (2014) 89:064406. 10.1103/physrevb.89.064406
15.
MoskvinEGrigorievSDyadkinVEckerlebeHBaenitzMSchmidtMet alComplex Chiral Modulations in Fege Close to Magnetic Ordering. Phys Rev Lett (2013) 110:077207. 10.1103/PhysRevLett.110.077207
16.
SekiSKimJ-HInosovDSGeorgiiRKeimerBIshiwataSet alFormation and Rotation of Skyrmion crystal in the Chiral-Lattice Insulator Cu2oseo3. Phys Rev B (2012) 85:220406. 10.1103/physrevb.85.220406
17.
BakPJensenMH. Theory of Helical Magnetic Structures and Phase Transitions in Mnsi and Fege. J Phys C: Solid State Phys (1980) 13:L881–L885. 10.1088/0022-3719/13/31/002
18.
ParkJ-HHanJH. Zero-temperature Phases for Chiral Magnets in Three Dimensions. Phys Rev B (2011) 83:184406. 10.1103/physrevb.83.184406
19.
HalderMChaconABauerASimethWMühlbauerSBergerHet alThermodynamic Evidence of a Second Skyrmion Lattice Phase and Tilted Conical Phase in Cu2oseo3. Phys Rev B (2018) 98:144429. 10.1103/physrevb.98.144429
20.
ChaconAHeinenLHalderMBauerASimethWMühlbauerSet alObservation of Two Independent Skyrmion Phases in a Chiral Magnetic Material. Nat Phys (2018) 14:936–41. 10.1038/s41567-018-0184-y
21.
LuoYLinS-ZFobesDMLiuZBauerEDBettsJBet alAnisotropic Magnetocrystalline Coupling of the Skyrmion Lattice in Mnsi. Phys Rev B (2018) 97:104423. 10.1103/physrevb.97.104423
22.
ShibataKIwasakiJKanazawaNAizawaSTanigakiTShiraiMet alLarge Anisotropic Deformation of Skyrmions in Strained crystal. Nat Nanotech (2015) 10:589–92. 10.1038/nnano.2015.113
23.
WangCDuHZhaoXJinCTianMZhangYet alEnhanced Stability of the Magnetic Skyrmion Lattice Phase under a Tilted Magnetic Field in a Two-Dimensional Chiral Magnet. Nano Lett (2017) 17:2921–7. 10.1021/acs.nanolett.7b00135
24.
MakinoKReimJDHigashiDOkuyamaDSatoTJNambuYet alThermal Stability and Irreversibility of Skyrmion-Lattice Phases in Cu2oseo3. Phys Rev B (2017) 95:134412. 10.1103/physrevb.95.134412
25.
GrigorievSVPotapovaNMMoskvinEVDyadkinVADewhurstCMaleyevSV. Hexagonal Spin Structure of A-phase in Mnsi: Densely Packed Skyrmion Quasiparticles or Two-Dimensionally Modulated Spin Superlattice?. Jetp Lett (2014) 100:216–21. 10.1134/s0021364014150065
26.
BannenbergLJQianFDalglieshRMMartinNChaboussantGSchmidtMet alReorientations, Relaxations, Metastabilities, and Multidomains of Skyrmion Lattices. Phys Rev B (2017a) 96:184416. 10.1103/physrevb.96.184416
27.
AdamsTMühlbauerSNeubauerAMünzerWJonietzFGeorgiiRet alSkyrmion Lattice Domains in Fe1−xCoxSi. J Phys Conf Ser (2010) 200:032001. 10.1088/1742-6596/200/3/032001
28.
BannenbergLJKakuraiKQianFLelièvre-BernaEDewhurstCDOnoseYet alExtended Skyrmion Lattice Scattering and Long-Time Memory in the Chiral magnetFe1−xCoxSi. Phys Rev B (2016) 94:104406. 10.1103/physrevb.94.104406
29.
BannenbergLJKakuraiKFalusPLelièvre-BernaEDalglieshRDewhurstCDet alUniversality of the Helimagnetic Transition in Cubic Chiral Magnets: Small Angle Neutron Scattering and Neutron Spin echo Spectroscopy Studies of Fecosi. Phys Rev B (2017b) 95:144433. 10.1103/physrevb.95.144433
30.
WanXHuYWangB. Exchange-anisotropy-induced Intrinsic Distortion, Structural Transition, and Rotational Transition in Skyrmion Crystals. Phys Rev B (2018) 98:174427. 10.1103/physrevb.98.174427
31.
HuY. Wave Nature and Metastability of Emergent Crystals in Chiral Magnets. Commun Phys (2018) 1:82. 10.1038/s42005-018-0071-y
32.
HuYWanX. Thermodynamics and Elasticity of Emergent Crystals. arXiv:1905.02165 (2019).
33.
ReichhardtCReichhardtCMilosevicM. Statics and Dynamics of Skyrmions Interacting with Pinning: A Review. arXiv:2102.10464 (2021).
34.
HuYWangB. Unified Theory of Magnetoelastic Effects in B20 Chiral Magnets. New J Phys (2017) 19:123002. 10.1088/1367-2630/aa9507
35.
EhlersDStasinopoulosIKézsmárkiIFehérTTsurkanVvon NiddaH-AKet alExchange Anisotropy in the Skyrmion Host Gav4s8. J Phys Condens Matter (2016) 29:065803. 10.1088/1361-648x/aa4e7e
36.
AdamsTMühlbauerSPfleidererCJonietzFBauerANeubauerAet alLong-range Crystalline Nature of the Skyrmion Lattice in Mnsi. Phys Rev Lett (2011) 107:217206. 10.1103/physrevlett.107.217206
37.
GrigorievSMaleyevSDyadkinVMenzelDSchoenesJEckerlebeH. Principal Interactions in the Magnetic System Fe1-Xcoxsi: Magnetic Structure and Critical Temperature by Neutron Diffraction and Squid Measurements. Phys Rev B (2007) 76:092407. 10.1103/physrevb.76.092407
38.
TakedaMEndohYKakuraiKOnoseYSuzukiJTokuraY. Nematic-to-smectic Transition of Magnetic Texture in Conical State. J Phys Soc Jpn (2009) 78:093704. 10.1143/jpsj.78.093704
39.
BannenbergLJWilhelmHCubittRLabhASchmidtMPLelièvre-BernaEet alMultiple Low-Temperature Skyrmionic States in a Bulk Chiral Magnet. Npj Quan Mater. (2019) 4:11. 10.1038/s41535-019-0150-7
40.
QianFBannenbergLJWilhelmHChaboussantGDebeer-SchmittLMSchmidtMPet alNew Magnetic Phase of the Chiral Skyrmion Material Cu2oseo3. Sci Adv (2018) 4:eaat7323. 10.1126/sciadv.aat7323
41.
WanXHuYWangB. First and Second Order Rotational Transitions of Skyrmion crystal in Multiferroic Cu2oseo3 under Electric Field. Appl Phys Lett (2020) 116:182403. 10.1063/5.0003880
42.
WhiteJSPršaKHuangPOmraniAAŽivkovićIBartkowiakMet alElectric-Field-Induced Skyrmion Distortion and Giant Lattice Rotation in the Magnetoelectric InsulatorCu2OSeO3. Phys Rev Lett (2014) 113:107203. 10.1103/physrevlett.113.107203
43.
LeonovAOKézsmárkiI. Skyrmion Robustness in Noncentrosymmetric Magnets with Axial Symmetry: The Role of Anisotropy and Tilted Magnetic fields. Phys Rev B (2017) 96:214413. 10.1103/physrevb.96.214413
44.
KézsmárkiIBordácsSMildePNeuberEEngLMWhiteJSet alNéel-type Skyrmion Lattice with Confined Orientation in the Polar Magnetic Semiconductor GaV4S8. Nat Mater (2015) 14:1116–22. 10.1038/nmat4402
45.
KurumajiTNakajimaTUkleevVFeoktystovAArimaT-h.KakuraiKet alNéel-Type Skyrmion Lattice in the Tetragonal Polar Magnet VOSe2O5. Phys Rev Lett (2017) 119:237201. 10.1103/physrevlett.119.237201
Summary
Keywords
skyrmion crystal, pinning effect, exchange anisotropy, magnetocrystalline anisotropy, helimagnet
Citation
Wan X, Hu Y, Hou Z and Wang B (2021) Pinning Effects of Exchange and Magnetocrystalline Anisotropies on Skyrmion Lattice. Front. Phys. 9:684346. doi: 10.3389/fphy.2021.684346
Received
23 March 2021
Accepted
27 May 2021
Published
21 June 2021
Volume
9 - 2021
Edited by
Cynthia Reichhardt, Los Alamos National Laboratory (DOE), United States
Reviewed by
Carles Navau, Universitat Autònoma de Barcelona, Spain
Anjan Soumyanarayanan, National University of Singapore, Singapore
Charles Reichhardt, Los Alamos National Laboratory (DOE), United States
Updates
Copyright
© 2021 Wan, Hu, Hou and Wang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yangfan Hu, huyf@dgut.edu.cn; Biao Wang, wangbiao@mail.sysu.edu.cn
This article was submitted to Condensed Matter Physics, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.