  We optimize the first and second intrinsic hyperpolarizabilities for a 1D
piecewise linear potential dressed with Dirac delta functions for $N$
non-interacting electrons. The optimized values fall rapidly for $N>1$, but
approach constant values of $\beta_{int}=0.40$, $\gamma_{int}^{+}=0.16$ and
$\gamma_{int}^{-}=-0.061$ above $N\gtrsim8$. These apparent bounds are achieved
with only 2 parameters with more general potentials achieving no better value.
In contrast to previous studies, analysis of the hessian matrices of
$\beta_{int}$ and $\gamma_{int}$ taken with respect to these parameters shows
that the eigenvectors are well aligned with the basis vectors of the parameter
space, indicating that the parametrization was well-chosen. The physical
significance of the important parameters is also discussed.
