The performance of the six second order linear response methods RPA(D), SOPPA, SOPPA(CCSD), CIS(D), CC2, and CCSD, which include either noniterative or iterative doubles contributions, has been studied in calculations of vertical excitation energies. The benchmark set consisted of 39 valence and 76 Rydberg states of benzene and five polycyclic aromatic hydrocarbons. As reference values we have used the results of the corresponding calculations with the third order method CCSDR(3), which includes noniterative triples contributions. In addition we have also carried out equivalent calculations at the level of the random phase approximation as well as with the configuration interaction singles and multireference configuration interaction singles and doubles methods.