Replication of "Optimal monetary policy in an estimated DSGE for the euro area"

[jpfeifer]

One problem you are almost surely facing is

An infinity of steady states with Taylor rules

Could it be that you are trying to use the steady state operator to compute something that cannot be set endogenously? Take inflation: every percentage point higher steady state inflation will increase the steady state nominal interest accordingly. But any consistent combination of both (via the link the Fisher equation provides) is feasible. You need to pin down one of them by setting the value. The model cannot do this for you.