Scalar relativistic corrections were not included in the results to reduce computational cost. However, some preliminary investigations were performed. For instance, incorporating such effects in determining vertical excitation features at the TDDFT level resulted in a systematic shift in state energies but did not affect the observed trends or underlying characteristics. Similarly, for 2-TP aceto impl and 4-TP aceto impl models, 50 trajectories were simulated using the Douglas-Kroll-Hess Hamiltonian 66 and the DKH-DEF2-TZVP basis set with SARC/J auxiliary basis. The results remained consistent with previous discussions, showing no significant variations in the excited-state processes.