A Monte Carlo study of symmetry breaking of I3− in aqueous solution using a multistate diabatic Hamiltonian