Valence electron contributions to the static molecular third-order polarizabilities (γ) are calculated for C 60 , C 70 , and two stable structures of C 84 (D 2 and D 2d ). The method utilized is based on the finite-field approach coupled with semiempirical polarization calculations on all-valence electrons. An increase in the third-order polarizability contributions is observed for molecular structures with a reduction in group symmetry, in agreement with recent experimental observations for these fullerenes. This increase is attributed mainly to the appearance of aromatic structures within the molecules as well as to the increase in molecular volume.Fullerenes are the subject of intensive research by many groups on several contents. There has also been a great deal of speculation about potential uses for fullerenes and their derivatives, ranging from electrochemistry and electronics (the latter based in part on the superconducting properties of alkalidoped fullerides) to catalysis, hydrogen storage, antiviral agents, and all-optical devices. 1 Associated with the delocalized π-conjugated electrons in fullerenes, their large nonlinear optical responses have been the subject of several recent experimental studies. 1-10 A wealth of experimental data concerning the thirdorder optical nonlinearities have been accumulated for the fullerenes C 60 , C 70 , and C 84 . The two most abundant fullerenes, C 60 and C 70 , have configurations with I h and D 5h symmetries, respectively. 11-14 C 84 has been recently synthesized, 15,16 and the most stable configurations have been determined to be D 2 and D 2d . 1,17 An interesting feature discovered by the experimental study is the increase of the third-order nonlinearity from C 60 to C 70 and to C 84 . 1 With the advances in the synthesis of macroscopic quantities of fullerenes, a theoretical understanding of the dependence of the nonlinear optical responses on the molecular symmetry and electronic properties is thus of particular interest.Theoretical calculations of third-order nonlinear effects have been carried out using various semiempirical and ab initio methods. π-electron calculations on unsaturated hydrocarbon chains have been used to predict third-order effects. 19 Ab initio methods applied to molecules with more than 15 atoms are extremely computer intensive. An alternative method using allvalence electrons and a semiempirical perturbation (sum-overstates approach) has generally been applied to molecules with less than 30 atoms. 20,21 In the present work, the procedure X