If you want to compare simulations paths as opposed to just moments, i.e. averages over a sample, then you need to use the same exogenous shocks. For that, you either need a seed or directly feed in the same exogenous shocks. So using the simult_-function is the way to go.