The minimum energy path on the S 1 potential energy surface of azobenzene connects the trans (and cis) geometries to the S 1 /S 0 conical intersection, located at approximately 95 • of torsion of the CNNC dihedral. Consequently, nonadiabatic transitions to S 0 in isolated trans-azobenzene occur predominantly in this region or slightly before, within a CNNC range of 90 to 120 • . During the decay to S 0 , this molecule can isomerize to cis-azobenzene. The partial 238 photoisomerization quantum yields Φ nπ * calculated using all four CTTSH methods, TSH, and TSH-ODC for the trajectories starting in the nπ * state are shown in the second column of Table . The calculated Φ nπ * values are consistent across methods, are in agreement with each other and in good agreement with experimental results, which fall within the range of 20-32% . It is worth noting that trajectories with faster decay predominantly remain in the trans configuration after relaxing to the ground state, whereas slower trajectories have a higher probability of undergoing photoisomerization. In the third column method Φ nπ * Φ ππ * CTTSH-FGAC 0.22 ± 0.04 0.11 ± 0.04 CTTSH-FGLA 0.19 ± 0.03 0.19 ± 0.05 CTTSH-TGAC 0.23 ± 0.04 0.13 ± 0.04 CTTSH-TGLA 0.20 ± 0.04 0.38 ± 0.08 TSH 0.18 ± 0.04 0.17 ± 0.05 TSH-0DC 0.22 ± 0.03 0.13 ± 0.04 TABLE II. Partial photoisomerization quantum yields (Φ) calculated with our four implementations of CTTSH, with TSH and TSH-ODC. The second column reports the results for the dynamics initiated in the nπ * state and the third column reports the results for the ππ * excitation.