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

Abstract: We numerically investigate the stabilization mechanisms of skyrmion crystals under thermal fluctuations and external magnetic field in itinerant centrosymmetric tetragonal magnets. By adopting an efficient steepest descent method with a small computational cost, we systematically construct the magnetic field-temperature phase diagrams of the effective spin model derived from the itinerant electron model on a two-dimensional square lattice. As a result, we find that a square-type skyrmion crystal is stabilized by either or both of the high-harmonic wave-vector interaction and the biquadratic interaction under an external magnetic field. Especially, we discover that the former high-harmonic wave-vector interaction can stabilize the skyrmion crystal only at finite temperatures when its magnitude is small. In addition to the skyrmion crystal, we also find other stable multiple-$Q$ states in the phase diagram. Lastly, we discuss the correspondence of the phase diagrams between the effective spin model and the skyrmion-hosting material GdRu$_2$Si$_2$. The present results suggest a variety of multiple-$Q$ states could be driven by thermal fluctuations and external magnetic fields in centrosymmetric itinerant magnets.