B-B̄ mixing in the HQET

We present a high statistics, quenched lattice calculation of the B-parameters BBd and BBs, 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 243 × 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 B̂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 obtain B̂Bs/B̂Bd = 1.01 ± 0.01 and f2BsB̂Bs/f2BdB̂Bd = 1.38 ± 0.07. © 1997 Elsevier Science B.V.