We use the bottom quark mass \(m_Q = 4.977 \, \text{GeV}\)  from quark models and the spin-averaged masses \(m_B = 5.313 \, \text{GeV}\)  and \(m_{B_s} = 5.403 \, \text{GeV}\)  from experiments. \(E_{\text{threshold}} = 10.790 \, \text{GeV}\)  in the Schrödinger equation (REF ) is closer to \(2 m_{B_s}\)  than to \(2 m_B\)  and reflects that the lattice QCD results from Ref. [1]} were obtained with a light \(u/d\)  quark mass rather close to the physical mass of the \(s\)  quark.
We propagate the uncertainties of the lattice potentials by resampling, i.e. we generate 1000 statistically independent samples and repeat all computations on each the 1000 samples. Statistical errors are defined via the 16th and 84th percentile.
<FIGURE>