Coupled-cluster theory with single and double excitations is applied to the calculation of optical properties of large polyaromatic hydrocarbons. Dipole polarizabilities are reported for benzene, pyrene, and the oligoacenes sequence n=2-6. Dynamic polarizabilities were calculated on polyacences as large as pentacene for a single frequency and for benzene and pyrene at many frequencies. The basis set effect was studied for benzene using a variety of basis sets in the Pople [Theor. Chim. Acta 28, 213 (1973)] and Dunning [J. Chem. Phys. 90, 1007 (1989)] families up to aug-cc-pVQZ and the Sadlej pVTZ basis [Collect. Czech. Chem. Commun. 53, 1995 (1998)], which was used exclusively for the largest molecules. Geometries were optimized using HF, B3LYP, PBE0, and MP2 and compared to experiment to measure method dependence and the possible role of bond-length alternation. Finally, the polarizability results were compared to four common density functionals (B3LYP, BLYP, PBE0, PBE).