Initializing variables for steady state

Hi Everyone! I am currently trying to run a model based on the.

The authors calibrated some values and used Bayesian estimation for the others, I am trying to initialize them just to arrive to a random steady state. I am getting NaN errors in all the residuals when running the steady function.
When I run resid 1, the Nan appear in a selection of variables but I could not find mistaked with my values. How do you initialize values in this case?

rbc.mod (5.6 KB)

a_t=a_t(-1)+ e_at;

implies that there are four unit roots and consequently infinitely many steady states. You cannot hope to simply find one numerically. You need to provide one analytical steady state.