Location of conical intersections in solution using a sequential quantum mechanics/molecular dynamics method