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.