Interplay of Anisotropy, Dzyaloshinskii Moriya Interaction and Symmetry breaking Fields in a 2D XY Ferromagnet

Abstract: A two dimensional ferromagnetic XY model with its bound vortex-antivortex dominated quasi long range ordered phase at low temperatures is a long standing as well as well studied problem of interest in the field of condensed matter. We conduct a detailed Monte Carlo study of such model with rather unexplored extensions where additional anisotropic exchange coupling and Dzyaloshinskii-Moriya interactions (DMI) together affect the Kosterlitz-Thoulass (KT) transition in presence/absence of symmetry breaking fields. Without DMI, the exchange term promotes collinear (ferromagnetic) order, whereas the DMI term induces spin cantings. By tuning anisotropy upto Ising limit, we document energy, specific-heat, magnetizations as well as helicity modulus and vortex densities for different tempeatures and DMI strength. We also compute the 2nd moment of correlation lengths in order to probe the spatial correlation of the spins. Furthermore, the effect of U(1) symmetry breaking 4-fold and 8-fold symmetric h4 and h8 fields are explored which shows how the double-peaked specific heat profiles changes in presence of DMI. Overall, our findings append many important updates in the low temperature phases of a topological XY ferromagnet when additional DMI and isotropy-breaking exchange and/or field terms are considered and thus providing a practical blueprint for suitably engineering topological spin systems.

• The paper demonstrates through Monte Carlo simulations how exchange anisotropy, DMI, and crystalline fields alter KT transitions and vortex dynamics in 2D XY ferromagnets.
• It reveals that increasing anisotropy sharpens specific heat peaks and drives a crossover from continuous U(1) to discrete Z2 order, impacting ferromagnetic behavior.
• The inclusion of DMI stabilizes chiral phases and shifts phase boundaries, offering insights for engineering ultrathin films and spintronic devices.

Anisotropy, Dzyaloshinskii–Moriya Interaction, and Symmetry Breaking in the 2D XY Ferromagnet

Introduction

This paper presents a comprehensive Monte Carlo analysis of the two-dimensional XY ferromagnet (XYFM) incorporating the effects of exchange anisotropy, bulk Dzyaloshinskii–Moriya interactions (DMI), and crystalline symmetry-breaking fields. The main objective is to elucidate the intertwined roles these perturbations play in the thermodynamic and topological properties of planar magnets. By systematically tuning the strength of anisotropic coupling, DMI, and various $Z_p$-symmetric field terms, the study investigates alterations in the Kosterlitz–Thouless (KT) transition, quasi-long-range order (qLRO), phase diagrams, vortex dynamics, and order parameter behavior—especially relevant for quantum materials with tunable spin–orbit coupling and engineered chiral textures.

XY Model Baseline: Thermodynamics and Topological Transition

The isotropic XY model on a 2D lattice serves as the fundamental starting point, showcasing the archetypal Berezinskii–Kosterlitz–Thouless transition driven by the binding and unbinding of vortex–antivortex pairs. Metropolis Monte Carlo simulations confirm the established phenomenology: a broad specific heat ($C_V$) hump marks the onset of topological defect proliferation at $T\sim0.895J/k_B$, while the spin stiffness (helicity modulus, $\rho_S$) exhibits a universal jump at $T_{KT}$, intersecting at $\rho_S=2T/\pi$, consistent for all simulated lattice sizes. The correlation function exponent, $\eta$, and finite-size scaling are quantitatively consistent with prior high-precision studies, affirming the robustness of the simulation protocol.

Figure 1: (a) Specific heat $C_V$ and (b) spin stiffness $\rho_S$ versus inverse temperature $J\beta$ for various lattice sizes; $C_V$0 exhibits a universal jump at the KT transition.

Exchange Anisotropy Effects

Introducing an anisotropy $C_V$1 in the XY exchange drastically alters the nature of the underlying transition. With increasing $C_V$2, the specific heat becomes sharper and transitions shift to higher temperatures, revealing crossover behavior toward Ising-like order for large anisotropy. The emergent phase exhibits true long-range ferromagnetic order at low temperatures, further corroborated by a nonzero magnetization. However, even in the strong anisotropic regime, vortex-like excitations persist, as evidenced by non-vanishing vortex densities at high temperatures. Exchange anisotropy thus breaks continuous $C_V$3 symmetry, making the system increasingly susceptible to discrete $C_V$4-type order.

Figure 2: (a) $C_V$5, (b) energy per spin $C_V$6, and (c) spin stiffness $C_V$7 versus $C_V$8 for different anisotropy parameters $C_V$9.

Inclusion of Dzyaloshinskii–Moriya Interaction

DM interaction, treated as a uniform chiral bulk term with strength $T\sim0.895J/k_B$0, fundamentally modifies the spin texture and phase diagram. Increasing DMI augments the critical temperature and stabilizes a canted qLRO, clearly observable from the shift in the $T\sim0.895J/k_B$1 peak and phase boundaries. Typical ground-state configurations display spiral modulations rather than collinear alignment, in excellent agreement with analytical expectations and prior MC studies. Notably, for commensurate DMI, the system supports nontrivial diagonal arrangements of spins at low temperature, which thermal fluctuations disrupt at elevated temperatures.

Figure 3: $T\sim0.895J/k_B$2 as a function of $T\sim0.895J/k_B$3 for different DMI strength $T\sim0.895J/k_B$4; inset: phase diagram for QLRO–DO phase.

Figure 4: Representative low-$T\sim0.895J/k_B$5 and high-$T\sim0.895J/k_B$6 spin configurations ($T\sim0.895J/k_B$7) for various $T\sim0.895J/k_B$8 pairs, displaying the vortex-antivortex structure and canted ordering induced by DMI.

Competition and Interplay: Anisotropy vs. DMI

When both anisotropy $T\sim0.895J/k_B$9 and DMI $\rho_S$0 are present, strong competition arises between collinear and chiral tendencies. DMI drives the pseudo-critical temperature upward and suppresses magnetization, whereas anisotropy fosters ferromagnetic order. The phase boundary's response to $\rho_S$1 in the presence of DMI is non-monotonic, with $\rho_S$2 profiles, magnetization curves, and spin stiffness reflecting the shifting balance and rich crossover phenomena. The second-moment correlation length, computed for a range of system sizes, provides key insights into the spatial scale of magnetic coherence and accurately identifies crossover or transition regions.

Figure 5: (a) $\rho_S$3 vs $\rho_S$4 for $\rho_S$5 at several $\rho_S$6; (b,c) phase boundaries from $\rho_S$7 peaks for $\rho_S$8 and $\rho_S$9 respectively.

Figure 6: (a) $T_{KT}$0, (b) $T_{KT}$1, (c) $T_{KT}$2, and (d) vortex density $T_{KT}$3 versus $T_{KT}$4 for several $T_{KT}$5 at $T_{KT}$6.

Figure 7: Second-moment correlation length as a function of temperature for different values of DMI and anisotropy.

Symmetry-Breaking Fields: Multipeak Specific Heat and Competing Orders

The inclusion of $T_{KT}$7 and $T_{KT}$8 symmetry-breaking crystalline fields (parameterized as $T_{KT}$9 and $\rho_S=2T/\pi$0) introduces further complexity to the phase structure. The principal observation is the appearance of multi-peak $\rho_S=2T/\pi$1 signatures indicating multiple phase transitions or crossovers. For compatible fields ($\rho_S=2T/\pi$2), a single transition plus low-$\rho_S=2T/\pi$3 crossover is observed; in the competing regime ($\rho_S=2T/\pi$4, $\rho_S=2T/\pi$5), a sharp low-$\rho_S=2T/\pi$6 peak and a broader KT-like peak manifest, signaling successive transitions between FM, KT, and paramagnetic phases. The presence of DMI shifts and sometimes suppresses these features, demonstrating the nontrivial interplay between chiral and symmetry-breaking perturbations.

Figure 8: $\rho_S=2T/\pi$7 versus $\rho_S=2T/\pi$8 for several field configurations and system sizes, with and without DMI, highlighting double-peak structure characteristic of competing symmetry-breaking fields.

Implications and Future Directions

The study establishes that DMI and exchange anisotropy act as powerful control parameters for tuning the thermal and topological behavior of 2D magnetic systems—relevant to the engineering of ultrathin films, interfacial magnets, and van der Waals materials with large spin–orbit coupling. The controlled introduction of DMI or anisotropy can stabilize novel chiral or Ising-like phases and manipulate vortex-driven transitions, while crystalline fields further enrich the landscape with multi-transition regimes. These results lay a rigorous foundation for subsequent studies on quantum and 3D generalizations, with implications for emergent phenomena such as skyrmion nucleation and transport (e.g., skyrmion Hall effect), as well as for understanding field-driven dynamics in topological magnonic materials.

Conclusion

By systematically mapping the combined effects of anisotropy, DMI, and symmetry-breaking fields, this paper advances a nuanced and quantitative understanding of the 2D XYFM. The equilibrium phase diagrams and thermodynamic observables reveal intricate competition and coupling effects, consistent with predictions from renormalization group and large-scale MC studies. The findings not only quantify how each perturbation modifies KT-like criticality but also provide specific predictions for experimental platforms where these couplings are tunable. The general framework and numerical strategies employed offer a versatile template for future investigations into topological matter, quantum fluctuations, and chiral spin structures in higher dimensional systems.