Valence band electronic structure of the van der Waals ferromagnetic insulators: VI3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_3$$\end{document} and CrI3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_3$$\end{document}

Ferromagnetic van der Waals (vdW) insulators are of great scientific interest for their promising applications in spintronics. It has been indicated that in the two materials within this class, CrI3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_3$$\end{document} and VI3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_3$$\end{document}, the magnetic ground state, the band gap, and the Fermi level could be manipulated by varying the layer thickness, strain or doping. To understand how these factors impact the properties, a detailed understanding of the electronic structure would be required. However, the experimental studies of the electronic structure of these materials are still very sparse. Here, we present the detailed electronic structure of CrI3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_3$$\end{document} and VI3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_3$$\end{document} measured by angle-resolved photoemission spectroscopy (ARPES). Our results show a band-gap of the order of 1 eV, sharply contrasting some theoretical predictions such as Dirac half-metallicity and metallic phases, indicating that the intra-atomic interaction parameter (U) and spin-orbit coupling (SOC) were not properly accounted for in the calculations. We also find significant differences in the electronic properties of these two materials, in spite of similarities in their crystal structure. In CrI3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_3$$\end{document}, the valence band maximum is dominated by the I 5p, whereas in VI3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_3$$\end{document} it is dominated by the V 3d derived states. Our results represent valuable input for further improvements in the theoretical modeling of these systems.

www.nature.com/scientificreports/ Magnetic measurements show the presence of magnetocrystalline anisotropy (MCA) in these systems, favoring an out-of-plane orientation of magnetic moments 18,22,25,32 . It is known that MCA in these systems originates due to the strong SOC induced by the iodine ligand 19,20,30 . The experimental studies of the electronic structure in this class of materials were restricted to Cr 2 Ge 2 Te 6 and Cr 2 Si 2 Te 6 and to the best of our knowledge, no reports are available on CrI 3 and VI 3 12,13,15 . Even though the magnetic and optical measurements indicate that the strong SOC and U play an important role, the band structure calculations often do not take into account these parameters, erroneously predicting properties such as Dirac half-metallicity or metallic phases 17,29,42 . The calculations that do account for U predict Mott-insulating ground state, both in monolayer as well as bulk VI 3 , in agreement with the experiments 23,25,43 .
Webster et al. 17 have calculated the band structure of a monolayer CrI 3 for non-magnetic, FM and AFM phases using local density approximation (LDA) and generalized-gradient approximation (GGA). These calculations show that the introduction of magnetism leads to metal-to-semiconductor transition. In addition, they indicate that the nature of band gaps (direct or indirect) in magnetic phases shows a delicate dependence on the magnetic ordering and spin-orbit coupling. Zhang et al. 44 have performed calculations for bulk and monolayer CrX 3 (X=F, Cl, Br and I) using two different methods. Their results show that a monolayer and bulk have very similar partial-and total DOS but the relative contributions within the DOS and energy positions differ significantly between the two methods. Orbital resolved band dispersions for a monolayer CrI 3 is reported by Jiang et al. 45 and Kim et al. 46 . They show that the the valence band maximum (VBM) is dominated by I 5p orbitals (mainly p x and p y ) and have also discussed the effects of SOC in the band dispersions. The density functional theory (DFT) calculations 19 using DFT+U also show the dominance of I 5p orbitals in the VBM of CrI 3 and highlight the importance of SOC for the observed magnetic anisotropy. Moreover, depending on the calculation method and choice of parameters, the resulting electronic structure of the same material may differ significantly, leaving huge discrepancies in predicted material's properties [43][44][45]47,48 . These uncertainties emphasize the importance of experimental studies of the electronic structure of these materials.
Here, we present the first experimental studies of the valence band electronic structure of CrI 3 and VI 3 using ARPES. Even though both materials have a similar crystal structure, our results indicate significant differences in their electronic structure. In CrI 3 we find that the VBM is dominated by I 5p contribution, while in VI 3 the VBM is mainly of V 3d character. A comparison of the experimental band dispersions with the theoretical ones indicates that even the monolayer band structure captures most of the features seen experimentally. Our results illustrate the importance of SOC and U in these materials, ruling out various metallic phases predicted by theories in which these parameters were neglected. Our study clearly shows the semiconducting/insulating behavior, with the band gap in the 1 eV range, and provides useful input for the improvements in the theoretical understanding of these systems. Figure 1 shows electronic band dispersions of CrI 3 at room temperature (300 K) along the high symmetry directions of surface Brillouin zone (SBZ). Band dispersions along the Ŵ -M-K-Ŵ path is shown in Fig. 1a. To understand the orbital character of these bands, the calculated band structure of a CrI 3 monolayer 45 considering FM ground state and SOC is shown on top of the experimental one in Fig. 1b. The details of orbital character can be found elsewhere in the supplimentary material of ref. 45 . The experimental VBM and the theoretical one were aligned to be at the same energy. Even though our data are from the cleaved bulk crystal, a qualitative comparison with a monolayer should be appropriate because the layers in the crystal are weakly coupled and their electronic structure does not depend significantly on the thickness 44 . Indeed, the overall features in the experimental electronic structure moderately agree with the theory. While the experimental bands are not sharp enough to identify all the fine features existing in the calculations, the measured band dispersions generally agree with the calculations 45 . Some discrepancies may arise due to the absence of U in the calculations. The ARPES spectra are intrinsically broad due to the existence of a large number of overlapping bands 44 and a relatively high temperature of the measurements. The broad, less dispersive features at ∼ 2.5 eV are mainly of Cr 3d character, according to calculations 45,46 . Around the Ŵ point, within the energy range 1.3-2.1 eV, a strongly dispersive feature is attributed to iodine 5p orbital 45,46 . In iso-energy contour, this band forms a hole-like circular feature, centered at Ŵ , as shown in Fig. 1c. This feature resembles the theoretical Fermi surface of the hole-doped CrI 3 monolayer with an out-of-plane magnetization 45 . Very similar features are reported for Cr 2 Ge 2 Te 6 12,12,13 , a material in the same class, which has predominantly the Te 5p character at the VBM. The iso-energy contour at 2.9 eV binding energy shows hexagonal symmetry, representing the symmetry of the CrI 3 surface (Fig. 1d). Figure 1f shows the second derivative of the ARPES intensity ( d 2 I/dE 2 ) along the Ŵ -M , measured using the He-Iα photons. The data were acquired at even higher temperature, 370 K, to minimize the charging. Similar to Fig. 1a, we see the strong bands at ∼ 2.5 eV originating from the Cr 3d. In addition, above those, the three more bands can now be identified more clearly. These three bands are mainly of I 5p character 17,19,45,46 around Ŵ and agree moderately with the bands calculated using the generalized gradient approximation (GGA) and accounting for the SOC and the ferromagnetic ground state with the out-of-plane moments 17,45 . The calculated bands 17,45 with antiferromagnetic and non-magnetic ground states do not agree with the experimental results. The ferromagnetic ground state with the out-of-plane moments is also in line with the magnetic measurements 18 . Direct comparison of experimental data with the calculated bands 17 for a monolayer and bulk (FM state) CrI 3 is shown in Fig. 1g. It appears that the positions, as well as the shapes of these bands agree better with the calculations for a monolayer CrI 3 . It is somewhat surprising that our data match moderately well with the calculations for the FM-ground state, even though the experiments were performed far above the Curie temperature, T C . This could indicate that short-range FM fluctuations might already be at play, while the long-range order is not established, as expected for a low-dimensional system 49  www.nature.com/scientificreports/ et al. 39 have also pointed out to the possibility of the presence of short-range magnetic correlations at much higher temperatures than T C . Recently, electron spin resonance (EPR) and magnetic measurements also suggest similar possibilities in other vdW magnetic materials including VI 3 3,12,23,50-52 . However, we cannot completely rule out the non-magnetic state at high temperature because at present there are no nonmagnetic calculations that account for SOC and U. Without accounting for these, the calculations result in a metallic ground state, contrasting the observed insulating character of CrI 3 . Further theoretical and experimental studies are needed for the proper understanding of this issue. We also note that the band structure calculations without SOC do not agree with our results 17,44,45 , indicating that the SOC is essential for the proper description of this system. This is also in line with the observed magnetic anisotropy, as the two-dimensional long-range ferromagnetic order in this system is stabilized by the magnetic anisotropy originating from the SOC [18][19][20] . We also remark that the calculations based on the hybrid-function (HSE06) place the VBM at Ŵ , in agreement with our results, but the position of the Cr 3d bands (w.r.t. VBM) is much deeper ( ∼ 3.0 eV) in these calculations 44 than in the experiment ( ∼ 0.85 eV, see Fig. 1e) or in the GGA calculations ( ∼ 0.7 eV) 44,45 , LDA+U ( ∼ 0.8 eV) 47 . However, the predicted band gap of ∼ 1.22-1.28 eV 44,47 is much closer to our results and the optical studies ∼ 1.24 eV 38 than the one obtained from GGA ∼ 0.7 eV 44,45 . It appears that the HSE06 calculations better predict the band gap, but underestimate the position of the Cr 3d bands w.r.t. the VBM. Now, we turn to VI 3 . Figure 2 shows the electronic structure of VI 3 along the high symmetry directions of SBZ in the FM phase, at 35 K. The band dispersions along the Ŵ -M-K-Ŵ path are shown in Fig. 2a. The iso-energy contours at different energies and their second derivatives are shown in Fig. 2b,c, respectively. Figure 2d,e represent the LEED patterns obtained from the VI 3 surface at 200 K. The DOS obtained by integrating the photoemission intensity over the whole Brillouin zone and its derivative are plotted in Fig. 2f. The bands with a predominant V 3d and I 5p character are marked in Fig. 2a. Below the V 3d-dominated flat bands, three more dispersing bands, similar to those observed near the VBM in CrI 3 , are clearly resolved. The white arrows point to the two nearly degenerate bands marked by I 5p, seen in both the Ŵ -M and Ŵ -K directions. These two bands are better resolved in the second derivative plots shown in Fig. 4a,b. The iso-energy contour of d 2 I/dE 2 at ∼1.6 eV shows these two bands forming two sheets, as indicated in Fig. 2c by the white arrows. Theoretically, very similar isoenergy features were predicted to form by Te 5p states for the prototype system, Cr 2 Ge 2 Te 6 15 . At higher binding energies, 1.85 and 1.97 eV, these two bands are getting closer and form an electron-like contour with the local minimum at K . Another hole-like contour centered at Ŵ is also identified. Interestingly, all iso-energy contours show a 3-fold symmetry, reflecting the symmetry of the surface. This is also observed in our LEED results. From Fig. 2d it is clear that two sets of spots (marked with yellow and pink circles) display different intensities. For www.nature.com/scientificreports/ the electron energy of 228 eV, only one set of spots are seen (Fig. 2e), suggesting the trigonal symmetry of the crystal (space group R 3 ) 26,27 .

Results and discussions
Here, we point out that the recent theoretical studies have demonstrated that AFM ordering in a bilayer CrI 3 breaks both time-reversal (T) and inversion (P) symmetries, leading to a giant nonreciprocal second-harmonic generation 35 and a magnetic photogalvanic effect 36 . In the underlying crystallographic structure, the inversion symmetry is preserved, but the AF ordering of the two FM layers within a bilayer breaks it, eliminating the need for the non-centrosymmetric crystal. As a consequence of this symmetry breaking, the electronic states at k and −k in the momentum space are not equivalent. The ARPES would be an ideal experimental probe to directly detect such an effect in an AF bilayer. In the present case, the FM order breaks the time-reversal symmetry, but the inversion symmetry is preserved in bulk in both materials. However, the inversion symmetry is always broken at any surface. Therefore, we might expect that in the FM state, both symmetries are broken in the surface region in both CrI 3 and VI 3 , leading to a situation similar to the one predicted for a CrI 3 bilayer. An exciting possibility is that the observed 3-fold symmetry of photoemission intensity in Fig. 2b,c might be partially caused by P − T symmetry breaking near the surface.
We note that the overall band dispersions in VI 3 are in general agreement with the reported calculations 23 as shown in the supplementary material (Fig. S1). However, as most of the calculations do not account for SOC, a detailed comparison would be somewhat inappropriate. A closer look of the VBM region in Fig. 2a indicates that the position of the VBM is not at Ŵ , but at k � ≃ 0.23 Å −1 (marked by red arrows) in both Ŵ -M and Ŵ -K directions. For better visualization of the VBM position, we show the set of corresponding energy distribution curves (EDC) in the supplementary material (Fig. S2). As a consequence of this displacement, the iso-energy surface at ∼ 0.85 eV shows a nearly circular feature around Ŵ , as can be seen in Fig. 2b,c. This is in sharp contrast with the reported theoretical calculations 23,26 , where VBM is located at Ŵ . Furthermore, the absence of on-site Coulomb interactions (U) erroneously predicts a metallic ground state 26,43 . Our results show a clear band gap, in line with the optical gap of ∼ 0.6-0.67 eV observed in this system 25,26 . These results clearly indicate the importance of on-site Coulomb interaction for a proper description of this system.
Another important observation contradicting the theoretical reports is that the VBM in VI 3 is dominated by the V 3d character and not by I 5p as was the case in CrI 3 43,48 . A qualitative assessment of the orbital character can be made by comparing the electronic structure of these two materials (Figs. 1a, 2a) where the more localized 3d orbitals of a transition metal display less dispersion. In the case of CrI 3 these correspond to the features observed around ∼ 2.5 eV, whereas in VI 3 a flat band at ∼ 1.0 eV, most probably corresponds to V 3d related states. The dominant V 3d character of the VBM can be also affirmed by comparing the photoemission intensities of the spectra taken with He-I α and He-II α . It is well-known that photoemission cross-sections 53 of V 3d states are comparable at these two photon energies, whereas I 5p states should be an order of magnitude more intense at He-I α than at He-II α . This would result in a significantly stronger intensity from the V 3d dominated states in the spectra taken at He-II α . Figure 3 shows that this is indeed the case. Figure 3a,b are the ARPES spectra around the Ŵ point recorded by using He I α and He II α photons, respectively. The corresponding energy distribution curves (EDCs) at the marked momenta are plotted in Fig. 3c. If normalized to the state closer to the Fermi energy, the spectral intensity of the second peak ( ∼ 1.75 eV) is quite suppressed at He-II α , indicating that the former peak is dominated by V 3d, while the latter is dominated by I 5p orbitals. We note that the similar dependence of spectral intensities on the light excitation can be also observed in CrI 3 , where the I 5p-derived states are better visible in the He-I α spectra, Fig. 1f,g. In contrast, the calculations show that the VBM is dominated by iodine www.nature.com/scientificreports/ states in both materials 43,48 . We note that in some other vanadium based vdW systems, such as VSe 2 and VTe 2 , the less-dispersive V 3d states are also at the VBM, while the strongly dispersive Se 4p and Te 5p states are just below it, similar to our results from VI 3 54,55 . Similarly, the iodine-derived bands associated with VI 3 in the VI 3 /CrI 3 heterostructure 33 have their maxima at Ŵ and then disperse downward in both Ŵ-M and Ŵ-K directions, resembling the states just below the flat band in our experiment. We believe that the observed discrepancy between the current theory and experiment in VI 3 could be eliminated by the proper optimization of theoretical parameters, including the exact symmetry of the crystal structure, SOC, U, etc.
The electronic band gap in VI 3 in our study is ≥ 0.85 eV (as the position of the Fermi level is ∼ 0.85 eV above the VBM) which is higher than the optical band gap, 0.6-0.67 eV 25,26 . The theoretical band gap varies from 0.9 eV for bulk to the 0.97 eV for monolayer 23 . Our results would place the Fermi level very close to the conduction band minimum (CBM), suggesting a slightly electron-doped system. This may reflect a slight excess of vanadium in the crystal that is thought to partially occupy the vacant interstitial sites within the honeycomb V layers 26 that act as electron donors in this system. Furthermore, to study the effects of structural phase transition and magnetization, we have performed ARPES measurements at room temperature (300 K) and 35 K. Figure 4a,b represent the second derivative of photoemission intensity, indicating band dispersions below and above magnetic transition along the cut indicated by dottedline (red) in Fig. 2b. Momentum distribution curves (MDC) and EDCs at binding energy of 1.6 eV and at k =0 Å −1 are presented in Fig. 4c,d, respectively. We do not observe major changes in the electronic structure between these two temperatures, except for the fact that the two nearly degenerate bands originating from I 5p states are better resolved at 35 K than at 300 K and a slight increase in their separation at low temperature ( Fig. 4c and  supplementary material, Fig. S3). In Fig. 4d, the EDCs show that the peak at 1.7 eV also gets sharper and shifts towards the higher binding energy in the FM phase. This is similar to what has been recently observed in another vdW ferromagnet, Cr 2 Ge 2 Te 6 , where the FM transition induces a band-width enhancement of Te 5p 13 and Cr t 2g 56 related states, along with the energy lowering of 5p-e g hybridized states. We note, however, that relatively minor effects of ferromagnetic transition on the electronic structure of these compounds might not be surprising. Without the applied magnetic field, the magnetic domains should be randomly oriented and their zero net magnetization would not produce a significant effect on the energy bands. As already noted for CrI 3 , this could also mean that the short-range FM correlations might be present above T C . We also do not detect any change related to the structural transition near 80 K. This is also not surprising as the structural transition involves only the changes in the stacking order along the c-axis while the intra-layer structure remains intact 23,25,43 . The weak vdW interaction between the layers and the high surface sensitivity of ARPES would naturally make the changes related to the different stacking of layers unnoticeable.
In summary, we have explored the valence band electronic structure of ferromagnetic vdW materials, CrI 3 and VI 3 using ARPES and identified the bands with dominant V 3d, Cr 3d and I 5p orbital characters. Our results highlight the importance of correlation effects and SOC in these systems and point out to the drawbacks of some calculation methods. Both the optical and transport measurements 8,26 show a semiconducting character of these materials, in good agreement with our results. Moreover, our results show the quasi-two dimensional behavior of the electronic structure in these materials. The future experimental studies of the electronic structure of thin flakes or films would offer a better insight into the development of dimensionality in these materials. Also, this would allow to study some truly exotic phenomena, predicted to occur at some thicknesses, including a giant nonreciprocal second-harmonic generation 35 and the creation of dc current by a linearly polarised light 36 in an AFM bilayer of CrI 3 . www.nature.com/scientificreports/

Methods
The experiments within this study were performed in an experimental facility that integrates oxide-molecular beam epitaxy (OMBE), ARPES, and scanning tunneling microscopy (STM) in a common ultra-high vacuum (UHV) system 57 . Bulk CrI 3 and VI 3 single crystals were grown by the chemical vapor transport method (CVT). The details of sample growth and characterization can be found elsewhere 21,27 . The as-grown samples were clamped to the sample holder under argon atmosphere, cleaved with kapton tape inside UHV, and studied by ARPES. The photoemission experiments were carried out on a Scienta SES-R4000 electron spectrometer with the monochromatized He-I α (21.2 eV) and He-II α (40.8 eV) radiation (VUV-5k). The total instrumental energy resolution was ∼ 5 meV and ∼ 20 meV at He-I α and He-II α , respectively. Angular resolution was better than ∼ 0.15 • and 0.4 • along and perpendicular to the slit of the analyzer, respectively. In-situ low-energy electron diffraction (LEED) measurements were conducted after ARPES measurements to establish the orientation and surface quality. Due to the insulating behavior of the measured samples, the charging was present in the photoemission experiments at the low temperature for VI 3 and even at 300 K for CrI 3 , resulting in energy shifts of photoemission intensity. In the case of VI 3 the charging was not present at 300 K, neither at He-I α nor at He-II α . At 35 K, there was only a ∼ 20 meV shift. Therefore, the VBM recorded at 300 K was taken as a reference and the 35 K data were aligned with that. Similarly, for CrI 3 , where charging was more severe, the VBM measured at 370 K was taken as a reference and the 300 K data were aligned accordingly. When charging was present, it was steady with time.