Three-photon absorption probabilities delta(3PA) have been calculated through application of a recently derived method for cubic response functions within density functional theory (DFT). Calculations are compared with Hartree-Fock (HF) and with a coupled cluster hierarchy of models in a benchmarking procedure. Except for cases having intermediate states near resonance, density functional theory is demonstrated to be in sufficient agreement with the highly correlated methods in order to qualify for predictions of delta(3PA). For the larger systems addressed, a set of acceptor A and donor D substituted pi-conjugated systems formed by trans-stilbene and dithienothiophene (DTT), we find noticeable differences in the magnitude of delta(3PA) between HF and DFT, although similar trends are followed. It is shown that the dipolar structures, TS-AD and DTT-AD, have substantially larger delta(3PA) than other types of modifications which is in accordance with observations for two-photon absorption. This is the first application of density functional theory to three-photon absorption beyond the use of few-state models. (C) 2004 American Institute of Physics.