As 2D IR spectroscopy is commonly applied to probe anharmonic vibrations of coupled local oscillators, methods of anharmonic theoretical vibrational spectroscopy that are based on localized vibrational modes ideally suited for computational 2D IR spectroscopy. Here, we have demonstrated that the L-VSCF/L-VCI method can be used for predicting 2D IR spectra.