B-B mixing in the HQET

We present a high statistics, quenched lattice calculation of the B-parametersB_Bd and B_Bs, 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 β = 6.0 on a 24³ × 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 isB_Bd = 1.08 +/- 0.06 +/- 0.08. B_Bd 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 +/-0.08 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 obtainB_Bs /B_Bd = 1.01 +/- 0.01 andf²_BsB_Bs /f²_BdB_Bd = 1.38 +/- 0.07.