Using three typical pi-conjugated molecules (1,3,5-hexatriene, 1-formyl-6-hydroxyhexa-1,3,5- triene, and 1,1-diamino-6,6-dinitrohexa-1,3,5-triene) we investigate the level of ab initio theory necessary to produce reliable values for linear and nonlinear optical properties, with emphasis on the vibrational contributions that are known to be important or potentially important. These calculations are made feasible by employing field-induced coordinates in combination with a finite field procedure. For many, but not all, purposes the MP2/6-31+G(d) level is adequate. Based on our results the convergence of the usual perturbation treatment for vibrational anharmonicity was examined. Although this treatment is initially convergent in most circumstances, a problematic situation has been identified.