# \<math display\="inline"\>\<msub\>\<mi\>B\</mi\>\<mi\>s\</mi\>\</msub\>\<mo stretchy\="false"\>→\</mo\>\<msup\>\<mi\>μ\</mi\>\<mo\>\+\</mo\>\</msup\>\<msup\>\<mi\>μ\</mi\>\<mo\>\-\</mo\>\</msup\>\<mi\>γ\</mi\>\</math\> decay rate at large \<math display\="inline"\>\<msup\>\<mi\>q\</mi\>\<mn\>2\</mn\>\</msup\>\</math\> from lattice QCD

**INSPIRE:** [2755393](https://inspirehep.net/literature/2755393)
**arXiv:** [2402\.03262](https://arxiv.org/abs/2402.03262)
**DOI:** [10\.1103/PhysRevD\.109\.114506](https://doi.org/10.1103/PhysRevD.109.114506)

**Authors:** Frezzotti, R\., Tantalo, N\., Gagliardi, G\., Sanfilippo, F\., Simula, S\., Lubicz, V\., Martinelli, G\., Sachrajda, C\.T\.

**Submitted:** 5 February 2024

**Subjects:**
- hep\-lat
- hep\-ph
- Lattice

**Journal reference:** Phys\.Rev\.D 109 \(11\) 114506 \(2024\)

**Citations:** 25

## Abstract

We determine, by means of lattice QCD calculations, the local form factors describing the \<math display\="inline"\>\<msub\>\<mi\>B\</mi\>\<mi\>s\</mi\>\</msub\>\<mo stretchy\="false"\>→\</mo\>\<msup\>\<mi\>μ\</mi\>\<mo\>\+\</mo\>\</msup\>\<msup\>\<mi\>μ\</mi\>\<mo\>\-\</mo\>\</msup\>\<mi\>γ\</mi\>\</math\> decay, in the so\-called electroquenched approximation\. For this analysis we make use of the gauge configurations produced by the ETM Collaboration with \<math display\="inline"\>\<msub\>\<mi\>N\</mi\>\<mi\>f\</mi\>\</msub\>\<mo\>\=\</mo\>\<mn\>2\</mn\>\<mo\>\+\</mo\>\<mn\>1\</mn\>\<mo\>\+\</mo\>\<mn\>1\</mn\>\</math\> flavor of Wilson\-Clover twisted\-mass fermions at maximal twist\. To obtain the \<math display\="inline"\>\<msub\>\<mi\>B\</mi\>\<mi\>s\</mi\>\</msub\>\</math\> meson form factors, we perform simulations for several heavy\-strange meson masses \<math display\="inline"\>\<msub\>\<mi\>m\</mi\>\<msub\>\<mi\>H\</mi\>\<mi\>s\</mi\>\</msub\>\</msub\>\</math\> in the range \<math display\="inline"\>\<msub\>\<mi\>m\</mi\>\<msub\>\<mi\>H\</mi\>\<mi\>s\</mi\>\</msub\>\</msub\>\<mo\>∈\</mo\>\<mo stretchy\="false"\>\[\</mo\>\<msub\>\<mi\>m\</mi\>\<msub\>\<mi\>D\</mi\>\<mi\>s\</mi\>\</msub\>\</msub\>\<mo\>,\</mo\>\<mn\>2\</mn\>\<msub\>\<mi\>m\</mi\>\<msub\>\<mi\>D\</mi\>\<mi\>s\</mi\>\</msub\>\</msub\>\<mo stretchy\="false"\>\]\</mo\>\</math\>, and extrapolate to the physical \<math display\="inline"\>\<msub\>\<mi\>B\</mi\>\<mi\>s\</mi\>\</msub\>\</math\> meson point \<math display\="inline"\>\<msub\>\<mi\>m\</mi\>\<msub\>\<mi\>B\</mi\>\<mi\>s\</mi\>\</msub\>\</msub\>\<mo\>≃\</mo\>\<mn\>5\.367\</mn\>\<mtext\> \</mtext\>\<mtext\> \</mtext\>\<mi\>GeV\</mi\>\</math\> making use of the HQET scaling laws\. We cover the region of large dimuon invariant masses \<math display\="inline"\>\<msqrt\>\<msup\>\<mi\>q\</mi\>\<mn\>2\</mn\>\</msup\>\</msqrt\>\<mo\>&gt;\</mo\>\<mn\>4\.16\</mn\>\<mtext\> \</mtext\>\<mtext\> \</mtext\>\<mi\>GeV\</mi\>\</math\>, and use our results to determine the branching fraction for \<math display\="inline"\>\<msub\>\<mi\>B\</mi\>\<mi\>s\</mi\>\</msub\>\<mo stretchy\="false"\>→\</mo\>\<msup\>\<mi\>μ\</mi\>\<mo\>\+\</mo\>\</msup\>\<msup\>\<mi\>μ\</mi\>\<mo\>\-\</mo\>\</msup\>\<mi\>γ\</mi\>\</math\>, which has been recently measured by LHCb in the region \<math display\="inline"\>\<msqrt\>\<msup\>\<mi\>q\</mi\>\<mn\>2\</mn\>\</msup\>\</msqrt\>\<mo\>&gt;\</mo\>\<mn\>4\.9\</mn\>\<mtext\> \</mtext\>\<mtext\> \</mtext\>\<mi\>GeV\</mi\>\</math\>\. The largest contribution to the uncertainty in the partial branching fractions at values of \<math display\="inline"\>\<msqrt\>\<msup\>\<mi\>q\</mi\>\<mn\>2\</mn\>\</msup\>\</msqrt\>\<mo\>&lt;\</mo\>\<mn\>4\.8\</mn\>\<mtext\> \</mtext\>\<mtext\> \</mtext\>\<mi\>GeV\</mi\>\</math\> is now due to resonance and other long\-distance effects, including those from “charming penguins,” which we estimate by summing over the contributions from the \<math display\="inline"\>\<msup\>\<mi\>J\</mi\>\<mi\>P\</mi\>\</msup\>\<mo\>\=\</mo\>\<msup\>\<mn\>1\</mn\>\<mo\>\-\</mo\>\</msup\>\</math\> charmonium resonances\.