In this article we report results for the electronic and vibrational hyperpolarizabilities of ten molecules: Lipbenzene, Lip pyridine, Li@pyrimidine, and Lippyrazine; Lipnaphthalene, Li@quinoline, Li@isoquinoline, Li@cinnolin, Li@quinazoline, and Li@quinoxaline. An electron correlation study shows that second-order many-body perturbation theory and density functional theory with CAM-B3LYP and MOS-2X functionals give results for the electronic hyperpolarizabilities in good agreement with coupled cluster with single and doubles reference values. Static and dynamic vibrational corrections were computed through the perturbation theoretical method of Bishop and Kirtman and using a variational approach. In general, we obtained notable discrepancies between the results obtained by the two methods for the pure vibrational corrections because of the deficiency of the perturbation method to properly treat low-frequency normal modes present in the investigated systems. However, both methods give results similar to the zero-point vibrational average corrections.