Density Functional Theory calculations, in the ground state of 2-Phenylbenzofuran, were carried out using the GGA-PBE, PBV86 and meta-GGA-TPSS hybrid functionals with 6-31G (d, p) as a basis set. First, theoretical calculations were performed using these functionals to obtain the stable conformer of the molecule. In addition, Mulliken population natural population and natural bond orbital analyses were calculated. The molecular electrostatic potential, band gap energies, global, local chemical reactivity descriptors and non-linear optical (NLO) properties were studied. Additionally, the NLO properties of 2-Phenylbenzofuran and those of its derivatives were investigated by GGA-PBE/6-31G (d,p) level of theory. The first-order hyperpolarizability value of all 2-Phenylbenzofuran derivatives was found within the range extending from 4.00 × 10-30 to 43.57 × 10-30 (esu). It indicated that they possess remarkable NLO properties. In addition, a multiple linear regression procedure was used to envisage the relationships between molecular descriptors and the activity of 2-Phenylbenzofuran derivatives; the quantitative structure-activity relationship (QSAR) studies were performed on them using quantum descriptors. The QSAR was applied to determine a correlation between the various physico-chemical parameters of the studied compounds and their biological activities. The statistical quality of the QSAR models was assessed using statistical parameters, i.e. R2, R2adj and R2cv.