GaAs/AlGaAs quantum wells doped with donor atoms are investigated for nonlinear optical applications in the THz range. The electronic properties of the quantum wells are obtained numerically by applying an iterative shooting algorithm. Donor binding energies are computed through the evaluation of variational wavefunctions. The solution of the density matrix equations of motion for non-interacting two-level atoms within the rotating wave approximation is used to formulate the third-order optical nonlinearities. Transitions between the 1s and 2p+ impurity states because of an incident light polarized perpendicularly to the growth direction are considered as the origins of optical nonlinearity. Following a set of computations for a quantum well doped at the center, it is found that the nonlinear susceptibility decreases when the well becomes wider or the Al concentration increases. Additionally, when the doping center is shifted to the well edge, the nonlinear susceptibility decreases too. A large nonlinear figure of merit is obtainable in wider wells compared to the narrower wells although the latter delivers larger nonlinear susceptibilities.