A method for modeling electron dynamics in large chemical systems