We have measured the intermolecular dissociation energies D 0 of supersonically cooled 1-naphthol (1NpOH) complexes with the solvents S=furan, thiophene, 2,5-dimethylfuran and tetrahydrofuran. The naphthol OH forms non-classical H-bonds to the aromatic π-electrons of furan, thiophene and 2,5dimethylfuran, and a classical H-bond to the tetrahydrofuran O atom. Using the stimulated-emission pumping resonant two-photon ionization (SEP-R2PI) method, the ground-state D 0 (S 0 ) values were bracketed as 21.8 ± 0.3 kJ/mol for furan, 26.6 ± 0.6 kJ/mol for thiophene, 36.5 ± 2.3 kJ/mol for 2,5-dimethylfuran and 37.6 ± 1.3 kJ/mol for tetrahydrofuran. The dispersion-corrected density functional theory (DFT-D) methods B97-D3, B3LYP-D3 (using the def2-TZVPP basis set) and ωB97X-D [using the 6-311++G(d,p) basis set] predict that the H-bonded (Edge) isomers are more stable than the Face isomers bound by dispersion; experimentally, we only observe Edge isomers. We compare the calculated and experimental D 0 values and extend the comparison to the previously measured 1NpOH complexes with cyclopropane, benzene, water, alcohols and cyclic ethers. The dissociation energies of the nonclassically H-bonded complexes increase roughly linearly with the average polarizability of the solvent,ᾱ(S). In contrast, the D 0 values of the classically H-bonded complexes are larger, increase more rapidly at lowᾱ(S) but saturate for largeᾱ(S). The calculated D 0 (S 0 ) values for the cyclopropane, benzene, furan and tetrahydrofuran complexes agree with experiment to within 1 kJ/mol, those of thiophene and 2,5-dimethylfuran are ∼ 3 kJ/mol smaller than experiment. The B3LYP-D3 calculated D 0 values exhibit the lowest mean absolute deviation (MAD) relative to experiment (MAD=1.7 kJ/mol), the B97-D3 and ωB97X-D MADs are 2.2 and 2.6 kJ/mol, respectively.