running simulation with density profile and rate constant output