Ferromagnetism in 2D Vanadium Diselenide

Tunable Valley Pseudospin and Electron-Phonon Coupling in WSe2/1T-VSe2 Heterostructures

Heterostructure engineering provides versatile platforms for exploring exotic physics and enhancing the device performance through interface coupling. Despite the rich array of physical phenomena presented by heterostructures composed of semiconductor and metal van der Waals materials, significant gaps remain in understanding their optical, thermal, and electronic properties. Here, we demonstrate that the valley pseudospin and electron-phonon coupling in monolayer WSe2 are significantly influenced by interface coupling with 1T-VSe2. The heterointerface alters the relaxation process of valley excitons, leading to a transition in magnetic-field-dependent valley polarization from a linear to a "V" shape. Furthermore, we uncover that enhanced electron-phonon coupling exacerbates variations in exciton and valley exciton behavior with temperature, involving higher phonon energies and a shift from acoustic to optical phonons. These findings highlight a promising pathway to manipulate valley excitons and investigate electron-phonon coupling through van der Waals interface interactions.

Controlled growth of 3R phase niobium diselenide and its properties

Inversion symmetry broken 3R phase transition metal dichalcogenides (TMDs) show fascinating prospects in spintronics, valleytronics, and nonlinear optics. However, the controlled synthesis of 3R phase TMDs is still a great challenge. In this work, two-dimensional 3R-NbSe2 single crystals up to 0.2 mm were synthesized for the first time through chemical vapor deposition method by designing a space-confined system. The crystal size and morphology can be controlled by the location of the stacked substrates and the amount of the Nb2O5 precursor. Scanning transmission electron microscopy and Raman measurements reveal the NbSe2 exhibits a pure 3R stacking mode with relatively weak interlayer van der Waals interactions. Importantly, 3R-NbSe2 shows obvious second harmonic generation signal which intensity intensified as thickness increases. Density functional theory calculations and optical absorption demonstrate the coexistence of metallic and semiconducting optical properties of 3R-NbSe2. We designed a NbSe2/WS2/NbSe2 photodetector utilizing the metallicity of 3R-NbSe2, which shows good performance especially an ultrafast response (6-7 μs, 0.5 ms - 7.9 s for Au electrodes in literature). The proposed strategy and findings are of great significance for the growth of many other 3R-TMDs and applications of nonlinear optical and ultrafast devices.

Ferroelectrically controlled electromagnetic and transport properties of VN2H2/Al2O3 van der Waals multiferroic heterostructures

The vertical integration of a ferromagnetic monolayer and a ferroelectric monolayer into van der Waals heterostructures offers a promising route to achieve two-dimensional multiferroic semiconductors owing to the lack of intrinsic single-phase multiferroic materials in nature. In this study, we propose a VN2H2/Al2O3 van der Waals magnetoelectric multiferroic heterostructure and investigate its electronic, magnetic, and transport properties using density functional theory combined with the Boltzmann transport theory. The VN2H2 monolayer is a room-temperature ferromagnetic semiconductor with a band gap of 0.24 eV and a Curie temperature of 411 K, while the Al2O3 monolayer is a ferroelectric semiconductor with a polarization value of 0.11 C m-2. In the VN2H2/Al2O3 van der Waals heterostructures, the conversion between the metal and the semiconductor can be controlled by altering the polarization of the Al2O3 layer. The VN2H2/Al2O3 van der Waals heterostructure retains room-temperature ferromagnetism, and the reverse of polarization is accompanied with a change in the direction of the easy magnetization axis. In addition, electrostatic doping can significantly improve the conductivity of the downward polarization state and transform the upward polarization state from a metal to a half-metal, achieving 100% spin polarization. Our results thus pave the way for achieving highly tunable electromagnetic and transport properties in van der Waals magnetoelectric heterostructures, which have potential applications in next-generation low-power logic and memory devices.

Phase-Dependent Magnetic Proximity Modulations on Valley Polarization and Splitting

Proximate-induced magnetic interactions present a promising strategy for precise manipulation of valley degrees of freedom. Taking advantage of the splendid valleytronic platform of transition metal dichalcogenides, magnetic two-dimensional VSe2 with different phases are introduced to intervene in the spin of electrons and modulate their valleytronic properties. When constructing the heterostructures, 1T-VSe2/WX2 (X = S and Se) showcases significant improvement in the valley polarizations at room temperature, while 2H-VSe2/WX2 exhibits superior performance at low temperatures and demonstrates heightened sensitivity to the external magnetic field. Simultaneously, considerable valley splitting with a large geff factor up to -29.0 is observed in 2H-VSe2/WS2, while it is negligible in 1T-VSe2/WX2. First-principles calculations reveal a phase-dependent magnetic proximity mechanism on the valleytronic modulations, which is dominated by interfacial charge transfer in 1T-VSe2/WX2 and the proximity exchange field in 2H-VSe2/WX2 heterostructures. The effective control over valley degrees of freedom will bridge the valleytronic physics and devices, rendering enormous potential in the field of valley quantum applications.

Near-room temperature ferromagnetism and a tunable anomalous Hall effect in atomically thin Fe4CoGeTe2

Itinerant ferromagnetism at room temperature is a key factor for spin transport and manipulation. Here, we report the realization of near-room temperature itinerant ferromagnetism in Co doped Fe5GeTe2 thin flakes. The ferromagnetic transition temperature TC (∼323 K-337 K) is almost unchanged when the thickness is as low as 12 nm and is still about 284 K at 2 nm (bilayer thickness). Theoretical calculations further indicate that the ferromagnetism persists in monolayer Fe4CoGeTe2. In addition to the robust ferromagnetism down to the ultrathin limit, Fe4CoGeTe2 exhibits an unusual temperature- and thickness-dependent intrinsic anomalous Hall effect. We propose that it could be ascribed to the dependence of the band structure on thickness that changes the Berry curvature near the Fermi energy level subtly. The near-room temperature ferromagnetism and tunable anomalous Hall effect in atomically thin Fe4CoGeTe2 provide opportunities to understand the exotic transport properties of two-dimensional van der Waals magnetic materials and explore their potential applications in spintronics.

Tunable polarization properties of charge, spin, and valley in Janus VSiGeZ4 (Z = N, P, As) monolayers

Polarization, as an important characterization of the symmetry breaking systems, has attracted tremendous attention in two-dimensional (2D) materials. Due to their significant symmetry breaking, Janus 2D ferrovalley materials provide a desirable platform to investigate the charge, spin, and valley polarization, as well as their coupling effects. Herein, using first-principles calculations, the polarization properties of charge, spin, and valley in Janus VSiGeZ4 (Z = N, P, and As) monolayers are systematically studied. The mirror symmetry breaking leads to a non-zero dipole moment and surface work function difference, indicating the presence of out-of-plane charge polarization. Magnetic properties calculations demonstrate that VSiGeN4 is a 2D-XY magnet with a Berezinskii-Kosterlitz-Thouless temperature of 342 K, while VSiGeP4 and VSiGeAs4 have an out-of-plane magnetization with a Curie temperature below room temperature. The magnetization can be rotated by applying biaxial strain, allowing manipulation of the spin polarization via nonmagnetic means. The spontaneous valley polarization is predicted to be 46, 49, and 70 meV for VSiGeN4, VSiGeP4, and VSiGeAs4, respectively, whose physical origin can be elucidated by employing the model analysis. In particular, the biaxial strain can induce the valley polarization switching from the valence (conduction) band to conduction (valence) band, but it hardly changes the valley polarization strength. Meanwhile, the valley extremum is transformed from the K' (K) to K (K') points. The present work not only provides an underlying insight into the polarization properties of Janus VSiGeZ4 but also offers a class of promising materials for spintronic and valleytronic devices.

Strong Room-Temperature Ferromagnetism in Ultrathin NiOOH Nanosheets through Surfactant Manipulation

Two-dimensional (2D) ferromagnetic (FM) materials with nanoscale thickness and spontaneous net magnetization have emerged as a promising class of functional materials for applications in next-generation spintronics, quantum processing, and data storage devices. However, most 2D materials exhibit weak FM even at low temperatures, limiting their potential applications in many technological fields. The fabrication of strong room-temperature FM 2D materials is highly desirable for the development of practical applications. Here, we demonstrate an ionic layer epitaxy strategy to synthesize few-layered NiOOH nanosheets with strong room-temperature FM and a saturation magnetization up to 409.86 emu cm-3 at 300 K. The results are consistent with the ab initio predictions of a stable FM NiOOH nanolayer structure with an FM configuration. The FM strength of the NiOOH nanosheets can be tuned by controlling the surfactant monolayer density and annealing. This work offers a promising strategy for achieving strong high-temperature FM in 2D materials for spintronic applications.

Magnetic semiconducting borophenes and their derivatives

Two semiconducting borophenes with layer-dependent magnetism are predicted based on spin-polarized density functional theory. Both monolayer borophenes are ferromagnetic. One is composed of B3 and B15 triangular motifs, exhibiting bipolar spin polarization and a magnetic moment of 1.00 μB per primitive cell. The other consists of B15 triangular motifs, possessing a Curie temperature of about 437 K and a magnetic moment of 3.00 μB per primitive cell. B atoms located between the triangular motifs are essential for inducing ferromagnetism in monolayer borophenes. However, bilayer borophenes with high-symmetry stacking orders are nonmagnetic. Furthermore, magnetic boron nanotubes and fullerenes could be made of monolayer borophenes. Finally, we propose to fabricate these magnetic semiconducting borophenes from the buckled triangular structure of borophenes via selective electron beam ionization of B atoms by scanning transmission electron microscopy.

Electron transport, ferroelectric, piezoelectric and optical properties of two-dimensional In2Te3: a first-principles study

Two-dimensional (2D) materials have garnered significant interest in the fields of optoelectronics and electronics due to their unique and diverse properties. In this work, the electron transport, ferroelectric, piezoelectric, and optical properties of 2D In2Te3 were systematically investigated using first-principles based on density functional theory. The analysis of the phonon spectrum and elastic modulus of the Born effective criterion indicates that the structure of the novel 2D In2Te3 is dynamically stable. The calculation results show that 2D In2Te3 exhibits a carrier mobility as high as 3680.99 cm2 V-1 s-1 (y direction), a high in-plane polarization of 2.428 × 10-10 C m-1, and an excellent ferroelectric phase transition barrier (52.847 meV) and piezoelectric properties (e11 = 1.52 × 10-10 C m-1). The higher carrier mobility is attributed to the band degeneracy and small carrier effective mass. In addition, biaxial strain is an effective way to modulate the band gap and optical properties of 2D In2Te3. These properties indicate that 2D In2Te3 is a promising candidate material for flexible electronic devices and ferroelectric photovoltaic devices.

Inducing itinerant ferromagnetism by manipulating van Hove singularity in epitaxial monolayer 1T-VSe2

The itinerant ferromagnetism can be induced by a van Hove singularity (VHS) with a divergent density of states at Fermi level. Utilizing the giant magnified dielectric constant εr of SrTiO3(111) substrate with cooling, here we successfully manipulated the VHS in the epitaxial monolayer (ML) 1T-VSe2 film approaching to Fermi level via the large interfacial charge transfer, and thus induced a two-dimensional (2D) itinerant ferromagnetic state below 3.3 K. Combining the direct characterization of the VHS structure via angle-resolved photoemission spectroscopy (ARPES), together with the theoretical analysis, we ascribe the manipulation of VHS to the physical origin of the itinerant ferromagnetic state in ML 1T-VSe2. Therefore, we further demonstrated that the ferromagnetic state in the 2D system can be controlled through manipulating the VHS by engineering the film thickness or replacing the substrate. Our findings clearly evidence that the VHS can serve as an effective manipulating degree of freedom for the itinerant ferromagnetic state, expanding the application potentials of 2D magnets for the next-generation information technology.

Spontaneous valley polarization and valley-nonequilibrium quantum anomalous Hall effect in Janus monolayer ScBrI

Topology and ferrovalley (FV) are two essential concepts in emerging device applications and the fundamental research field. To date, relevant reports are extremely rare about the coupling of FV and topology in a single system. By Monte Carlo (MC) simulations and first-principles calculations, a stable intrinsic FV ScBrI semiconductor with high Curie temperature (T_C) is predicted. Because of the combination of spin-orbital coupling (SOC) and exchange interaction, the Janus monolayer ScBrI shows a spontaneous valley polarization of 90 meV, which is located in the top valence band. For the magnetization direction perpendicular to the plane, the changes from FV to half-valley-metal (HVM), to valley-nonequilibrium quantum anomalous Hall effect (VQAHE), to HVM, and to FV can be induced by strain engineering. It is worth noting that there are no particular valley polarization and VQAHE states for in-plane (IP) magnetic anisotropy. By obtaining the real magnetic anisotropy energy (MAE) under different strains, due to spontaneous valley polarization, intrinsic out-of-plane (OOP) magnetic anisotropy, a chiral edge state, and a unit Chern number, the VQAHE can reliably appear between two HVM states. The increasing strains can induce VQAHE, which can be clarified by a band inversion between d_x²-y²/d_xy and d_z² orbitals, and a sign-reversible Berry curvature. Once synthesized, the Janus monolayer ScBrI would find more significant applications in topological electronic, valleytronic, and spintronic nanodevices.

A Quantum Monte Carlo Study of the Structural, Energetic, and Magnetic Properties of Two-Dimensional H and T Phase VSe2

Previous works have controversially claimed near-room-temperature ferromagnetism in two-dimensional (2D) VSe₂, with conflicting results throughout the literature. These discrepancies in magnetic properties between both phases (T and H) of 2D VSe₂ are most likely due to the structural parameters being coupled to the magnetic properties. Specifically, both phases have a close lattice match and similar total energies, which makes it difficult to determine which phase is being observed experimentally. In this study, we used a combination of density functional theory, highly accurate diffusion Monte Carlo (DMC), and a surrogate Hessian line-search optimization technique to resolve the previously reported discrepancy in structural parameters and relative phase stability. With DMC accuracy, we determined the free-standing geometry of both phases and constructed a phase diagram. Our findings demonstrate the successes of the DMC method coupled with the surrogate Hessian structural optimization technique when applied to a 2D magnetic system.

Collective Magnetic Behavior in Vanadium Telluride Induced by Self-Intercalation

Self-intercalation of native magnetic atoms within the van der Waals (vdW) gap of layered two-dimensional (2D) materials provides a degree of freedom to manipulate magnetism in low-dimensional systems. Among various vdW magnets, the vanadium telluride is an interesting system to explore the interlayer order-disorder transition of magnetic impurities due to its flexibility in taking nonstoichiometric compositions. In this work, we combine high-resolution scanning transmission electron microscopy (STEM) analysis with density functional theory (DFT) calculations and magnetometry measurements, to unveil the local atomic structure and magnetic behavior of V-rich V_1+xTe₂ nanoplates with embedded V₃Te₄ nanoclusters grown by chemical vapor deposition (CVD). The segregation of V intercalations locally stabilizes the self-intercalated V₃Te₄ magnetic phase, which possesses a distorted 1T'-like monoclinic structure. This phase transition is controlled by the electron doping from the intercalant V ions. The magnetic hysteresis loops show that the nanoplates exhibit superparamagnetism, while the temperature-dependent magnetization curves evidence a collective superspin-glass magnetic behavior of the nanoclusters at low temperature. Using four-dimensional (4D) STEM diffraction imaging, we reveal the formation of collective diffuse magnetic domain structures within the sample under the high magnetic fields inside the electron microscope. Our results shed light on the studies of dilute magnetism at the 2D limit and on strategies for the manipulation of magnetism for spintronic applications.

Giant thermal switching in ferromagnetic VSe2 with programmable switching temperature

Active and reversible modulation of thermal conductivity can realize efficient heat energy management in many applications such as thermoelectrics. Using first-principles calculations, this study reports a giant thermal switching ratio of 12, much higher than previously reported values, in monolayer 2H-VSe₂ above room temperature. Detailed analysis indicates that the high thermal switching ratio is dominated by the ferromagnetic ordering induced phonon bandgap, which significantly suppresses the phonon-phonon scattering phase space across the entire vibration spectrum. The thermal switching in bulk 2H-VSe₂ is also investigated and the thermal switching ratio reaches 9.2 at the magnetic transition temperature. Both the phonon-phonon scattering space phase and phonon anharmonicity are responsible for the 9.2-fold thermal switching. This study advances the understanding of heat energy transport in two-dimensional ferromagnets, and also provides new insight into heat energy control and conversion.

Above-Room-Temperature Ferromagnetism in Thin van der Waals Flakes of Cobalt-Substituted Fe5GeTe2

Two-dimensional (2D) magnetic van der Waals materials provide a powerful platform for studying the fundamental physics of low-dimensional magnetism, engineering novel magnetic phases, and enabling thin and highly tunable spintronic devices. To realize high-quality and practical devices for such applications, there is a critical need for robust 2D magnets with ordering temperatures above room temperature that can be created via exfoliation. Here, the study of exfoliated flakes of cobalt-substituted Fe₅GeTe₂ (CFGT) exhibiting magnetism above room temperature is reported. Via quantum magnetic imaging with nitrogen-vacancy centers in diamond, ferromagnetism at room temperature was observed in CFGT flakes as thin as 16 nm corresponding to 16 layers. This result expands the portfolio of thin room-temperature 2D magnet flakes exfoliated from robust single crystals that reach a thickness regime relevant to practical spintronic applications. The Curie temperature T_c of CFGT ranges from 310 K in the thinnest flake studied to 328 K in the bulk. To investigate the prospect of high-temperature monolayer ferromagnetism, Monte Carlo calculations were performed, which predicted a high value of T_c of ∼270 K in CFGT monolayers. Pathways toward further enhancing monolayer T_c are discussed. These results support CFGT as a promising platform for realizing high-quality room-temperature 2D magnet devices.

Salt-Assisted Selective Growth of H-phase Monolayer VSe2 with Apparent Hole Transport Behavior

Vanadium diselenide (VSe₂) exhibits versatile electronic and magnetic properties in the trigonal prismatic (H-) and octahedral (T-) phases. Compared to the metallic T-phase, the H-phase with a tunable semiconductor property is predicted to be a ferrovalley material with spontaneous valley polarization. Herein we report an epitaxial growth of the monolayer 2D VSe₂ on a mica substrate via the chemical vapor deposition (CVD) method by introducing salt in the precursor. Our first-principles calculations suggest that the monolayer H-phase VSe₂ with a large lateral size is thermodynamically favorable. The honeycomb-like structure and the broken symmetry are directly observed by spherical aberration-corrected scanning transmission electron microscopy (STEM) and confirmed by giant second harmonic generation (SHG) intensity. The p-type transport behavior is further evidenced by the temperature-dependent resistance and field-effect device study. The present work introduces a new phase-stable 2D transition metal dichalcogenide, opening the prospect of novel electronic and spintronics device design.

Two-dimensional chalcogenide-based ferromagnetic semiconductors

Two-dimensional (2D) magnetic materials draw an enormous amount of attention due to their novel physical properties and potential spintronics device applications. Room-temperature ferromagnetic (FM) semiconductors have long been pursued in 2D magnetic materials, which show a long range magnetic order down to atomic-layer thickness. The intrinsic ferromagnetism has been predicted in a series of 2D materials and verified in experiments and the magnetism can be modulated by multiple physical fields, exhibiting promising application prospects. In this review, we overview several types of 2D chalcogenide-based FM semiconductors discovered in recent years. We summary and compare their basic physical properties, including the crystal structures, electronic structures, and mechanical stability. The 2D magnetism can be described by several physical models. We also focus on the recent progresses about theoretical prediction of FM semiconductors and experimental observation of external-field regulation. Most of investigations have shown that 2D chalcogenide-based FM semiconductors have relatively high Curie temperature (Tc) and structural stability. These materials are promising to realize the room-temperature ferromagnetism in atomic-layer thickness, which is significant to design spintronics devices.

Massive Monte Carlo simulations-guided interpretable learning of two-dimensional Curie temperature

Monte Carlo (MC) simulation of the classical Heisenberg model has become the de facto tool to estimate the Curie temperature (T _C) of two-dimensional (2D) magnets. As an alternative, here we develop data-driven models for the five most common crystal types, considering the isotropic and anisotropic exchange of up to four nearest neighbors and the single-ion anisotropy. We sample the 20-dimensional Heisenberg spin Hamiltonian and conceive a bisection-based MC technique to simulate a quarter of a million materials for training deep neural networks, which yield testing R ² scores of nearly 0.99. Since 2D magnetism has a natural tendency toward low T _C, learning-from-data is combined with data-from-learning to ensure a nearly uniform final data distribution over a wide range of T _C (10-1,000 K). Global and local analysis of the features confirms the models' interpretability. We also demonstrate that the T _C can be accurately estimated by a purely first-principles-based approach, free from any empirical corrections.

Stress-driven structural and bond reconstruction in 2D ferromagnetic semiconductor VSe2

Two-dimensional (2D) semiconducting transition metal dichalcogenides can be used to make high-performance electronic, spintronic, and optoelectronic devices. Recently, room-temperature ferromagnetism and semiconduction in 2D VSe₂nanoflakes were attributed to the stable 2H-phase of VSe₂in the 2D limit. Here, our first-principles investigation shows that a metastable semiconducting H' phase can be formed from the H VSe2 monolayer through uniaxial stress or uniaxial strain. The calculated phonon spectra indicate the dynamical stability of the metastable H' VSe₂and the path of phase switching between the H and H' VSe₂phases is calculated. For the uniaxial stress (or strain) scheme, the H' phase can become lower in total energy than the H phase at a transition point. The H' phase has stronger ferromagnetism and its Curier temperature can be enhanced by applying uniaxial stress or strain. Applying uniaxial stress or strain can substantially change spin-resolved electronic structures, energy band edges, and effective carrier masses for both of the H and H' phases, and can cause some flat bands near the band edges in the strained H' phase. Further analysis indicates that one of the Se-Se bonds in the H' phase can be shortened by 19% and the related Se-V-Se bond angles are reduced by 23% with respect to those of the H phase, which is believed to increase the Se-Se covalence feature and reduce the valence of the nearby V atoms. Therefore, structural and bond reconstruction can be realized by applying uniaxial stress in such 2D ferromagnetic semiconductors for potential spintronic and optoelectronic applications.

Ferroelectric control of band alignments and magnetic properties in the two-dimensional multiferroic VSe2/In2Se3

Our first-principles evidence shows that the two-dimensional (2D) multiferroic VSe₂/In₂Se₃experiences continuous change of electronic structures, i.e. with the change of the ferroelectric (FE) polarization of In₂Se₃, the heterostructure can possess type-I, -II, and -III band alignments. When the FE polarization points from In₂Se₃to VSe₂, the heterostructure has a type-III band alignment, and the charge transfer from In₂Se₃into VSe₂induces half-metallicity. With reversal of the FE polarization, the heterostructure enters the type-I band alignment, and the spin-polarized current is turned off. When the In₂Se₃is depolarized, the heterostructure has a type-II band alignment. In addition, influence of the FE polarization on magnetism and magnetic anisotropy energy of VSe₂was also analyzed, through which we reveal the interfacial magnetoelectric coupling effects. Our investigation about VSe₂/In₂Se₃predicts its wide applications in the fields of both 2D spintronics and multiferroics.

The Magnetic Genome of Two-Dimensional van der Waals Materials

Magnetism in two-dimensional (2D) van der Waals (vdW) materials has recently emerged as one of the most promising areas in condensed matter research, with many exciting emerging properties and significant potential for applications ranging from topological magnonics to low-power spintronics, quantum computing, and optical communications. In the brief time after their discovery, 2D magnets have blossomed into a rich area for investigation, where fundamental concepts in magnetism are challenged by the behavior of spins that can develop at the single layer limit. However, much effort is still needed in multiple fronts before 2D magnets can be routinely used for practical implementations. In this comprehensive review, prominent authors with expertise in complementary fields of 2D magnetism (i.e., synthesis, device engineering, magneto-optics, imaging, transport, mechanics, spin excitations, and theory and simulations) have joined together to provide a genome of current knowledge and a guideline for future developments in 2D magnetic materials research.

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

The first magnetic 2D material discovered, monolayer (ML) CrI₃, is particularly fascinating due to its ground state ferromagnetism. However, because ML materials are difficult to probe experimentally, much remains unresolved about ML CrI₃'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 CrI₃. We exploit a recently developed surrogate Hessian DMC line search technique to determine CrI₃'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 CrI₃'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 CrI₃ 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 CrI₃'s I-p orbital suggest a ligand superexchange-dominated magnetic anisotropy in ML CrI₃, corroborating recent observations of magnons in its 2D limit. We also find that ML CrI₃ exhibits a substantial spin-phonon coupling of ∼3.32 cm^-1. Our work, thus, establishes many of ML CrI₃'s key properties, while also continuing to demonstrate the pivotal role that DMC can assume in the study of magnetic and other 2D materials.