I've tried to implement NK-SOE model by Alpanda et al. (2010). First, I rewrote the equations and calibrated parameters to their posterior mean from the paper. IRFs look the same, however Dynare 4.6 says my model contains unit root, whether original did not as authors calculated asymptotic variance decomposition. I tried to run IRFs in horizon 1000 and all variables returned to the steady state.
I attach paper and the code.