For both pathways, energy profiles computed using SF-TDDFT with the BH&HLYP functional (called "SF-BH&HLYP" hereafter) agree well with the MRCISD+Q results, except that SF-BH&HLYP predicts a larger energy gap in the CX regions, due to the fact that the intersection seams obtained by these two methods are located in somewhat different regions of coordinate space. This is not surprising as MRCISD+Q fails to give the correct topology of the CX, because the Davidson correction (+Q) modifies the energy but not the wave function. In the present study, the ωB97X functional is used for the LR-TDDFT calculations and we abbreviate this as LR-ωB97X. Previously, Casida and co-workers compared the performance of different density functionals for the out-of-plane relaxations of the retinal protonated Schiff base models PSB4 and PSB5. Taking the CASPT2 results as the reference, they found that a particular longrange corrected (LRC) functional, LRC-µBLYP, outperformed other functionals including ωB97X. For the PSB3 reaction pathways considered here, however, we find that the performance of ωB97X and LRC-µBLYP is almost identical. This lends some credence to the idea that the LR-TDDFT results reported here are unlikely to change qualitatively if other functionals commonly used for LR-TDDFT are substituted in place of ωB97X.