Langrange multipliers are an integral part of dsge models. Theoretically, they are always positive. But in dynare model section, we do not introduce anything to restrict the values of these multipliers to positive values. Does that imply that while we can ensure that their steady state values are positive, they can squander away into the negative territory when the model is away from the steady state? I have often see the langrange multipliers moving from negative to positive in the IRFs, hence the question.

That comes from doing a polynomial approximation. See e.g.