The dash-dotted curve in Figure 3b is a fit to |σ shg zzz | based on the two-photon term in equation (3). Given the simplicity of the model relative to the complexity of the full 3D bandstructure of TaAs, the 3D version of the RM model describes the data remarkably well. The peak position, low, and high energy tails of the spectrum are best described using parameterst = 1.5,δ = 1.4,t AB = 0.02, and ∆ = 0.428. The lattice parameters for TaAs are c = 1.165 nm, and b = a = 0.344 nm.The theory discussed above helps to identify the factors that contribute to a large SHG response. The first is a polar crystal structure such as 4mm (TaAs, NbAs, etc.) or mm2 (single layer monochalcogenides) consistent with a phenomenological description in terms of coupled ferroelectric chains. However, even when this description is possible there is a global bound on the frequency integrated second-order response, Σ shift z , of a system described by H RM that is given by,where G is also a dimensionless function of the dimensionless parameters, with a global maximum 0.604 [17]. Given this bound on Σ shift z , it is clear that the single most important factor in optimizing the second-order response of a 3D crystal is the effective aspect ratio parameter c 2 /ab.We conclude by addressing the question of ultimate bounds on Σ shift z , which is relevant not only to SHG, but to sensitivity of photodetectors and efficiency of solar cells based on Berry curvature generated intrinsic photogalvanic effects [21,28] as well. We have discovered a new connection between Σ shift z and the modern theory of polarization that generalizes the formulae for spectral weight beyond the RM Hamiltonian. Specifically we show that Σ shift z is unbounded and potentially divergent when the possibility of next-neighbor hopping is included.As demonstrated in the Supplementary Information, Σ shift z in two band systems, regardless of dimensionality d and range of electron hopping amplitude is given by,