$B - \bar B$ Mixing in the HQET
  We present a high statistics, quenched lattice calculation of the B-parameters $B_{B_d}$ and $B_{B_s}$, computed at lowest order in the HQET. The results were obtained using a sample of 600 quenched gauge field configurations, generated by Monte Carlo simulation at $\beta=6.0$ on a $24^{3}\times 40$ lattice. For the light quarks the SW-Clover action was used; the propagator of the lattice HQET was also tree-level improved. Our best estimate of the renormalization scale independent B-parameter is $\hat{B}_{B_d} = 1.03 \pm 0.06 \pm 0.18$. $\hat{B}_{B_d}$ has been obtained by using ``boosted'' perturbation theory to calculate the renormalization constants which relate the matrix elements of the lattice operators to the corresponding amplitudes in the continuum. Due to the large statistics, the errors in the extraction of the matrix elements of the relevant bare operators are rather small. The main systematic error, corresponding to $\pm 0.18$ in the above result, comes from the uncertainty in the evaluation of the renormalization constants, for which the one-loop corrections are rather large. The non-perturbative evaluation of these constants will help to reduce the final error. We also obtain $\hat{B}_{B_s}/\hat{B}_{B_d} = 1.01 \pm 0.01$ and $f^2_{B_s}\hat{B}_{B_s}/f^2_{B_d}\hat{B}_{B_d} = 1.38 \pm 0.07$. 
V. Gimenez and G. Martinelli