Mid-infrared spectra of an O-H center dot center dot center dot pi hydrogen-bonded 1:1 complex between formic acid and benzene were measured by isolating the complex in an argon matrix at a temperature of 8 K. The O-H stretching fundamental of formic acid (nu(O-H)) undergoes a red shift of 120 cm(-1), which is the largest among the known pi-hydrogen bonded complexes of an O-H donor with respect to benzene as acceptor. Electronic structure theory methods were used extensively to suggest a suitable geometry of the complex that is consistent with a recent study performed at CCSD(T)/CBS level by Zhao et al. (J Chem. Theory Comput. 2009, 5, 2726-2733), as well as with the measured IR spectral shifts of the present study. It has been determined that density functional theory (DFT) D functionals as well as parametrized DFT functionals like M06-2X, in conjunction with modestly sized basis sets like 6-31G (d, p), are sufficient for correct predictions of the spectral shifts observed in our measurement and also for reproducing the value of the binding energy reported by Zhao et al. We also verified that these low-cost methods are sufficient in predicting the nu(O-H) spectral shifts of an analogous O-H center dot center dot center dot pi hydrogen-bonded complex between phenol and benzene. However, some inconsistencies with respect to shifts of nu(O-H) arise when diffuse functions are included in the basis sets, and the origin of this anomaly is shown to lie in the predicted geometry of the complex. Natural bond orbital (NBO) and atoms-in-molecule (AIM) analyses were performed to correlate the spectral behavior of the complex with its geometric parameters.