Analytic expressions of corrections for the breakdown of the Born−Oppenheimer approximation to Dunham's Y ij with optimal parameters, i.e., determinable clusters of expansion coefficients, are applied to a data analysis of the rotational and vibrational-rotational transitions of HF reported in the literature. All the available spectral lines of the two isotopologues, DF and HF, are simultaneously fitted to a single set of molecular parameters of HF within experimental errors. Fitting of a data set of 595 spectral transitions for DF and HF has generated only 20 minimal independent parameter values, i.e., "irreducible" molecular constants of HF, that are sufficient to precisely generate 82 Y ij coefficients and 144 band constants in total: 41 Y ij and 72 band constants each for DF and HF.