We present a semiempirical Hamiltonian that provides an accurate description of the first singlet and triplet potential energy surfaces of azobenzene for use in direct simulations of the excited-state dynamics. The parameterization made use of spectroscopic and thermochernical data and the best ab initio results available to date. Two-dimensional potential energy surfaces based on constrained geometry optimizations are presented for the states that are most relevant for the photochemistry of azobenzene, namely, S(0), S(1), and S(2). In order. to run simulations of the photodynamics of azobenzene in hydrocarbons or hydroxylic solvents, we determined the interactions of methane and methanol with the azo group by ab initio calculations and fitted the interactions with a QM/MM interaction Hamiltonian.