Magnetic Field Induced Competing Phases in Spin-Orbital Entangled Kitaev Magnets

There has been a great interest in magnetic field induced quantum spin liquids in Kitaev magnets after the discovery of neutron scattering continuum and half quantized thermal Hall conductivity in the material $\alpha$-RuCl$_3$. In this work, we provide a semiclassical analysis of the relevant theoretical models on large system sizes, and compare the results to previous studies on quantum models with small system sizes. We find a series of competing magnetic orders with fairly large unit cells at intermediate magnetic fields, which are most likely missed by previous approaches. We show that quantum fluctuations are typically strong in these large unit cell orders, while their magnetic excitations may resemble a scattering continuum and give rise to a large thermal Hall conductivity. Our work provides an important basis for a thorough investigation of emergent spin liquids and competing phases in Kitaev magnets.


Introduction
Discovery of quantum spin liquids 1,2 with emergent quasiparticles has been an important subject in modern condensed matter physics. This serves as an ultimate test of our understanding of highly quantum entangled phases in interacting electron systems. Recent research has invested tremendous effort on a number of materials with strong spin-orbit coupling 3,4 , which leads to intriguing bond dependent exchange interactions between spin-orbital entangled pseudospin-1/2 moments. These studies are largely motivated by the exact solution of the Kitaev honeycomb model 5 . The Kitaev interaction is naturally present in the systems with 4d/5d transition metal elements 6 , such as honeycomb/hyperhoneycomb iridates [7][8][9] and α-RuCl 3 10 . However, other exchange interactions are present too 11 , which often lead to magnetically ordered ground states instead of the desired quantum spin liquid [12][13][14][15] . Hence much effort has been spent to suppress the magnetic orders and gain access to the possible spin liquid phases.
Over the past few years, great experimental progress has been achieved in α-RuCl 3 . At zero magnetic field, this material orders magnetically in the zigzag (ZZ) order 14,15 . Upon the application of an external field, neutron scattering experiments [16][17][18] find an intermediate window of fields before the system enters the polarized state, where sharp magnon modes are absent, but a scattering continuum appears instead. Under a [111] field (perpendicular to the honeycomb plane), the measured thermal Hall conductivity above the ordering temperature T N ≈ 7 K follows the predicted trend of itinerant Majorana fermions in the pure Kitaev model 19 . When the field is tilted away from the [111] direction by 45°and 60°, half quantized thermal Hall conductivity is observed 20 . These observations raise the hope that the paramagnetic state in the intermediate field regime may be the sought-after chiral spin liquid with Majorana edge modes.
Theoretical models for α-RuCl 3 include substantial Kitaev and symmetric anisotropic Γ interactions, both strongly dependent on the bond directions, with additional small exchanges such as the neareast neighbor Heisenberg J, the third nearest-neighbor Heisenberg J 3 , and the anisotropic Γ 21-23 , on the honeycomb lattice. Previous analyses are largely done on quantum models with small system sizes (typically a 24-site cluster) via exact diagonalization (ED) [24][25][26][27][28] or in quasi-one dimensional limit via density matrix renormalization group (DMRG) 26,27,29 , with varying degree of complexity. For example, a recent work 27 on the KΓΓ model in an external magnetic field suggests that it allows an intermediate spin liquid phase continuously connected to the pure Kitaev model between the low field ZZ order and high field polarized state.
In this article, we investigate the possible competing phases in the classical KΓΓ model under a [111] magnetic field for large system sizes. The purpose is to critically examine what kind of competing phases may be present and how these phases may be related to potential spin liquids in the quantum model. Rather surprisingly, we find a series of competing magnetic orders with large unit cells in the intermediate field regime. In particular, in the KΓΓ model with small Γ , the ground state in the zero field limit is the ZZ order, which is consistent with previous experiments and theoretical calculations. Upon increasing the field, the ZZ order is replaced by a series of magnetically ordered phases with 8, 18, 32, 50, 70 and 98-site unit cells before the system enters the polarized state (see Fig. 2). Hence the magnetic field reveals a series of competing orders, which form an intermediate region in the phase diagram. Most of these large unit cell orders had not been identified in previous works.
We compute the zero point quantum fluctuations for these magnetic orders and estimate the reduction of the size of the local moments. We find that quantum fluctuations are strong in the large unit cell orders so that the renormalized local moment is only about 50% of the full magnitude on average. The flat and dense spin wave spectra in the large unit cell orders, in particular the 70 and 98-site orders, essentially look like continua of spin excitations. Furthermore, we calculate the thermal Hall conductivity due to magnons in some of the large unit cell orders and find that it is as large as that observed experimentally at low temperatures. While strong quantum fluctuations are present and hence it is likely for the series of competing phases to turn into spin liquids in the quantum limit, it is also evident that previous theoretical studies on quantum models with small system sizes [24][25][26][27][28][29] cannot resolve many of these large unit cell orders. Therefore, in future analyses of such quantum models, it will be important to understand the role of quantum fluctuations in the large unit cell orders unveiled in the current work. Our findings demonstrate the possibility of novel and exotic ordering patterns in spin-orbital entangled Kitaev magnets, which provide an important basis for further investigations of the origin and the nature of quantum spin liquids that they may host.

Results
Model. We investigate the nearest neighbor KΓΓ model on the honeycomb lattice in a [111] magnetic field h, where K is the Kitaev interaction, Γ and Γ are off diagonal spin exchanges, (λ, µ, ν) is a cyclic permutation of (x, y, z) and the field h = h(1, 1, 1)/ √ 3. We have also assumed an isotropic g tensor. In (1), h actually carries a factor of S but for notational simplicity we will just write h in units of the Kitaev interaction, for instance h = 0.1|K | instead of h = 0.1|K |S, in the rest of this article.
In the experimentally relevant parameter regime, K < 0 and Γ > 0 are large while Γ < 0 is small. In contrast to many of the previous studies 24-30 , we investigate the classical limit of this model, that is, by treating the spins (1) as three dimensional vectors of fixed magnitude |S i | = S for all i. We use simulated annealing to determine the ground state spin configuration of the system. Details of the simulated annealing calculation can be found in the Methods section.
We explore Γ ∈ [0, 0.5], h ∈ [0, 1.2] and map out the phase diagram, as shown in Fig. 1. Apart from the extensively degenerate Kitaev limit (Γ, h) = (0, 0), we find that the vast majority of the parameter space favors particular magnetic orders. All these ordered phases, except the zigzag (ZZ) order and the ferromagnet (FM), are labeled by the number of sites contained in their respective magnetic unit cells.
In the zero field limit, the degeneracy of the Kitaev manifold is lifted as the ZZ order and a 12-site order (a 6-site order and an 18-site order) are selected at small (intermediate) Γ. These two phases have exactly the same energy at h = 0, but the ZZ order or the 6-site order is prefered once h 0. However, the 18-site order reemerges at higher fields and replaces the 6-site order as the ground state. Tracing back to the parameter region with small Γ and h, we see that the ZZ, 6-site, 12-site and 18-site orders are continuously connected to the Kitaev limit. The 6-site order (the 18-site order) was first reported in Ref. 31 (Ref. 32) and termed the X phase (the diluted star phase). At sufficiently large values of Γ and h, even larger cluster ordering patterns like the 30-site and 50-site orders are stabilized. There is also an 18-site order with C 3 symmetry, which we label by 18-C 3 to distinguish it from the previous 18-site order as they are described by different arrangements of spins on the honeycomb lattice.
Next, we set Γ = −0.02 and map out the phase diagram within the same ranges of Γ and h, as shown in Fig. 2. The addition of such a small Γ term to the KΓ model alters the phase diagram quite significantly. The degenerate manifold in the Kitaev limit and its neighborhood are replaced by the FM phase. The ZZ order is stabilized over a large portion of the parameter space at zero 33 and low fields. Once again, we find at intermediate fields several large cluster ordering patterns, a 32-site order, a 70-site order, a 98-site order and a 50-site order with C 3 symmetry which we label by 50-C 3 . Finally, the strong Γ high field regime of the phase diagram displays some similarities to the Γ = 0 case, where the same 50-site, 18-site and ZZ order are the lowest energy spin configurations, before the system becomes a FM.
Details of the magnetic orders (the real space spin configurations and the static spin structure factors, etc.) that show up in the phase diagrams Figs. 1 and 2, from the four sublattice ZZ order to the 98-site order, can be found in the Supplementary Materials. We make some qualitative observations as follows. Firstly, stronger Γ interaction stabilizes magnetic orders with larger unit cells. This is true for both zero and finite Γ . We expect that ordering patterns with even larger unit cells than those mentioned above may appear if Γ is further increased beyond 0.5. Secondly, the large unit cell orders, like the 70-site and 98-site orders, are closely competing in the parameter region where they are stabilized. The difference in energy is typically 10 −3 to 10 −4 of the energies of these orders. Thirdly, the magnetic orders can be classified into two categories, one with an inversion symmetry and the other with a three fold rotational symmetry. The ZZ order and the magnetic orders labeled by numbers fall into the former, while the magnetic orders labeled by numbers appended with -C 3 fall into the latter. More details can be found in the Supplementary Materials. Magnetization. The proposed spin model for the material α-RuCl 3 is parametrized by dominant K and Γ exchanges, with K < 0 and Γ ≈ −K/2, plus some small additional interactions like Γ , J and J 3 , where J (J 3 ) is the (third) nearest neighbor Heisenberg exchange [21][22][23] . Therefore, in the phase diagram Fig. 2 of the KΓΓ model in a [111] magnetic field, we take a cut along Γ = 0.5 and plot the magnetization m = S ·ĥ as a function of the field h, as shown in Fig. 3a. The magnetization increases monotonically with the field, and jumps at the phase transitions. The discontinuities are not very obvious at the transitions between the large unit cell orders, but are significant when the system enters to (exit from) a large unit cell order from (to) ZZ, and from ZZ to FM. This suggests the difficulty of detecting phase transitions at intermediate fields by inspecting the magnetization, if they exist at all in the quantum model. Linear spin wave theory. As a first approach to study the effect of quantum fluctuations on the classical orders, we apply the linear spin wave theory 34,35 to calculate the reduction of ordered moments in the zero temperature limit. For simplicity, we assume the same underlying magnetic orders, and do not consider how the classical phase diagram may be changed due to quantum correction to the energy because there are too many competing phases. Details of the linear spin wave calculation can be found in the Methods section. In Fig. 3b, we plot the average fraction of spins achieved in the linear spin wave theory with S = 1/2 as a function of the field, where b † i is the magnon creation operator at site i and N site is the total number of sites in the system. Blank regions indicate that the spin wave Hamiltonian is not positive definite at one or more momenta, i.e. the lowest magnon band becomes gapless. At low and intermediate fields h 1, the average reduction of ordered moments is about 50% of the full magnitude S, hinting at strong quantum fluctuations. At high fields h 1, S LSWT /S increase monotonically with h in the ZZ phase, but the spin wave solution becomes unstable in the region h ∈ [1.29, 1.37], where the system is in the FM phase with the spins not completely aligning with the [111] field (see Fig. 3a). Not only is this region likely to host quantum spin liquids, but it is also interesting from the classical aspect, which we will discuss in details later. Finally, for h > 1.37, the system enters the fully polarized state and S LSWT /S = 1 achieves saturation.
The spin wave dispersion of a (very) large unit cell order typically appears flat and dense. As an example, we show the spin wave dispersion of the 98-site order along certain high symmetry directions in the first Brillouin zone of the honeycomb lattice in Fig. 4.
Thermal Hall conductivity. We calculate the thermal Hall conductivity due to magnons [36][37][38] , Details of the calculation can be found in the Methods section. Expressing the field in (1) as h = gµ B µ 0 H, assuming the g factor g = 2.3 35,39 and the magnitude of the Kitaev interaction |K | ≈ 80 K 19,21 , µ 0 H = 12 T roughly corresponds to h = 0.23. At this field and with the parametrization We plot the thermal Hall conductivity κ xy as a function of temperature T, as shown in Fig. 5a. We show only data below T c ≈ 11 K, defined as the temperature at which S LSWT /S drops to zero, i.e. the magnetic order is destroyed by thermal fluctuations (see Fig. 5c). We find that κ xy is close to zero but slightly negative at 10 K. It gradually develops a positive value as T decreases, and peaks at 5 K before diminishing again as T −→ 0. The maximum value of κ xy /T is about 0.3 × 10 −3 W/K 2 m, which is of the same order of magnitude as the half quantized value 0.826 × 10 −3 W/K 2 m measured in Ref. 20.
We also calculate the thermal Hall conductivity for another large unit cell order, the 32-site order, at the field h = 0.82 (which would roughly correspond to µ 0 H = 43 T) and with the same parametrization, as shown in Fig. 5b. This time T c ≈ 6 K (see Fig. 5c) and κ xy is negative. Starting from zero temperature, κ xy grows in magnitude as T increases, and reaches −0.9 × 10 −3 W/Km at 6 K. The trend and the magnitude of κ xy are similar to those reported in Ref. 19 at lower fields (µ 0 H = 6, 12 and 15 T). Hence the opposite signs of κ xy may indicate the presence of different magnetic orders.
Frustrated ferromagnet. We notice that there is a window of h where the system is a FM but not fully polarized, i.e. the spins align uniformly but not in the direction of the [111] field. Such a phase is also stabilized in the high field regime at other parametrizations (K, Γ, Γ ) including the KΓ model, and the width of the window is usually larger for stronger Γ. In the following, we attempt to derive some analytical understanding of why this situation occurs.
We start from the KΓ model with K < 0, Γ > 0, and Γ = 0 in (1). Assuming a FM state, that is, S i = S for all sites i, the Hamiltonian reduces to with the matrices and N being the total number of unit cells. The Kitaev interaction becomes "isotropic" in the FM state, behaving like the Heisenberg interaction. The Γ interaction still appears quite anisotropic at this stage, but a change of basis will bring it to a simpler and more illuminating form.
In the abc basis, the Kitaev interactionH K = RH K R T = H K remains the same, while the Γ interaction assumes the form of an XXZ modelH We can then analyze (4) in the abc basis term by term. It can be shown analytically that the energy of the classical Kitaev model is K |S| 2 per unit cell 41 . Thus any FM phase will minimize the energy of the K term. On the other hand, the Γ term attains maximum (minimum) when the spin points along (lies on) the c-axis (ab-plane). The energy profile of the Γ term is shown in Fig. 6.
Suppose that the fieldh = (0, 0, h) is along the c-axis or the [111] direction. The h term wants to align the spin with the c-axis, but this will be costly in energy for the Γ term. The competition between Γ and h tilts the spin away from the c-axis. Therefore, such a FM phase can be stabilized between the fully polarized state and some other orders, e.g. ZZ and 18, in the high field regime. In contrast, if the field is along any of the in-plane directions, then all the K, Γ and h terms in (4) can be minimized simultaneously. Now let us consider the KΓΓ model. A finite Γ term acts similarly as Γ. One can easily show that, assuming a FM state, H Γ has the same structure as H Γ in (5). Thus a small Γ < 0 (Γ > 0) weakens (enhances) the effect of Γ. A similar FM but not fully polarized phase due to the presence of a large Γ interation in the JKΓJ 3 model under a [001] field was also found and discussed in Ref. 42.

Discussion
Classcially, the pure Kitaev model is extremely sensitive to an external magnetic field. It is polarized whenever the field h 0. From the result of our simulated annealing calculation, a finite Γ interaction on top of K gives rise to a multitude of ordered phases, many of them possess fairly large magnetic unit cells, at finite fields. As Γ increases, the window of these nontrivial magnetic orders becomes wider, and the system becomes polarized at greater value of h. Thus the combination of Γ and h effects like a prism that produce a rich and colorful phase diagram. Adding a small Γ term stabilizes even larger cluster magnetic orders at intermediate fields. We successfully demonstrate that the KΓΓ honeycomb model is a playground for many exotic field induced magnetic orders, not simply the zigzag (ZZ) order and the polarized state as largely perceived in the past.
We discuss the implications of these large unit cell orders. First of all, the size of the system has to be sufficiently large to host them. If the system is smaller than or incommensurate with the magnetic order, the ground state spin configuration may appear like a disordered state. This calls for a serious reconsideration of the results from quantum calculations on small systems where finite size effect can be important, such as ED on the 24-site cluster 27 and iDMRG on the cylinder geometry 29 , which report a quantum spin liquid ground state. The large unit cell magnetic orders found in this work cannot be captured by these and similar computations 24-26,28 on quantum models with small system sizes. Nevertheless, the possibility of a quantum spin liquid still exists, especially in the vicinity of the large unit cell orders where quantum fluctuations are strong. The large unit cell orders are very close in energy in the parameter region where they are stabilized. In addition, the average spin wave correction to the ordered moments in the large unit cell orders for a representative parametrization of α-RuCl 3 is found to be more than 50%. One can imagine that quantum fluctuations may melt these competing orders and promote a spin liquid state, but we will not know whether this is true until the magnetic orders are explicitly taken into account in the quantum model. If the large unit cell orders (partially) survive under quantum fluctuations, the magnon bands are typically flat and very close to each other such that they appear like the excitation continuum seen in inelastic neutron scattering experiments, which is often interpreted as fractionalized excitations in a quantum spin liquid [16][17][18]43,44 . The resulting two magnon excitations will also form a very broad continuum at low energies. Moreover, we calculate the magnon thermal Hall conductivity for two of the large unit cell orders and show that it resembles the trend and/or the magnitude as that measured in experiments 19,20 below the ordering temperature. In contrast, as computed in Ref. 35, the magnon thermal Hall conductivity in the ZZ order is in general quite small in magnitude.
We also discover the existence of a ferromagnetic (FM) but not fully polarized state at high fields in the KΓΓ model with zero or small Γ , which can be understood through the competition between the Γ and h terms. Here K < 0 and Γ > 0 are assumed. While the field always wants to orient the spins in its direction, the Γ interaction is only minimized (maximized) when the spins are all lying on the ab-plane (pointing along the c-axis). This may explain why the system is more prone to polarization when the tilting angle of the field from the [111] direction is larger. This also suggests that frustration is stronger (weaker) when the field is along (in) the c-axis (ab-plane). For instance, the simulated annealing calculation on the classical JKΓJ 3 honeycomb model in an in-plane field 31 with the parametrization (J, K, Γ, J 3 ) = (−0.035, −1, 0.5, 0.035) only yield the 6-site order (termed the X phase) at intermediate fields, between the ZZ order at low fields and the polarized state at high fields, thus leading to a relatively simple phase diagram.

Methods
Simulated annealing. The simulated annealing calculation is performed on a honeycomb lattice with L × L unit cells (or L × L × 2 sites) with periodic boundary conditions. Most of the computations are done with L = 12, 15, 20, but sometimes other L is used when the ground state spin configuration is not obvious. The procedure of simulated annealing is outlined as follows. In the beginning, we generate a totally random spin configuration on the honeycomb lattice and define a "temperature" parameter T. We randomly select a site on the honeycomb lattice and propose a random orientation for the spin on that site. Next, we calculate the difference in energy and accept the change with the probability min{1, exp(−∆E/T)}. This step is repeated for ∼ 10 7 times at a fixed T, which is then decreased gradually. Once T < T c for some critical temperature ∼ 10 −8 |K |, we update the spin at site i deterministically by aligning it in the direction of the local field 32 defined as where H i j is the three dimensional matrix that encodes the interaction between the spins at i and j.
If the sublattice structure of a magnetic order is known, we can carry out the above procedure for a small number of spins and calculate the energy to very high precision. This allows us to better determine the phase boundary between competing orders.
Linear spin wave theory. The content in this section is mainly derived from Ref. 34. For each sublattice i in the magnetic unit cell, we first choose a local coordinates system in which the spin S i aligns in the z-direction. The amounts to a change of basis characterized by the rotation matrix where θ i and φ i are the two angles parametrizing the orientation of S i in the cubic x yz coordinates, The third column of R i is precisely S i up to the factor S, while the first and second columns are chosen such that the three columns are mutually orthonormal and satisfy the right hand rule. We define S i = R iSi . Classically, we haveS i = (0, 0, S). Quantum effects on the ordered moments are introduced through spin wave excitations (magnons), where we have used the linear spin wave approximation that neglects the third and higher order terms in b i in the series expansion of (10b) and (10c). Next, we rewrite the spin Hamiltonian as RepresentingS i using (10a)-(10c), keeping only terms quadratic in b i , and performing a Fourier transform where, from now on, i denotes the position in the Bravais lattice defined by the translational symmetries of the magnetic order, s denotes the sublattice in the magnetic unit cell, and N is the total number of magnetic unit cells. We then obtain the spin wave Hamiltonian in momentum space where Ψ k = (b k,1 , . . . , b k, N , b † −k,1 , . . . , b † −k, N ) and N is the total number of sublattices in the magnetic unit cell. D k is a 2N dimensional matrix of the form where A k and B k are N dimensional matrices. To obtain the spin wave dispersion, we diagonalize D k by a Bogoliubov transformation in order to to preserve the canonical commutation relation of the bosons, where E k = diag(ε k,1 , . . . , ε k, N , ε −k,1 , . . . , ε −k, N ) and σ 3 is a diagonal matrix with the first N entries equal to 1 and the last N entries equal to −1. The average reduction of ordered moments (10a) at temperature T can be calculated from the matrix elements of the Bogoliubov transformation, Thermal Hall conductivity. We explain the various symbols that appear in the formula (3) for the calculation of the thermal Hall conductivity 37,38 . n is the magnon band index that runs from 1 to N . The function c 2 is given by where Li 2 is the dilogarithm. f BE is the Bose-Einstein distribution as defined in (17). Ω k,n is the Berry curvature defined as where σ 3 and T k are defined as in (15). For the calculation of the total volume V of the system, we use the inter-plane distance 5.72 Å between the honeycomb layers in α-RuCl 3 19,20 . The exact value of the in-plane lattice constant does not enter the calculation explicitly because, while 1/V contribute two inverse factors of it, Ω k,n contribute two factors, so they cancel out. When performing the summation over momenta in (3), we partition the first Brillouin zone (of the magnetic order) evenly such that it contains a total number of L × L k points. We check the convergence of κ xy with increasing L up to L = 800. We also ensure that the Chern number of each magnon band, where A is the total area of the system, converges to an integer with increasing L.

Data availability.
The data that support the findings of this study are available from the corresponding author upon reasonable request.

DETAILS OF THE MAGNETIC ORDERS
For each of the nontrivial magnetic orders, we show in Figs. S1-S14 the sublattice structure in the magnetic unit cell, the real space spin configuration (except for the 70-site and 98-site orders) and the static spin structure factor where r i = R i + sa/ √ 3ŷ, R i is the coordinates of the Bravais lattice, a is the lattice constant, s = 0 (s = 1) if r i belongs to the even (odd) sublattice of the honeycomb lattice and N site is the total number of sites in the system. We choose the direct lattice vectors to be a 1 = ax, a 2 = a(cos π 3x + sin π 3ŷ ), so that the reciprocal lattice vectors are . The structure factors are plotted within the second Brillouin zone, and the first Brillouin zone is enclosed by dashed lines. Whenever our discussion involves more than one of the three inequivalent M points in the momentum space, we use the notation M i , i = 1, 2, 3 to distinguish them.
We observe that the magnetic orders respect either an inversion I symmetry or a three fold rotational C 3 symmetry. Therefore, the magnetic orders can be divided into two classes. In Class I, if two sites i and j are related by inversion, then the spins on these sites align exactly in the same direction. i and j belong to distinct sublattices of the honeycomb lattice, either i even and j odd, or i odd and j even. The inversion only acts in real space, i.e. it does not flip the direction of the spin. In Class C 3 , the spin configuration remains invariant under a C 3 rotation about the [111] direction, with the axis of rotation piercing a site. The rotation takes place in both the real space and the spin space, In the spin space, C 3 permutes the x, y and z components of the spin. We illustrate the I and C 3 symmetries using the 18-site and 18-C 3 orders as examples in Figs. S5b and S7b. In addition, each magnetic order in Class I has three different domains with the same energy, which are related by a C 3 rotation in both the real space and the spin space. Consequently, the profile of the structure factors corresponding to these three domains differs by a C 3 rotation. For clarity, however, there is no C 3 symmetry in each domain. Each magnetic order in Class C 3 has only one domain by definition.
Furthermore, we notice that there is perhaps some kind of number rule governing the size of the magnetic unit cells, which has yet to be understood. The 98-site order appears like an augmented version of the 50-site order, which in turn appears like an augmented version of the 18-site order (see Figs. S5b and S10b, for instance). Their respective magnetic unit cells contain 7 × 7, 5 × 5 and 3 × 3 unit cells of the honeycomb lattice. Besides, the 6-site, 30-site and 70-site order, whose magnetic unit cells contain 3 × 1, 5 × 3 and 7 × 5 physical unit cells, appear like parts of the 18-site, 50-site and 98-site orders, respectively.     The static spin structure factor of the 50-site order at high fields, which peaks at k = Γ, 4 5 M i , with the same parameters in (b). The intensities along the three fold directions are same. Comparing to Fig. S10c, the intensities at 4 5 M i become lower, while the intensity at Γ becomes higher, as more spins are aligning towards the field direction.