The Selberg zeta function for a hyperbolic space X links the geodesic flow on X with the Laplace spectrum of X. On hyperbolic surfaces, the Selberg zeta function has been studied for both unitary and non-unitary twists. In this talk, I will extend this investigation to higher-dimensional hyperbolic spaces for non-unitary twists. I will introduce the hyperbolic spaces and present the current state of our research, focusing on the challenges inherent in higher dimensions.