New reduced three-dimensional (3D) diabatic potential energy surfaces (PESs) involving the ππ,ππ*, and πσ* states for the nonadiabatic photodissociation CHSH(S) + hv → CHSH(ππ*/πσ*) → H + CHS•(A/X) were constructed at a high computational level, namely explicitly correlated multireference configuration interaction (MRCI-F12) method with the cc-pVTZ-F12 basis. The diabatization of the PESs was achieved by a simple, efficient, and reliable "regularized diabatization" method [Köppel, H.; Gronki, J.; Mahapatra, S. J. Chem. Phys. 2001, 115, 2377-2388]. The dissociation energy of the S state and the excitation energies of the excited S and S states were found to be in reasonably good agreement with the experimental values. The vibronic energy levels of the thiophenol (PhSH) and deuterated thiophenol (PhSD) for S and S states were calculated using a three-dimensional model, and they are in reasonably good agreement with the available experimental results, which validate the high accuracy of the adiabatic PESs and the reasonability of the diabatic couplings.