A linear variational recipe is used to compute the second non-linear optical response (gamma) of single carrier 2D quantum dots and determine conditions required for obtaining high values of (gamma). The effects of anharmonicities in the confinement potential are analyzed. (c) 2005 Elsevier B.V. All rights reserved.