with SO coupling has been employed using a range of functionals, including the statistical average of orbital potential (SAOP), CAM-B3LYP, and local density approximation (LDA). These studies have shown that TDDFT calculations with functionals other than SAOP often significantly deviate from results obtained through relativistic wavefunction approaches, suggesting that the SAOP functional provides a more robust and reliable description of excited states in the presence of significant SO coupling. Furthermore, Bast et al. demonstrated that when the SAOP functional is consistently applied, the impact of different relativistic treatments -such as 4c, X2C, and ZORA approach-is minimal.