Skyrmion Crystal in Bilinear–Biquadratic–Bicubic Model on a Centrosymmetric Triangular Lattice

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.