Experimental excitation energies (EEs) of low-lying states for uranyl halides are relatively independent of the type of the equatorial ligands, while many theoretical results indicate that the ligands have a significant effect on the electronic structure of uranyl. In this work, vertical and adiabatic EEs for [UO 2 X 4 ] 2À and UO 2 X 2 (X = F, Cl) are compared at the same computational level, that is, the equation-of-motion coupled-cluster method, to investigate effects of the number and type of ligands on EEs of uranyl. Furthermore, accuracy of a perturbative treatment of spin-orbit coupling (SOC) on EEs of these molecules is also studied. Our results show that reliable EEs are achieved when a sufficient number of spin-free excited states are considered in dealing with SOC perturbatively. Difference between the lowest adiabatic EEs for these molecules is calculated to be within 0.2 eV, which are in reasonable agreement with experimental results. On the other hand, difference between EEs for other higher states of [UO 2 X 4 ] 2À and those of UO 2 X 2 is sizeable. In addition, reasonable EEs can still be obtained when the ligands are replaced by negative point charges.Our results indicate that the number of ligands has a more pronounced effect on EEs of these uranyl chlorides and fluorides than the type of the ligands.