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.
Read full abstract