Introduction

Angle-resolved photoemission spectroscopy (ARPES) is the most direct experimental technique to measure the momentum-resolved electronic eigenvalues (energy bands) of crystalline solids1,2. The technique is based on the photoelectric effect3,4, in which electrons inside solids absorb a photon with energy larger than the work function and escape into the vacuum. A key aspect of ARPES is that the energy (E) and momentum (k)-resolved photoemission signal is proportional to the single-particle spectral function A(k, E)5, a fundamental quantity containing key information about many-body interactions inside solids.

However, more than two decades ago, it was realized that the direct relationship between ARPES spectral lineshapes and the spectral function was obstructed by additional intensity modulation coming from the transition dipole matrix element connecting the initial Bloch wave and the final photoelectron states6. Because of its semimetallic 2D nature and its speculated textbook Fermi-Liquid behavior, the layered transition metal dichalcogenide (TMDC) 1T-TiTe2 was the central material around the debate on directly linking the spectral function and the ARPES lineshape7,8,9,10. Indeed, 1T-TiTe2 served as a benchmark material to demonstrate the signatures of many-body effects in ARPES lineshape11, the possibility to retrieve the photohole lifetime12, as well as the potential to dissect the role of the different quasiparticle scattering processes from ARPES measurements10.

On the other hand, very little attention has been devoted to the other key ingredient governing ARPES spectra of 1T-TiTe2: the transition dipole matrix element. The effect of the transition dipole matrix element is often viewed as detrimental, since it breaks the quasi-one-to-one correspondence between the ARPES spectrum and the occupied electronic band structure, by imprinting a complex and often strong energy-dependent and momentum-dependent modulation onto the photoemission intensity13,14. Matrix element effects can sometimes be confused with signatures of many-body interactions6,15 in ARPES spectra, and are often interpreted by hand-waving arguments. While computational schemes to remove matrix element effects from measured ARPES spectra are an attractive route to focus on the band structure16, it has become clear that additional information on the electronic wavefunction can be obtained by understanding matrix element effects. While the matrix element has some extrinsic contributions coming from the experimental geometry (e.g. photon momentum, light polarization, and crystal orientation), it also encodes intrinsic information on the nature of the initial state from which the electrons are emitted. Therefore, disentangling extrinsic and intrinsic contributions is of capital importance to access the additional information encoded in the photoemission matrix elements. Important examples include gaining insights into the orbital character of the electronic band structure in solids17,18,19,20,21, the chirality of charge carriers22,23, quantum geometric properties, such as orbital pseudospin texture24,25 as well as Berry curvature26,27,28,29. The orbital texture in TMDCs is very sensitive to the lattice structure and collective effects such as charge density waves (CDW)30,31,32,33,34 and the debated excitonic insulator state35,36,37. As demonstrated recently, CDW phases and orbital texture can be intertwined38. Tracing the orbital texture from matrix element effects in photoemission will also greatly benefit discerning structural and collective phases in TMDCs. More specifically, 1T-TiTe2 has been recognized to be a very interesting prototypical 2D material: it exhibits a strong thickness33,39 and substrate-dependent34 emergence of CDW and pseudogap phases33.

To disentangle between extrinsic and intrinsic contributions to dichroism in ARPES, we recently introduced a measurement methodology to extract the modulation of the photoemission intensity upon a crystal rotation by a well-chosen angle θ in the laboratory frame (\({{{{{{\mathcal{R}}}}}}}_{\theta }\))24. We benchmarked this approach by investigating the modulation of photoemission intensity from 2H-WSe2 upon a 60° azimuthal crystal rotation. In 2H-TMDCs, due to the broken inversion symmetry within each monolayer and the strong spin-orbit splitting, a 60° crystal rotation (\({{{{{{\mathcal{R}}}}}}}_{6{0}^{\circ }}\)) is mimicking a time-reversal operation (\(\hat{T}\)), i.e. \({{{{{{\mathcal{R}}}}}}}_{6{0}^{\circ }}\equiv \hat{T}\), which is equivalent to swapping valley-pseudospin indices, i.e., K → K′, and vice versa. We have shown that since the extrinsic contribution to dichroism was invariant upon \({{{{{{\mathcal{R}}}}}}}_{6{0}^{\circ }}\), the antisymmetric components of the dichroism upon \({{{{{{\mathcal{R}}}}}}}_{6{0}^{\circ }}\) (i.e. the differential signal which undergoes a sign change, called time-reversal dichroism in the photoelectron angular distribution—TRDAD) is a background-free intrinsic dichroic signal sensitive to the orbital pseudospin texture in 2H-WSe224.

In this work, we use the same experimental methodology—i.e. we measure the photoemission intensity modulation upon sample rotation—to unveil the orbital texture of 1T-TiTe2. We introduce the intrinsic linear dichroism in photoelectron angular distributions (iLDAD)—a generalized version of TRDAD, in samples featuring inversion symmetry. Indeed, TRDAD can be seen as a special case of iLDAD, when a sample possesses broken inversion symmetry and when the crystal rotation is thus equivalent to a time-reversal operation. The iLDAD from 1T-TiTe2, extracted from ARPES measurements for two different crystal orientations, is fully reproduced by first-principle simulations based on the Korringa–Kohn–Rostoker (KKR) method. The excellent agreement of experiment and theory allows us to disentangle different contributions to the ARPES signal and to construct a tight-binding (TB) model capturing the orbital texture. Combining the TB model with qualitative modeling of the photoemission matrix elements, we trace the measured dichroism to the tilted real-space character of the relevant orbitals. The intrinsic dichroism is found to be particularly sensitive to the orbital character, which is beyond standard LDAD based on the anisotropy of the photoemission intensity upon rotating the light polarization direction by 90°17,19,20. Hence, iLDAD is a powerful addition to the toolbox of dichroism in ARPES, in the spirit of going beyond band structure mapping and gaining information about energy- and momentum-resolved electronic wavefunction of solids.

Results

Band structure mapping

The experimental apparatus is composed of a table-top femtosecond (fs) XUV source (centered around 21.7 eV, spectral bandwidth 110 meV FWHM)40 coupled to a time-of-flight momentum microscope spectrometer (METIS 1000, SPECS GmbH) (see Fig. 1a). This detector allows measuring all electrons emitted above the surface of the sample and recording a stream of individually detected photoelectrons, resolved in the three intrinsic dimensions of the detector: kinetic energy (E) and both parallel momentum components (kx, ky). One of the major advantages of this type of detector is the parallel detection of the 3D photoemission intensity I(E, kx, ky) in a single measurement, without having to rotate the sample in the laboratory frame41. More details about the experimental setup can be found elsewhere40,42 and in the “Methods” section. Giving the energy-dependence of the inelastic mean free path of the outgoing electrons43,44, working in the XUV spectral range ensure that most of the detected electrons are coming from the topmost 1T-TiTe2 layer (first atomic trilayer)24. More details about this statement, as well as its theoretical validation, will be given later.

Fig. 1: Scheme of the multidimensional photoemission spectroscopy experiments.
figure 1

a A p-polarized fs-XUV (21.7 eV) pulse is focused onto a freshly cleaved bulk 1T-TiTe2 crystal at an angle of incidence of 65° with respect to the surface normal, emitting photoelectrons which are detected by a time-of-flight momentum microscope (METIS1000, SPECS GmbH). The crystal is orientated in the laboratory frame such that the light polarization axis is along the Γ−M/M' direction. b Schematic of the equivalence between a crystal rotation by 60° (\({{{{{{\mathcal{R}}}}}}}_{6{0}^{\circ }}\)) in the laboratory frame and a mirror operation (z → −z), in 1T-TiTe2. Indeed, the 1T- polytype is strictly equivalent to a mirror operation across the z plane of each layer (z → −z). c Three-dimensional photoemission intensity I(E, kx, ky), and associated constant energy contour (CEC) at the Fermi energy (E−EF = 0 eV). d Same as (c), but following a crystal rotation by 60° in the laboratory frame.

The experimentally measured 3D photoemission intensities and associated constant energy contours (CECs) at the Fermi energy (E−EF = 0 eV) are shown in Fig. 1c and d. Looking at the 3D photoemission intensity (Fig. 1c), we recognize the well-known semimetallic nature of 1T-TiTe2, featuring spectrally overlapping Te-5p hole pockets at Γ and Ti-3d electron pockets around two inequivalent M points (denoted by M/M′)45. Looking at the CECs in Fig. 1c, the most striking observation is the strong photoemission intensity anisotropy between Ti-3d electron pockets located at M and M′. Indeed, while the photoemission intensity from M′ pockets is strong, it is almost completely suppressed for adjacent M pockets. The emergence of this strong linear dichroism in photoemission a priori can originate from extrinsic (e.g. experimental geometry) and/or intrinsic (i.e. orbital/Bloch wavefunction symmetries) effects.

Motivated by the recently introduced approach for disentangling extrinsic and intrinsic contributions to dichroism in ARPES (see ref. 24 and the section “Introduction”), we also performed band structure mapping of 1T-TiTe2 for a second crystal orientation rotated by 60°. In inversion symmetric 1T-TiTe2, this crystal rotation (\({{{{{{\mathcal{R}}}}}}}_{6{0}^{\circ }}\)) is strictly equivalent to a mirror-image symmetry of each layer (see Fig. 1b), i.e. \({{{{{{\mathcal{R}}}}}}}_{6{0}^{\circ }}\equiv z\to -z\). Although the effect of a 60° crystal rotation in 1T-TMDC (\({{{{{{\mathcal{R}}}}}}}_{6{0}^{\circ }}\equiv z\to -z\)) is different than in 2H-TMDC (\({{{{{{\mathcal{R}}}}}}}_{6{0}^{\circ }}\equiv \hat{T}\)), in both cases it corresponds to a well-defined symmetry operation on the electronic wavefunction of the crystal. Linking a given symmetry operation on the electronic wavefunction to the modulation of the photoemission intensity, complemented by proper theoretical analysis, is a powerful tool to investigate the energy-dependent and momentum-dependent orbital texture, a quantity of fundamental importance but hardly accessible with other techniques.

Intrinsic linear dichroism in photoelectron angular distributions

In Fig. 2, we show the procedure to extract the antisymmetric components of the dichroism (referred to as iLDAD) from photoemission data obtained from two crystal orientations rotated by 60° with respect to each other. In Fig. 2a, one can see the CEC at the Fermi energy featuring the triangular doughnut-shaped Te-5p hole pockets around Γ, which shows subtle but non-vanishing angular anisotropies. More remarkably, the strongly anisotropic photoemission from Ti-3d electron pockets located around M and M′ is dominating the signal. As already discussed in the context of Fig. 1, while the photoemission intensity is strong for the M′ pockets, it is almost vanishing for the adjacent M pockets. Upon a 60° azimuthal rotation of the crystal (see Fig. 2b), this asymmetric photoemission yield from M and M′ pockets remained mostly unchanged, leading to a 60° rotation of the momentum-resolved photoemission patterns. This is a strong indication that this asymmetric photoemission yield is of intrinsic origin and should encode information about the symmetry properties of the orbitals.

Fig. 2: Extraction of the intrisic linear dichroism in photoelectron angular distributions (iLDAD).
figure 2

a and b \({I}_{{0}^{\circ }}\) and \({I}_{6{0}^{\circ }}\), the CECs at the Fermi energy, measured for two crystal orientations rotated by 60° with respect to each others. c ΔI/I, the raw normalized difference, i.e. (\({I}_{{0}^{\circ }}\) − \({I}_{6{0}^{\circ }}\))/(\({I}_{{0}^{\circ }}\) + \({I}_{6{0}^{\circ }}\)) between CECs shown in (a) and (b). d and e \({A}_{{{{\mathrm{LDAD}}}}}^{{0}^{\circ }}\) and \({A}_{{{{\mathrm{LDAD}}}}}^{6{0}^{\circ }}\), the “left–right asymmetries'', reflect the photoemission intensity asymmetry between the kx < 0 and kx > 0 hemispheres, for two crystal orientations, respectively. f iLDAD represents the component of \({A}_{{{{\mathrm{LDAD}}}}}^{{0}^{\circ }/6{0}^{\circ }}\) which is antisymmetric upon 60∘ azimuthal rotation of the crystal.

While Fig. 2b shows the raw normalized photoemission intensity difference between the two crystal orientations, using our multidimensional detection scheme with the p-polarized fs-XUV pulses impinging in the x-z plane (along with Γ−M/M′), the normalized intensity differences between the forward Iα(E, kx, ky) and backward Iα(E, − kx, ky) hemisphere, i.e. the linear dichroism asymmetry in the photoelectron angular distribution \({A}_{{\mathrm {LDAD}}}^{\alpha }(E,{k}_{x},{k}_{y})\), can be extracted (see Eq. (1)), without the need to rearrange the sample geometry nor the light-polarization state46,47.

$${A}_{{{{\mathrm{LDAD}}}}}^{\alpha }=\frac{{I}_{\alpha }(E,{k}_{x},{k}_{y})-{I}_{\alpha }(E,-{k}_{x},{k}_{y})}{{I}_{\alpha }(E,{k}_{x},{k}_{y})+{I}_{\alpha }(E,-{k}_{x},{k}_{y})}$$
(1)

Looking at \({A}_{{{{\mathrm{LDAD}}}}}^{{0}^{\circ }}\) (Fig. 2d), we notice very strong dichroism (up to ±97%) around M/M′ pockets with an opposite sign for adjacent valleys, and weaker dichroism (up to ±35%) with the same threefold symmetry and sign texture for the Te-5p hole pockets around Γ. The dichroic signal is not perfectly up/down (ky > 0/ky < 0) symmetric, suggesting small experimental imperfection in the crystal orientation, XUV beam alignment, or photoelectron imaging, since the up/down symmetry, should not be broken within this experimental configuration (i.e. p-polarization along Γ−M/M′). Moreover, since we are dealing with normalized signals (see Eq. (1)), we only show the dichroism for momentum regions where the right-left (kx > 0/kx < 0) symmetrized signal is > 0.1Imax, where Imax is the maximum photoemission intensity for this CEC. This thresholding procedure ensures that we do not over-interpret the potentially strong dichroism coming from momentum regions of vanishing photoemission intensity. Figure 2e shows \({A}_{{{{\mathrm{LDAD}}}}}^{6{0}^{\circ }}\), i.e. the dichroism extracted after rotating the crystal orientation by 60∘. The measured \({A}_{{{{\mathrm{LDAD}}}}}^{6{0}^{\circ }}\) dichroic signal is almost perfectly antisymmetric with respect to \({A}_{{{{\mathrm{LDAD}}}}}^{{0}^{\circ }}\), again indicating its intrinsic origin since an extrinsic (experimental geometry) origin would have led to a dichroic signal invariant upon an operation.

In order to quantitatively isolate the antisymmetric part of the \({A}_{{{{\mathrm{LDAD}}}}}^{\alpha }\) dichroism upon crystal rotation (\({{{{{{\mathcal{R}}}}}}}_{6{0}^{\circ }}\)), i.e. to remove any spurious contributions from experimental geometry, we introduce the iLDAD:

$${i}{{{\mathrm{LDAD}}}}=\frac{{A}_{{{{\mathrm{LDAD}}}}}^{\alpha }-{A}_{{{{\mathrm{LDAD}}}}}^{\alpha ^{\prime} }}{2}$$
(2)

where the crystal rotation is hereby an angle α−\(\alpha ^{\prime}\) = 60∘. The resulting iLDAD is shown in Fig. 2f. The strong similarity between \({A}_{{{{\mathrm{LDAD}}}}}^{{0}^{\circ }}\) and iLDAD is another observation revealing that the dichroism in 1T-TiTe2 is almost exclusively originating from intrinsic origins (orbital symmetry/orientation). We will investigate the microscopic origin of the emergence of this strong intrinsic dichroism and how it can be used to reveal key information about the orbital texture of this layered semimetallic TMDC.

Comparison with DFT band structure and one-step model of photoemission

To extract quantitative information about orbital texture from such differential (iLDAD) ARPES measurements, a systematic comparison with quantitative theories is highly beneficial. Here, we use inputs from two different theoretical approaches in a synergistic manner. First, we performed one-step model photoemission calculations based on fully relativistic density functional theory (DFT). The one-step model of photoemission is implemented in the fully relativistic KKR method48,49. The calculated photoemission signal includes all matrix element effects such as experimental geometry, photon energy, polarization state, and final state effects (for more details, see the “Methods” section). Second, we constructed a first-principle TB model by computing projective Wannier functions (see the “Methods” section). The TB model leads to an excellent description of the electronic structure of 1T-TiTe2 and provides an intuitive picture of the link between orbital degrees of freedom and dichroism in photoemission.

An important step, before calculating iLDAD, is to ensure that both theoretical approaches can accurately capture the electronic structure of 1T-TiTe2. The experimentally measured photoemission intensity, the simulated KKR photoemission intensity, and the calculated DFT band structure along the M–Γ–M′ cut are shown in Fig. 3. The experimental measured photoemission intensity are shown in Fig. 3a, from which one can clearly recognize the well-known semimetallic nature of 1T-TiTe2, featuring Te-5p hole pockets at Γ and Ti-3d electron pockets around M/M′ points. The experimentally obtained band structure is well captured by the calculated DFT band structure. Within the local density approximation (LDA)—including spin–orbit coupling (SOC)—used for the KKR method (Fig. 3b), the measured band structure near the Fermi level is well reproduced, both at the Γ point and at the M/M′ pockets. At lower energies (E−EF ~ −0.9 eV) the measured spectrum shows bands that are absent in the calculated spectrum. We also performed DFT calculations using the Perdew–Burke–Ernzerhof (PBE) functional (Fig. 3c), which captures these bands. The size of the pockets at M/M’ is slightly smaller than in the experiments and within the LDA, which is consistent with calculations including exchange effects39. We have also calculated the band structure of a free-standing monolayer (Fig. 3d) for comparison. Albeit the band structure is generally modified with respect to the bulk (especially near Γ), the electron pockets at M/M′ are still present. Constructing a TB model from projective Wannier functions (red lines in Fig. 3d) provides an excellent fit of PBE band structure. The TB description provides us with a simple model for understanding the measured photoemission features at M/M′, which is discussed below.

Fig. 3: Band structure of 1T-TiTe2 along the M-Γ-M′ cut.
figure 3

a Experimentally measured photoemission intensity using the same experimental geometry as described in Fig. 1. b Photoemission intensity along the same cut obtained using the one-step model (KKR) method. The electronic structure has obtained from the local density approximation (LDA). c Band structure of bulk TiTe2 obtained from nonrelativistic DFT calculations. We used the Perdew–Burke–Ernzerhof (PBE) functional. d Band structure of a free-standing monolayer from DFT calculations (gray lines). The thin red lines represent the bands obtained from the ab initio TB model, which has been constructed from projective Wannier functions.

By looking more carefully at Fig. 3a, and as seen in Fermi level CECs (Fig. 1c and d), one can notice strong photoemission anisotropies between positive and negative kx (corresponding to Γ−M′ and Γ−M direction, respectively). These photoemission anisotropies, originating from transition dipole matrix element effects, are well-captured by the calculated KKR photoemission intensity, shown in Fig. 3b. As mentioned earlier, the most striking asymmetry in the photoemission signal is coming from Ti-3d electron pockets around M/M’ points, near the Fermi level. Moreover, because the role of orbital texture near the Fermi level is of capital importance in determining materials’ optical and transport properties, we will focus the subsequent analysis on the photoemission anisotropies near the Fermi energy.

Calculated iLDAD and its photon-energy dependence using KKR

In Fig. 4a and b we compare the iLDAD calculated using the one-step model (KKR) for two different photon energies (18.7 and 21.7 eV, respectively). For both photon energies, the calculated iLDAD at M/M′ pockets is very similar to the experimentally measured iLDAD. However, surprisingly, using 18.7 eV instead of the experimentally used 21.7 eV photon energy captures better the iLDAD features around Γ. This can be interpreted as a limitation of DFT: it is well known that the calculated time-reversed low-energy electron diffraction (TRLEED) final state can be shifted in energy with respect to the experimental value (for more details see Fig. 4). However, the dispersion and the orbital character of the time-reversed LEED state are typically well captured. Moreover, it was pointed out by Strocov et al. 45 that for 1T-TiTe2 the photoemission final state in the UV photon energy regime features dramatic final state effects such as multiband composition and non-free electron like dispersions around the center of the Brillouin zone. In addition, it was shown that the states around Γ at the Fermi level show three-dimensional character due to the interactions between 1T-TiTe2 monolayers (atomic trilayers)45.

Fig. 4: Photon energy dependence of iLDAD calculated using the KKR method.
figure 4

a and b iLDAD calculated using the KKR and a photon energy of 18.7 eV. The regions of interest label 1–4 are selected to investigate the photon energy dependence of iLDAD. b iLDAD as a function of photon energy, in the selected region of interest shown in (a).

To investigate the role of the final states, we performed a detailed theoretical photon-energy-dependent iLDAD scan (kz-scan). In Fig. 4c, we show the photon energy dependence of the iLDAD in the selected regions of interest (marked as 1–4). While the iLDAD of Ti-3d derived electron pockets around M/M′ points (features 1 and 2) are extremely stable against variation of the photon energy, the hole pockets around Γ point (features 3 and 4) exhibits strong modulations including even sign reversal upon modification of the photon energy. We attribute this to the strong 3D character of Te-5p states. A similar trend was observed for the soft-X-ray photoemission, where the final state becomes truly free-electron-like (see SI), ruling out the predominance of final-state effects in the emergence of iLDAD.

Performing such a photon-energy scan of the iLDAD is currently not possible using our HHG-driven XUV source. This is specifically for this reason that we took advantage of our quantitative KKR theoretical framework to investigate the photon-energy-dependence of the iLDAD. In the future, synchrotron light source could be used for photon-energy-dependent iLDAD studies.

The role of orbital texture in photoemission: tight-binding analysis

While the iLDAD exhibits a pronounced photon energy dependence near the Γ point, the robustness of the iLDAD near the M/M′ pockets gives us confidence that we can link its origin to microscopic properties of the electronic structure of 1T-TiTe2. To find an intuitive link between the orbital texture and this intrinsic dichroic signal, we employ an ab initio TB model. Within this model, the Bloch wavefunction is represented by

$${\psi }_{{{{{{\bf{k}}}}}}n}({{{{{\bf{r}}}}}})=\frac{1}{N}\mathop{\sum}\limits_{{{{{{\bf{R}}}}}}}{e}^{i{{{{{\bf{k}}}}}}\cdot {{{{{\bf{R}}}}}}}\mathop{\sum}\limits_{j}{C}_{jn}({{{{{\bf{k}}}}}}){\phi }_{j}({{{{{\bf{r}}}}}}-{{{{{\bf{R}}}}}}),$$
(3)

where ϕj(r) denote the Wannier functions localized at particular sites R. The coefficients Cjn(k) are the coefficients connection orbital (index j) and band (index n) space; they are the eigenvectors of the TB Hamiltonian \({H}_{j{j}^{\prime}}^{{{{{{\rm{TB}}}}}}}({{{{{\bf{k}}}}}})\). The Wannier representation (3) provides a direct of real and reciprocal space. For each momentum k the Bloch wavefunction is a coherent superposition of the Wannier functions ϕj(r), whose orbital symmetry is directly reflected in the ARPES signal.

From a layer-resolved KKR calculation, we found that the characteristic intensity variation can be captured by including the top layer TiTe2 only (see SI). Therefore, we constructed the TB model for a free-standing monolayer, including five Ti-d orbitals three p orbitals for each of the two Te chalcogen atoms in Eq. (3). The model provides a systematic way to examine the orbital character of the bands; in particular, we find that the contribution of the Te-p orbitals is negligible near EF around M/M′ pockets (see SI). Figure 5 shows the corresponding orbital-projected band structure for the dominant Ti-d orbitals along the M–Γ–M′ (upper row) and a M′–Γ–M cut rotated by 60° (lower row in Fig. 5), respectively. Focusing on the electron pockets crossing the Fermi level near M/M′, it is clear that the out-of-plane Ti-\({d}_{{z}^{2}}\) orbitals (Fig. 5b, e) contribute equally to the density of state along both paths. Inspecting the Ti-dxz and Ti-dyz orbital-projected band structure (Fig. 5c, d and f, g, respectively), the situation is drastically different. Along the path parallel to kx the weight of the Ti-dxz orbital is negligible at the M/M’ pockets, while the Ti-dyz orbital contributes significantly. Along the rotated path the band at the M/M′ pocket is comprised of the Ti-\({d}_{{z}^{2}}\) and the Ti-dxz orbital with almost equal weight, with a smaller contribution from the Ti-dyz orbital. It is important to note that while the band structure possesses six-fold rotational symmetry in momentum space, the orbital texture exhibits a three-fold symmetry corresponding to the lattice symmetry.

Fig. 5: Orbital-projected band structure of monolayer 1T-TiTe2.
figure 5

a First surface Brillouin zone (black) and momentum cut (dashed lines). b–d Fat-band representation of the band structure projected onto the \({d}_{{z}^{2}}\), dxz and dyz orbitals. The line thickness represents the weight of the respective orbital. e–g Fat-band representation of the orbital weight for the lower cut in (a).

This momentum-dependent orbital content of the band structure underlines the important role of the orbital degree of freedom in photoemission from 1T-TiTe2. For exploring the manifestations in the ARPES signal, we combined the TB model with modeling the Wannier functions and the plane-wave (PW) approximation to the final photoelectron states (see the “Methods” section). We calculated the photoemission spectrum within the resulting TB + PW model for both sample orientations (see the “Methods” section) and extracted the iLDAD. Including specific orbitals only then allows for isolating the individual channels giving rise to the measured iLDAD.

First, in Fig. 6a, we show the iLDAD obtained from the TB + PW model including the Ti-\({d}_{{z}^{2}}\), dxz and dyz orbitals. The calculated iLDAD is in excellent qualitative agreement with the measured signal near the M/M′ pockets, which confirms that these three orbitals are the main ingredients of the iLDAD signal. Moreover, the coherent superposition of these orbitals is important: when incoherently summing \(i{{{{{\rm{LDAD}}}}}}=i{{{{{{\rm{LDAD}}}}}}}_{{z}^{2}}+i{{{{{{\rm{LDAD}}}}}}}_{xz}+i{{{{{{\rm{LDAD}}}}}}}_{yz}\), iLDAD is extremely weak and does not exhibit the expected symmetry (as experimentally measured), as shown in Fig. 6b. Including only the \({d}_{{z}^{2}}\) and dyz orbitals (Fig. 6c) captures the full signal along the kx axis, which is consistent with the corresponding orbital weight (see Fig. 5). In contrast, including only the \({d}_{{z}^{2}}\) and dxz orbitals (Fig. 6d) accounts for most of the iLDAD for the upper and lower two pockets, respectively. Note that the Ti-d orbitals are defined in the laboratory coordinate system here, which is the natural choice when computing ARPES signals. Earlier experiments50,51 attempted to connect the photoemission signal to the symmetry-adapted crystal-field orbitals. The crystal axes defined by the Ti–Te bonds are not strictly orthogonal for the presented experimental geometry, which complicates this picture. In the SI, we present a similar orbital analysis as above in the crystal-field basis for a complementary picture.

Fig. 6: Role of orbital mixing in photoemission.
figure 6

a iLDAD calculated using the TB + PW model where the initial state is composed of a coherent superposition of \({d}_{{z}^{2}}\), dxz and dyz orbitals. b Same as (a), for the initial state being an incoherent sum the signal including the \({d}_{{z}^{2}}\), dxz or dyz orbitals. c Same as (a) but including only \({d}_{{z}^{2}}\) and dyz orbitals. d Same as (a) for the initial state being a coherent combination of \({d}_{{z}^{2}}\), and dxz orbitals. The binding energy in all plots is fixed at EB = −0.05 eV; the photon energy is ℏω = 21.7 eV as in the experiments.

The TB + PW model allows for an even more intuitive picture for the characteristic iLDAD by inspecting the real-space dependence of the wavefunction. As explained in the “Methods” section, the photoemission matrix elements can be expressed in terms of the Wannier functions (for the Ti-d orbitals) and eigenvectors Cjn(k). Directly at the six M/M’ points on the edge of the Brillouin zone (labeled by km, m = 1, …, 6; see Fig. 7a), the matrix element within the TB + PW model reduces to

$$M({{{{{{\bf{k}}}}}}}_{m},E)=\int {{{\mathrm{d}}}}{{{{{\bf{r}}}}}}\ {{{{\mathrm{e}}}}}^{-i{{{{{\bf{p}}}}}}\cdot {{{{{\bf{r}}}}}}}{{{{{\bf{e}}}}}}\cdot {{{{{\bf{r}}}}}}{\widetilde{\psi }}_{m}({{{{{\bf{r}}}}}}),$$
(4)

where e denotes the polarization vector, while p = (km,x, km,y, p⊥) is the three-dimensional momentum of the photoelectron. The hybridized initial state \({\widetilde{\psi }}_{m}({{{{{\bf{r}}}}}})\) is given by the superposition

$${\widetilde{\psi }}_{m}({{{{{\bf{r}}}}}})=\mathop{\sum}\limits_{j}{C}_{jn}({{{{{{\bf{k}}}}}}}_{m}){\phi }_{j}({{{{{\bf{r}}}}}})\,,$$
(5)

where the band index n corresponds to the band near EF. The magnitude of the matrix element (4) is thus determined by the hybrid orbital (5), in particular its orientation with respect to the photoelectrom momentum p. The hydrid orbital (5) corresponds to taking the R = 0 contribution in Eq. (3). Including only the well-localized Ti-d orbitals, we find \({\widetilde{\psi }}_{m}({{{{{\bf{r}}}}}})\approx N{\psi }_{{{{{{\bf{{k}}}}}_{m}}}n}({{{{{\bf{r}}}}}})\) (up to a trivial phase) in the unit cell; the hybrid orbital (5) is thus a good representation of the Bloch wavefunction.

Fig. 7: Physical origin iLDAD: momentum-dependent orbital orientation.
figure 7

a Sketch of the calculated (TB + PW) iLDAD with labeled M/M' pockets (1–6). b Sketch of the hybrid orbital \({\widetilde{\psi }}_{m}({{{{{\bf{r}}}}}})\), for each M/M' points shown in (a). The orbital is represented by the constant-value contour \(| {\widetilde{\psi }}_{m}({{{{{\bf{r}}}}}})| =0.02\), while the red (blue) color shading indicates \({{{{{\rm{Re}}}}}}[{\widetilde{\psi }}_{m}({{{{{\bf{r}}}}}})] \, > \, 0\) (\({{{{{\rm{Re}}}}}}[{\widetilde{\psi }}_{m}({{{{{\bf{r}}}}}})] \, < \, 0\)). The coordinate systems corresponds to the laboratory frame. c projection of the hybrid orbitals onto the x–z plane. The photoelectron wave vector p is indicated by the black arrow.

In Fig. 7b we plot \({\widetilde{\psi }}_{m}({{{{{\bf{r}}}}}})\) for the α = 0∘ sample orientation. Inspecting Fig. 7b lead us to the following observation: the electronic wavefunction can be qualitatively described by tilted \({d}_{{z}^{2}}\) orbital, where the tilt angle is determined by the coefficients Cjn(km) and is thus momentum-dependent. This rotated \({d}_{{z}^{2}}\) is mainly composed of the \({d}_{{z}^{2}}\), dxz and dyz orbitals in the laboratory frame. By projecting the hybrid wavefunction at each M/M’ point onto the x–z plane and comparing their relative orientation to the outgoing electron momentum p, the origin of the strong iLDAD becomes clear. For m = 1, 3, 5 the lobes of \({\widetilde{\psi }}_{m}({{{{{\bf{r}}}}}})\) align with p, leading to enhanced photoemission intensity, while for m = 2, 4, 6 the photoelectron momentum is roughly pointing along the nodal direction of the hybrid orbitals. Rotating by 60° or equivalently transforming z → −z inverts this behavior, as the alignment of p with the orientation of the lobes (nodal direction) transforms into alignment with the nodal direction (orientation of the lobes).

The iLDAD is thus sensitive to the orbital orientation. In contrast to regular linear dichroism, the sample rotation \({{{{{{\mathcal{R}}}}}}}_{6{0}^{\circ }}\) changes the relative orientation of photoelectron momentum and orbitals while keeping the light polarization fixed. Revisiting Eq. (4), the photoemission matrix element (within the PW approximation) in essence is a Fourier transform of the hybrid orbital. From the momentum-space dependece of iLDAD, we gain insights into the real-space orientation of the orbitals. This information is complementary to varying the light polarization, where the matrix element is determined by the selection rules. In particular, orbital tilting as in 1T-TiTe2 is hard to discern directly from standard linear dichroism—when comparing the photoemission intensity modulation upon swapping between s- and p-polarized light.

Discussion

We have extended the recently introduced differential measurement methodology in multidimensional photoemission spectroscopy to investigate the role of orbital texture in photoemission from a semimetallic layered TMDC (1T-TiTe2). Indeed, we introduced the intrinsic linear dichroism in photoelectron angular distributions (iLDAD)—a generalized version of TRDAD, in samples featuring inversion symmetry. By tracking the photoemission intensity modulation upon specific crystal rotation in the laboratory frame, we experimentally extracted the intrinsic linear dichroism in the photoelectron angular distribution. Comparing our results with two different theoretical calculations based on a quantitative one-step model of photoemission and an intuitive tight-binding model, we provide an intuitive interpretation for the link between orbital texture and the pronounced photoemission intensity anisotropies from momentum regions with equivalent energy eigenvalues but with nonequivalent electronic wavefunctions. Our joint experimental and theoretical work demonstrated the direct relation between momentum-dependent hybridized orbital orientations and photoemission intensities, unambiguously ruling out final state effects as the primary origin of the observed dichroism50,51.

The methodology that we used is general and can be exploited to investigate a plethora of phenomena governed by orbital physics. Indeed, one expects very similar types of orbital texture to be common in governing properties of similar TMDC exhibiting multi-component charge order (e.g. 1T-TiSe252 and 1T-VSe253), so that our methodology can be used to investigate the link between the orbital texture of the emergence of CDW38,54. For example, we envision that performing a temperature-dependent iLDAD investigation for monolayer 1T-TiTe2 across its CDW transition might allows to shed light on the possible intricate nature of the orbital order and the emergence of CDW phases. In particular, the proposed excitonic insulator state55,56 gives rise to the admixture of the orbital character from the Γ point at the M/M’ pockets. Hence, the iLDAD is expected to exhibit pronounced changes, thus serving as a hallmark of the emergence of the excitonic insulator phase. Our approach could potentially also be used to investigate the orbital texture of exotic and stable hidden quantum states (e.g. in 1T-TaS257) which are not accessible using equilibrium tuning parameters.

Moreover, applying stress to the material investigated here (1T-TiTe2) has been predicted to drive strain-induced topological phase transitions58 associated with a change of the orbital texture. Strain-dependent measurement of the iLDAD could thus allow providing energy- and momentum-resolved view on the modification of orbital texture underlying topological phase transitions. This approach could also be used to investigate the role of orbital texture in the emergence of spontaneous electronic nematic ordering giving rise to the momentum-anisotropy of superconducting gap, for example in FeSe59.

One additional route that we envision is the ultrafast time-resolved investigation of dynamical modification of orbital texture upon impulsive photoexcitation of solids. Angle-resolved photoemission spectroscopy in general, and more specifically the experimental setup used in this study, is compatible with time-resolved studies with femtosecond temporal resolution40,42,60. Thus, combining this differential observable with pump-probe schemes could allow following in real-time modification of orbital texture during non-equilibrium topological phase transitions61,62, photoinduced orbital order63 and coherent phonon excitation breaking fundamental symmetries of the materials64.

Methods

Angle-resolved photoemission spectroscopy

The experimental apparatus features a table-top femtosecond XUV beamline coupled to a photoemission end-station. Briefly, a home-built optical parametric chirped-pulse amplifier (OPCPA) delivering 15 W (800 nm, 30 fs) at 500 kHz repetition rate65 is used to drive high-order harmonic generation (HHG) by tightly focusing the second harmonic of the laser pulses (400 nm) onto a thin and dense Argon gas jet. The nonperturbative nonlinear interaction between the laser pulses and the Argon atoms leads to the generation of a comb of odd harmonics of the driving laser, extending up to the 11th order. A single harmonic (7th order, 21.7 eV) is isolated by reflection on a focusing multilayer XUV mirror and propagation through a 400 nm-thick Sn metallic filter. A photon flux of up to 2 × 1011 photons/s at the sample position is obtained (110 meV FWHM). The bulk 1T-TiTe2 samples (HQ Graphene) were cleaved at room temperature and base pressure of 5 × 10−11 mbar, and handled by a six-axis manipulator (SPECS GmbH). The photoemission data is acquired using a time-of-flight momentum microscope (METIS1000, SPECS GmbH). This detector allows for simultaneous detection of the full surface Brillouin zone, over an extended binding energy range, without the need to rearrange the sample geometry41. For the data post-processing, we use a recently developed open-source workflow66 to efficiently convert the raw single-event-based datasets into binned calibrated data hypervolumes of the desired dimension, including axes calibration and artifact corrections (including symmetry distortion corrections67). The resulting 3D photoemission intensity data have the coordinates I(E, kx, ky).

DFT and tight-binding model

We have performed DFT calculations for both the bulk and the monolayer 1T-TiTe2 using the QUANTUM ESPRESSO code68. All calculations were performed using the PBE functional and norm-conserving pseudopotentials from the PSEUDODOJO project69. Self-consistent calculations were converged on a 15 × 15 Monkhorst–Pack grid for the Brillouin zone. For the monolayer, we constructed projective Wannier on the same grid with the WANNIER90 code70. We included the Ti-d orbitals and the p orbitals on each of the two Te atoms, yielding an 11-band TB model. This procedure yields the TB Hamiltonian \({H}_{j{j}^{\prime}}^{{{{{{\rm{TB}}}}}}}({{{{{\bf{k}}}}}})\) and thus the coefficients defining the Bloch wavefunction in Eq. (3). The real-space dependence of the Wannier functions is approximated by

$${\phi }_{j}({{{{{\bf{r}}}}}})={R}_{j}(r){X}_{{\ell }_{j}{m}_{j}}(\hat{{{{{{\bf{r}}}}}}})\,,$$
(6)

where \({X}_{\ell m}(\hat{{{{{{\bf{r}}}}}}})\) denotes the real spherical harmonics. Due to the small weight of the p orbitals we only construct the Ti d orbitals explicitly; the radial dependence entering Rj(r) is taken from the \({d}_{{z}^{2}}\) Wannier orbital and fixed for all d orbitals. We use atomic units in what follows.

From the Wannier functions (6) and the eigenvectors Cjn(k), we construct the Bloch states via Eq. (3), which allows us to calculate the ARPES intensity according to Fermi’s Golden rule

$$I({{{{{\bf{k}}}}}},E)=\mathop{\sum}\limits_{n}{f}_{n}({{{{{\bf{k}}}}}})| {M}_{n}({{{{{\bf{k}}}}}},E){| }^{2}\delta ({\varepsilon }_{n}({{{{{\bf{k}}}}}})+\omega -E)\,,$$
(7)
$${M}_{n}({{{{{\bf{k}}}}}},E)=\langle {{{{{\bf{k}}}}}},E| {{{{{\bf{e}}}}}}\cdot \hat{{{{{{\bf{D}}}}}}}| {\psi }_{{{{{{\bf{k}}}}}}n}\rangle \,.$$
(8)

Here, fn(k) is the occupation of band n with energy εn(k), ω is the photon energy, and E denotes the energy of the final state \(\left|{{{{{\bf{k}}}}}},E\right\rangle\). The band structure is shifted to reproduce the experimentally measured work function. The dipole operator is evaluate in the dipole gauge (\(\hat{{{{{{\bf{D}}}}}}}={{{{{\bf{r}}}}}}\)), which is defined in terms of the Berry connection. We ignore nonlocal dipole contributions since photoemission at the M/M’ pockets is predominantly originating from the Ti d orbitals. The final states are approximated by plane waves. In this set-up, the dipole gauge has proven to be more accurate than the velocity gauge (\(\hat{{{{{{\bf{D}}}}}}}=-i\nabla\)) for a number of related systems14,24,71. Within this TB + PW model, the photoemission matrix element is then given by

$${M}_{n}({{{{{\bf{k}}}}}},E)=\mathop{\sum}\limits_{j}{C}_{jn}({{{{{\bf{k}}}}}})\int {{{\mathrm{d}}}}{{{{{\bf{r}}}}}}{{{{\mathrm{e}}}}}^{-i{{{{{\bf{p}}}}}}\cdot {{{{{\bf{r}}}}}}}{{{{{\bf{e}}}}}}\cdot {{{{{\bf{r}}}}}}{\phi }_{j}({{{{{\bf{r}}}}}})\,.$$
(9)

Here, p = (kx, ky, p⊥) is the three-dimensional momentum of the photoelectron; the out-of-plane component p⊥ is fixed by the kinetic energy E = p2/2. The orbital-resolved photoemission analysis in the section “The role of orbital texture in photoemission: tight-binding analysis” is performed by including a subset of \(j\in \{{d}_{{z}^{2}},{d}_{xz},{d}_{yz},{d}_{{x}^{2}-{y}^{2}},{d}_{xy}\}\) orbitals in Eq. (9).

One-step model of photoemission: KKR method

The one-step model of photoemission is implemented in the spin-polarized relativistic Korringa–Kohn–Rostoker (SPR-KKR) scheme of the Munich bandstructure software package, based on Green’s function and multiple scattering spin-density matrix formalisms48,49. Relativistic phenomena such as SOC are treated by the Dirac equation. The LDA has been chosen to approximate the exchange-correlation part of the potential. The bulk potential is converged within the atomic sphere approximation (ASA) geometry where the empty spheres were used. We have used a lattice constant of 3.779 Å which corresponds to the interlayer Te–Te distance of 4.015 Å. During the tuning procedure of the theoretical model, it was necessary to modify the Wigner–Seitz radius of individual atomic types to the following ratio: Ti = 1.03, Te = 1.08 and the vacuum type = 1. We used an angular momentum expansion up to a lmax = 3. After the self-consistency was reached, the ARPES calculations have been performed and are based on the one-step model of photoemission in its spin density matrix formulation using the same geometry as used in the experiments. This theory accounts for effects induced by the light polarization, matrix-element effects, final state effects and surface effects. The final states are described by TRLEED states. These states are calculated by multiple scattering formalism and it are considered to give the best single particle description of final states in photoemission. In particular, this approach goes far beyond free-electron-like approximation for the final state. Lifetime effects in the final states have been simulated via a constant imaginary part V0f = 1 eV in the inner potential and the lifetime of the initial state was simulated by an imaginary part of V0i = 0.01 eV. The photoemission signal factors in all the matrix element effects such as experimental geometry, photon energy, polarization state, and surface termination or final state effects.