--- title: 'Trilayer WTe₂: Structure, Semimetallicity & More' url: https://www.emergentmind.com/topics/trilayer-wte2 type: topic --- # Trilayer WTe₂: Structure, Semimetallicity & More Searching arXiv for relevant papers on trilayer WTe₂ to ground the article in the literature. Trilayer WTe\(_2\) is a three-monolayer film of tungsten ditelluride in the distorted orthorhombic \(T_d\) phase whose properties occupy an intermediate regime between atomically thin insulating layers and bulk semimetallic or Weyl behavior. In this thickness range, the material retains a distinct layer-resolved Raman fingerprint, exhibits ambient degradation that is slower than in monolayer and bilayer flakes because oxidation is confined to a self-limiting surface region, and undergoes an electronic reconstruction in which the indirect gap present in thinner films is suppressed so that valence and conduction states overlap at the Fermi level [1508.03244] [1608.00097] [2103.11885] [2604.16860]. Trilayer WTe\(_2\) is therefore a focal thickness for studies of dimensional crossover, interlayer-coupling-driven topology, Berry-curvature-related transport, and environmental stability. ## 1. Crystal structure and layer-specific identity Trilayer WTe\(_2\) consists of three \(T_d\)-phase WTe\(_2\) sheets stacked along the \(c\)-axis by van-der-Waals forces. Each monolayer is a W plane sandwiched by two Te planes in a distorted octahedral coordination [1608.00097]. In the structural description used for few-layer electronic studies, each WTe\(_2\) monolayer in the \(T_d\) phase is a distorted \(1T'\) unit with local inversion symmetry, and adjacent layers stack with a \(180^\circ\) (\(\pi\)) rotation about the \(z\)-axis and no net in-plane shift [2103.11885]. Bulk \(T_d\)-WTe\(_2\) has space group \(Pmn2_1\); for \(n=3\), the real structure has point group \(C_s\) with a single \(yz\) mirror, while the effective in-plane potential gradient alternates in sign layer by layer [2103.11885]. The lattice is strongly anisotropic. W atoms displace from ideal octahedral sites to form zigzag chains along the \(a\)-axis, yielding quasi-1D bonding and weaker coupling along the \(b\)- and \(c\)-axes [1508.03244]. This anisotropy is central to the lattice dynamics of the trilayer because one prominent Raman mode corresponds to vibrations along the tungsten-chain direction and is nearly insensitive to thickness reduction [1508.03244]. A practical consequence is that trilayer WTe\(_2\) can be distinguished spectroscopically from thinner films. In the Raman study of Kim et al., the four principal first-order Raman peaks observed in bulk WTe\(_2\) are fully conserved down to three layers, whereas bilayer and monolayer progressively lose one or two of these peaks [1508.03244]. This makes trilayer a spectroscopically well-defined thickness even before considering its electronic structure. ## 2. Raman and lattice-dynamical signatures Two complementary Raman descriptions appear in the literature. In one, the trilayer vacuum Raman fingerprint at room temperature shows four robust first-order peaks labeled \(P1\)–\(P4\): \(P1 \simeq 105\ \mathrm{cm^{-1}}\), \(P2 \simeq 161\ \mathrm{cm^{-1}}\), \(P3 \simeq 208\ \mathrm{cm^{-1}}\), and \(P4 \simeq 224\ \mathrm{cm^{-1}}\), assigned respectively to an in-plane W–Te vibration, an out-of-plane Te–Te/W–Te mixed mode, an in-plane W–Te bending mode, and an out-of-plane breathing mode [1608.00097]. In the layer-resolved Raman work of Kim et al., the four principal first-order peaks for the trilayer are reported as \(A_1^3\) at \(121\ \mathrm{cm^{-1}}\), \(A_2^4\) at \(138\ \mathrm{cm^{-1}}\), \(A_1^7\) at \(166\ \mathrm{cm^{-1}}\), and \(A_1^9\) at \(214\ \mathrm{cm^{-1}}\), with relative intensities of approximately \(90\%\), \(92\%\), \(145\%\), and \(145\%\) of the bulk values, respectively [1508.03244]. Both descriptions agree that trilayer WTe\(_2\) retains four clear first-order Raman modes. The layer dependence is anomalous relative to hexagonal dichalcogenides. As the layer number decreases from bulk to trilayer to monolayer, the \(A_3\), \(A_4\), and \(A_9\) modes blueshift, while the \(A_7\) mode remains almost fixed in frequency [1508.03244]. Specifically, the measured bulk-to-3L-to-2L-to-1L evolution is \(\sim 118 \to 121 \to 122\ \mathrm{cm^{-1}}\) for \(A_3\), \(\sim 134 \to 138 \to 138\ \mathrm{cm^{-1}}\) for \(A_4\), \(\sim 164 \to 166 \to 166\ \mathrm{cm^{-1}}\) for \(A_7\), and \(\sim 212 \to 214 \to 215\ \mathrm{cm^{-1}}\) for \(A_9\) [1508.03244]. The negligible shift of \(A_7\) is attributed to its atomic-displacement pattern along the one-dimensional tungsten chains, so it is almost unaffected by interlayer dielectric screening [1508.03244]. First-principles phonon calculations support this interpretation. Bulk, bilayer, and trilayer phonons were computed using density functional theory as implemented in VASP + PHONON, with LDA for exchange-correlation, PAW pseudopotentials, cutoff energy \(340\ \mathrm{eV}\), and a Monkhorst–Pack \(24\times14\times6\) mesh for dielectric and Raman calculations [1508.03244]. The phonon eigenvalue problem at \(\Gamma\) is written as $$ \sum_j D_{ij}(0)\,\xi_{j,\nu}=\omega_\nu^2\,\xi_{i,\nu}, $$ and Raman intensities are extracted from derivatives of the dielectric tensor: $$ I_\nu \propto \left|\, e_s^i \left(e_i^j \frac{\partial \varepsilon_{ij}}{\partial Q_\nu}\right)\right|^2. $$ The calculations reproduce the observed trend that all modes blueshift except \(A_7\) [1508.03244]. For identification purposes, trilayer WTe\(_2\) is characterized by the simultaneous presence of four peaks near \(121\), \(138\), \(166\), and \(214\ \mathrm{cm^{-1}}\), with \(A_7\) nearly unshifted from bulk and \(A_7\) and \(A_9\) reaching maximum integrated intensities at \(N=3\) [1508.03244]. ## 3. Electronic structure and the insulator–semimetal transition Trilayer WTe\(_2\) is the critical thickness at which few-layer WTe\(_2\) becomes semimetallic. Laser-\(\mu\)-ARPES on exfoliated, graphene-encapsulated flakes strictly sorted by layer number shows that, along \(\Gamma\)–\(X\), the highest valence band and lowest conduction band in 3-layer WTe\(_2\) approach and overlap at \(E_F\) [2103.11885]. The momentum-dependent overlap or gap is defined as $$ \Delta E(k)=E_c(k)-E_v(k). $$ Experimentally, \(\Delta E(k)\) changes sign as a function of \(k_x\): near \(k_x \approx 0\) (\(\Gamma\)), the valence-band maximum lies just below \(E_F\) and \(\Delta E(k)>0\), whereas at \(k_x \approx -0.32\) to \(-0.36\ \mathrm{\AA^{-1}}\), weak ARPES intensity from the conduction-band bottom crosses \(E_F\) and \(\Delta E(k)<0\) [2103.11885]. The numerical scale of the overlap is small. At \(\Gamma\), \(E_v(0)\) is \(\sim -20\ \mathrm{meV}\) below \(E_F\), with no conduction band observed up to \(+50\ \mathrm{meV}\), so \(\Delta E(0)\approx +20\) to \(+50\ \mathrm{meV}\) [2103.11885]. At \(k_x=-0.32\ \mathrm{\AA^{-1}}\), $$ \Delta E(-0.32)\approx (-0\ \text{to}\ +5\,\mathrm{meV})-(-10\,\mathrm{meV})\approx -5\text{ to }+5\,\mathrm{meV}, $$ that is, effectively zero or slightly negative overlap of order \(5\)–\(10\ \mathrm{meV}\), with a maximum semimetallic overlap of \(\sim 10\)–\(20\ \mathrm{meV}\) where both bands are best resolved [2103.11885]. This establishes trilayer WTe\(_2\) as a compensated semimetal at the onset of the \(2 \to 3\) layer transition. ARPES with \(h\nu=21.2\ \mathrm{eV}\) and \(T=10\ \mathrm{K}\) reports the same qualitative electronic reconstruction. Along the high-symmetry \(X\)–\(\Gamma\)–\(X\) and \(F\)–\(X\)–\(F\) directions, the top of the valence band crosses \(E_F\) in the trilayer, with no obvious gap, and second-derivative spectra confirm that the indirect gap closes between valence-band maximum and conduction-band minimum [2604.16860]. Energy-distribution-curve fitting yields no resolvable gap at \(\Gamma\) or along \(F\), in contrast to the \(\sim 55\ \mathrm{meV}\) indirect gap in the monolayer and \(\sim 30\ \mathrm{meV}\) gap in the bilayer [2604.16860]. First-principles calculations reproduce the measured semimetallicity. HSE06 in the \(\mu\)-ARPES study confirms the small overlap in \(n=3\) [2103.11885], while GGA+\(U\) with \(U_W=2\ \mathrm{eV}\) and HSE06 in the later thickness-evolution study find that the valence and conduction bands overlap by \(\sim 30\ \mathrm{meV}\) near \(\Gamma\), with a finite density of states \(N(E_F)\simeq 0.5\ \mathrm{states/(eV\cdot cell)}\), characteristic of a semimetal [2604.16860]. The valence band is composed primarily of W \(5d_{xy}\) and Te \(p_x\) orbitals, and the conduction-band bottom derives from W \(5d_{xz/yz}\) states [2604.16860]. ## 4. Interlayer coupling, symmetry, and topological characterization The thickness evolution of WTe\(_2\) is non-monotonic. Monolayer WTe\(_2\) is reported as insulating with indirect QSH gap \(\simeq 55\ \mathrm{meV}\) and \(Z_2=1\), bilayer WTe\(_2\) as a reduced-gap trivial insulator with \(Z_2=0\), and trilayer WTe\(_2\) as semimetallic with \(Z_2=1\) in the sense that valence and conduction bands remain separated in \(k\)-space near the time-reversal invariant momenta even though there is no global full gap [2604.16860]. In the bulk limit, the system becomes a 3D type-II Weyl semimetal with chiral Weyl nodes of \(C=\pm 1\) and Fermi-arc surface states [2604.16860]. The mechanism is interlayer-coupling-driven band reconfiguration. The van der Waals interlayer hopping \(t_\perp\) between W \(d\)-orbitals in adjacent layers is reported as \(\approx 50\ \mathrm{meV}\) [2604.16860]. As the layer number \(N\) increases, quantum-well subbands form, driving the valence band up and the conduction band down; a critical coupling occurs when the \(k\)-dependent mass term $$ M(k)=M_0+\delta k^2-t_\perp $$ changes sign, closing the gap at \(k\approx 0.1\ \mathrm{\AA^{-1}}\) [2604.16860]. In the effective Hamiltonian near \(\Gamma\), $$ H(k)=[M(k)\sigma_z + A k_x \sigma_x + A k_y \sigma_y]\otimes \tau_0 + t_\perp(N)\sigma_0\otimes \tau_x, $$ with \(\sigma_i\) acting on orbital parity and \(\tau_i\) on layer pseudospin, the mass term \(M(k)=M_0-B(k_x^2+k_y^2)\) changes sign when \(t_\perp(N=3)>M_0\), leading to gap closing at finite \(k\) [2604.16860]. For the trilayer, the \(Z_2\) invariant is defined through parity products at the four time-reversal invariant momenta \(\Lambda_i\in\{\Gamma,X,Y,M\}\): $$ Z_2=\prod_{i=1}^4 \delta(\Lambda_i)\ \mathrm{mod}\ 2, $$ with $$ \delta(\Lambda_i)=\prod_m^{N_{occ}}\xi_{2m}(\Lambda_i), $$ where \(\xi_{2m}(\Lambda_i)=\pm 1\) is the parity eigenvalue of the \(2m\)-th occupied band [2604.16860]. The reported DFT result is \(\delta(\Gamma)=-1\), \(\delta(X)=+1\), \(\delta(Y)=+1\), and \(\delta(M)=+1\), implying \(Z_2=1\), and hybrid Wannier-center evolution along \(\Gamma\)–\(X\) crosses the reference line an odd number of times, confirming the nontrivial invariant [2604.16860]. Because the trilayer is globally gapless, global Chern numbers are ill-defined, but local Berry curvature can still be computed around avoided crossings to track the approach to a Weyl transition [2604.16860]. The even–odd effect in layer number is a recurrent theme. The reported interpretation is that odd \(N\) yields net inverted band ordering and \(Z_2=1\), while even \(N\) yields \(Z_2=0\), before semimetallicity emerges when \(N\ge 3\) and eventually connects to bulk Weyl physics [2604.16860]. This places trilayer WTe\(_2\) at the first odd-layer semimetallic member of the sequence. ## 5. Spin splitting, Berry curvature, and nonlinear transport context Trilayer WTe\(_2\) also occupies a distinct symmetry regime in spin-resolved and Berry-curvature-related phenomenology. In the \(\mu\)-ARPES measurements, the spin-resolved splitting at momentum \(k\) is defined as $$ \Delta E_s(k)=E_{\uparrow}(k)-E_{\downarrow}(k). $$ For \(n=3\), the valence bands show only very small or unresolvable splittings; fits of energy-distribution curves at \(k_x\approx -0.1\ \mathrm{\AA^{-1}}\) give $$ \Delta E_s(-0.1)\lesssim 30\,\mathrm{meV}, $$ typically in the \(0\)–\(30\ \mathrm{meV}\) range, with most bands \(\lesssim 15\ \mathrm{meV}\) [2103.11885]. This contrasts with even-\(n\) films, where splittings up to \(70\ \mathrm{meV}\) are observed [2103.11885]. The structural interpretation is that incomplete cancellation of local dipoles in even-\(n\) layers gives large spin–orbit splittings, whereas in odd-\(n\) trilayers much of the triangular-Te dipole arrangement partially cancels, leaving only small residual \(\Delta E_s\) [2103.11885]. This is consistent with the recovered partial symmetry of the middle layer in trilayer stacking and with the weaker effective structural asymmetry relative to even-\(n\) films [2103.11885]. Few-layer WTe\(_2\) permits nonzero Berry curvature because inversion symmetry is broken. The band Berry curvature is written as $$ \Omega_n(\mathbf{k})=-2\,\mathrm{Im}\sum_{m\neq n}\frac{\langle u_{n\mathbf{k}}|\partial_{k_x}H|u_{m\mathbf{k}}\rangle \langle u_{m\mathbf{k}}|\partial_{k_y}H|u_{n\mathbf{k}}\rangle}{(\varepsilon_{m\mathbf{k}}-\varepsilon_{n\mathbf{k}})^2}, $$ and the Berry-curvature dipole as $$ D_a=\sum_{n,\mathbf{k}} f_{n,\mathbf{k}}\,\partial_{k_a}\Omega_n(\mathbf{k}). $$ Although \(\Omega_n(\mathbf{k})\) and \(D_a\) are not directly reported for trilayer in the ARPES work, previous transport measurements on 3-layer WTe\(_2\) are stated to have revealed nonlinear anomalous Hall currents and ferroelectric switching, both interpreted there as fingerprints of a nonzero \(D_y\) [2103.11885]. This suggests that trilayer WTe\(_2\) couples its near-compensated semimetallicity to inversion-breaking Berry-curvature physics, albeit with weaker spin splitting than in even-layer films. A plausible implication is that trilayer WTe\(_2\) is experimentally important not because it maximizes a single order parameter, but because it combines semimetallic overlap, odd-layer topological character, and finite inversion-breaking responses within the same thickness window [2103.11885] [2604.16860]. ## 6. Environmental instability, oxidation chemistry, and passivation Trilayer WTe\(_2\) is substantially more stable in air than monolayer or bilayer WTe\(_2\), but it still degrades by oxidation. Raman monitoring after exfoliation shows that 3L WTe\(_2\) is stable in vacuum, whereas in ambient air each Raman mode decays approximately exponentially toward a nonzero plateau [1608.00097]: $$ I(t)=I_1\exp[-(t-t_0)/\tau]+I_2, $$ where \(t\) is elapsed time since exfoliation, \(t_0\approx 0\)–\(5\ \mathrm{min}\) accounts for measurement delay, \(I_1\) is the Raman intensity that decays, \(I_2\) is the residual saturation intensity, and \(\tau\) is the characteristic decay time [1608.00097]. For trilayer WTe\(_2\), \(\tau_{3L}\approx 7{,}000\)–\(9{,}000\ \mathrm{min}\) (\(\approx 5\)–\(6\) days), and residual intensities \(I_2/I_1\approx 20\)–\(30\%\) remain after two weeks [1608.00097]. This decay is much slower than in thinner flakes. Monolayer WTe\(_2\) shows complete disappearance of Raman modes within \(\sim 13\ \mathrm{min}\) in air, bilayer WTe\(_2\) has \(\tau_{2L}\approx 140\)–\(160\ \mathrm{min}\) with saturation after \(\sim 3\ \mathrm{h}\), and trilayer requires \(\sim 2\) weeks for full saturation [1608.00097]. The thickness dependence is explained by the fact that only the topmost layer oxidizes and the oxidized surface constitutes a smaller fraction of the total volume as the layer number increases [1608.00097]. Surface analysis identifies the chemistry. In 3L WTe\(_2\) aged for \(\sim 15\) days, XPS before etching shows new Te \(3d_{5/2,3/2}\) oxide components at \(576.1\ \mathrm{eV}\) and \(587.1\ \mathrm{eV}\), assigned to TeO\(_2\), and W \(4d_{5/2,3/2}\) oxide components at \(247.4\ \mathrm{eV}\) and \(260.6\ \mathrm{eV}\), assigned to WO\(_x\) with \(2