In this manuscript, a relatively simple and inexpensive INDO/SCI finite-field (FF) method for calculating polarizabilities (R) is demonstrated to give good agreement with results obtained by both the INDO/MRD/ SDCI sum-over-states routine and published results using the RPA method. The FF method is as effective as the other techniques in predicting both ground and excited state R's in all substituted and unsubstituted polyenes studied. We observe the correlation described by Marder and co-workers between bond-order-alternation (BOA) and dipolar properties, such as the change in R between the ground and excited states (∆R). In addition, qualitative, but not quantitative, agreement is seen between the calculated ∆R's of polar polyenes and those measured by Stark-effect spectroscopy.Recent interest in substituted polyenes which exhibit nonlinear optical (NLO) behavior has made accurate calculations of their properties an important predictive tool in materials development. For example, the hyperpolarizabilities ( ) of molecules are of interest because this parameter is directly related to efficiency in optical frequency doubling and to the Pockel effect. Experimentally, there are two common strategies for optimizing : [1][2][3][4][5][6][7] (1) changing the donor (D) and/or acceptor (A) strength of substituents on the polyene and (2) changing the solvent polarity and/or polarizability. Both alter the local field of the polyene chain, leading to a change in bond-order alternation (BOA) 8 and in NLO properties, including . In the two-state model, 9,10 is proportional to the change in dipole moment between the ground and excited states (∆µ), so an accurate calculation of solvated ∆µ is needed to correctly model the effects of solvation on . 3,6,[11][12][13][14][15][16] In dielectric cavity models, 17 the solvent-corrected ∆µ also depends on the change in polarizability between the ground and excited states (∆R). Therefore, quantitatively accurate calculations of both ground and excited state R are needed to correctly model the solvent-effects on and ∆µ, particularly in highly polarizable systems such as the substituted polyenes. This paper compares the two most common routines for calculating ground state R, the sum-over-states and finite-field methods, and applies them to the calculation of excited state R.While techniques for the calculations of first ( ) and second (γ) order hyperpolarizabilities and ground-state polarizabilities (R g ) of polyenes are well developed, 12,14,15,[18][19][20][21][22][23] we are aware of only one example of excited-state polarizability (R e ) calculations for such systems in the literature. 24 In this published work, the random phase approximation (RPA) was utilized to investigate a series of linear, unsubstituted polyenes. Marder and coworkers have investigated trends in ∆R 25 and R g 1 with respect to BOA, but to our knowledge, no quantitative calculations of R e or ∆R have been published on substituted systems similar to those being developed for NLO applications. Because linear unsub...