Geometry and interaction energy in complexes of the Ph-L type (L ϭAr, N 2 , CO, H 2 O, NH 3 , CH 4 , CH 3 OH, CH 3 F) involving neutral or cationic phenol were determined using the density functional theory formalism based on the minimization of the total energy bifunctional and gradient-dependent approximations for its exchangecorrelation and nonadditive kinetic-energy parts. For the neutral complexes the calculated interaction energies range from 1 kcal/mol for the Ph-Ar complex to about 10 kcal/mol for Ph-NH 3 . The interactions are stronger if the cationic phenol is involved (up to 25 kcal/mol). It is found, except for neutral Ph-Ar, that the hydrogen-bonded structure is more stable than the -bound one. Calculated interaction energies (D e ) correlate well with the experimental dissociation energies (D 0 ).