Sextuple-Q Spin States in Centrosymmetric Hexagonal Magnets

We theoretically investigate topological transitions between coplanar and noncoplanar magnetic states in centrosymmetric itinerant magnets on a square lattice. A canonical effective spin model incorporating bilinear and biquadratic exchange interactions at finite wave vectors is analyzed to elucidate the emergence of multiple-Q magnetic orders. By taking into account high-harmonic wave-vector interactions, we demonstrate that a coplanar double-Q spin texture continuously evolves into a noncoplanar triple-Q state carrying a finite scalar spin chirality. The stability of these multiple-Q states is examined using simulated annealing as a function of the relative strengths of the high-harmonic coupling, the biquadratic interaction, and the external magnetic field. The resulting phase diagrams reveal a competition between double-Q and triple-Q states, where the noncoplanar triple-Q phase is stabilized through the cooperative effect of the high-harmonic and biquadratic interactions. Real-space spin textures, spin structure factors, and scalar spin chirality distributions are analyzed to characterize the distinct magnetic phases and the topological transitions connecting them. These findings provide a microscopic framework for understanding the emergence of noncoplanar magnetic textures driven by the interplay between two- and four-spin interactions in centrosymmetric itinerant magnets.

In this study, we investigate the stability of a magnetic skyrmion crystal with short-period magnetic modulations in a centrosymmetric body-centered tetragonal system. By performing the simulated annealing for the spin model, incorporating the effects of the biquadratic interaction and high-harmonic wave–vector interaction in momentum space, we find that the double-Q square skyrmion crystal consisting of two spin density waves is stabilized in an external magnetic field. We also show that double-Q states appear in both low- and high-field regions; the low-field spin configuration is characterized by an anisotropic double-Q modulation consisting of a superposition of the spiral wave and sinusoidal wave, while the high-field spin configuration is characterized by an isotropic double-Q modulation consisting of a superposition of two sinusoidal waves. Furthermore, we show that the obtained multiple-Q instabilities can be realized for various ordering wave vectors. The results provide the possibility of realizing the short-period skyrmion crystals under the body-centered tetragonal lattice structure.

We theoretically investigate multiple-$Q$ spin textures, which are composed of superpositions of spin density waves with different wave numbers, for an effective spin model of centrosymmetric itinerant magnets. Our focus is on the interplay between biquadratic interactions arising from the spin-charge coupling and magnetic anisotropy caused by the spin-orbit coupling. Taking into account two types of the magnetic anisotropy, single-ion anisotropy and bond-dependent anisotropy, we elucidate magnetic phase diagrams for an archetypal triangular-lattice system in the absence and presence of an external magnetic field. In the case of the single-ion anisotropy, we find a plethora of multiple-$Q$ instabilities depending on the strength and the sign of the anisotropy (easy plane or easy axis), including a noncoplanar triple-$Q$ state regarded as a skyrmion crystal with topological number of two, and coplanar and noncoplanar double-$Q$ states. In an external magnetic field, we find that another noncoplanar triple-$Q$ state, a skyrmion crystal with topological number of one, is stabilized by the in-plane (out-of-plane) magnetic field under the easy-plane (easy-axis) anisotropy. A part of the results, especially for the relatively large biquadratic interaction, qualitatively reproduce those in the Kondo lattice model, which explicitly includes itinerant electrons [S. Hayami and Y. Motome, Phys. Rev. B 99, 094420 (2019)]. We also examine the stability of the field-induced skyrmion crystal by rotating the field direction. As a by-product, we show that a different triple-$Q$ state with nonzero chirality appears in the presence of the biquadratic interaction and the easy-axis anisotropy. Meanwhile, we find that the bond-dependent anisotropy also stabilizes both types of skyrmion crystals. We show that, however, for the skyrmion crystal with topological number of one, Bloch- and N\'eel-type skyrmion crystals are selectively realized depending on the sign of the bond-dependent anisotropy, since this anisotropy selects a particular set of the helicity and vorticity. Moreover, we find yet another multiple-$Q$ states with nonzero spin scalar chirality, including two types of meron crystals with the skyrmion numbers of one and two. The systematic investigation of multiple-$Q$ instabilities in triangular itinerant magnets will provide a reference to complex magnetic textures in centrosymmetric magnetic metals.

Widely-sweeping magnetic field–temperature phase diagrams for skyrmion-hosting centrosymmetric tetragonal magnets

Equilibrium properties of a spin-1 Ising system with bilinear, biquadratic and odd interactions

A renormalization-group study into the effects of competing biquadratic and crystal-field interactions in the Blume–Emery–Griffiths model

Dark solitonic interaction and conservation laws for a higher-order [formula omitted]-dimensional nonlinear Schrödinger-type equation in a Heisenberg ferromagnetic spin chain with bilinear and biquadratic interaction

Dynamic magnetic properties the spin-1 Ising model with bilinear and biquadratic interactions within the path probability method

The lowest energy configurations of short odd open chains with classical spins are determined for antiferromagnetic bilinear and biquadratic nearest-neighbor exchange interactions. The zero field residual magnetization generates differences with the magnetic behavior of even chains, as the odd chain is like a small magnet for weak magnetic fields. The lowest energy configuration is calculated as a function of the total magnetization M, even for M less than the zero field residual magnetization. Analytic expressions and their proofs are provided for the threshold magnetic field needed to drive the system away from the antiferromagnetic configuration and the spin polar angles in its vicinity, when the biquadratic interaction is relatively weak. They are also given for the saturation magnetic field and the spin polar angles close to it. Finally, an analytic expression along with its proof is given for the maximum magnetization in zero magnetic field for stronger biquadratic interaction, where the lowest energy configuration is highly degenerate.

Based on group theoretical arguments we derive the most general Hamiltonian for the Bi2Se3-class of materials including terms to third order in the wave vector, first order in electric and magnetic fields, first order in strain and first order in both strain and wave vector. We determine analytically the effects of strain on the electronic structure of Bi2Se3. For the most experimentally relevant surface termination we analytically derive the surface state (SS) spectrum, revealing an anisotropic Dirac cone with elliptical constant energy contours giving rise to a direction-dependent group velocity. The spin-momentum locking of strained Bi2Se3 is shown to be modified. Hence, strain control can be used to manipulate the spin degree of freedom via the spin–orbit coupling. We show that for a thin film of Bi2Se3 the SS band gap induced by coupling between the opposite surfaces changes opposite to the bulk band gap under strain. Tuning the SS band gap by strain, gives new possibilities for the experimental investigation of the thickness dependent gap and optimization of optical properties relevant for, e.g., photodetector and energy harvesting applications. We finally derive analytical expressions for the effective mass tensor of the Bi2Se3 class of materials as a function of strain and electric field.

Evidence for biquadratic exchange interactions in GdAg1_xZnx

We numerically investigate the effect of multi-spin interactions on the stability of skyrmion crystals and other multiple-Q magnetic states, with a particular emphasis on the momentum-resolved bicubic interaction. By performing simulated annealing for an effective spin model that incorporates bilinear, biquadratic, and bicubic interactions on a two-dimensional triangular lattice, we construct the corresponding low-temperature phase diagram. Our results reveal that a positive bicubic interaction stabilizes a skyrmion crystal with a skyrmion number of two, whereas a negative bicubic interaction favors a single-Q spiral state. Moreover, we demonstrate that the stability region of the field-induced skyrmion crystal with the skyrmion number of one is largely enlarged in the presence of a positive bicubic interaction.

Noncollinear and noncoplanar magnetic textures including skyrmions and vortices act as emergent electromagnetic fields and give rise to novel electronic and transport properties. We here report a unified understanding of noncoplanar magnetic orderings emergent from the spin-charge coupling in itinerant magnets. The mechanism has its roots in effective multiple spin interactions beyond the conventional Ruderman-Kittel-Kasuya-Yosida (RKKY) mechanism, which are ubiquitously generated in itinerant electron systems with local magnetic moments. Carefully examining the higher-order perturbations in terms of the spin-charge coupling, we construct a minimal effective spin model composed of the bilinear and biquadratic interactions with particular wave numbers dictated by the Fermi surface. We find that our effective model captures the underlying physics of the instability toward noncoplanar multiple-$Q$ states discovered recently in itinerant magnets: a single-$Q$ helical state expected from the RKKY theory is replaced by a double-$Q$ vortex crystal with chirality density waves even for an infinitely small spin-charge coupling on generic lattices [R. Ozawa $\mathit{et \ al.}$, J. Phys. Soc. Jpn. $\mathbf{85}$, 103703 (2016)], and a triple-$Q$ skyrmion crystal with a high topological number of two appears with increasing the spin-charge coupling on a triangular lattice [R. Ozawa, S. Hayami, and Y. Motome, to appear in Phys. Rev. Lett. (arXiv: 1703.03227)]. We find that, by introducing an external magnetic field, our effective model exhibits a plethora of multiple-$Q$ states. Our findings will serve as a guide for exploring further exotic magnetic orderings in itinerant magnets.

We investigate an instability toward a square-lattice formation of magnetic skyrmions in centrosymmetric layered systems. By focusing on a bilayer square-lattice structure with the inversion center at the interlayer bond instead of the atomic site, we numerically examine the stability of the square skyrmion crystal (SkX) based on an effective spin model with the momentum-resolved interaction in the ground state through the simulated annealing. As a result, we find that a layer-dependent staggered Dzyaloshinskii–Moriya (DM) interaction built in the lattice structure becomes the origin of the square SkX in an external magnetic field irrespective of the sign of the interlayer exchange interaction. The obtained square SkX is constituted of the SkXs with different helicities in each layer due to the staggered DM interaction. Furthermore, we show that the interplay between the staggered DM interaction and the interlayer exchange interaction gives rise to a double-Q state with a uniform component of the scalar chirality in the low-field region. The present results provide another way of stabilizing the square SkX in centrosymmetric magnets, which will be useful to explore further exotic topological spin textures.

We conduct a numerical investigation into the stability of a quadruple-Q skyrmion crystal, a structure generated by the superposition of four spin density waves traveling in distinct directions within three-dimensional space, hosted on a centrosymmetric body-centered tetragonal lattice. Using simulated annealing applied to an effective spin model that includes momentum-resolved bilinear and biquadratic interactions, we construct a magnetic phase diagram spanning a broad range of model parameters. Our study finds that a quadruple-Q skyrmion crystal does not emerge within the phase diagram when varying the biquadratic interaction and external magnetic field. Instead, three distinct quadruple-Q states with topologically trivial spin textures are stabilized. However, we demonstrate that the quadruple-Q skyrmion crystal can become the ground state when an additional high-harmonic wave–vector interaction is considered. Depending on the magnitude of this interaction, we obtain two types of quadruple-Q skyrmion crystals exhibiting the skyrmion numbers of one and two. These findings highlight the emergence of diverse three-dimensional multiple-Q spin states in centrosymmetric body-centered tetragonal magnets.