The structure, conformational behavior, and spectroscopic parameters of bibenzyl have been investigated by a computational protocol including proper treatment of anharmonic and hindered rotor contributions. Conventional hybrid functionals overstabilize the anti conformer while low-order post-Hartree-Fock (MP2) approaches strongly favor the gauche conformer. However, inclusion of semiempirical dispersion effects in density functionals or coupled cluster post-Hartree-Fock models agree in forecasting the simultaneous presence of both conformers in the gas phase with a slightly larger stability (0.7 kcal·mol(-1)) of the gauche conformer. Addition of thermal and entropic effects finally leads to very close Gibbs free energies for both conformers and, thus, to a slight preference for the gauche form due to statistical factors (2 vs 1). The situation remains essentially the same in solution. On these grounds, perturbative vibrational computations including both electrical and mechanical anharmonicities lead to IR and Raman spectra in remarkable agreement with experiment. Full assignment of the IR spectra explains the presence of peaks from gauche or anti conformers. Comparison between computed and experimental Raman spectra confirms that both conformers are present in liquid phase, whereas the anti conformer seems to be preponderant in the solid state. Also computed NMR parameters are in good agreement with experiment.