Quartz tuning forks (QTFs), which were coated with gold and with self-assembled monolayers (SAM) of a lower-rim functionalized calix[4]arene methoxy ester (CME), were used for the detection of benzene, toluene, and ethylbenzene in water samples. The QTF device was tested by measuring the respective frequency shifts obtained using small (100 µL) samples of aqueous benzene, toluene, and ethylbenzene at four different concentrations (10−12, 10−10, 10−8, and 10−6 M). The QTFs had lower limits of detection for all three aromatic hydrocarbons in the 10−14 M range, with the highest resonance frequency shifts (±5%) being shown for the corresponding 10−6 M solutions in the following order: benzene (199 Hz) > toluene (191 Hz) > ethylbenzene (149 Hz). The frequency shifts measured with the QTFs relative to that in deionized water were inversely proportional to the concentration/mass of the analytes. Insights into the effects of the alkyl groups of the aromatic hydrocarbons on the electronic interaction energies for their hypothetical 1:1 supramolecular host–guest binding with the CME sensing layer were obtained through density functional theory (DFT) calculations of the electronic interaction energies (ΔIEs) using B3LYP-D3/GenECP with a mixed basis set: LANL2DZ and 6-311++g(d,p), CAM-B3LYP/LANL2DZ, and PBE/LANL2DZ. The magnitudes of the ΔIEs were in the following order: [Au4-CME⊃[benzene] > [Au4-CME]⊃[toluene] > [Au4-CME]⊃[ethylbenzene]. The gas-phase BSSE-uncorrected ΔIE values for these complexes were higher, with values of −96.86, −87.80, and −79.33 kJ mol−1, respectively, and −86.39, −77.23, and −67.63 kJ mol−1, respectively, for the corresponding BSSE-corrected values using B3LYP-D3/GenECP with LANL2dZ and 6-311++g(d,p). The computational findings strongly support the experimental results, revealing the same trend in the ΔIEs for the proposed hypothetical binding modes between the tested analytes with the CME SAMs on the Au-QTF sensing surfaces.