A set of 13 bis(difluoroboron)-1,2-bis((pyrrol-2-yl)methylene)hydrazine (BOPHY) dyes is studied through a hybrid time-dependent density functional theory(TD-DFT)-scaled opposite spin-configuration interaction singles with a double correction [SOS-CIS(D)] approach accounting for solvent effects, to shed light onto the structure property relationships of these recently developed chromophores. In the first step, we calculate the absorption fluorescence crossing points with refined TD-DFT models considering the influences of both vibrational and solvent contributions. We found that the systematic overestimation of the 0-0 energies is effectively reduced by combining polarizable continuum model-TD-DFT with a scaled opposite spin-configuration interaction singles with a double correction [SOS-CIS(D)]. Next, for a representative system, the vibrationally resolved spectrum within the harmonic approximation is computed on the basis of TD-DFT vibrational signatures and an excellent match with experiment is found. Finally, the influence of different lateral groups on the spectroscopic properties is rationalized by investigating charge transfer density difference maps. It is found that one can tune the position of the absorption/emission maxima by a judicious choice of the lateral substituents or by using pi-extended segments. The largest absorption and emission wavelengths as well as the largest Stokes shifts are obtained for BOPHYs containing strong electron-donor dimethylaminophenyl groups attached to the a positions of the pyrrole units through vinyl linkers, making these chromophores promising candidates for bioluminescence applications.