  We show that the icosahedral packings of protein capsomeres proposed by
Caspar and Klug for spherical viruses become unstable to faceting for
sufficiently large virus size, in analogy with the buckling instability of
disclinations in two-dimensional crystals. Our model, based on the nonlinear
physics of thin elastic shells, produces excellent one parameter fits in real
space to the full three-dimensional shape of large spherical viruses. The
faceted shape depends only on the dimensionless Foppl-von Karman number
\gamma=YR^2/\kappa, where Y is the two-dimensional Young's modulus of the
protein shell, \kappa is its bending rigidity and R is the mean virus radius.
The shape can be parameterized more quantitatively in terms of a spherical
harmonic expansion. We also investigate elastic shell theory for extremely
large \gamma, 10^3 < \gamma < 10^8, and find results applicable to icosahedral
shapes of large vesicles studied with freeze fracture and electron microscopy.
