Introduction

Since the elaboration of the way of its preparation1, manganese-doped gallium arsenide, (Ga,Mn)As, has been the most intensively studied dilute magnetic semiconductors for two interrelated reasons2,3. First, Mn ions substituted for Ga act as an acceptor and provide holes, which mediate the ferromagnetic interaction among localized Mn moments. The presence of randomly distributed acceptors leads to a substantial disorder and to metal-to-insulator transition (MIT) in the range of hole densities relevant to ferromagnetism. Interplay of carrier-mediated exchange interactions and quantum Anderson-Mott localization results in striking properties, often difficult to describe quantitatively either analytically or by the state-of-the-art ab initio methods. This has opened a room for extensive qualitative debates on the character of states in the vicinity of the Fermi energy and on the mechanism accounting for ferromagnetism in this system4,5.

Second, because of the novel carrier-induced nature of ferromagnetism as well as the compatibility with well-established GaAs-based devices, (Ga,Mn)As and related compounds have become a platform on which a number of new concepts of spintronic devices has been demonstrated2,3. In one kind of these devices carrier-concentration dependent ferromagnetism (e.g., spin-FET6 or spin p-i-n diode7) or the presence of spin current (e.g., spin-LED8,9 and related spin-injection devices10, domain-wall displacement11) was exploited. In another type of functional structures the direction of the magnetization vector was altered12,13,14,15 via spin-orbit coupling originating from the electric field brought about either by intra-atomic potentials16,17,18,19 [and giving rise to splitting between heavy/light hole (Γ8) and split-off (Γ7) valence-band subbands in zinc-blende semiconductors] or by inversion asymmetry leading to the presence of Dresselhaus terms in the band dispersion E(k)15,18,20. Tunnelling anisotropic magnetoresistance21 is one more example of devices demonstrated first for (Ga,Mn)As.

Photoemission spectroscopy and its angle-resolved version that we employ in this work have already provided a number of important information, particularly on hybridization of Mn d-levels and valence-band p-states and the associated participation of d-states in the wave function of carriers at the Fermi level EF in (Ga,Mn)As22,23,24,25,26. These studies have also brought into light issues associated with surface contamination.

In this work, we at first provide information on samples preparation, protocols implemented to avoid contamination as well as on experimental procedure employed in our studies of (Ga,Mn)As by angle-resolved photoemission spectroscopy (ARPES). The obtained spectra are presented, together with results of tight-binding computations. This comparison allows us to assign particular features to individual valence-band subbands and to demonstrate that EF is located within the heavy/light hole band in ferromagnetic (Ga,Mn)As. At the same time, density of states (DOS) is considerably depleted at EF, which we take as new evidence for the presence of the Coulomb gap, driven by disorder-enhanced carrier correlations, as predicted by Altshuler and Aronov27. By studying a series of samples we have found that the lower Curie temperature TC, the larger DOS depletion, i.e., TC decreases with approaching the MIT. It is emphasized that such a Coulomb gap cannot be reproduced by the current ab initio method, since it results from quantum interference of carrier-carrier interaction amplitudes at the mesoscopic scale; therefore the experimental output on the Coulomb gap is of particular importance. Another new effect revealed by our studies is magnetic linear dichroism (MLD), particularly well resolved for the spilt-off subband. By direct computations of light absorption for transitions from this band to free-electron states we demonstrate that MLD originates from the Dresselhaus spin-orbit terms and is present for both s and p light polarization. This is intriguing since the MLD observed in ferromagnets by core-level photoemission28 is active only for p-polarization29.

Results

Samples and experimental

(Ga,Mn)As films with 25-nm thickness were grown by the molecular beam epitaxy (MBE) method1 and then transferred to the ARPES chamber without being exposed to the air by using a portable ultrahigh vacuum (UHV) chamber, which keeps 10−10 Torr during the transportation (see Methods). This procedure has proved very effective to obtain reliable ARPES data because ARPES is very sensitive to the condition of sample surface22.

The properties of (Ga,Mn)As are strongly affected by the growth conditions during MBE30. To avoid the formation of MnAs precipitates, it is necessary to grow (Ga,Mn)As below 250 °C (ref. 1), while the low-temperature growth likely produces anti-site As defects acting as double-donors31. It is also known that Mn atoms in GaAs occupy either the Ga sites or the interstitial positions32. Ga-substitutional Mn provides a hole by acting as an acceptor, while the interstitial Mn compensates holes by acting as a double-donor and its magnetic moment couples antiferromagnetically with that of the substitutional Mn33. To obtain higher metallicity of films, one needs to grow (Ga,Mn)As with less amount of interstitial Mn atoms, whose formation energy decreases with increasing Mn content34. Although the post-annealing of films effectively removes the interstitial Mn atoms, the annealing alternatively causes precipitation of Mn atoms and/or Mn-related compounds at the surface35. Because ARPES is very surface-sensitive technique, it is crucial to exclude these extrinsic factors that distort ARPES results.

We have grown 25-nm-thick Ga1−xMnxAs layers with nominal Mn compositions x between 0.05 and 0.07 on an n+-GaAs (001) substrate through a 100-nm thick n+-GaAs buffer layer at the growth temperature between 230 and 245 °C. Magnetic characterization has been performed after the ARPES measurements and a typical temperature dependence of in-plane magnetization of (Ga,Mn)As with a relatively high TC above 100 K is presented in Fig. 1a. We find that the higher growth temperature tends to result in the higher values of TC. We prepared a series of (Ga,Mn)As films with TC between 35 and 112 K in the as-grown state. A reference n-type GaAs film doped with Si concentration of 1019 cm−3 has been grown on an n+-GaAs (001) substrate in the same MBE chamber.

Figure 1
figure 1

Valence-band structure of heavily Mn-doped GaAs.

(a) Field-cooled magnetization M along the in-plane < 100 > orientation of (Ga,Mn)As with x = 0.05 grown at 245 °C, measured on heating the sample in the virtual absence of an external magnetic field (magnetic field μ0H of 0.8 mT was applied in order to compensate the remanent field in the magnetometer). These data indicate that the easy axis is along the < 100 > crystal direction56. (b) Bulk and surface Brillouin zones of (Ga,Mn)As, together with the emission plane in ARPES measurements along the cut (green shade). (c) Experimental band structure along the cut for (Ga,Mn)As (TC = 101 K) at T = 30 K measured with the He-Iα resonance line, obtained by plotting the ARPES intensity in the linear scale as a function of binding energy and wave vector. Black dashed curves are a guide for eyes to trace the band dispersions. Inset shows the same plot with an enhanced colour contrast in the area enclosed by red rectangle. (d) Experimental band structure obtained by plotting the second derivative intensity of MDCs. Bands are labeled with A-D. (e) Comparison of experimental band dispersions of (Ga,Mn)As with the tight-binding calculation for pristine GaAs at kz = 0 (ΓKX plane; solid curves) and kz = 2π/a0 (XKΓ plane; dashed curves). HH, LH and SO denote the heavy-hole, light-hole and split-off bands, respectively. The experimental band dispersions are extracted by tracing the peak position in the second-derivative intensity of MDCs (open circles) and ARPES spectra (filled circles).

Location of the Fermi level

Figure 1c shows the experimental band structure of Ga0.95Mn0.05As (TC = 101 K) along the cut obtained by plotting the ARPES-spectral intensity as a function of binding energy and wave vector. We have measured ARPES spectra with the He-Iα resonance line (hν = 21.218 eV) with the main polarization vector along the cut (Fig. 1b) at temperature T = 30 K. We observe three dispersive hole-like bands (A-C) centred at the point, all of which are ascribed mostly to the As 4p orbitals in GaAs36 and also in (Ga,Mn)As24. Band A has the top of dispersion at the binding energy of 2.9 eV. Bands B and C tend to converge on approaching the point and their spectral intensity is markedly suppressed in the region within 1 eV from EF. In the inset with enhanced colour contrast, one can see a signature of band B nearly touching EF. This feature is more clearly seen in the experimental band structure obtained by the second derivative of the momentum distribution curves (MDCs) presented in Fig. 1d. It is obvious that band B almost reaches EF at the point, while band C is hardly traced around EF because of the very faint feature. We also observe another band (band D) directing toward EF on approaching the second point.

We compare in Fig. 1e the experimental band structure with our tight-binding band calculations for GaAs (see Methods) along the ΓKX (kz = 0, solid curves) and XKΓ (kz = 2π/a0 dashed curves) high-symmetry lines. There are good correspondences between the experiment and the calculation. Comparison shows that bands B and C are assigned to the bands at kz = 0, while bands A and D are at kz = 2π/a0. This suggests that the ARPES data reflect the electronic states averaged over a wide kz region in the bulk Brillouin zone, so that the electronic states along the high-symmetry lines, i.e. at kz = 0 and 2π/a0, have a dominant contribution to the total “one-dimensional” DOS averaged over kz36,37. Therefore, the overall agreement between the experimental and the calculated dispersions as seen in Fig. 1e leads us to conclude that bands B and C are the split-off (SO) and light-hole (LH) bands, respectively.

We have investigated the detailed band structure near EF at the Γ point. As displayed in Fig. 2a, band B shows a “Λ”-shaped spectral-intensity distribution. To see the dispersive feature more clearly, we have subtracted the featureless background from the raw spectrum (Fig. 2b) and show the result and its intensity plot in Fig. 2c,d, respectively. Intriguingly, a clear Fermi-edge structure is recognized in the spectrum when the band approaches EF (Fig. 2c). It is noted that a similar Fermi-edge-like structure is already seen even in the raw spectrum (Fig. 2b), indicating its intrinsic nature. It is stressed here that the present first observation of the Fermi edge in (Ga,Mn) As is a consequence of the careful surface treatment using the UHV sample transfer chamber as well as the high-energy resolution in the ARPES measurement. As seen in Fig. 2c, the spectral intensity at the Fermi edge systematically increases when band B approaches EF, as highlighted by the MDC at EF (white curve in Fig. 2d). We also observed a similar spectral characteristics in another (Ga,Mn)As sample with lower TC of 62 K (Fig. 2e,f). As shown in Fig. 2g, band B in both samples with TC = 62 K or 101 K shows a good correspondence to the SO band lying ~0.35 eV below the LH/HH (heavy hole) band at the Γ point in GaAs. Thus, the present ARPES results unambiguously demonstrate that the Fermi level resides deeply inside the As-4p valence band. This conclusion is in accord with recent on-Mn resonance photoemission studies23,24 which confirm that the main spectral weight of Mn d states is at 3.5 eV below EF. Accordingly, taking also into account a large difference in the concentration of As and Mn, a local d-weight maximum below EF accounts for only a few percent of the total one-particle DOS in this region.

Figure 2
figure 2

Evidence for the metallic valence band in (Ga,Mn)As.

(a) Plot of near-EF ARPES intensity around the Γ point for (Ga,Mn)As with TC = 101 K as a function of binding energy and wave vector k. (b) Raw ARPES spectrum (red curve) measured at a cut shown by solid green line in (a) and the corresponding background spectrum (black curve) obtained by integrating ARPES-spectral intensity over a wide k region (0.4 Å−1) around the Γ point. Background-subtracted ARPES spectrum is shown with blue curve. (c,d) A set of background-subtracted ARPES spectra and corresponding intensity plot, respectively, for (Ga,Mn)As with TC = 101 K. White curve is MDC at EF. (e,f) Same as (c,d) but for (Ga,Mn)As with TC = 62 K. (g) Comparison of the band dispersions among (Ga,Mn)As with TC = 101 K and 62 K and n-GaAs, extracted from the peak position of MDCs, together with the band-structure calculation (same as Fig. 1e). Representative fitting result to the MDCs (red curve) with two Lorentzians (black curve) is shown in the bottom. (h) ARPES intensity around the Γ point of n-GaAs. (i) Intensity contour plots for (Ga,Mn)As with TC = 62 K as a function of in-plane wave vector at several energy slices from EF to 0.6 eV.

The hole-doped nature of the valence-band states is also corroborated by a comparison of the experimental band structure between (Ga,Mn)As with TC = 62 K (Fig. 2f) and nonmagnetic n-type GaAs (n-GaAs) doped with Si at concentration of 1019 cm−3 (Fig. 2h). The doped electron concentration of 1019 cm−3 results in the surface depletion layer with a thickness of several tens of nm, which is much thicker than the probing depth with the He Iα line in ARPES measurements. Hence, when discussing the binding energy in n-GaAs one needs to take into account the band bending effect at the surface, which is ineffective for metallic (Ga,Mn)As with a depletion layer less than 1 nm. By shifting the chemical potential of n-GaAs upward by ~0.9 eV, the band dispersion of n-GaAs fits well with that of (Ga,Mn)As (Fig. 2g). Taking into account that the sum of bandgap and SO splitting at the Γ point is ~1.85 eV, it is inferred that there is an upward surface band bending of ~0.95 eV (=1.85-0.9 eV) in n-GaAs, in agreement with the reported Fermi-level pinning position of ~0.8 eV38. The effective mass of band B obtained by ARPES is 0.08 ± 0.02 m0 (m0: the free-electron mass), in good agreement with the calculated effective mass of the SO band (LH and HH bands should have a much larger effective mass)39.

Figure 2i displays the contour plot of ARPES-spectral intensity at several binding energies from EF to 0.6 eV. We find that the intersection of band B has a circular shape indicative of the SO-band nature, supporting the above interpretation that the Fermi level of (Ga,Mn)As is deeply embedded in the valence band. On the other hand, we do not observe the HH and LH bands in Fig. 2, which is likely due to the intensity reduction by the disorder in sample as well as the photoemission matrix-element effect23, as inferred from the strong orbital dependence of the calculated photoelectron intensity in GaAs40 (see Supplementary Information). A part of the loss in the spectral intensity may be accounted for by the presence of a Coulomb gap discussed below.

Coulomb gap

Having established the position of the Fermi energy in the As 4p valence band in (Ga,Mn)As, the next question is how the metallic state evolves from pristine GaAs by Mn doping. To address this question, we performed systematic ARPES measurements on a variety of samples with different TC’s. Figure 3a shows a set of ARPES spectra at the Γ point revealing a rounded shape in the vicinity of EF unlike typical metals such as gold (Au) which shows a sharp Fermi-edge cut-off as seen in Fig. 3a. This rounded feature of spectrum provides a new evidence for a depression of one-particle DOS at EF in (Ga,Mn)As, so far seen by tunnelling spectroscopy and assigned to the electron-electron correlation in the vicinity to a disorder-driven MIT41,42,43.

Figure 3
figure 3

Evolution of metallic states and soft Coulomb gap.

(a) Comparison of raw ARPES spectra around the Γ point (integrated over ±0.08 Å−1 centered at the Γ point) for (Ga,Mn)As with various TC’s (35, 60, 62, 101 and 112 K) and nonmagnetic n-GaAs. The intensity is normalized to the spectral weight integrated over 0–1 eV. Numerical simulation of the ARPES spectrum (blue solid curve) for (Ga,Mn)As with TC = 101 K by taking the spectral DOS with a soft Coulomb gap (dotted curve). The simulation was performed with a linearly decreasing DOS with square root energy dependence of the gap27, multiplied by the Fermi-Dirac distribution function at T = 30 K and convoluted with the resolution function (energy resolution 15 meV). Shaded area highlights the energy region of the gap. ARPES spectrum of gold (Au) is also shown for comparison. (b) Spectral DOS at EF plotted as a function of TC. The DOS is obtained by integrating the ARPES-spectral intensity within ± 40 meV of EF. Insets show schematics of the DOS in (Ga,Mn)As with high and low TC.

According to the Altshuler-Aronov theory27,44, the interaction-induced correction to single-particle DOS in a three-dimensional metal assumes the form δν(ε) = a + b|ε|1/2 at kBT < |ε| < ħ/τ < EF, where ε is the quasiparticle energy with respect to EF and τ is the momentum relaxation time. The sign and magnitude of a and b are determined by mutually compensating interference effects in the singlet and triplet particle-hole diffusion channels45,46,47. However, spin-disorder scattering in a paramagnetic phase or large band spin-splitting in a ferromagnetic phase destroy the interference in the triplet channel, which according to the perturbation46,47 and dynamic renormalization group theory45,48, should result in a sizable decrease of DOS at EF, i.e., a < 0, b > 0 and ν(0) → 0 on approaching the MIT from the metal side.

As shown in Fig. 3a, a numerically simulated ARPES spectrum with taking into account the expected form of the DOS depression, δν(ε) = a + b|ε|1/2, where a < 0 and b > 0 reproduces satisfactorily the experimental spectrum up to ~0.1 eV relative to EF (Fig. 3b). This value coincides well with the energy at which the tunnelling DOS starts to be depleted43 and also the optical conductivity is peaked49, suggesting that the suppression of DOS near EF is an inherent feature of ferromagnetic (Ga,Mn)As. We observe a monotonic decrease of the near-EF spectral weight upon decreasing TC. This is more clearly visible by plotting the DOS at EF, i.e., ν(0) as a function of TC (Fig. 3b), demonstrating a close relationship between the proximity to MIT, TC and ν(0). Accordingly, an effect from the soft Coulomb gap should be considered when discussing the one-particle excitations. At the same time, many important physical parameters are determined by the thermodynamic DOS which in terms of the hole concentration p assumes the form ρ(EF) = −∂p/∂EF, so that its magnitude is only weakly perturbed by disorder and carrier correlations on the metallic side of the MIT. According to the theoretical predictions46,47, the temperature dependence of electric conductance at low temperatures belongs to this category. In agreement with these expectations, experimental and theoretical studies of millikelvin conductance in various dimensionality structures of (Ga,Mn)As provided the magnitude of ρ(EF) close to that expected for holes in an unperturbed valence band of GaAs (ref. 50). The present ARPES results thus unambiguously demonstrate that (Ga,Mn)As with high TC possesses the metallic electronic structure with the Fermi level in the valence band and is strongly influenced by the conjunction of correlation and disorder effects.

Magnetic linear dichroism

We now turn our attention to the “magnetic linear dichroism” (MLD) effect in the ARPES intensity of the valance band in (Ga,Mn)As (for details of the MLD measurements, see Methods). According to Fig. 4a,b the symmetry of the ARPES intensity with respect to the Γ point shows a magnetization-direction dependence in (Ga,Mn)As with TC = 101 K. As shown in Fig. 4a, the near-EF photoelectron intensity of the SO band in the ferromagnetic phase (T = 30 K) in the k region (left-hand side, negative k values) is stronger than that in the k[110] region (right-hand side) when we magnetize the sample along the [100] direction as illustrated in Fig. 4c (see also Fig. 2d,f). On the other hand, when we reverse the magnetization direction to the photoelectron intensity in the positive k region becomes more dominant (Fig. 4b). The intensity difference between two opposite magnetization cases is estimated to be 2–8% (depending on the binding energy) by taking account of background contribution of ~70%. When we overlaid the band dispersion for two magnetization cases, we recognize a finite difference in their energy positions as shown in Fig. 4d. Interestingly, when we cool down the sample without external magnetic field, the intensity asymmetry becomes less clear (Fig. 4e). In addition, the asymmetry is not seen in nonmagnetic n-GaAs (Fig. 4f).

Figure 4
figure 4

Magnetization dependence of ARPES intensity in (Ga,Mn)As.

(a,b) Near-EF ARPES intensity at T = 30 K of (Ga,Mn)As with TC = 101 K, magnetized along the [100] and directions, respectively. Calculated band dispersion of (Ga,Mn)As with incorporating the exchange splitting57 is shown with yellow curves for comparison. (c) Experimental geometry of sample axes, magnetization direction and incident photons. Solid and dashed rectangles indicate the plane of incidence (110) and the emission plane of photoelectrons , respectively. (d) Location of the energy band extracted from the peak positions of the momentum distribution curves (MDCs) in (a,b). (e) Same as (a,b) but with zero-field (z.f.) cooling. (f) Near-EF ARPES intensity at T = 30 K of GaAs (n-type; Si-doped).

In order to determine the origin of this MLD we note that for the valence-band photoemission in question, photon-induced transitions occur between the extended valence band states (SO band in our case) and the free-electron-like states which propagate towards the detector in the vacuum while damped inside the crystal (the one-step model). Since there is no periodicity in the direction normal to the surface, kz is not conserved. Therefore valence states with various kz can contribute to the magnitude of emitted photocurrent at given Ekin and k// values. Hence, following the method elaborated previously51, we have calculated for a given orientation of k// and energy ESO in the SO band the absorption coefficient for linearly polarized light considering electric-dipole transitions from the SO band split by the p-d exchange and spin-orbit interactions to spin-degenerate free-electron states. Since only k// is conserved, the absorption coefficient for four possible transitions is a product of corresponding DOS ν(ESO) = k//2(2π2dESO/dk//)−1 at kz = 0 and the square of the absolute value of the matrix element averaged over kz values. No MLD is found if the valence band is described within the six band Luttinger model51. In contrast, allowing for the presence in the 6 × 6 kp Hamiltonian of terms brought about by the inversion asymmetry52,53, MLD of characteristics observed experimentally shows up. There are three such contributions in the case under consideration: terms proportional to k3, originally introduced by Dresselhaus for zinc-blende crystals and two components linear in k, appearing under biaxial and shear strain, respectively. In the absence of the exchange splitting, time reversal symmetry leads to the same magnitude of photoemission for k and -k if transitions from both spin subbands are taken into account. However, this is no longer the case if time reversal symmetry is broken by non-zero magnetization.

Figure 5 shows the computed magnitude of MLD = (Ak[110] − Ak) /Ak[110] at M || [100] for the SO band as well as for kp parameters specified in Methods11,15,18,52 and for the expected values of strain2 εxx = −0.4% and εxy = 0.1% as well as for the magnitude exchange splitting of the valence band at Γ8 point Δv = 180 meV, corresponding to the saturation magnetization for x = 0.05. Since the character of light polarization is uncertain under our experimental conditions, the calculations have been performed for both s and p polarization. We have found that for the employed parameters, the k3 contribution dominates. As seen, the computed magnitude of the average value of MLD over the two polarization is in accord with experimental results, i.e., the theory predicts properly the sign, magnitude and energy width of the effect as well as its symmetry, Ak[110](M || 100) = Ak (M || ). It is worth noting that to resolve directly spin splitting and k asymmetry of the valence-band subbands ultrahigh-resolution spin-resolved ARPES measurements would be necessary.

Figure 5
figure 5

Theoretical magnitudes of magnetic linear dichroism generated by inversion asymmetry in ferromagnetic (Ga,Mn)As epilayers.

Here, MLD = (Ak[110] −Ak ) /Ak[110], where A is a product of the square of the absolute value of the matrix element for electric dipole transitions average over |kz| ≤ 0.2π/a0 and density of states at k || [110] and , respectively corresponding to a given energy in the split-off band at kz = 0. The magnitude of the exchange splitting of the valence band Δv = 180 meV (at Γ8 point) corresponds to the value of saturation magnetization M for Mn content x = 0.05. Empty symbols are for M || [100] and for s and p light polarization (see Fig. 4c). The full symbols depict average values weighted by relative magnitudes of A for s and p polarization. The sign of MLD is reversed for M || .

Conclusions

To conclude, our ARPES studies of (Ga,Mn)As resolve particularly well the band split off by spin-orbit coupling (i.e., Γ7 valence band subband at k = 0). Since the magnitude of spin-orbit splitting is known quite precisely and varies little across the arsenide family of III-V compounds, we have been able to evaluate quite accurately the position of the Fermi level with respect to the top of the valence band. The obtained Fermi energy of about −0.3 eV in (Ga,Mn)As with TC of the order of 100 K is consistent with the p-d Zener model16,51. At the same time our data point to the presence of a substantial depression in DOS at EF. We assign this observation, in accord with tunnelling studies41,42,43, to electron correlation in disordered systems, which result in a depletion of DOS for the universality class in question45,46,47,48. We have also detected magnetic linear dichroism that results from inversion asymmetry of zinc-blende structure. Our computations within kp theory shows that a dominant contribution to the magnitude of MLD in the SO band comes from cubic k3 terms. It would be interesting to check their contribution to spin-orbit torque analysed so-far considering only linear terms in k (refs 13, 14, 15).

Methods

Sample preparation and ARPES experiments

25-nm thick Ga1−xMnxAs films with x between 0.05 and 0.07 were grown at substrate temperature between 230 and 245 °C on n+-GaAs (001) substrate through a 100-nm thick n+-GaAs buffer layer grown at ~560 °C by molecular beam epitaxy (MBE). After the growth, we immediately transferred the films to a portable ultrahigh-vacuum (UHV) chamber equipped with non-evaporating getter and ion pumps, disconnected the portable chamber from the MBE system and then connected to the ARPES system with keeping the UHV condition better than 5 × 10−10 Torr. All the sample-transfer procedures were carried out within one hour. ARPES measurements were performed with the MBS-A1 electron analyzer equipped with a high-intensity He plasma discharge lamp. We use the He-Iα resonance line (photon energy: hν = 21.218 eV) to excite photoelectrons. The energy resolution for the ARPES measurements was set at 15–40 meV. The sample temperature was kept at 30 K during the measurements.

For the MLD measurements, first we cooled down the sample to T = 30 K without magnetic field in the UHV chamber directly connected to the ARPES chamber with μ-metal shielding (~0.5 μT). Next, we magnetized the (Ga,Mn)As thin film (TC = 101 K) by approaching the bar-shaped neodymium magnet to align its magnetic field parallel to the magnetization easy axis of the sample, [100] or . It is noted that this easy axis is indicated by the observation of sizable magnetization and smooth temperature dependence of the magnetization curve at the remanent magnetic state in Fig. 1a. After the magnetization, we transferred the sample to the ARPES chamber with keeping the sample temperature at 30 K and measured the MLD data without external magnetic field (Fig. 4a,b). For the zero-field cooling measurement (Fig. 4e), we demagnetized the sample in the ARPES chamber by increasing the temperature up to 150 K (above TC) and then slowly cooled down to T = 30 K.

Calculations

A tight-binding method described previously54 was employed to determine the band structure E(k) of GaAs and Ga0.95Mn0.05As in the whole Brillouin zone. Density of states and matrix elements for electric-dipole optical transitions between split-off band (Γ7 at k = 0) and free electron states were computed from the six band Luttinger-Kohn kp theory with the standard values51 of the Luttinger parameters, the p-d exchange integral β and the elastic moduli cij, taking additionally into account odd in k terms resulting from inversion asymmetry of the zinc-blende structure. The numerical values of kp parameters ( = −58.71 eVÅ2,  = −81.93 eVÅ2 and  = −101.9 eVÅ2), describing the dominant k3 terms, are taken from ref. 53. There are additional two contributions to the six band kp Hamiltonian linear in k, contributing also to the sector of the kp Hamiltonian describing the split-off band Hv7v7. These terms are generated by shear εij and diagonal εii strain tensor components that enters into Hv7v7 via σ⋅φ and σ⋅ψ, respectively52, where σi are Pauli matrices, φx = kyεxy − kzεxz and ψx = kx(εyy − εzz) (and cyclic permutations). The form of the six band Hamiltonian involving φi and the corresponding deformation potential constant C4/ħ = −2.2 × 106 m/s were given previously18. The magnitude of the deformation potential describing the term with ψi is unknown. Following a previous approach15, the same value of the deformation potential and also the same form of the Hamiltonian (with φi replaced by ψi) are adopted. In the case under consideration εxy ≠ 0, εyz = εzx = 0 and εxx = εyy = −εzzc11/2c12, where z axis is taken along the growth direction51,55.

Additional Information

How to cite this article: Souma, S. et al. Fermi level position, Coulomb gap and Dresselhaus splitting in (Ga,Mn)As. Sci. Rep. 6, 27266; doi: 10.1038/srep27266 (2016).