Bug fix in vacancy concentration calculation in dilute_solution_model