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 CrI3, 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

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 NiX2 (X=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 NiCl2 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.

The magnetic properties of the two-dimensional VI3 bilayer are the focus of our first-principles analysis, highlighting the role of t2g orbital splitting and carried out in comparison with the CrI3 prototypical case, where the splitting is negligible. In VI3 bilayers, the empty a1g 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 VI3 stackings (labeled AB and AB′), 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 AB′ stacking can be switched from antiferromagnetic to ferromagnetic, suggesting the VI3 bilayer as an appealing candidate for electric-field-driven miniaturized spintronic devices.