GdRu2Si2 has recently drawn significant attention as a centrosymmetric magnet capable of hosting a short-period skyrmion square lattice (SkL) in the absence of Dzyaloshinskii-Moriya interaction (DMI). In this system, Gd atoms are arranged on a square lattice, forming 2D layers separated by the Ru-Si network in the out-of-plane direction. RKKY-type exchange between the Gd moments results in an exchange frustration, which is the main source of non-collinearity in the spiral phases of the system. So far, the literature has emphasized the importance of in-plane Q⃗ vectors in understanding the magnetic phases of the system as they appear and have been observed in the 2D Gd layers. In this work, we calculate the Gd-Gd magnetic exchange interactions (Jij) and perform atomic spin dynamics (ASD) simulations, providing new insights about GdRu2Si2. Our calculated Jij shows that the strongest magnetic interaction occurs between Gd atoms along the 111 body-diagonal direction of the unit cell. This, along with the body-centered tetragonal structure of the Gd sublattice, points to the presence of a hitherto ignored modulation vector, Q⃗111, along the 111 direction in the spiral phases of the system. We confirm this with ASD simulations and establish that the magnetic phases in GdRu2Si2 are far more complex than they appear to be in the 2D Gd layers. This interlayer modulation vector Q⃗111, along with the intralayer Q⃗100 and Q⃗010, determines the total magnetic ordering of the system. Despite its layered appearance, both structurally and magnetically, it is in reality a strong 3D magnet, requiring additional considerations for its proper theoretical modeling. Taking these into account, our ASD simulations provide an excellent description of the experimentally observed magnetic phase transitions. Furthermore, our work reveals the critical influence of dipolar interactions on the ground-state magnetic properties of GdRu2Si2, a direct consequence of the large magnetic moment of Gd. We observed a competition between dipolar interactions and uniaxial anisotropy, with each favoring a different magnetic ordering. However, the dominance of the dipolar interaction is evident in the observed ground state. This suggests a weak uniaxial anisotropy in the system, which aligns with the spherical symmetry of the Gd 4f states. Our findings offer a deeper understanding of GdRu2Si2's complex magnetism and highlight the potential for similar important interlayer effects in other layered magnetic systems.
Sarkar et al. (2025) studied this question.