Charge doping into spin minority states mediates doubling of TC in ferromagnetic CrGeTe3

Introduction
Dimensionality has a profound impact on the physical properties of materials. Control of the effective system dimensionality can be used to tune band gaps in semiconductors1, control catalytic activity2, and manipulate the strength of electronic interactions3. As such, it is of core importance both for fundamental understanding and technological application. This is true no more so than in the field of magnetism. While three-dimensional magnets are commonplace, long-range order is strictly forbidden to occur in one-dimensional systems4. Layered materials hosting magnetic ions present a novel environment in which to study the critical dimensionality between these two extremes5: finite inter-plane coupling allows long-range order to develop, while the quasi-two-dimensionality of the system can be expected to have a strong influence on the magnetic anisotropy and the role of fluctuations.
Here, we study the quasi-2D layered magnet CrGeTe3. In this compound, magnetic layers are separated by van der Waals gaps, leading to weak inter-layer interactions6,7,8. Nonetheless, recent studies employing electrostatic gating have reported both a dramatic increase in the magnetic ordering temperature, TC, and a striking switch in the magnetic anisotropy with increasing gate voltage9, pointing to a subtle interplay between dimensionality and doping effects. We utilise alkali-metal surface doping in ultra-high vacuum to mimic the conditions of field-effect doping, and probe the resulting changes in the surface electronic structure by angle-resolved photoemission spectroscopy (ARPES). Combined with density-functional theory (DFT) calculations, we demonstrate how electrostatic-type doping in this system drives the formation of a half-metallic 2D electron gas of minority-spin carriers, and show how this opens new exchange pathways which underpin the increase in TC.
Results
Doping-dependent electronic structure of CrGeTe3
Figure 1 shows an overview of the electronic structure of pristine CrGeTe3, and the measured electronic structure following the deposition of sub-monolayer coverage of Rb atoms at the surface. In agreement with our DFT calculations (Fig. 1e), ARPES of the pristine compound (Fig. 1a, b) indicates several dispersive states near the top of the valence band, which are dominantly comprised of Te orbital character (see also Supplementary Fig. 1 for orbitally-resolved calculations). These approach, but do not reach, the Fermi level, with negligible spectral weight at EF consistent with the semiconducting nature of this compound. The rather flat band observed at E − EF ≈ − 1.5 eV is evident in our calculations as the spin-majority Cr-derived states. Consistent with previous calculations10, we find very little weight of the Ge states in the valence band here (see Supplementary Fig. 1). While Ge is crucial to obtain the chemically stable Cr3+ charge state in CrGeTe3, it can be considered as electrically inert in the parent compound.

a ARPES dispersions measured (hν = 80 eV) in the magnetic ground state along the ({{rm{K}}}^{{prime} }-Gamma -{rm{K}}-{rm{M}}-{{rm{K}}}^{{prime} }) direction of the Brillouin zone for undoped CrGeTe3, reproduced from Ref. 16. b Corresponding Fermi surface map. Consistent with the semiconducting nature of CrGeTe3, we find no states at the Fermi level. c, d Equivalent to (b, a) for the Rb-doped surface. A clear set of Fermi pockets arises around the K points of the surface Brillouin zone. DFT band structure calculations for (e) undoped and (g) electron-doped CrGeTe3. f Crystal structure of CrGeTe3, showing a schematic of the Rb surface doping.
Upon surface Rb deposition (Fig. 1(d)), the dispersive Te-derived states and the flatter Cr-derived bands are both shifted downwards by around 500 meV (1(d)). This points to an efficient electron doping at the surface, where the deposited Rb atoms become readily ionised, donating their electrons into the near-surface region. We note that similar population of the conduction band has previously been achieved by a similar method via the deposition of K at the surface of CrGeTe311. Consistent with this, we find a finite population of the conduction band states around the K-M-K’ path of the surface Brillouin zone. Our measurements indicate the conduction band minimum to be located at the K points (see also Fig. 2), indicating an indirect band gap in this system. Small triangular electron-like Fermi pockets are created around these zone-corner points, clearly visible in our Fermi surface measurements in Fig. 1c. From photon energy-dependent measurements (Supplementary Fig. 2), we find that these states are non-dispersive along the out-of-plane direction, and are thus two-dimensional. We therefore conclude that the alkali metal dosing here leads to a near-surface band bending which, in turn, drives the formation of a 2D electron gas localised at the sample surface12.

a–c ARPES dispersions measured along the (Gamma -{rm{K}}-{rm{M}}-{{rm{K}}}^{{prime} }) direction of the Brillouin zone for a surface carrier density of (a) N = 3.1 × 1013 cm−2 (hν = 48 eV), (b) N = 4.8 × 1013 cm−2 (hν = 80 eV), and (c) N = 7.9 × 1013 cm−2 (hν = 80 eV). Fermi surface measurements around the M point of the Brillouin zone are shown as insets. The topology of the Fermi surface changes, indicating a doping-dependent Lifshitz transition. d DFT calculated band structure, indicating the spin-minority nature of the conduction band. Dashed lines represent EF for the ARPES measurements shown in (a–c).
Interestingly, our DFT calculations for an electron-doped CrGeTe3 layer (see Methods) show that the lowest-energy conduction band, which becomes populated upon carrier doping (Fig. 1g), is of spin-minority character. We find a similar conduction band structure in the undoped system (Fig. 1e), a calculation where we manually change the electron density (Fig. 1g) and in supercell calculations where we explicitly include a Rb atom (Supplementary Fig. 3), indicating that this is a robust feature of the electronic structure, intrinsic to CrGeTe3. In fact, previous calculations13,14 have come to differing conclusions regarding the nature of the lowest conduction band in CrGeTe3. We find that the exact band ordering, and the spin and orbital decomposition of the lowest conduction band is sensitive to the value of the onsite energy U used for the calculations (see Supplementary Fig. 4), explaining the previous variability in calculated electronic structures. We note that only the spin-minority band has the correct momentum-space dispersion to reproduce our experimental measurements, with band minima located at the zone-corner K points. We thus conclude that the low-U regime provides the most accurate calculation scheme for the electron-doped system here, rather than U = 3 eV as has been previously reported in the literature. Importantly, the calculations indicate how the low-U ultimately facilitates the formation, on doping, of a two-dimensional half-metal of spin-minority character at the surface here.
Figure 2 shows how this half-metallic 2DEG becomes populated with increasing surface carrier doping. At a surface doping density of 3.1 × 1013 cm−2 (0.06 electrons/Cr, Fig. 2a), nearly V-shaped conduction bands are observed centred at the K and K’ points of the Brillouin zone, with the conduction band crossing through the Fermi level along the Γ-K and K-M-K’ lines. The resulting Fermi surface (Fig. 2a, inset) is comprised of closed triangular pockets centred at the K and K’ points. With increasing deposition of Rb on the surface, the conduction band pockets become more occupied, with a concomitant increase in the size of the triangular Fermi pockets (Fig. 2b). At a critical electron doping of around 5 × 1013 cm−2 (0.1 electrons/Cr), the Fermi level moves above a saddle point located at the M point of the Brillouin zone (Fig. 2c). This drives a Lifshitz transition15, where the triangular electron-like Fermi surfaces at low doping merge to form a large hole-like pocket centred at the Brillouin zone centre. These features are well described by our DFT calculations of the electronic structure of a monolayer of CrGeTe3 (Fig. 2d, see Methods, and also Supplementary Fig. 5), indicating that the quantum-confined 2DEG here is well represented within a monolayer description.
Increase of T
C with carrier doping
We now turn to the magnetic order of this surface-doped layer. Previously, some of us have identified a particular state in the valence band of the undoped compound – derived from hybridised Te 5p and Cr eg states – which shows a strong temperature-dependent shift associated with the development of long-range magnetic order (see Fig. 3c, inset)16. This band feature is also evident in our doped samples, for example in Fig. 3a as the bright spectral weight marked by a red arrow at the Γ point at E − EF ≈ − 2.9 eV. As is evident in the raw data in Fig. 3a, and in energy distribution curves (EDCs) extracted at the Γ point (Fig. 3b), this hybridised state shifts away from the Fermi level with decreasing temperature. We monitor the temperature-dependent separation of this hybridised state from the top of the Te 5p-derived valence bands as a fingerprint of the onset of long-range magnetic order here, where the cross-referencing to the valence band top ensures that temperature-dependent chemical potential shifts or residual sample charging due to the semiconducting bulk cannot affect our conclusions. From this (Fig. 3c), we find that the onset of the shift to higher binding energies of the hybridised state starts at around 120 K for the sample shown here with an electron doping of N = 3.1 × 1013 cm−2, considerably higher than for the undoped bulk where such energy shifts onset at the bulk TC of ≈ 65 K. We also note that the valence bands are located at slightly higher binding energy for the undoped compound than the doped compound, reflecting a potential role of quantum confinement from the surface doping here for the top-most Te-derived states.

a ARPES measurements around the Γ point (hν = 48eV) for a CrGeTe3 sample with carrier density N = 3.1 × 1013 cm−2, as a function of temperature, with the relevant states highlighted by red arrows. b Corresponding temperature-dependant EDCs extracted at Γ. c Extracted temperature-dependent separation of the Cr-Te hybridised state from the Te 5p valence band maximum (see arrows in (a) and (b)). Grey and black dashed lines represent an estimated baseline and a 10% rise from that baseline. Approximate values for TC are given by the intersection of each data set with the black dashed line, indicated by the coloured arrows. The magnetic TC for pristine CrGeTe3 is also shown, and the corresponding temperature-dependent band shifts are shown inset, reproduced from Ref. 16. d Experimental phase diagram, plotting the extracted TC estimates as a function of doping along with estimates from transport measurements from Ref. 9.
Irrespective, we find that the onset of the energetic shifts of the Cr-Te hybridised state here is systematically moved to higher temperatures with increasing surface doping level (Fig. 3(c)). This indicates that the induced carriers are mediating an increased TC here, consistent with both studies of the bulk where carriers are generated by organic ion intercalation with tetrabutyl ammonium17, and with the observation of an enhanced TC in electrostatically-doped CrGeTe39. While TC is clearly enhanced for the doped surfaces here, extracting a quantitative measure of TC is challenging as the transition is not particularly sharp in our measurements, reflecting the complexity of temperature-dependent measurements for an alkali-doped surface. To allow making robust comparisons between the samples, we therefore estimate a baseline of the binding energy at high temperature in our Evs.T curves (grey dashed line from Fig. 3c) and determine a TC value as the temperature at which our measured band shifts cross a threshold at 10% of the total binding energy shift (black dashed line from Fig. 3c). From this, we extract the doping-dependent phase diagram shown in Fig. 3d. We observe a rapid increase in TC from the bulk value already with our lowest electron doping. TC then continues increasing with increasing doping density, also in line with the results from electrostatic gating which we include in our phase diagram9, but with a much slower rate. It thus appears that the introduction of free carriers in the relatively low-doping regime, at and around the point where the spin-minority states are populated and the van Hove singularity in the conduction band is traversed, are the most important for mediating a TC-enhancement in this system.
Mechanisms of enhanced magnetic order
To explore the origins of this further, we study the evolution of the exchange parameters, Ji, as a function of doping from our DFT calculations, considering up to third-nearest neighbour terms (i = 3, Fig. 4a, b). Interestingly, we find modest increases in both the nearest and next-nearest neighbour terms with increasing electron doping, and a rather dramatic increase in J3 (Fig. 4c). We can understand this from the nature of the conduction band states which become populated upon electron doping. In an octahedral configuration, the naive choice for the next orbital to become occupied upon electron doping would be a spin-majority eg state. Indeed, there is a significant amount of unoccupied eg spectral weight peaked at around 0.7 eV above the band gap (see Supplementary Fig. 6), and virtual hopping processes into these states are thought to be relevant for the superexchange in the undoped system16. However, our calculated charge density distribution about the Cr sites for the lowest conduction band (Fig. 4a) exhibits lobes pointing between the neighbouring Te sites, indicating a t2g-like character (where we are considering within the octahedral basis, i.e., where the x, y, z axes are oriented along the octahedral axes defined by the Cr-Te bonds). This is consistent with the spin-minority nature of the states at the bottom of the conduction band, as well as a direct projection of its orbital character from our calculations (Supplementary Fig. 6 and Supplementary Fig. 7). This band, however, sits well below the majority of the unoccupied t2g weight18 (see also Supplementary Fig. 6), and is surprisingly dispersive, dropping below the eg states at the K and M points. In fact, our calculations reveal a notable charge density of this band also around the Ge sites. This points to a significant hybridisation of the states at the bottom of the conduction band with Ge, which had typically been thought to play a relatively inert role in the magnetic properties of (undoped) CrGeTe3. Considering the J3 exchange pathway, however, (Fig. 4(b)), it is clear how in the electron-doped system, hopping via the Ge sites could mediate a significant increase in J3. We note that for doping into the spin-majority eg-like states, a much more modest involvement of the Ge states can be expected, as evident from our calculated charge density for these states shown in Supplementary Fig. 7.

a Calculated charge density (yellow surface), integrated over an energy range of 0.2 eV around the spin-minority conduction band minimum and superimposed onto the crystal structure of CrGeTe3 and (b) corresponding top-down view of the surface layer. The lobes of the charge density around each Cr site point in between the Cr-Te bond axis, a signature of t2g character, while the charge density around the Ge dimers is dramatically enhanced, pointing to a new conduction path involving the Ge sites. J1, J2, and J3 paths are labeled with colour coded arrows. c Calculated values for J1, J2, and J3, plotted as a function of doping concentration. The absolute value of J3, in particular, shows a dramatic enhancement with doping. d TC extracted from Monte Carlo simulations, considering different combinations of J1,2,3, plotted alongside the experimental data. Here, J3 is required (and has the equivalent impact of J1 and J2 combined) to match the shape of the experimental phase diagram. The dashed line through the experimental data points is the calculation for J1 + J2 + J3 offset by 45 K, as a guide for the reader in showing the relative trends between experiment and simulation. e Schematic representation of the mechanism to enhance TC enabled by the doped carriers in the minority-spin t2g manifold. Hopping paths from Cr to Te include direct Cr-Te-Cr hopping (J1), and higher-order paths (J2) as well as effective Cr-Cr exchange pathways that are mediated via the Ge sites (i.e. J3).
We show in Fig. 4d an estimate from Monte Carlo calculations (see also Supplementary Fig. 8) of the doping-dependent TC that would be expected, given such a change in exchange parameters as calculated here. TC increases monotonically and approximately linearly in the low-doping regime. While modest changes in TC can be expected for very low dopings, in the doped-semiconductor regime, when the doping becomes sufficient that the system is metallic, a significant increase in TC can be expected, consistent with our experimental measurements. Indeed, considering the evolution of all calculated exchange constants (purple line in Fig. 4d), we find that the calculated TC slightly underestimates the experimental values both for the bulk and the doped systems, but with a doping-dependent trend which closely follows our experimental results (see also offset calculation, purple dashed line in Fig. 4d). Excluding either the increase in J3 (green line in Fig. 4d) or in J1 and J2 (red line in Fig. 4d), it is clear that approximately half of the increase observed here originates from the new hopping pathway which is opened up via the Ge dimers, shown in Fig. 4a and b.
The effect of magnetic anisotropy changes vs doping was included in the spin-Hamiltonian used in the Monte Carlo simulations to estimate the Curie temperature, indicating a doping-dependent magneto-crystalline anisotropy (MCA) which changes sign with increasing carrier density. While the precise doping density at which this occurs will be sensitive to the detailed parameters used in the calculations19, this general doping-dependent trend is in good agreement with previous calculations9 (see Methods and Supplementary Fig. 9). However, despite generally being a very important microscopic ingredient in the magnetism of 2D systems20, in this specific case, the magnetic anisotropy seems not to play a major role. The trend of the Tc (estimated including magnetic anisotropy terms) closely follows that of J1 + J2 + J3 isotropic Heisenberg-only terms. Moreover, we performed Monte Carlo calculations where we fix the magnetic exchange parameters as calculated for the undoped system, but change the MCA in accordance with its doping-dependent evolution: as shown in Supplementary Fig. 10, we find no significant change in the predicted TC.
Discussion
We attribute the resulting increase in TC here to an unusual form of double exchange, mediated via a dispersive branch of the spin-minority t2g carriers, as summarised schematically in Fig. 4(e). In a simple picture, we can consider the undoped electronic structure of CrGeTe3 as having a half-filled t2g manifold in a d3 electron configuration. (We note that, in reality, even for the undoped compound, finite mixing with the ligand states plays a crucial role, ultimately underpinning the ferromagnetic order via a superexchange mechanism16.) The doped carriers must therefore fill either the majority-spin eg states or the minority-spin t2g states. Our DFT calculations above indicate the latter occurs here at low doping, thus leading to an additional Cr spin on some sites which is anti-aligned with the majority spin moment. Given the dispersive nature of the t2g-derived state at the bottom of the conduction band here, this doped carrier can readily delocalise through the crystal, hopping via the Te sites (including the Cr-Te-Ge-Te-Cr hopping pathway identified above). This, however, requires the creation of a hole on the Te site that has equal spin as the doped electron, which thus acts to mediate a net ferromagnetic coupling here (Fig. 4e).
Since the doped carriers occupy the minority-spin t2g states, the net magnetic moment will be partially cancelled. Indeed, as shown in Supplementary Fig. 11, we find that the magnetic moment decreases linearly with the increase of doping concentration at low doping levels. As the doping concentration increases further, however, the spin-majority eg states start to become populated (Fig. 2d), resulting in fluctuations in the magnetic moment (Supplementary Fig. 11). Consistent with this, we find that the increase in J and the associated increase in calculated TC starts to saturate above a critical doping of around 0.3 electrons/Cr (Fig. 4c, d). These trends are in good agreement with experiments (Fig. 3d), validating our DFT approach. It will be an interesting subject for further study to consider whether beyond-DFT methods can obtain even better quantitative accuracy with our measurements. Already here, our study highlights how the population of spin-minority conduction band states with electron doping in CrGeTe3 mediates a more than doubling of its ferromagnetic TC, explaining recent findings of gate-tuned magnetism in this system. Our findings also highlight an unusual situation that occurs here under such electrostatic doping, with a half-metal forming at the surface which has opposite spin-polarisation to the bulk, thus suggesting potential for novel spintronic devices such as intrinsic, and gate-tunable, spin valve structures, magnetic tunnel junctions, or spin-filters21.
Methods
Sample growth
CrGeTe3 single crystals were prepared following the procedure in Ref. 7. The crystals were separated from the flux by centrifuging at 500∘C and then cleaned of any excess Te. The crystals were found to exhibit a sharp magnetic transition at approximately 65 K16, in good agreement with previous literature reports for pristine bulk CrGeTe38.
Angle-resolved photoemission spectroscopy
ARPES measurements were performed at the high-resolution branch of the I05 beamline at Diamond Light Source, UK. Bulk single crystal samples were cleaved in situ after being cooled to approximately 40 K, and then electron-doped by depositing Rb atoms on the surface via evaporation from a SAES getter source. To obtain variable electron density, the samples were dosed for varying time, between 30 s and 360 s. The chamber pressure did not exceed 7 × 10−10 mbar during the deposition. The alkali surface doping leads to an increase in carrier concentration only close to the surface, with a spatial extent governed by the Thomas-Fermi screening length, and an associated near-surface band bending. We find that the doped surface layer is in reasonable approximation to a monolayer-like system. Note, we would expect qualitatively similar findings as presented in the main text if other alkali metals were used, with the primary role of the alkali metal being to electron dope into the near-surface region.
ARPES measurements were performed using linearly polarised photons (with photon energies as specified in the figure captions) at sample temperatures ranging between 40 K and 200 K. Since CrGeTe3 is a semiconductor, sample charging due to photoemission can become problematic at low temperatures22. To minimise the impact of this, we limited our measurement range to above 40 K and used rather low photon fluxes. To further rule out any influence of sample charging on our conclusions, we extract TC based on analysis of energy differences in our measured spectra, as outlined in the text.
The doping levels were estimated from our measured ARPES data using Luttinger’s Theorem23, where the volume enclosed by the Fermi surface is proportional to the particle density by,
where N is the number of charge carriers, D is the dimensionality of the system, and g is the spin and valley degeneracy, characterising how many equivalent bands exist within the first Brillouin zone, with D = 2 here. The area enclosed by the Fermi surface was extracted from our fitted Fermi surface data by implementing a Green’s Theorem algorithm.
Density functional theory and Monte Carlo calculations
Our calculations were carried out within the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation24, a commonly-used approximation to treat the exchange-correlation potential term within DFT, here as implemented in the Vienna ab Initio Simulation Package (VASP)25. To better consider the vdW forces, we performed calculations with a non-local vdW density functional in the form of optB86b-vdW26,27. The interactions of electrons with ionic cores were described using projector-augmented wave (PAW) pseudopotentials28,29. For all the calculations, we chose the energy cutoff to be 600 eV, the maximum force of all atoms was optimized until it fell to less than 0.01 eV/Å, and the criterion for the total energy convergence was set as 10−6 eV. The set of 12 × 12 × 1 k-point samplings30, in the reciprocal space by the Gamma-centered grid, was used to model the Brillouin zone. A vacuum space exceeding 20 Å was included to avoid interaction between adjacent periodic layers for the monolayer calculations.
In this work, all of the geometry structures were fully relaxed, according to the maximum force and total energy convergence, listed previously. The lattice constant a and b of optimized monolayer CrGeTe3 are 6.90 Å, which is close to 6.91 Å reported by Fang et al.31. We calculated the energy of different dopant positions and found that the configuration in which the Rb is located directly above the Cr sites has the lowest energy (Supplementary Fig. 12). This superstructure was then used to calculate the band structure. Our calculations indicated that the effects of electron doping are dominated by a rigid shift of the band structure. For the calculation of the exchange parameters and magnetic moments as a function of electron doping, we thus artificially changed the number of electrons to simulate the non-integer doping concentration. This allows us to calculate a continuously varying doping concentration, including the regime of low doping levels, which would be unfeasible for direct supercell calculations with chemical doping. We have validated this approach against supercell calculations for high carrier densities, where the supercell with chemical dopants becomes of a feasible size. Throughout, we refer to the supercell calculations as Rb-doped in the caption, while the DFT results calculated by artificially changing the electron number are referred to as electron-doped. We note that these doping calculations were performed without spin-orbit coupling. Whilst spin-orbit coupling leads to significant changes for the Te-derived valence bands16 see Supplementary Fig. 13, we confirmed that they have minimal impact on the conduction band electronic structure and on the calculated exchange parameters, as shown in Supplementary Fig. 14.
We used the VAMPIRE32 code to perform Monte Carlo (MC) simulations. The magnetic properties can be modelled via the following spin Hamiltonian for classical spins:
assuming isotropic exchange constants up to the third nearest neighbour. Anisotropic effects were modelled via an effective magneto-crystalline anisotropy (MCA) term, with coupling constant Ai, obtained upon including spin-orbit coupling in the DFT calculations. The J terms were extracted by the energy mapping method33, with a negative value indicating ferromagnetic interactions. The energy mapping method is a method to solve the magnetic exchange parameters within a Heisenberg model by combining the energies of the magnetic configurations. Here we calculated the energy of the four configurations shown in Supplementary Fig. 15 through DFT, and substituted them into the left side of the energy expression shown. Then one can obtain J1 − J3 by combining the four equations. The unit cell we use is rectangular, which is consistent with periodic boundary conditions of the honeycomb lattice and is able to host – within the same unit cell – four different spin configurations. The standard Metropolis algorithm34 has been used for MC simulations with 5 × 104 MC steps for equilibration and 3 × 105 MC steps for averaging. Since the VAMPIRE software can only simulate cuboid cells, we used cuboid cells, which contain 4 Cr sites. We performed calculations for a 40 × 40 × 1 supercell with Ns = 6400 spins. The transition temperature can be estimated from the peaks that appear in the temperature evolution of specific heat (Supplementary Fig. 8), evaluated as:
where E is the energy calculated using VAMPIRE, kB is the Boltzmann constant and β = 1/kBT, while 〈. . . 〉 indicates statistical averages.
Responses