Hydrogen bonding in liquids of the constitution isomers of heptan-1-ol mixed with 1-alkyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ionic liquids (ILs), [C_nmim][NTf₂], is investigated using both computational and experimental techniques. All-atom non-polarizable molecular-dynamics (MD) simulations predict that the hydrogen bonds gradually decay with increasing temperature. This decay is more pronounced for the branched alcohols and in the presence of the ionic liquids. The primary and linear isomer, heptan-1-ol, and its tertiary and bulky analogue 3-ethylpentan-3-ol are identified as the opposite extremes of the spectrum of hydrogen bonding stability in the bulk liquid. While neat heptan-1-ol exhibits strong hydrogen bonding at 350 K, 3-ethylpentan-3-ol is prone to hydrogen bonding decay already at 300 K. The presence of ionic liquids is found to affect the hydrogen bonding comparably as a 50 K temperature increase. Since the heat capacities of the associating liquids are very sensitive to any variation in hydrogen bonding strength and to the character of the hydrogen-bonded clusters in the bulk liquid, the calorimetric effort provides useful experimental data to confirm the results predicted by MD simulations. In this work, excess heat capacity is measured for equimolar single-phase mixtures of alcohols and ILs, and it differs largely in its sign and magnitude for individual heptanol isomers. Temperature trends of the excess heat capacities suggest that the stability of hydrogen bonding for individual heptanol isomers is temperature-shifted, based on their capability of hydrogen bonding. The predicted hierarchy of hydrogen bonding in individual alcohols and its impact on the excess heat capacity trends are qualitatively confirmed via thermodynamic modelling of the associative contribution to the excess heat capacities. These terms are found to predetermine the observed non-monotonous excess heat capacity trends.