A combined first principles study of the structural, magnetic, and phonon properties of monolayer CrI3.

The first magnetic 2D material discovered, monolayer (ML) CrI3, is particularly fascinating due to its ground state ferromagnetism. However, because ML materials are difficult to probe experimentally, much remains unresolved about ML CrI3's structural, electronic, and magnetic properties. Here, we leverage Density Functional Theory (DFT) and high-accuracy Diffusion Monte Carlo (DMC) simulations to predict lattice parameters, magnetic moments, and spin-phonon and spin-lattice coupling of ML CrI3. We exploit a recently developed surrogate Hessian DMC line search technique to determine CrI3's ML geometry with DMC accuracy, yielding lattice parameters in good agreement with recently published STM measurements-an accomplishment given the ∼10% variability in previous DFT-derived estimates depending upon the functional. Strikingly, we find that previous DFT predictions of ML CrI3's magnetic spin moments are correct on average across a unit cell but miss critical local spatial fluctuations in the spin density revealed by more accurate DMC. DMC predicts that magnetic moments in ML CrI3 are 3.62 µB per chromium and -0.145 µB per iodine, both larger than previous DFT predictions. The large disparate moments together with the large spin-orbit coupling of CrI3's I-p orbital suggest a ligand superexchange-dominated magnetic anisotropy in ML CrI3, corroborating recent observations of magnons in its 2D limit. We also find that ML CrI3 exhibits a substantial spin-phonon coupling of ∼3.32 cm-1. Our work, thus, establishes many of ML CrI3's key properties, while also continuing to demonstrate the pivotal role that DMC can assume in the study of magnetic and other 2D materials.

First principles calculations of the electric field gradient tensors of Ba2NaOsO6, a Mott insulator with strong spin orbit coupling.

We present first principles calculations of the electrostatic properties of Ba2NaOsO6 (BNOO), a 5d1 Mott insulator with strong spin orbit coupling (SOC) in its low temperature quantum phases. In light of recent NMR experiments showing that BNOO develops a local octahedral distortion that is accompanied by the emergence of an electric field gradient (EFG) and precedes the formation of long range magnetic order (Lu et al 2017 Nat. Commun. 8 14407, Liu et al 2018 Phys. Rev. B 97 224103; Liu et al 2018 Physica B 536 863), we calculated BNOO's EFG tensor for several different model distortions. The local orthorhombic distortion that we identified as most strongly agreeing with experiment corresponds to a Q2 distortion mode of the Na-O octahedra, in agreement with conclusions given in (Liu et al 2018 Phys. Rev. B 97 224103). Furthermore, we found that the EFG is insensitive to the type of underlying magnetic order. By combining NMR results with first principles modeling, we have thus forged a more complete understanding of BNOO's structural and magnetic properties, which could not be achieved based upon experiment or theory alone.

How electronic dynamics with Pauli exclusion produces Fermi-Dirac statistics.

It is important that any dynamics method approaches the correct population distribution at long times. In this paper, we derive a one-body reduced density matrix dynamics for electrons in energetic contact with a bath. We obtain a remarkable equation of motion which shows that in order to reach equilibrium properly, rates of electron transitions depend on the density matrix. Even though the bath drives the electrons towards a Boltzmann distribution, hole blocking factors in our equation of motion cause the electronic populations to relax to a Fermi-Dirac distribution. These factors are an old concept, but we show how they can be derived with a combination of time-dependent perturbation theory and the extended normal ordering of Mukherjee and Kutzelnigg for a general electronic state. The resulting non-equilibrium kinetic equations generalize the usual Redfield theory to many-electron systems, while ensuring that the orbital occupations remain between zero and one. In numerical applications of our equations, we show that relaxation rates of molecules are not constant because of the blocking effect. Other applications to model atomic chains are also presented which highlight the importance of treating both dephasing and relaxation. Finally, we show how the bath localizes the electron density matrix.

Calculation of multipolar exchange interactions in spin-orbital coupled systems.

A new method of computing multipolar exchange interaction in spin-orbit coupled systems is developed using multipolar tensor expansion of the density matrix in local density approximation+U electronic structure calculation. Within the mean field approximation, exchange constants can be mapped into a series of total energy calculations by the pair-flip approximation technique. The application to uranium dioxide shows an antiferromagnetic superexchange coupling in dipoles but a ferromagnetic one in quadrupoles which is very different from past studies. Further calculation of the spin-lattice interaction indicates it is of the same order with the superexchange and characterizes the overall behavior of the quadrupolar part as a competition between them.