I want to calculate the spin orbit coupling matrix elements (SOCME) between first singlet excited state (S1) and triplet excited states (T1, T2, T3, ....). I have done the TDDFT calculation using Gaussian16 using td (50:50) command. Now, I want to determine the SOCME from the output file. Is there any program/script to do so?