The perturbative contribution to the correlator of two HQET heavy-light currents expanded in light-quark masses up to quadratic terms is calculated up to 4 loops. The leading large-${\beta }_0$ limit is also considered, so that terms with the highest degrees of $n_f$ are calculated to all orders in ${\alpha }_s$. Borel images of coefficient functions in this limit contain renormalon poles. Naive nonabelianization works surprisingly poorly for the coefficient functions considered here.