General formalism of absorption and emission spectra, and of radiative and nonradiative decay rates are derived using a thermal vibration correlation function formalism for the transition between two adiabatic electronic states in polyatomic molecules. Displacements, distortions, and Duschinsky rotation of potential energy surfaces are included within the framework of a multidimensional harmonic oscillator model. The Herzberg-Teller (HT) effect is also taken into account. This formalism gives a reliable description of the Q(x) spectral band of free-base porphyrin with weakly electric dipole-allowed transitions. For the strongly dipole-allowed transitions, e.g., S-1 -> S-0 and S-0 -> S-1 of linear polyacenes, anthracene, tetracene, and pentacene, the HT effect is found to enhance the radiative decay rates by -> 10% compared to those without the HT effect. For nonradiative transition processes, a general formalism is presented to extend the application scope of the internal conversion theory by going beyond the promoting-mode approximation. Numerical calculations for the nonradiative S-1 -> S-0 decay rate of azulene well explain the origin of the violation of Kasha's rule. When coupled with first-principles density functional theory (DFT) calculations, the present approach appears to be an effective tool to obtain a quantitative description and detailed understanding of spectra and photophysical processes in polyatomic molecules.