Unified Treatment of Magnons and Excitons in Monolayer CrI_{3} from Many-Body Perturbation Theory

Linear response theories for interatomic exchange interactions

The linear response is a perturbation theory establishing the relationship between given physical variable and the external field inducing this variable. A well-known example of the linear response theory in magnetism is the susceptibility relating the magnetization with the magnetic field. In 1987, Liechtensteinet alcame up with the idea to formulate the problem of interatomic exchange interactions, which would describe the energy change caused by the infinitesimal rotations of spins, in terms of this susceptibility. The formulation appears to be very generic and, for isotropic systems, expresses the energy change in the form of the Heisenberg model, irrespectively on which microscopic mechanism stands behind the interaction parameters. Moreover, this approach establishes the relationship between the exchange interactions and the electronic structure obtained, for instance, in the first-principles calculations based on the density functional theory. The purpose of this review is to elaborate basic ideas of the linear response theories for the exchange interactions as well as more recent developments. The special attention is paid to the approximations underlying the original method of Liechtensteinet alin comparison with its more recent and more rigorous extensions, the roles of the on-site Coulomb interactions and the ligand states, and calculations of antisymmetric Dzyaloshinskii-Moriya interactions, which can be performed alongside with the isotropic exchange, within one computational scheme. The abilities of the linear response theories as well as many theoretical nuances, which may arise in the analysis of interatomic exchange interactions, are illustrated on magnetic van der Walls materials CrX3(X=Cl, I), half-metallic ferromagnet CrO2, ferromagnetic Weyl semimetal Co3Sn2S2, and orthorhombic manganitesAMnO3(A=La, Ho), known for the peculiar interplay of the lattice distortion, spin, and orbital ordering.

GPAW: An open Python package for electronic structure calculations

We review the GPAW open-source Python package for electronic structure calculations. GPAW is based on the projector-augmented wave method and can solve the self-consistent density functional theory (DFT) equations using three different wave-function representations, namely real-space grids, plane waves, and numerical atomic orbitals. The three representations are complementary and mutually independent and can be connected by transformations via the real-space grid. This multi-basis feature renders GPAW highly versatile and unique among similar codes. By virtue of its modular structure, the GPAW code constitutes an ideal platform for the implementation of new features and methodologies. Moreover, it is well integrated with the Atomic Simulation Environment (ASE), providing a flexible and dynamic user interface. In addition to ground-state DFT calculations, GPAW supports many-body GW band structures, optical excitations from the Bethe-Salpeter Equation, variational calculations of excited states in molecules and solids via direct optimization, and real-time propagation of the Kohn-Sham equations within time-dependent DFT. A range of more advanced methods to describe magnetic excitations and non-collinear magnetism in solids are also now available. In addition, GPAW can calculate non-linear optical tensors of solids, charged crystal point defects, and much more. Recently, support for graphics processing unit (GPU) acceleration has been achieved with minor modifications to the GPAW code thanks to the CuPy library. We end the review with an outlook, describing some future plans for GPAW.

Control of magnetic states and spin interactions in bilayer CrCl3 with strain and electric fields: an ab initio study

Using ab initio density functional theory, we demonstrated the possibility of controlling the magnetic ground-state properties of bilayer CrCl[Formula: see text] by means of mechanical strains and electric fields. In principle, we investigated the influence of these two fields on parameters describing the spin Hamiltonian of the system. The obtained results show that biaxial strains change the magnetic ground state between ferromagnetic and antiferromagnetic phases. The mechanical strain also affects the direction and amplitude of the magnetic anisotropy energy (MAE). Importantly, the direction and amplitude of the Dzyaloshinskii-Moriya vectors are also highly tunable under external strain and electric fields. The competition between nearest-neighbor exchange interactions, MAE, and Dzyaloshinskii-Moriya interactions can lead to the stabilization of various exotic spin textures and novel magnetic excitations. The high tunability of magnetic properties by external fields makes bilayer CrCl[Formula: see text] a promising candidate for application in the emerging field of two-dimensional quantum spintronics and magnonics.

Plane wave implementation of the magnetic force theorem for magnetic exchange constants: application to bulk Fe, Co and Ni

We present a plane wave implementation of the magnetic force theorem, which provides a first principles framework for extracting exchange constants parameterizing a classical Heisenberg model description of magnetic materials. It is shown that the full microscopic exchange tensor may be expressed in terms of the static Kohn-Sham susceptibility tensor and the exchange-correlation magnetic field. This formulation allows one to define arbitrary magnetic sites localized to predefined spatial regions, hence rendering the problem of finding Heisenberg parameters independent of any orbital decomposition of the problem. The susceptibility is calculated in a plane wave basis, which allows for systematic convergence with respect to unoccupied bands and spatial representation. We then apply the method to the well-studied problem of calculating adiabatic spin wave spectra for bulk Fe, Co and Ni, finding good agreement with previous calculations. In particular, we utilize the freedom of defining magnetic sites to show that the calculated Heisenberg parameters are robust towards changes in the definition of magnetic sites. This demonstrates that the magnetic sites can be regarded as well-defined and thus asserts the relevance of the Heisenberg model description despite the itinerant nature of the magnetic state.

Systematic DFT+U and Quantum Monte Carlo Benchmark of Magnetic Two-Dimensional (2D) CrX3 (X = I, Br, Cl, F)

The search for two-dimensional (2D) magnetic materials has attracted a great deal of attention because of the experimental synthesis of 2D CrI₃, which has a measured Curie temperature of 45 K. Often times, these monolayers have a higher degree of electron correlation and require more sophisticated methods beyond density functional theory (DFT). Diffusion Monte Carlo (DMC) is a correlated electronic structure method that has been demonstrated to be successful for calculating the electronic and magnetic properties of a wide variety of 2D and bulk systems, since it has a weaker dependence on the Hubbard parameter (U) and density functional. In this study, we designed a workflow that combines DFT +U and DMC in order to treat 2D correlated magnetic systems. We chose monolayer CrX₃ (X = I, Br, Cl, F), with a stronger focus on CrI₃ and CrBr₃, as a case study due to the fact that they have been experimentally realized and have a finite critical temperature. With this DFT+U and DMC workflow and the analytical method of Torelli and Olsen, we estimated a maximum value of 43.56 K for the T_c of CrI₃ and 20.78 K for the T_c of CrBr₃, in addition to analyzing the spin densities and magnetic properties with DMC and DFT+U. We expect that running this workflow for a well-known material class will aid in the future discovery and characterization of lesser known and more complex correlated 2D magnetic materials.

Computational methods for 2D materials modelling

Materials with thickness ranging from a few nanometers to a single atomic layer present unprecedented opportunities to investigate new phases of matter constrained to the two-dimensional plane. Particle-particle Coulomb interaction is dramatically affected and shaped by the dimensionality reduction, driving well-established solid state theoretical approaches to their limit of applicability. Methodological developments in theoretical modelling and computational algorithms, in close interaction with experiments, led to the discovery of the extraordinary properties of two-dimensional materials, such as high carrier mobility, Dirac cone dispersion and bright exciton luminescence, and inspired new device design paradigms. This review aims to describe the computational techniques used to simulate and predict the optical, electronic and mechanical properties of two-dimensional materials, and to interpret experimental observations. In particular, we discuss in detail the particular challenges arising in the simulation of two-dimensional constrained fermions and quasiparticles, and we offer our perspective on the future directions in this field.