Atomic-scale descriptions of linear-optical properties such as reflection are nearly a century old, but surprisingly, analogous models describing nonlinear-optical (NLO) properties as the natural dynamic response of bond charges driven by an external field are a recent development. These bond-charge models have proven to be particularly useful in describing the relevant physics of second-, third-, and fourth-harmonic generation, identifying previously unrecognized contributions to NLO responses, and uncovering well-disguised correlations in tensor parameters determined phenomenologically from symmetry conditions. Current capabilities are discussed, and opportunities for improved understanding noted. 1 Introduction The availability of continuum sources of radiation covering wide spectral ranges has not only made linear-optical spectroscopy possible, but also linear-optical techniques indispensable for determining information nondestructively about materials, interfaces, and structures. Owing to the higher-order tensors involved, in principle nonlinear-optical (NLO) probes are much more powerful, but spectroscopic data [1-5] have been difficult to acquire. However, commercial femtosecond systems that can tune over wide wavelength ranges are becoming available, and we are rapidly reaching the point where NLO spectroscopy will be routine.The potentially greater diagnostic power of NLO comes at a price, specifically the need to use more complex models and methods to interpret these more complex data. For most linear-optical applications, measurements can be analyzed using macroscopic parameters such as the dielectric function e, which for isotropic materials is a scalar and anisotropic materials a second-rank tensor. Here a return to fundamentals is generally unnecessary, and the last contact most linear-optical spectroscopists have had with atomic scale properties is the textbook derivation of the ClausiusMossotti relation [6], which relates the dielectric constant to local polarizabilities in the static limit.For NLO, quantities equivalent to e are the susceptibility tensors x ijk , x ijkl , etc. [7]. Symmetry-allowed NLO tensor components have been determined for all crystal classes, and