Emergent Kitaev materials in synthetic Fermi-Hubbard bilayers

We investigate the emergence of bond-directional spin-spin interactions in a synthetic Fermi-Hubbard bilayer that can be realized with ultracold fermions in Raman optical lattices. The model exploits synthetic dimensions to couple two honeycomb layers, each corresponding to a different hyperfine atomic state, via Raman-assisted tunneling and, moreover, via an inter-layer Hubbard repulsion due to the cold-atom scattering. In the strong-coupling regime at half filling, we derive effective spin Hamiltonians for the kinetic exchange featuring Kitaev, Heisenberg, off-diagonal exchange (-couplings), as well as tunable Dzyaloshinskii-Moriya interactions. We identify specific configurations that generate both ferromagnetic and antiferromagnetic Kitaev couplings with various perturbations of relevance to Kitaev materials, providing a tunable platform that can explore how quantum spin liquids emerge from itinerant fermion systems. We analyze the Fermi-liquid and Mott-insulating phases, highlighting a correspondence between Dirac and Majorana quasi-particles, with possible phase transitions thereof. In an extreme anisotropic limit, we show that the model reduces to an inter-layer ribbon in a quasi-1D ladder, allowing for a numerical study of the correlated ground state using matrix product states. We find a transition from a symmetry-protected topological insulator to a Kitaev-like regime characterized by nonlocal string order. Our results establish that cold-atom quantum simulators based on Raman optical lattices can be a playground for extended Kitaev models, bridging itinerant fermionic systems and spin-liquid physics.