Electron-vibration coupling (or vibronic coupling) of polyatomic molecules in the condensed phase has received considerable attention experimentally and theoretically using linear spectroscopy and four-wave mixing techniques in an effort to probe the structure and dynamics of the system at hand. For this reason, a detailed study of vibronic coupling in harmonic polyatomic molecules featuring a greater degree of computational efficiency is presented. A full treatment of non-Condon systems whereupon linear and nonlinear Herzberg-Teller vibronic coupling and Duschinsky mixing effects in harmonic systems taking place upon electronic excitation is provided. The utilization of an exponential function to express the nuclear dependence of the electronic transition dipole moment, thereby avoiding the finite sum and eigenstate representation that is normally used in computing a non-Condon interaction, leads to a simpler electronic transition dipole moment time correlation function with rapid convergence and better numerical stability than previously reported works. A closed-form expression is obtained for the electronic transition dipole moment time correlation function of polyatomic molecules in which linear and nonlinear Herzberg-Teller vibronic coupling and Duschinsky mixing effects are accounted for. An important numerical observation regarding dealing with branch cuts, which manifest themselves as discontinuities in the function itself or its first derivative, and are often exhibited by complex-valued correlation functions, is noted and treated using the Riemann surface approach. The resultant dipole moment correlation function is in turn employed to calculate linear absorption and hole-burning signals, accounting for the aforementioned spectroscopic effects. A link between wavelets and the electronic transition dipole moment time correlation function in multidimensional harmonic systems is made in the concluding remarks.