Penguin Contributions to CP Phases in B(d,s) Decays to Charmonium

The precision of the CP phases 2β and 2β(s) determined from the mixing-induced CP asymmetries in B(d)→J/ψK(S) and B(s)→J/ψϕ, respectively, is limited by the unknown long-distance contribution of a penguin diagram involving up quarks. The penguin contribution is expected to be comparable in size to the precision of the LHCb and Belle II experiments and, therefore, limits the sensitivity of the measured quantities to new physics. We analyze the infrared QCD structure of this contribution and find that all soft and collinear divergences either cancel between different diagrams or factorize into matrix elements of local four-quark operators up to terms suppressed by Λ(QCD)/m(ψ), where m(ψ) denotes the J/ψ mass. Our results, which are based on an operator product expansion, allow us to calculate the penguin-to-tree ratio P/T in terms of the matrix elements of these operators and to constrain the penguin contribution to the phase 2β as |Δϕ(d)|≤0.68°. The penguin contribution to 2β(s) is bounded as |Δϕ(s)(0)|≤0.97°, |Δϕ(s)(∥)|≤1.22°, and |Δϕ(s)(⊥)|≤0.99° for the case of longitudinal, parallel, and perpendicular ϕ and J/ψ polarizations, respectively. Further, we place bounds on |Δϕ(d)| for B(d)→ψ(2S)K(S) and the polarization amplitudes in B(d)→J/ψK(*). In our approach, it is further possible to constrain P/T for decays in which P/T is Cabibbo unsuppressed, and we derive upper limits on the penguin contribution to the mixing-induced CP asymmetries in B(d)→J/ψπ(0), B(d)→J/ψρ(0), B(s)→J/ψK(S), and B(s)→J/ψK(*). For all studied decay modes, we also constrain the sizes of the direct CP asymmetries.