We propose that resonant inelastic X-ray scattering (RIXS) is an effective probe of the fractionalized excitations in three-dimensional (3D) Kitaev spin liquids. While the non-spin-conserving RIXS responses are dominated by the gauge-flux excitations and reproduce the inelastic-neutron-scattering response, the spin-conserving (SC) RIXS response picks up the Majorana-fermion excitations and detects whether they are gapless at Weyl points, nodal lines, or Fermi surfaces. As a signature of symmetry fractionalization, the SC RIXS response is suppressed around the Γ point. On a technical level, we calculate the exact SC RIXS responses of the Kitaev models on the hyperhoneycomb, stripyhoneycomb, hyperhexagon, and hyperoctagon lattices, arguing that our main results also apply to generic 3D Kitaev spin liquids beyond these exactly solvable models.Quantum spin liquids (QSLs) are exotic and entirely quantum phases of matter [1,2] From a theoretical point of view, KSLs are particularly appealing because each of them has an exactly solvable limit governed by a Kitaev model [3]. In general, the Kitaev model is defined on a tricoordinated lattice with S = 1/2 spins σ x,y,z r at the sites r, which are coupled to their neighbors via bonddependent Ising interactions. The Hamiltonian readswhere J x,y,z are the coupling constants for the three types of bonds x, y, and z. Remarkably, this model is exactly solvable whenever there is precisely one bond of each type around each site of the tricoordinated lattice.These exactly solvable Kitaev models have been defined on a wide range of tricoordinated 3D lattices [4][5][6][7][8], including the hyperhoneycomb, stripyhoneycomb, hyperhexagon, and hyperoctagon lattices (see Fig. 1). In the experimentally relevant isotropic regime (J x ≈ J y ≈ J z ), the ground state is a gapless Z 2 QSL, while the (fractionalized) excitations are gapless Majorana fermions and gapped Z 2 gauge fluxes. Importantly, the Majorana fermions (spinons) exhibit a rich variety of nodal structures due to the different (projective) ways symmetries can act on them [5][6][7]. Indeed, they are gapless along nodal lines for the hyperhoneycomb and the stripyhoneycomb models [4], on Fermi surfaces for the hyperoctagon model [5], and at Weyl points for the hyperhexagon model [7].From an experimental point of view, however, it is difficult to identify and characterize QSLs due to the lack of any local order parameters that could be used as "smoking-gun" signatures. In recent years, a remarkable theoretical and experimental progress has been achieved in understanding that fractionalization is one of the most promising hallmarks of a QSL. Indeed, it has been demonstrated that fractionalized excitations, which are Majorana fermions and Z 2 gauge fluxes for KSLs, can be probed by conventional spectroscopic techniques, such as inelastic neutron scattering (INS) [26,[31][32][33][34], Raman scattering with visible light [21,25,[35][36][37][38][39], and resonant inelastic X-ray scattering (RIXS) [40][41][42].In this Letter, we propo...