- You must use the natural logarithm, not the base 10 one.
- Have you tried to compare the means and standard deviations of the model simulated with what you consider a sensible calibration to the moments from the data you are feeding in? That might give you a clue why a unit root is favored by the data.