Heterocyclic compounds have potential applications in many fields of life. We synthesized novel tetra substituted imidazoles by four-component condensation of benzil, substituted aldehydes, substituted anilines and ammonium acetate as a source of ammonia and acetic acid as the solvent. Their chemical structures were resolved through X-ray crystallographic and spectroscopic (Fourier transform IR and UV–vis) techniques. In addition to experimental analysis, density functional theory (DFT) calculations at the B3LYP/6-311 + G(d,p) level were performed on 4-bromo-2-(1-(4-methoxyphenyl)-4,5-diphenyl-1H-imidazole-2-yl)phenol (1), 4-bromo-2-(1-(1-naphthalen-yl)-4,5-diphenyl-1H-imidazole-2-yl)phenol (2), and 2-(1-(2-chlorophenyl)-4,5-diphenyl-1-H-imidazole-2-yl)-6-methoxyphenol (3) to obtain the optimized geometry and spectroscopic (Fourier transform IR and UV–vis) and non-linear optical properties. Frontier molecular orbital analysis was performed at the Hartee-Fock/6-311+g(d,p) and DFT/B3LYP/6-311+G(d,p) levels of theory. Natural bond orbital (NBO) and UV–vis spectral analyses were performed at the M06-2X/6-31+G(d,p) and time-dependent DFT/B3LYP/6-311+G(d,p) levels, respectively. Overall, the DFT findings show good agreement with the experimental data. The hyper conjugative interaction network, which is responsible for the stability of compounds 1, 2 and 3 was explored by the NBO approach. The global reactivity parameters were explored with use of the energy of the frontier molecular orbitals. DFT calculations predict the first-order hyperpolarizabilities of compounds 1, 2 and 3 are 294.89 × 10−30, 219.45 × 10−30 and 146.77 × 10−30 esu, respectively. A two-state model was used to describe the non-linear optical properties of the compounds investigated.