MATLAB Code: Raman Scattering Intensity for Quadratic Hamiltonians

These are codes were used to generate the Raman scattering intensity spectra of Kitaev Spin Liquid models using the Loudon-Fleury approach. In its most basic form this code diagonalizes a quadratic fermionic Hamiltonian and computes the spectra by constructing the Raman operator, using the eigenfunctions to compute matrix elements, and the eigenvalues to plot the spectrum as a function of energy. Variants are included that consider 2D and 3D lattices, finite systems, as well as the resonant Raman scattering. Neither interactions nor the bosonic case are considered.