Topological spin structures, such as magnetic skyrmions, hold great promises for data storage applications, thanks to their inherent stability. In most cases, skyrmions are stabilized by magnetic fields in non-centrosymmetric systems displaying the chiral Dzyaloshinskii-Moriya exchange interaction, while spontaneous skyrmion lattices have been reported in centrosymmetric itinerant magnets with long-range interactions. Here, a spontaneous anti-biskyrmion lattice with unique topology and chirality is predicted in the monolayer of a semiconducting and centrosymmetric metal halide, NiI₂. Our first-principles and Monte Carlo simulations reveal that the anisotropies of the short-range symmetric exchange, when combined with magnetic frustration, can lead to an emergent chiral interaction that is responsible for the predicted topological spin structures. The proposed mechanism finds a prototypical manifestation in two-dimensional magnets, thus broadening the class of materials that can host spontaneous skyrmionic states.

Two-dimensional (2D) van der Waals (vdW) magnets provide an ideal platform for exploring, on the fundamental side, new microscopic mechanisms and for developing, on the technological side, ultracompact spintronic applications. So far, bilinear spin Hamiltonians have been commonly adopted to investigate the magnetic properties of 2D magnets, neglecting higher order magnetic interactions. However, we here provide quantitative evidence of giant biquadratic exchange interactions in monolayer ${\mathrm{NiX}}_2$ ($\mathrm{X}=\mathrm{Cl}$, Br and I), by combining first-principles calculations and the newly developed machine learning method for constructing Hamiltonian. Interestingly, we show that the ferromagnetic ground state within ${\mathrm{NiCl}}_2$ single layers cannot be explained by means of the bilinear Heisenberg Hamiltonian; rather, the nearest-neighbor biquadratic interaction is found to be crucial. Furthermore, using a three-orbitals Hubbard model, we propose that the giant biquadratic exchange interaction originates from large hopping between unoccupied and occupied orbitals on neighboring magnetic ions. On a general framework, our work suggests biquadratic exchange interactions to be important in 2D magnets with edge-shared octahedra.

Curved magnets attract considerable interest for their unusually rich phase diagram, often encompassing exotic (e.g., topological or chiral) spin states. Micromagnetic simulations are playing a central role in the theoretical understanding of such phenomena; their predictive power, however, rests on the availability of reliable model parameters to describe a given material or nanostructure. Here we demonstrate how noncollinear-spin polarized density-functional theory can be used to determine the flexomagnetic coupling coefficients in real systems. By focusing on monolayer ${\mathrm{CrI}}_3$, we find a crossover as a function of curvature between a magnetization normal to the surface to a cycloidal state, which we rationalize in terms of effective anisotropy and Dzyaloshinskii-Moriya contributions to the magnetic energy. Our results reveal an unexpectedly large impact of spin-orbit interactions on the curvature-induced anisotropy, which we discuss in the context of existing phenomenological models.

Motivated by the recent wealth of exotic magnetic phases emerging in two-dimensional frustrated lattices, we investigate the origin of possible magnetism in the monolayer family of triangular lattice materials $M{X}_2$ ($M$=V, Mn, Ni and $X$=Cl, Br, I). We first show that consideration of general properties such as filling and hybridization enables to formulate the trends for the most relevant magnetic interaction parameters. In particular, we observe that the effects of spin-orbit coupling (SOC) can be effectively tuned through the ligand elements as the considered $3d$ transition metal ions do not strongly contribute to the anisotropic component of the intersite exchange interaction. Consequently, we find that the corresponding SOC matrix elements differ significantly from the atomic limit. In the next step and by using two ab initio based complementary approaches, we extract realistic effective spin models and find that in the case of heavy ligand elements, SOC effects manifest in anisotropic exchange and single-ion anisotropy only for specific fillings.

The effects of competing magnetic interactions in stabilizing different spin configurations are drawing renewed attention in order to unveil emerging topological spin textures and to highlight microscopic mechanisms leading to their stabilization. The possible key role of the two-site exchange anisotropy in selecting specific helicity and vorticity of skyrmionic lattices has only recently been proposed. In this work, we explore the phase diagram of a frustrated localized magnet characterized by a two-dimensional centrosymmetric triangular lattice, focusing on the interplay between the two-ion anisotropy and the single-ion anisotropy. The effects of an external magnetic field applied perpendicularly to the magnetic layer, are also investigated. By means of Monte Carlo simulations, we find an abundance of different spin configurations, going from trivial to high-order Q skyrmionic and meronic lattices. In closer detail, we find that a dominant role is played by the two-ion over the single-ion anisotropy in determining the planar spin texture; the strength and the sign of single ion anisotropy, together with the magnitude of the magnetic field, tune the perpendicular spin components, mostly affecting the polarity (and, in turn, the topology) of the spin texture. Our analysis confirms the crucial role of the anisotropic symmetric exchange in systems with dominant short-range interactions; at the same time, we predict a rich variety of complex magnetic textures, which may arise from a fine tuning of competing anisotropic mechanisms.

The magnetic properties of the two-dimensional ${\mathrm{VI}}_3$ bilayer are the focus of our first-principles analysis, highlighting the role of $t_{2g}$ orbital splitting and carried out in comparison with the ${\mathrm{CrI}}_3$ prototypical case, where the splitting is negligible. In ${\mathrm{VI}}_3$ bilayers, the empty $a_{1g}$ state is found to play a crucial role in both stabilizing the insulating state and in determining the interlayer magnetic interaction. Indeed, an analysis based on maximally localized Wannier functions allows one to evaluate the interlayer exchange interactions in two different ${\mathrm{VI}}_3$ stackings (labeled AB and ${\mathrm{AB}}^{\prime }$), to interpret the results in terms of the virtual-hopping mechanism, and to highlight the strongest hopping channels underlying the magnetic interlayer coupling. Upon application of electric fields perpendicular to the slab, we find that the magnetic ground state in the ${\mathrm{AB}}^{\prime }$ stacking can be switched from antiferromagnetic to ferromagnetic, suggesting the ${\mathrm{VI}}_3$ bilayer as an appealing candidate for electric-field-driven miniaturized spintronic devices.