In an attempt to understand the phase behavior of aqueous hydrogen fluoride, the clustering in the mixture is investigated at the molecular level. The study is performed at the mPW1B95/6-31+G(d,p) level of theory. Several previous studies attempted to describe the dissociation of HF in water, but in this investigation, the focus is only on the association patterns that are present in this binary mixture. A total of 214 optimized geometries of (HF)(n)(H2O)(m) clusters, with m + n as high as 8, were investigated. For each cluster combination, several different conformations are investigated, and the preferred conformations are presented. Using multiple linear regressions, the average strengths of the four possible H-bonding interactions are obtained. The strongest H-bond interaction is reported to be the H2O center dot center dot center dot H-F interaction. The most probable distributions of mixed clusters as a function of composition are also deduced. It is found that the larger (HF)(n)(H2O)(m) clusters are favored both energetically and entropically compared to the ones that are of size m + n <= 3. Also, the clusters with equimolar contributions of HF and H2O are found to have the strongest interactions.