Unveiling the orbital texture of 1T-TiTe2 using intrinsic linear dichroism in multidimensional photoemission spectroscopy

The momentum-dependent orbital character in crystalline solids, referred to as orbital texture, is of capital importance in the emergence of symmetry-broken collective phases, such as charge density waves as well as superconducting and topological states of matter. By performing extreme ultraviolet multidimensional angle-resolved photoemission spectroscopy for two different crystal orientations linked to each other by mirror symmetry, we isolate and identify the role of orbital texture in photoemission from the transition metal dichalcogenide 1T-TiTe2. By comparing our experimental results with theoretical calculations based on both a quantitative one-step model of photoemission and an intuitive tight-binding model, we unambiguously demonstrate the link between the momentum-dependent orbital orientation and the emergence of strong intrinsic linear dichroism in the photoelectron angular distributions. Our results represent an important step towards going beyond band structure (eigenvalues) mapping and learning about electronic wavefunction and orbital texture of solids by exploiting matrix element effects in photoemission spectroscopy.


INTRODUCTION
Angle-resolved photoemission spectroscopy (ARPES) is the most direct experimental technique to measure the momentumresolved electronic eigenvalues (energy bands) of crystalline solids 1,2 . The technique is based on the photoelectric effect 3,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 states 6 . Because of its semimetallic 2D nature and its speculated textbook Fermi-Liquid behavior, the layered transition metal dichalcogenide (TMDC) 1T-TiTe 2 was the central material around the debate on directly linking the spectral function and the ARPES lineshape [7][8][9][10] . Indeed, 1T-TiTe 2 served as a benchmark material to demonstrate the signatures of many-body effects in ARPES lineshape 11 , the possibility to retrieve the photohole lifetime 12 , as well as the potential to dissect the role of the different quasiparticle scattering processes from ARPES measurements 10 .
On the other hand, very little attention has been devoted to the other key ingredient governing ARPES spectra of 1T-TiTe 2 : 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 momentumdependent modulation onto the photoemission intensity 13,14 . Matrix element effects can sometimes be confused with signatures of many-body interactions 6,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 structure 16 , 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 solids [17][18][19][20][21] , the chirality of charge carriers 22,23 , quantum geometric properties, such as orbital pseudospin texture 24,25 as well as Berry curvature [26][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 state [35][36][37] . As demonstrated recently, CDW phases and orbital texture can be intertwined 38 . Tracing the orbital texture from matrix element effects in photoemission will also greatly benefit discerning structural and collective phases in TMDCs. More specifically, 1T-TiTe 2 has been recognized to be a very interesting prototypical 2D material: it exhibits a strong thickness 33,39 and substratedependent 34  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 (R θ ) 24 . We benchmarked this approach by investigating the modulation of photoemission intensity from 2H-WSe 2 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 (R 60 ) is mimicking a time-reversal operation (T), i.e. R 60 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 R 60 , the antisymmetric components of the dichroism upon R 60 (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-WSe 2 24 . 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-TiTe 2 . 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-TiTe 2 , 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°1 7,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 momentumresolved electronic wavefunction of solids.

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 (k x , k y ). One of the major advantages of this type of detector is the parallel detection of the 3D photoemission intensity I(E, k x , k y ) in a single measurement, without having to rotate the sample in the laboratory frame 41 . More details about the experimental setup can be found elsewhere 40,42 and in the "Methods" section. Giving the energydependence of the inelastic mean free path of the outgoing electrons 43,44 , working in the XUV spectral range ensure that most of the detected electrons are coming from the topmost 1T-TiTe 2 layer (first atomic trilayer) 24 . More details about this statement, as well as its theoretical validation, will be given later.
The experimentally measured 3D photoemission intensities and associated constant energy contours (CECs) at the Fermi energy (E−E F = 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-TiTe 2 , 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-TiTe 2 for a second crystal orientation rotated by 60°. In inversion symmetric 1T-TiTe 2 , this crystal rotation (R 60 ) is strictly equivalent to a mirror-image symmetry of each layer (see Fig. 1b), i.e. R 60 z ! Àz. Although the effect of a 60°crystal rotation in 1T-TMDC (R 60 z ! Àz) is different than in 2H-TMDC (R 60 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.
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, k x , k y ) and backward I α (E, − k x , k y ) hemisphere, i.e. the linear dichroism asymmetry in the photoelectron angular distribution A α LDAD ðE; k x ; k y Þ, can be extracted (see Eq. (1)), without the need to rearrange the sample geometry nor the light-polarization state 46,47 .
Looking at A 0 LDAD (Fig. 2d), we notice very strong dichroism (up to ±97%) around M/M′ pockets with an opposite sign for adjacent S. Beaulieu et al. 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 (k y > 0/k y < 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 (k x > 0/k x < 0) symmetrized signal is > 0.1I max , where I max 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 60 LDAD , i.e.
the dichroism extracted after rotating the crystal orientation by 60 ∘ . The measured A 60 LDAD dichroic signal is almost perfectly antisymmetric with respect to A 0 LDAD , 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 α LDAD dichroism upon crystal rotation (R 60 ), i.e. to remove any spurious contributions from experimental geometry, we introduce the iLDAD: where the crystal rotation is hereby an angle α−α 0 = 60 ∘ . The resulting iLDAD is shown in Fig. 2f. The strong similarity between 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, k x , k y ), and associated constant energy contour (CEC) at the Fermi energy (E−E F = 0 eV). d Same as (c), but following a crystal rotation by 60°in the laboratory frame.
A 0 LDAD and iLDAD is another observation revealing that the dichroism in 1T-TiTe 2 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 method 48,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-TiTe 2 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-TiTe 2 . 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 wellknown semimetallic nature of 1T-TiTe 2 , 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−E F~− 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 effects 39 . 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.
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 k x (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 onestep 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 timereversed 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-TiTe 2 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-TiTe 2 monolayers (atomic trilayers) 45 .
To investigate the role of the final states, we performed a detailed theoretical photon-energy-dependent iLDAD scan (k zscan). 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 freeelectron-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 Fig. 3 Band structure of 1T-TiTe 2 along the M-Γ-M′ cut. 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 TiTe 2 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.
for this reason that we took advantage of our quantitative KKR theoretical framework to investigate the photon-energydependence 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-TiTe 2 . 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 where ϕ j (r) denote the Wannier functions localized at particular sites R. The coefficients C jn (k) are the coefficients connection orbital (index j) and band (index n) space; they are the eigenvectors of the TB Hamiltonian H TB jj 0 ð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 TiTe 2 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 E F 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-d xz and Ti-d yz orbital-projected band structure (Fig. 5c, d and f, g,  respectively), the situation is drastically different. Along the path parallel to k x the weight of the Ti-d xz orbital is negligible at the M/M' pockets, while the Ti-d yz 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-d xz orbital with almost equal weight, with a smaller contribution from the Ti-d yz 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.
This momentum-dependent orbital content of the band structure underlines the important role of the orbital degree of freedom in photoemission from 1T-TiTe 2 . 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 , d xz and d yz 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 iLDAD ¼ iLDAD z 2 þ iLDAD xz þ iLDAD 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 d yz orbitals (Fig. 6c) captures the full signal along the k x axis, which is consistent with the corresponding orbital weight (see Fig. 5). In contrast, including only the d z 2 Fig. 6 Role of orbital mixing in photoemission. a iLDAD calculated using the TB + PW model where the initial state is composed of a coherent superposition of d z 2 , d xz and d yz orbitals. b Same as (a), for the initial state being an incoherent sum the signal including the d z 2 , d xz or d yz orbitals. c Same as (a) but including only d z 2 and d yz orbitals. d Same as (a) for the initial state being a coherent combination of d z 2 , and d xz orbitals. The binding energy in all plots is fixed at E B = −0.05 eV; the photon energy is ℏω = 21.7 eV as in the experiments. and d xz 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 experiments 50,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.
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 C jn (k).
where e denotes the polarization vector, while p = (k m,x , k m,y , p ⊥ ) is the three-dimensional momentum of the photoelectron. The hybridized initial state e ψ m ðrÞ is given by the superposition e ψ m ðrÞ ¼ X j C jn ðk m Þϕ j ðrÞ ; where the band index n corresponds to the band near E F . 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 e ψ m ðrÞ % Nψ kmn ðrÞ (up to a trivial phase) in the unit cell; the hybrid orbital (5) is thus a good representation of the Bloch wavefunction.
In Fig. 7b we plot e ψ m ð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 C jn (k m ) and is thus momentum-dependent. This rotated d z 2 is mainly composed of the d z 2 , d xz and d yz 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 e ψ m ð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 R 60 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-TiTe 2 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-TiTe 2 ). 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 dichroism 50,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-TiSe 2 52 and 1T-VSe 2 53 ), so that our methodology can be used to investigate the link between the orbital texture of the emergence of CDW 38,54 . For example, we envision that performing a temperature-dependent iLDAD investigation for monolayer 1T-TiTe 2 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 state 55,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-TaS2 57 ) which are not accessible using equilibrium tuning parameters.
Moreover, applying stress to the material investigated here (1T-TiTe 2 ) has been predicted to drive strain-induced topological phase transitions 58 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 FeSe 59 .
One additional route that we envision is the ultrafast timeresolved 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 timeresolved studies with femtosecond temporal resolution 40,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 transitions 61,62 , photoinduced orbital order 63 and coherent phonon excitation breaking fundamental symmetries of the materials 64 .

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 rate 65 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 × 10 11 photons/s at the sample position is obtained (110 meV FWHM). The bulk 1T-TiTe 2 samples (HQ Graphene) were cleaved at room temperature and base pressure of 5 × 10 −11 mbar, and handled by a sixaxis 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 geometry 41 . For the data post-processing, we use a recently developed open-source workflow 66 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 corrections 67 ). The resulting 3D photoemission intensity data have the coordinates I(E, k x , k y ).

DFT and tight-binding model
We have performed DFT calculations for both the bulk and the monolayer 1T-TiTe 2 using the QUANTUM ESPRESSO code 68 . All calculations were performed using the PBE functional and norm-conserving pseudopotentials from the PSEUDODOJO project 69 . 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 code 70 . 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 TB jj 0 ðkÞ and thus the coefficients defining the Bloch wavefunction in Eq. (3). The real-space dependence of the Wannier functions is approximated by where X ℓm ð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 R j (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 C jn (k), we construct the Bloch states via Eq. (3), which allows us to calculate the ARPES intensity according to Fermi's Golden rule Iðk; EÞ ¼ X n f n ðkÞjM n ðk; EÞj 2 δðε n ðkÞ þ ω À EÞ ; M n ðk; EÞ ¼ hk; Eje ÁDjψ kn i : Here, f n (k) is the occupation of band n with energy ε n (k), ω is the photon energy, and E denotes the energy of the final state k; E j i. The band structure is shifted to reproduce the experimentally measured work function. The dipole operator is evaluate in the dipole gauge (D ¼ 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 (D ¼ Ài∇) for a number of related systems 14,24,71 . Within this TB + PW model, the photoemission matrix element is then given by M n ðk; EÞ ¼ dre ÀipÁr e Á rϕ j ðrÞ : Here, p = (k x , k y , p ⊥ ) is the three-dimensional momentum of the photoelectron; the out-of-plane component p ⊥ is fixed by the kinetic energy E = p 2 /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 2 fd z 2 ; d xz ; d yz ; d x 2 Ày 2 ; d xy g orbitals in Eq. (9).

One-step model of photoemission: KKR method
The one-step model of photoemission is implemented in the spinpolarized relativistic Korringa-Kohn-Rostoker (SPR-KKR) scheme of the Munich bandstructure software package, based on Green's function and multiple scattering spin-density matrix formalisms 48,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 l max = 3. After the selfconsistency 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, matrixelement 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 V 0f = 1 eV in the inner potential and the lifetime of the initial state was simulated by an imaginary part of V 0i = 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.

DATA AVAILABILITY
All experimental photoemission data used for the presented analysis are publicly available on Zenodo open-access repository 72 .

CODE AVAILABILITY
The codes used to perform data analysis are available on request from the authors.