Classical continuum theories such as elasticity, electrostatics, and piezoelectricity are the long-wave and low-frequency limits of the equations of lattice dynamics. The partial differential equations of these continuum theories can be obtained from the finite difference equations of lattice dynamics by Taylor expansions and truncations. The equations of classical continuum theories are accurate for phenomena with a characteristic length much larger than microscopic characteristic lengths, for example, the distance between neighboring atoms in a lattice. When the characteristic length of a problem is not much larger than the microscopic characteristic length, classical continuum theories do not predict results consistent with lattice dynamics, and hence are no longer valid. For example, lattice waves are dispersive but the theory of elasticity only predicts nondispersive plane waves which are the long wave limit of lattice waves.There are different ways to modify the classical continuum theories so that their range of applicability can be extended to problems with smaller characteristic lengths, with results closer to lattice dynamics in a wider range of wave lengths.One is to generalize the constitutive relations in classical continuum theories from local to nonlocal [1][2][3][4][5]. In classical continuum theories, the constitutive relations are for stress and strain (or electric field and polarization) at the same material point, with a zero interaction distance (local theories). These constitutive relations neglect the microstructure of real materials and lead to macroscopic theories for large characteristic lengths. In nonlocal constitutive relations, a finite microscopic interaction distance is