The rigorous and transparent treatment of the effects of nuclear vibrational motion in two-photon absorption (TPA) was discussed. Perturbation formula for diatomic molecules were developed and applied to the X¹Σ+–A¹Π transition in CO. The analysis showed that the vibrations played an important role in TPA, just as their role in the calculation of conventional nonlinear optical (NLO) hyperpolarizabilities