To construct two-dimensional infrared (2D IR) spectra having all vibrational modes of a molecule included is still quite challenging, both experimentally and theoretically. Here we report an ab initio-based all-mode 2D IR spectra simulation approach. Using deuterated glycolaldehyde (CH2OHCDO), the smallest sugar as a model molecule, we have calculated correlation 2D IR spectrum of its entire 3N - 6 (N = 8) normal modes in the mid-to-far-IR region (4000-0 cm(-1)), using quantum chemical anharmonic frequency and anharmonicity computations in conjunction with time-domain third-order nonlinear response functions. The calculated 2D IR spectra were found to contain a network of structural and dynamical parameters of the molecule. It is found that certain spectral regions, once enlarged, show features that are in reasonable agreement with limited but already available single-and dual-frequency 2D IR experimental results. The extension of narrow-band 2D IR spectroscopy into the full mid-to-far-IR regime would allow us to characterize the structural distributions and dynamics of molecular complexes in condensed phases with sufficient number of parameters.