<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
arXiv: 2402.03262
DOI: 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.