arxiv: v4 [hep-ph] 19 Apr 2017

Similar documents
arxiv: v3 [hep-ph] 12 Apr 2018

arxiv: v1 [hep-ph] 5 Dec 2014

PoS(FPCP2017)018. CP violation in the baryon sector. Eluned Smith RWTH Aachen On behalf of the LHCb collaboration

arxiv:hep-ph/ v4 18 Nov 1999

Lecture 12 Weak Decays of Hadrons

Polarization Study in B u,d,s,c Decays to Vector Final States

CKM Matrix and CP Violation in Standard Model

arxiv: v3 [hep-ph] 14 Jan 2015

CP Violation at the LHC - workshop 2017

Advances in Open Charm Physics at CLEO-c

Rare Hadronic B Decays

CP violation in charmless hadronic B decays

Electroweak Theory: 5

The Cabibbo-Kobayashi-Maskawa (CKM) matrix

arxiv:hep-ph/ v1 26 Feb 1998

Hadronic Effects in B -Decays

P -wave Dilution, Penguin and Rescattering Effects

Recent CP violation measurements. Advanced topics in Particle Physics: LHC physics, 2011 Jeroen van Tilburg 1/38

Long distance weak annihilation contribution to

Measurements of CP violating phases in B decays at LHCb

B Factories. Alan Watson University of Birmingham, UK

Hiroyuki Sagawa KEK OHO 1-1, Tsukuba, Ibaraki, Japan

arxiv: v1 [hep-ex] 21 Jan 2012

arxiv: v1 [hep-ph] 14 Sep 2018

CP violation in the quark sector: What have we learned?

Absolute Branching Fractions of Λ + c at BESIII

arxiv:hep-ph/ v3 10 Jun 2005

Moriond QCD La Thuile, March 14 21, Flavour physics in the LHC era. An introduction. Clara Matteuzzi. INFN and Universita Milano-Bicocca

Hadronic B/D decays in factorization assisted topology diagram approach

Measurement of CP violation in B J/ψK 0 S. decays. Frank Meier TU Dortmund. XXIX Rencontres de Physique de la Vallée d Aoste March 1 7, 2015

PoS(FPCP 2010)017. Baryonic B Decays. Jing-Ge Shiu SLAC-PUB National Taiwan University, Taiwan

Elementary Particle Physics

CP Violation: A New Era

Theoretical study of rare B/D decays

La Fisica dei Sapori Pesanti

Department of Physics and Astronomy, University of California, Riverside, CA, USA

The weak interaction Part II

FLAVOR PHYSICS BEYOND THE STANDARD MODEL

Recent CP violation measurements

Angular Analysis of the Decay

CP Violation Beyond the Standard Model

Charmless b-meson and b-baryon decays at LHCb. Adam Morris

arxiv:hep-ph/ v1 30 Dec 1994

Fleischer Mannel analysis for direct CP asymmetry. Abstract

PoS(CKM2016)113. Measurement of γ from B meson decay to D ( ) K ( )

Theory of CP Violation

Parity violation. no left-handed ν$ are produced

A Summary of Ξ ο Physics at KTEV. Erik J. Ramberg Fermilab 2 December, 2005

ub + V cd V tb = 0, (1) cb + V td V

Flavor physics. Yuval Grossman. Cornell. Y. Grossman Flavor physics APS2009, Denver, 5/2/2009 p. 1

arxiv: v2 [hep-ph] 29 Oct 2015

CP Violation in B Decays and the CKM Matrix. Emmanuel Olaiya Rutherford Appleton Laboratory Chilton,Didcot,Oxon,OX11 0QX,UK

Remainder of Course. 4/22 Standard Model; Strong Interaction 4/24 Standard Model; Weak Interaction 4/27 Course review 5/01 Final Exam, 3:30 5:30 PM

Search for new physics in three-body charmless B mesons decays

Electroweak interactions of quarks. Benoit Clément, Université Joseph Fourier/LPSC Grenoble

Recent BaBar results on CP Violation in B decays

Amplitude analyses with charm decays at e + e machines

Particle Physics. Lecture 12: Hadron Decays.!Resonances!Heavy Meson and Baryons!Decays and Quantum numbers!ckm matrix

Theory of hadronic B decays

Adrian Bevan Department of Physics Liverpool University Liverpool, United Kingdom (from the BABAR Collaboration.)

Lessons for New Physics from CKM studies

HADRONIC B DECAYS. My Ph.D. work. Maxime Imbeault. Université de Montréal 24/04/2009

Standard Model of Particle Physics

Elementary Particles, Flavour Physics and all that...

Measurements of the phase φ s at LHCb

Flavour physics in the LHC era

Neutron Beta-Decay. Christopher B. Hayes. December 6, 2012

Theory and Phenomenology of CP Violation

The other window on New Physics: CP violation at the B factories

Charmless B decays in Factorization Assisted Topological amplitude approach

Recent results on CKM/CPV from Belle

Understanding the penguin amplitude in B φk decays

Quark Model of Hadrons

arxiv: v2 [hep-ph] 15 Jan 2018

Direct CPV. K.Trabelsi kek.jp

Introduction to Operator Product Expansion

An Introduction to the Standard Model of Particle Physics

Shahram Rahatlou University of Rome

Hadronic Charm Decays: Experimental Review

Precision measurement of

Unitary Triangle Analysis: Past, Present, Future

Flavour Physics at hadron machines

arxiv:hep-ph/ v1 17 Oct 2003

Theoretical Issues in B PV Decays

Measurement of ϕ 2 /α using

CP violation in B 0 π + π decays in the BABAR experiment. Muriel Pivk, CERN. 22 March 2004, Lausanne

Lecture III. Measurement of sin 2 in B K 0 S. Measurement of sin 2 in B + -, B + - Measurement of in B DK. Direct CP Violation

IoP Masterclass FLAVOUR PHYSICS. Tim Gershon, University of Warwick April 20 th 2007

arxiv: v2 [hep-ex] 8 Aug 2013

Recent V ub results from CLEO

Status of the Cabibbo-Kobayashi-Maskawa Quark Mixing Matrix

12 Best Reasons to Like CP Violation

Weak Decays, CKM, Anders Ryd Cornell University

Lecture 11. Weak interactions

Particules Élémentaires, Gravitation et Cosmologie Année Le Modèle Standard et ses extensions. The Flavour Sector

Measuring the Unitarity Triangle Angle α with BaBar. Adrian Bevan KEK, 15 th September 2006

A. J. Schwartz Physics Department, University of Cincinnati, Cincinnati, Ohio USA

Particle Physics II CP violation. Lecture 3. N. Tuning. (also known as Physics of Anti-matter ) Niels Tuning (1)

arxiv:hep-ph/ v1 22 Mar 1999

Cracking the Unitarity Triangle

Transcription:

Rates and CP asymmetries of Charmless Two-body Baryonic B uds Decays Chun-Khiang Chua Department of Physics and Center for High Energy Physics Chung Yuan Christian University Chung-Li Taiwan 30 Republic of China arxiv:161.0449v4 [hep-ph] 19 Apr 017 (Dated: September 4 018) Abstract With the experimental evidences of B 0 pp and B Λp decays it is now possible to extract both tree and penguin amplitudes of the charmless two-body baryonic B decays for the first time. The extracted penguin-tree ratio agrees with the expectation. Using the topological amplitude approach with the experimental results on B 0 pp and B Λ p decay rates as input predictions on all other B q BB BD DB and DD decay rates where B and D are the low lying octet and decuplet baryons respectively are given. It is non-trivial that the results do not violate any existing experimental upper limit. From the analysis it is understandable that why B 0 pp and B Λ p modes are the first two modes with experimental evidences. Relations on rates are verified using the numerical results. We note that the predicted B p ++ rate is close to the experimental bound which has not been updated in the last ten years. Direct CP asymmetries of all B q BB BD DB and DD modes are explored. Relations on CP asymmetries are examined using the numerical results. The direct CP asymmetry of B 0 pp decay can be as large as ±50%. Some of the CP asymmetries can serve as tests of the Standard Model. Most of them are pure penguin modes which are expected to be sensitive to New Physics contributions. In particular B 0 s Ξ Ξ B 0 Ξ Σ B 0 Ω Ξ B 0 s Σ Σ B 0 s Ω Ω B 0 s Ξ Ξ B 0 s Ξ Ξ B 0 Ξ Σ B 0 Ω Ξ and B 0 s Ξ Ξ decays are S = 1 pure penguin modes with unsuppressed rates which can be searched in near future. Their CP asymmetries are constrained to be of few % and are good candidates to be added to the list of the tests of the Standard Model. 1

I. INTRODUCTION Recently LHCb collaboration reported the evidence for the first penguin dominated charmless two-body baryonic mode B Λ p decay at 4.1σ level [1] giving B(B Λ p) = (.4 +1.0.8 ± 0.3) 10 7. (1) The result is consistent with the predictions given in a pole model calculation [] and a topological amplitude approach [3]. The latter work made use of the large m B mass the experimental rate of the tree-dominated B 0 p p decay [4] and a naïve estimation on tree-penguin ratio. B(B 0 pp) = (1.47 +0.6+0.35.51.14) 10 8 () In fact it was advocated in [3] that the B Λ p decay mode could be the second charmless two-body baryonic mode to be found experimentally. With both the experimental evidences on B 0 pp and B Λ p decays it is now possible to extract both tree and penguin amplitudes at the same time. Most of the decay amplitudes of the two-body baryonic decays are non-factorizable. 1 Various models such as pole model [ 6 8] sum rule [9] diquark model [10 11] and approaches employing flavor symmetry [1 15] were used (for recent reviews see [16 17]) to calculate the amplitudes. Nevertheless all predictions on the B 0 p p rate are off by several orders of magnitude comparing to the LHCb result [4 16 17]. Indeed the B 0 p p decay mode is more suppressed than expected. To see this one scales the B 0 Λ + c p rate by the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements V ub /V cb and with a possible dynamical suppression factor f dyn included [18] B(B 0 p p) = B(B 0 Λ + c p) V ub /V cb f dyn 10 7 f dyn. (3) The data demands f dyn 0.1. It was pointed out in Ref. [18] that for a given tree operator O i one needs to consider additional contribution which were missed in the literature from its Fiertz transformed operator O i. It was found that there are cancelations of Feynman diagrams induced by O i by that from O i and consequently the smallness of the treedominated charmless two-body baryonic B decays results from this partial cancellation. 1 For a study on the factorization contributions see Ref. [5].

Furthermore as pointed out in Ref. [18] the internal W -emission tree amplitude should be proportional to the Wilson coefficient combination c 1 + c rather than c 1 c where the latter is usually claimed in the literature. We shall adapt this point in the present work. We can make use of the newly measured B 0 p p and B Λ p rates to give information on other modes in a symmetry related approach. The quark diagram or the so-called topological approach has been extensively used in mesonic modes [19 4]. In fact the approach is closely related to the flavor SU(3) symmetry [19 5]. In [6] the approach was extended to the charmless two-body baryonic case. It was further developed in [3] where amplitudes for all low lying charmless two-body baryonic modes with full topological amplitudes are obtained. Note that in general a typical amplitude has more than one tree and one penguin amplitudes. Asymptotic relations in the large m B limit [7] can be used to relate various topological amplitudes [3 6]. In this work using the experimental results on the B 0 p p and B Λ p decay rates we will extract both tree and penguin amplitudes for the first time. Rates and direct CP asymmetries of all low lying charmless two body baryonic decays can be explored. Rates and CP asymmetries of some modes can be checked experimentally in the near future in LHCb and Belle-II. Note that CP asymmetries of some modes can be added to the list of the tests of the Standard Model. In particular S = 1 pure penguin modes have small CP asymmetries and they are expected to be sensitive to New Physics contributions. These modes are good candidates to be added to the lists of the tests of the Standard Model especially for those with unsuppressed rates. The layout of this paper is as following. In Sec. II we briefly review and update our formulation for charmless two-body baryonic decays modes. In Sec. III we present the numerical results on rates and direct CP asymmetries of all low lying baryon modes. Relations on rates and CP asymmetries will also be studied. Conclusion is given in Sec. IV which is followed by two appendices. II. TWO-BODY CHARMLESS BARYONIC B DECAY AMPLITUDES There are more than 160 B BB BD DB DD decay amplitudes with B and D the low lying octet and decuplet baryons respectively. Their decay amplitudes expressed in term of topological amplitudes can be found in [3] and are collected in Appendix A since they will 3

be used extensively in this work. We show a few examples here A(B Λp) = 1 (T + T ) + 1 (10P P ) 1 6 1BB 3BB 6 1BB BB 3 6 (P P 1EW BB EW BB 4P + 4P ) + 1 (10A 3EW BB 4EW BB 6 1BB A ) BB A( B 0 pp) = T BB + T 4BB + P BB + 3 P EW BB 5E 1BB + E BB 9P A BB A(B Ξ Ξ 0 ) = P BB + 1 3 P EW BB A BB A( B 0 s pp) = 5E 1BB + E BB 9P A BB (4) A(B p ++ ) = 6(T 1BD T BD ) + 6P BD + 3 P 1EW BD + 6A BD A( B 0 Ξ Σ ) = P BD 3 P 1EW BD (5) A(B 0 p) = T 1DB P DB + 3 (3P 1EW DB + P EW DB ) A DB A( B 0 Σ + p) = T DB + P DB + 3 P EW DB (6) A( B 0 0 0 ) = T DD + 4P DD + 3 P EW DD + E DD + 18P A DD A(B Σ + ++ ) = 3T DD + 3P DD + 4 3 P EW DD + 3A DD (7) where T ( ) P ( ) E ( ) A ( ) P A ( ) and P ( ) EW are tree penguin W -exchange annihilation penguin annihilation and electroweak penguin amplitudes respectively for S = 0( 1) decays (see Fig. 1). In most cases we needs more than one tree and one penguin amplitudes in the baryonic decay amplitudes. Note that from the above amplitudes B 0 p p decay is expected to be a tree dominated mode B Λ p decay a penguin dominated mode and B 0 s p p decay a suppressed mode. By considering the chirality nature of weak interaction and asymptotic relations [7] at the large m B limit (m q /m B m BD /m B 1) the number of independent amplitudes is significantly reduced [3 6]: T ( ) = T ( ) 1BBBB3BB4BB = T ( ) 1BDBD = T ( ) 1DBDB = T ( ) DD P ( ) = P ( ) 1BBBB = P ( ) BD = P ( ) DB = P ( ) DD P ( ) EW = P ( ) 1EW BBEW BB3EW BB4EW BB = P ( ) 1EW BDEW BD = P ( ) 1EW DBEW DB = P ( ) EW DD E 1BBBBBDDBDD A 1BBBBBDDBDD P A BBDD 0. (8) 4

b q i q k q l b q k q i q l q m q l q j q m q m q l q i q m (a) (b) b q i q l q m b q k q l q m q k q m q l q j q i q m q l q j (c) (d) b q l q m q n b q i q k q l q k q n q m q l q m q l q j q m (e) (f) FIG. 1: Flavor flow or topological diagrams of (a) T ( ) (tree) (b) P ( ) (penguin) (c) E ( ) (W - exchange) (d) A ( ) (annihilation) (e) P A ( ) (penguin annihilation) and (f) P ( ) EW penguin) amplitudes in B to baryon pair decays for S = 0( 1) decays. (electroweak and there are only one tree one penguin and one electroweak penguin amplitudes which are estimated to be T ( ) P ( ) P ( ) EW = V ub V = V tb V G f ud(s) = 3 V tbv (c 1 + c )χū (1 γ 5 )v G f td(s) [c 3 + c 4 + κ 1 c 5 + κ c 6 ]χū (1 γ 5 )v G f td(s) [c 9 + c 10 + κ 1 c 7 + κ c 8 ]χū (1 γ 5 )v (9) where c i are the next-to-leading order Wilson coefficients given by c 1 = 1.081 c =.190 c 3 = 0.014 c 4 =.036 c 5 = 0.009 c 6 =.04 5

c 7 =.011α EM c 8 = 0.060α EM c 9 = 1.54α EM c 10 = 0.3α EM (10) in the naive dimensional regularization scheme at scale µ = 4.GeV [8]. Note that the relative signs of c 139 and c 410 in Eq. (9) are determined according to Ref. [18]. Since the Fierz transformation of O 5678 are different from O 134 additional coefficients κ i in front of c 5(7) and c 6(8) in Eq. (9) are assigned. For simplicity we take κ 1 = κ = κ. The parameters χ and κ will be extracted from B 0 p p and B Λ p data. Note that κ is expected to be of order O(1). In reality the m B mass is finite. The decays amplitudes are not in the asymptotic forms. Some corrections are expected. The correction on T ( ) i P ( ) i and P ( ) EW i are estimated to be of order m B /m B (the baryon and B meson mass ratio). Hence we have with T ( ) i = (1 + r ( ) ti)t ( ) P ( ) i = (1 + r ( ) pi)p ( ) P ( ) EW i = (1 + r ( ) ewpi)p ( ) EW (11) r ( ) ti r ( ) pi r ( ) ewpi m p /m B (1) parametrizing the correction to the asymptotic relations Eq. (8). Note that the phases of these r ( ) can be any value. For P ( ) i we replace the CKM factor V tb Vtd(s) V ub V ud(s) and V cbv cd(s). The penguins with V ubv ud(s) and V cbv cd(s) by the sum of ( )u are u-penguin (P i ) and c-penguin (P ( )c i ) respectively. The r pi of the u-penguin and the c-penguin are independent. Furthermore although in the asymptotic limit the ū u and ū γ 5 u terms have the same coefficients see Eq. (9) this will no longer be true in the finite m B case. In other words the r ( ) for ū u and ū γ 5 u terms are independent. For subleading terms such as exchange penguin annihilation and annihilation amplitudes we have E ( ) i f B m B η i T ( ) A ( ) f B m B j η j T ( ) P A ( ) k m B m B m B m B f B η m B k P ( ) (13) m B m B where the ratio f B /m B is from the usual estimation the factor m B /m B is from the chirality structure. Note that η ijk are estimated to be of order 1 and explicitly we take 0 η ijk η = 1 with the bound η set to 1. With the above amplitudes and formulas given in Appendices A and B we are ready to study the two-body baryonic decay rates. Note that for B BD DB decays we need to add a correction factor of p cm /(m B /) to the topological amplitudes shown in Eq. (9) with p cm the momentum of the final state baryons in the center of mass frame. The factor will further correct the amplitudes from their asymptotic forms. 6

III. NUMERICAL RESULTS ON RATES AND DIRECT CP ASYMMETRIES A. Tree and penguin amplitudes Using the recent data on the B 0 p p rate and the B Λ p rate the unknown parameters χ and κ in asymptotic amplitudes Eq. (9) are fitted to be χ = (5.08 +1.1 1.0) 10 3 GeV κ = 1.9 +0.39.46. (14) Note that the value of κ is indeed of order 1. The tree-penguin ratio is close to the naïve estimation. The penguin-tree and tree-penguin ratios of the asymptotic amplitudes [see Eq. (8)] for S = 0 and 1 transitions respectively can be extracted for the first time to be P = 0.4 ± 0.04 T The above equation is one of the major results of this work. T P = 0.1+0.05.03. (15) B. Numerical Results on Rates In the following we make use the B 0 p p and B Λ p data as inputs and predict rates of all other B BB BD DB DD modes. Results are shown in Tables I VIII. They are major results of this work. A first glimpse on these tables reveals that all experimental upper bounds are satisfied. This is a non-trivial check. We will go through the discussions of the updated results on B BB BD DB DD decay rates in below and make suggestions on future experimental searches. We will give a summary of our suggestions at the end of this subsection. 1. Rates of B q BB decays Predictions on S = 0 B q BB decay rates are shown in Table I. The first uncertainty is from the uncertainties of the χ and κ parameters [see Eq. (14)] reflecting the uncertainties in the measurements of B(B 0 pp) and B(B Λ p). The second uncertainty is obtained 7

TABLE I: Decay rates of S = 0 B q BB modes. The first uncertainty is from the uncertainties of the χ and κ parameters reflecting the uncertainties in the measurements of B(B 0 pp) and B(B Λ p) the second uncertainty is obtained by varying the tree and penguin strong phase φ the third uncertainty is from relaxing the asymptotic relations by varying r ti r pi r ewpi (see Eqs. (11) and (1)) and the last uncertainty is from sub-leading contributions terms with η ijk (see Eq. (13)). Occasionally the last uncertainty is shown to larger decimal place. The latest experimental results are given in parentheses under the theoretical results. The experimental B 0 p p rate is one of the inputs. Mode B(10 8 ) Mode B(10 8 ) B np 3.45 +1.50 B Σ 0 Σ + 3.9 +1.56 B Σ Σ 0 0.6 +0.8 B Σ Λ 0.47 +0.1 1.3 +0.68 +1.99 1.18 +0.5 +.09 1.50 ± 0.09 B0 s pσ + 1.4 +0.69 1.58 ± 0.11 B0 s nσ 0 0.75 +0.36. ± 0+0.4.31 +0.006.17 ± 0+0.0.17 +0.003.004 B 0 s nλ.96 +1.36.00 B 0 s Σ 0 Ξ 0 10.85 +5.3.51 +0.13 +.01 1.1 ± 0.7 +0.05 +0.3.6 ± 0 1.06 +0.57 +1.88 1.41 ± 0 3.90 +1. +5.3 4.0 ± 0 B Ξ Ξ 0 0.07 +0.03.03 ± 0 ± 0.03+0.0005.0003 B 0 s Σ Ξ 1.76 +0.78.63 ± 0+0.76.63 ± 0 B ΛΣ + 0.47 +0.1 B 0 pp 1.47 +0.71.17 ± 0+0.38.17 +0.003.53 +0.14 +.07 1.47 +0.6+0.35.51.14 [4] B 0 nn 6.66 +3.15.39 +1.05 +4.5.00 B 0 s ΛΞ 0 0.11 +0.05.04 ± 0+0.63.08 ± 0 1.16 ± 0.1 B0 Σ + Σ + 0 ± 0 ± 0 ± 0 +0.003 3.0 ± 0.07 B0 Σ 0 Σ 0 1.5 +0.7.55 +0.4 +0.97.73 ± 0.07 B 0 Ξ 0 Ξ 0 0 ± 0 ± 0 ± 0 +0.0004 B 0 Σ Σ 1.15 +0.51.41 ± 0+0.78.58 ± 0.04 B 0 Ξ Ξ 0.07 +0.03.0 ± 0+0.03.0 ± 0.01 B0 Σ 0 Λ 4.10 +1.98 1.47 +0.39 +1.90 1.54 ± 0.05 B 0 ΛΛ 0 ± 0 ± 0 +0.3 +0.0005 B 0 ΛΣ 0 0. +0.10.08 ± 0+0.18.08 ± 0.001 (< 3) [9] by varying the tree and penguin strong phase φ where we assign to the penguin amplitude: P ( ) = V tb V G f td(s) [c 3 + c 4 + κc 5 + κc 6 ]χe iφ ū (1 γ 5 )v (16) and a similar expression for P ( ) EW. We use a common strong phase for simplicity. The third uncertainty is from relaxing the asymptotic relations by varying r ti r pi r ewpi [see Eqs. (11) The strong phase of the tree amplitude is factored out and set to zero. Therefore the strong phase φ is the relative strong phase of penguin and tree amplitudes. 8

and (1)] and the last uncertainty is from sub-leading contributions such as annihilation penguin annihilation and exchange amplitudes [see Eq. (13)]. From the tables we see that errors are reduced at least by a factor of two compared to the previous analysis in [3]. These errors can provide useful informations: (i) As noted before the first errors reflect the uncertainties in B(B 0 pp) and B(B Λ p). (ii) The second errors reflect the size of tree-penguin interferences. We can see from the table that in general the effects of tree-penguin interference on rates are not sizable. This is consistent with the tree-penguin ratios shown in Eq. (15). (iii) The third errors which correspond to corrections to amplitudes away from the asymptotic limit are usually the largest ones. (iv) Occasionally we show the last errors which are from sub-leading contributions to larger decimal place. These terms are with η ijk [see Eq. (13)]. Note that for modes only have sub-leading contributions the rates are proportional to η ijk while for modes having tree and/or penguin terms as well these (sub-leading) contributions are roughly proportional to η ijk. For S = 0 B q BB decays there are three modes that can decay cascadely to all charged final states namely B 0 p p B 0 Ξ Ξ and ΛΛ decays. 3 We see from Table I that among these modes the B 0 p p decay has the highest rate and highest detectability. The rates of the other two modes are in fact one or two orders of magnitude smaller the p p rate. In particular the predicted B 0 ΛΛ rate is much smaller than the present experimental upper limit [9]. More modes can be searched for with π 0 γ π 0 π 0 and π 0 γ in the future experiments such as in Belle II. For example with one π 0 or γ one can search for B 0 Σ 0 Λ and B 0 s pσ + decays with π 0 γ one can also search for B 0 s Σ 0 Ξ 0 and B Σ 0 Σ + decays while with γγ one can search for B 0 Σ 0 Σ 0 decays. In fact the B 0 s Σ 0 Ξ 0 decay rate is of the order of 10 7 which is the highest rate in the table but the mode is reconstructed through the cascade decay B 0 s Σ 0 Ξ 0 Λγ Λπ 0 which requires both γ and π 0 for the detection. It is understandable that why B 0 p p decay is the first mode with experimental evidence. 3 To study the accessibilities of searching of the charmless two-body baryonic modes we note that (i) ++0 Λ Ξ Σ ± Ξ 0 and Ω have non-suppressed decay modes of final states with all charged particles (ii) + Σ + Ξ 0 Σ 0 and Ξ can be detected by detecting a π 0 (iii) Σ 0 can be detected by detecting γ (iv) one needs n in detecting and Σ [3 30]. Note that although some particles such as Ξ and Ξ 0 can decay to final states with all charged particles they may suffer from low reconstruction efficiencies. 9

TABLE II: Same as Table I but with S = 1 B q BB modes. The latest experimental result is given in the parenthesis under the theoretical results. The experimental B Λ p rate is one of the inputs. Mode B(10 8 ) Mode B(10 8 ) B Σ 0 p 0.76 +0.35.7 +0.19 +0.50.37 ± 0.001 B0 Σ + p 1.83 +0.76 B Σ n 1.67 +0.74 B Ξ 0 Σ + 39.98 +17.53.59 ± 0+0.71.58 +0.00 14.3 +.14 +17.06.00 B 0 Σ 0 n 1.1 +0.47 14.0 ± 0.03 B0 Ξ 0 Σ 0 18.50 +8.11 B Ξ Σ 0 19.40 +8.59 6.90 ± 0+8.38 6.88 ± 0.0 B0 Ξ 0 Λ.59 +1.07.65 +0.47 +0.8.61 ± 0.40 +0.5 +0.47.36 ± 0 6.58 +0.99 +7.89.9 +0.67 +.80 6.48 ± 0 1.74 ± 0 B Ξ Λ.36 +1.05.84 ± 0+.65 1.67 ± 0.005 B0 Ξ Σ 35.88 +15.89 1.77 ± 0+15.50 1.73 ± 0 B Λp 4.00 +10.44 8.54 +.13 +1.48 4 +10 8 ± 3 [1] 9.85 ± 0.0 B0 Λn 3.48 +10.00 B 0 s pp 0 ± 0 ± 0 +0.007 B 0 s Σ + Σ + 1.76 +0.73 (.84 +.03+0.85 1.68.18 ) [4] B 0 s nn 0 ± 0 ± 0 +0.007 B 0 s Σ 0 Σ 0 1.61 +0.69 B 0 s Ξ 0 Ξ 0 4.46 +10.53 B 0 s Ξ Ξ.63 +10.0 B 0 s ΛΛ 14.90 +6.4 8.71 +3.4 +16.8 1.07 +0.75 8.05 ± 0+15.7 11.36 +0.7 5.31 +1.97 +7.58 5.99 +0.61.74 B 0 s Σ Σ 1.49 +0.66.71 B 0 s Σ 0 Λ 0.05 +0.03 8.36 +4.1 +1.13 9.50 ± 0.63 +0.45 +0.79.59 +0.1.0.57 +0.1 +0.69.55 +0.0.19.53 ± 0+0.63.5 +0.19.18.0 +0.04 +0.06.04 ± 0.001.60 B 0 s ΛΣ 0 0.05 +0.03.0 +0.04 ± 0.0 ± 0.001 From Eqs. (A1) (A) and (A3) we see that there are several modes without any tree (T ibb ) contribution. These include B Σ Σ 0 B Σ Λ B Ξ Ξ 0 B 0 Ξ Ξ B 0 Σ Σ B 0 Ξ 0 Ξ 0 B 0 Σ + Σ + and B 0 s Σ Ξ decays. As shown in Table I the second uncertainties of the rates of these modes are vanishing. Note that B 0 Ξ Ξ B 0 Σ Σ and B 0 s Σ Ξ decays are pure penguin modes which only have P ibb P iew BB and P A BB terms while B 0 Ξ 0 Ξ 0 and B 0 Σ + Σ + decays only have subleading contributions namely E ibb and P A BB. Note that although B ΛΣ + B 0 s ΛΞ 0 B 0 ΛΛ and B 0 ΛΣ 0 decays have tree amplitudes T ibb these tree amplitudes cancel out in the asymptotic limit [see Eqs (A1) (A) (A3) and (8)]. In particular the tree penguin electroweak penguin and exchange amplitudes in the B 0 ΛΛ amplitude all cancel out in the asymptotic limit. This mode is therefore sensitive to the correction to the asymptotic relations Eqs. (8) (11) and (1). 10

Predictions on S = 1 B q BB decay rates are shown in Table II. There are 9 modes that have rates of order 10 7 namely B Ξ 0 Σ + Ξ Σ 0 Λ p B 0 Ξ 0 Σ 0 Ξ Σ Λ n B 0 s Ξ 0 Ξ 0 Ξ Ξ and ΛΛ decays. On the other hand there are 5 modes that can cascadely decay to all charged final states namely B Ξ Λ Λ p B 0 s p p Ξ Ξ and ΛΛ. Comparing these two sets we see that B Λ p B 0 s ΛΛ and Ξ Ξ decays are the only three modes that have rates of order 10 7 and can cascadely decay to all charge final states. In fact the B Λ p has the highest rate among these three modes and has the best detectability as the others need both Λ and Λ for detections (B 0 s ΛΛ B 0 s Ξ Ξ Λπ Λπ + ). It is understandable that why B Λ p is the first penguin mode with experimental evidence. Nevertheless it is interesting to search for B 0 s ΛΛ and Ξ Ξ decays as well. One can also search for other modes. For example the B Ξ Λ mode can also cascadely decay to all charged final state and its rate is of order 10 8 but this mode suffers from the low reconstruction efficiency of Ξ. With γ one can search for B Ξ Σ 0 decay which has rate of order 10 7. With π 0 one can search for B 0 Ξ 0 Λ and B 0 Σ + p decays at 10 8 level. With π 0 π 0 one can search for B Ξ 0 Σ + B 0 s Ξ 0 Ξ 0 and B 0 s Σ + Σ + decays where the first two modes have rate of order 10 7 while with π 0 γ one can search for B 0 Ξ 0 Σ 0 which has rate of order 10 7 and finally with γγ one can search for B 0 s Σ 0 Σ 0 decay at 10 8 level. Note that with all charged final states γ π 0 π 0 and π 0 γ most modes having rates of the order of 10 7 can be searched for in the future. From Eqs. (A4) (A5) and (A6) we see that B Σ n B Ξ Σ 0 B Ξ Λ B 0 Ξ Σ B 0 s Σ Σ B 0 s Ξ Ξ B 0 s pp and B 0 s nn decays do not have any tree (T ) contribution. As shown in Table II the rates of these modes have vanishing second ibb uncertainties. In particular we note that B 0 Ξ Σ B 0 s Σ Σ and B 0 s Ξ Ξ are pure penguin modes which only have P P and P ibb iew BB A terms while BB B0 s pp and B 0 s nn decays only have subleading contributions the E and P ibb A terms. We will BB return to these modes later. The predicted B 0 s p p rate is several order smaller than the present experimental result which however has large uncertainty. As noted this mode only receives contributions from sub-leading terms [see Eq. (4)]. One may enhance the subleading contributions but will soon run into contradictions. For example after enhancing the subleading contributions in the B 0 p p mode the so-called subleading contributions will oversize the leading tree contribution. Note that some enhancement on subleading contributions is possible in 11

TABLE III: Same as Table I but with S = 0 B q BD modes. The latest experimental result is given in the parenthesis. Mode B(10 8 ) Mode B(10 8 ) B p ++ 6.1 +3.01.3 +0.58 +8.77 (< 14) [31] 4.89 ± 0.08 B0 s pσ + 1.84 +0.89 B n +.18 +1.05.79 +0.15 +0.93.75 ± 0.03 B0 s nσ 0 0.97 +0.47 B Σ 0 Σ + 3.14 +1.53 1.13 +0.17 +1.34 1.10 ± 0.0 B0 s Σ 0 Ξ 0.77 +1.35.66 +0.17 +.60 1.45 ± 0.35 +0.07 +0.41.33 ± 0 1.00 +0.15 +1.19.97 ± 0 B Σ Σ 0 0.04 ± 0.0 ± 0 ± 0.0 +0.0003.000 B 0 s Σ Ξ 0.08 ± 0.03 ± 0 ± 0.03 ± 0 B Ξ Ξ 0 0.07 +0.03.0 ± 0+0.03.0 +0.0004.0003 B 0 s Ξ Ω 0.18 +0.08.07 ± 0+0.08.06 ± 0 B ΛΣ + 0.14 +0.06.05 ± 0+0.1.05 ± +0.001 B0 s ΛΞ 0 0.1 +0.05.04 ± 0+0.19.05 ± 0 B 0 p + 1.9 +0.93.69 +0.18 +.71 1.51 ± 0.0 B0 Σ + Σ + 0 ± 0 ± 0 +0.0001 B 0 n 0.01 +0.97.73 +0.14 +0.86.69 ± 0.03 B0 Σ 0 Σ 0 1.45 +0.71.5 +0.08 +0.6.51 ± 0.01 B 0 Ξ 0 Ξ 0 0 ± 0 ± 0 +0.0001 B 0 Σ Σ 0.08 +0.04.03 ± 0 ± 0.03 ± 0 B 0 Ξ Ξ 0.06 +0.03.0 ± 0+0.03.0 ± 0 B0 ΛΣ 0 0.06 +0.03.0 ± 0+0.10.0 +0.0004.0003 the presence of final state rescattering (see for example [3 33]) but it is unlikely that the enhancement to be so significant. Note that in Ref. [5] when partial conservation of axial-vector current is relaxed the calculated B s p p rate can be close to data but the predicted B Λ p rate is of the order 10 8 and is in tension with the data [1]. We need a more precise measurement on the B s p p rate to settle the issue.. Rates of B q BD decays Predictions on rates of S = 0 B q BD decays are shown in Table III. Modes that can cascadely decay to all charged final states and with unsuppressed rates are B p ++ and B 0 s pσ + decays. In particular the predicted rate of B p ++ decay is the highest one in the table and is just roughly a factor of smaller than the present experimental upper bound [31] which however has not been updated for quite a while. This mode could be just around the corner. On the other hand four other modes that can cascadely decay to all charged final states namely B Ξ Ξ 0 ΛΣ + B 0 s Ξ Ω and ΛΞ 0 decays are one or two order of magnitude more suppressed. Note that with γ one can search for B Σ 0 Σ + 1

and B 0 s Σ 0 Ξ 0 decays with π 0 one can search for B 0 p + decay while with γπ 0 one can search for B 0 Σ 0 Σ 0 decay in the future. All of these modes have rates of orders 10 8. From Eqs. (A7) (A8) and (A9) we see that B Σ Σ 0 B Ξ Ξ 0 B 0 s Σ Ξ B 0 Ξ Ξ B 0 Σ Σ B 0 s Ξ Ω B 0 Σ + Σ + and B 0 Ξ 0 Ξ 0 decay amplitudes do not have any tree (T ibd ) contribution. As shown in Table III the rates of these modes have vanishing second uncertainties. We note that B 0 Ξ Ξ and B 0 Σ Σ decays are pure penguin modes with only P BD and P iew BD contributions while B 0 Σ + Σ + and B 0 Ξ 0 Ξ 0 decays are pure exchange modes with only E BD contributions. Although B ΛΣ + B 0 s ΛΞ 0 and B 0 ΛΣ 0 decay amplitudes have tree amplitudes these tree amplitudes cancel out in the asymptotic limit. There are relations among rates. Using formulas in Appendices A and B we have [3] B(B Ξ Ξ 0 ) = B(B Σ Σ 0 ) B(B Σ Σ 0 ) ( pcm (B Ξ ) Ξ 0 3 ) p cm (B Σ Σ 0 ) 3τ Bd B( B 0 Σ Σ ) 3τ Bd B( B 0 Ξ Ξ ) 3τ Bs B( B 0 s Σ Ξ ) τ Bs B( B 0 s Ξ Ω ) B( B 0 Σ + Σ + ) B( B 0 Ξ 0 Ξ 0 ). (17) One can check from Table III that the rates of these modes roughly satisfy the above relations and the agreement will be improved when the SU(3) breaking effects are taken into account. Note that these relations do not relay on the asymptotic relations. Predictions on S = 1 B q BD decay rates are shown in Table IV. Both B 0 Ξ Σ and Λ 0 modes can cascadely decay to all charge final states where only the B 0 Ξ Σ rate can reach 10 8 level. In principle this mode can be detected through B 0 Ξ Σ Λπ Λπ + decay but may be restricted by the low reconstruction efficiency of Ξ. Note that with π 0 one can search for B 0 s Σ + Σ + B Ξ 0 Σ + B Σ + ++ B 0 s Ξ Ξ B 0 s Ξ 0 Ξ 0 and B Ξ Σ 0 decays which have rates of orders 10 8 while with γ one can search for B 0 Σ 0 0 decay in the future. In particular the B Σ + ++ decay rate is the highest one in the table. With γπ 0 π 0 π 0 or γγ one can search for B Σ 0 + B 0 Σ + + B 0 s Σ 0 Σ 0 and B 0 Ξ 0 Σ 0 decays which have rates of orders (or close to) 10 8 in the future. Note that the predicted B 0 Λ 0 and B Λ + rates are both 13

TABLE IV: Same as Table I but with S = 1 B q BD modes. The latest experimental results are given in parentheses. Mode B(10 8 ) Mode B(10 8 ) B Σ + ++ 6.99 +.90 B Σ 0 + 4.5 +1.83.49 +1.80 +3.13 1.51 +0.56 +1.83.33 ± 0.00 B0 Σ + +.16 +0.90 1.46 ± 0.003 B0 Σ 0 0 3.93 +1.69.77 +0.56 +0.97.7 ± 0 1.40 +0.5 +1.69 1.35 ± 0 B Σ 0 1.96 +0.87.70 ± 0+0.83.69 ± 0.00 B0 Σ 5.46 +.4 1.94 ± 0+.3 1.91 ± 0 B Ξ 0 Σ + 1.49 +0.67.53 +0.38 +0.93.70 ± 0.00 B0 Ξ 0 Σ 0 0.69 +0.31.4 +0.17 +0.43.33 ± 0 B Ξ Σ 0 0.81 +0.36.9 ± 0+0.34.8 ± 0.001 B0 Ξ Σ 1.49 +0.66.53 ± 0+0.63.5 ± 0 B Λ + 0.15 +0.07.05 +0.11 +0.07.05 ± 0 B0 Λ 0 0.14 +0.07 (< 8) [34] (< 93) [34] B 0 s p + 0 ± 0 ± 0 +0.000005 B 0 s Σ + Σ +.07 +0.86 B 0 s n 0 0 ± 0 ± 0 +0.000005 B 0 s Σ 0 Σ 0 1.89 +0.81.05 +0.10 +0.06.05 ± 0.74 +0.53 +0.93.67 +0.5 +0.81.69 ± 0.001.65 ± 0.001 B 0 s Ξ 0 Ξ 0 1.31 +0.59.47 +0.33 +0.8.6 ± 0.00 B0 s Σ Σ 1.74 +0.77.6 ± 0+0.74.61 ± 0 B 0 s Ξ Ξ 1.4 +0.63.51 ± 0+0.60.50 ± 0 B0 s ΛΣ 0 0.07 +0.03.0 +0.05 +0.03.0 ± 0.001 roughly two orders of magnitude below the present experimental bounds [34]. From Eqs. (A10) (A11) and (A1) we see that B Σ 0 B Ξ Σ 0 B 0 Σ B 0 Ξ Σ B 0 s Σ Σ B 0 s Ξ Ξ B 0 s p + and B 0 s n 0 decays do not have any tree (T ) contribution. As shown in Table IV the rates of these modes ibd have vanishing second uncertainties. Note that B 0 Σ B 0 Ξ Σ B 0 s Σ Σ and B 0 s Ξ Ξ decays are pure penguin modes with only P BD and P iew BD contributions while B 0 s p + and B 0 s n 0 are pure exchange modes with only E BD contributions. The rates of some modes are related. From Appendices A and B we have [3] B(B Ξ Σ 0 ) = B(B Σ 0 ) 3τ Bs B( B 0 s Σ Σ ) = 3τ Bs B( B 0 s Ξ Ξ ) = 3τ Bd B( B 0 Ξ Σ ) = τ Bd B( B 0 Σ ) B( B 0 s p + ) = B( B 0 s n 0 ) (18) where these relations are subjected to corrections from SU(3) breaking in p cm 3. These relations do not relay on the asymptotic relations. As shown in Table III the rates of these modes roughly satisfy the above relations and the agreement will be improved when the 14

TABLE V: Same as Table I but with S = 0 B q DB modes. The latest experimental result is given in the parenthesis. Mode B(10 8 ) Mode B(10 8 ) B 0 p.19 +1.05.79 +0.15 +0.9.74 ± 0.03 B0 s + Σ + 1.70 +0.8 (< 138) [31] B n 0.33 +0.15.1 ± 0+0.14.1 +0.00.001 B 0 s 0 Σ 0 0.96 +0.46.61 +0.16 +0.79.64 ± 0.34 +0.15 +0.50.40 ± 0 B Σ 0 Σ + 0.85 +0.41.31 +0.06 +0.36.9 ± 0.01 B0 s Σ 0.7 +0.1.10 ± 0+0.11.09 ± 0 B Σ Σ 0 0.04 ± 0.0 ± 0 ± 0.0 +0.0003.000 B 0 s Σ 0 Ξ 0.73 +1.33 B Ξ Ξ 0 0.07 +0.03.0 ± 0+0.03.0 +0.0004.98 +0.15 +1.17.96 ± 0.0003 B 0 s Σ Ξ 0.07 ± 0.03 ± 0 ± 0.03 ± 0 B Σ Λ 0.14 +0.06.05 ± 0+0.06.05 ± 0.001 B0 s 0 Λ.59 +1.7 B 0 + p 1.9 +0.93 B 0 0 n 7.46 +3.64.69 +0.18 +0.89.69 +0.40 +3.19.93 +0.03 +1.01.84 ± 0.7 ± 0.0 B0 Σ + Σ + 0 ± 0 ± 0 ± 0 +0.0001.6 ± 0.05 B0 Σ 0 Σ 0 0.39 +0.19.14 +0.03 +0.17.13 ± 0.005 B 0 Ξ 0 Ξ 0 0 ± 0 ± 0 ± 0 +0.0001 B 0 Σ Σ 0.08 +0.04.03 ± 0 ± 0.03 ± 0 B 0 Ξ Ξ 0.06 +0.03.0 ± 0+0.03.0 ± 0 B0 Σ 0 Λ 1.18 +0.57.4 +0.11 +0.55.44 ± 0.0 SU(3) breaking effects are taken into account. 3. Rates of B q DB decays Predictions on S = 0 B q DB decay rates are shown in Table V. The mode with the highest decay rate is B 0 0 n decay which is unfortunately hard to detect. Both B 0 p and B 0 s 0 Λ decays can cascadely decay to all charge final states. These two modes both have rates at the order of 10 8 but the B 0 p decay has better detectability. Note that the predicted rate of this mode is roughly two orders of magnitude below the present experimental limit [31] which however has not been updated since 008. With π 0 one can search for B 0 + p and B 0 Σ 0 Λ decays while B 0 s 0 Σ 0 decay which has a slightly smaller rate can be searched for with π 0 in the future. With π 0 π 0 one can search for B 0 s Σ 0 Ξ 0 B 0 s + Σ + and B Σ 0 Σ + decays in the future. These modes all have rates at the level of 10 8. Note that B n B Σ Σ 0 B Ξ Ξ 0 B Σ Λ B 0 Ξ Ξ B 0 s Σ B 0 s Σ Ξ B 0 Σ Σ B 0 Ξ 0 Ξ 0 and B 0 Σ + Σ + decays do not 15

have any tree (T idb ) contribution [see Eqs. (A13) (A14) and (A15)] and consequently their rates in Table V have vanishing second uncertainties. In particular the B 0 Ξ Ξ B 0 s Σ B 0 s Σ Ξ and B 0 Σ Σ decays are pure penguin modes which only have P DB and P iew DB contributions while B 0 Ξ 0 Ξ 0 and B 0 Σ + Σ + decays are pure exchange (E DB ) modes. The rates of some modes are related. Using the formulas in Appendices A and B we have [3] B(B 0 p) = B(B Σ 0 Σ + ) B(B n) = 3B(B Ξ Ξ 0 ) = 6B(B Σ Σ 0 ) = B(B Σ Λ) 3τ Bd B( B 0 Σ Σ ) = 3τ Bd B( B 0 Ξ Ξ ) = τ Bs B( B s 0 Σ ) = 3τ Bs B( B s 0 Σ Ξ ) B( B 0 Σ + Σ + ) = B( B 0 Ξ 0 Ξ 0 ) (19) where these relations are subjected to corrections from SU(3) breaking in p cm 3. These relations do not relay on the asymptotic relations. Note that the relations on B decay rates are new compared to those in [3]. As shown in Table V the rates of these modes roughly satisfy the above relations and the agreement will be improved when the SU(3) breaking effects are taken into account. Predictions on S = 1 B q DB decay rates are shown in Table VI. Note that B 0 Ω Ξ B 0 Ξ 0 Λ and B 0 Σ + p decays are the only three modes that can cascadely decay to all charged final states. All of them have rates of order 10 8. We need more data to search for them as the reconstruction efficiencies of the final states of the first two modes are low while the predicted B 0 Σ + p rate is one order of magnitude below the experimental limit which has not been updated since 007 [34]. Note that with one π 0 one can search for B Ω Ξ 0 B Ξ Λ B 0 s Ξ 0 Ξ 0 B 0 s Σ + Σ + B Ξ 0 Σ + B 0 s Ξ Ξ and B Σ 0 p decays which have rates of order 10 8 in the future. In particular the B Ω Ξ 0 decay has the highest rate in the table and the predicted B Σ 0 p rate is more than one order of magnitude below the experimental limit which has not been updated since 007 [34]. With π 0 γ one can search for B 0 s Σ 0 Σ 0 and B Ξ Σ 0 decays which have rates of order (or close to) 10 8 in the future. 16

TABLE VI: Same as Table I but with S = 1 B q DB modes. Mode B(10 8 ) Mode B(10 8 ) B Σ 0 p 0.94 +0.43.33 +0.4 +0.51.40 ± 0.001 B0 Σ + p.5 +0.93 < 47 [34] < 6 [34] B Σ n.05 +0.91.73 ± 0+0.87.7 ± 0.00 B0 Σ 0 n 1.39 +0.58 B Ξ 0 Σ + 1.44 +0.65.51 +0.37 +0.78.61 ± 0.00 B0 Ξ 0 Σ 0 0.67 +0.30.80 +0.58 +0.93.75 ± 0.49 +0.65 +0.54.45 ± 0.4 +0.17 +0.36.8 ± 0 B Ξ Σ 0 0.78 +0.35.8 ± 0+0.33.7 ± 0.001 B0 Ξ Σ 1.45 +0.64.51 ± 0+0.61.51 ± 0 B Ω Ξ 0 3.77 +1.67 1.34 ± 0+1.60 1.3 ± 0.003 B0 Ω Ξ 3.47 +1.54 1.4 ± 0+1.47 1.1 ± 0 B Ξ Λ.48 +1.10.88 ± 0+1.05.87 ± 0.00 B0 Ξ 0 Λ.71 +1.13 B 0 s + p 0 ± 0 ± 0 ± 0 +0.000005 B 0 s Σ + Σ + 1.98 +0.8 B 0 s 0 n 0 ± 0 ± 0 ± 0 +0.000005 B 0 s Σ 0 Σ 0 1.81 +0.78 B 0 s Ξ 0 Ξ 0 1.99 +0.83.71 +0.93 +0.78.64 ± 0.000 B0 s Σ Σ 1.67 +0.74.97 +0.70 +1.1.90 ± 0.71 +0.51 +0.8.64 +0.4 +0.74.59 ± 0+0.71.66 ± 0.001.61 ± 0.001.58 ± 0 B 0 s Ξ Ξ 1.36 +0.60.48 ± 0+0.58.47 ± 0 B0 s Σ 0 Λ 0.06 +0.03.0 +0.05 ± 0.0 ± 0.001 Note that B Σ n B Ξ Σ 0 B Ω Ξ 0 B Ξ Λ B 0 Ω Ξ B 0 s Σ Σ B 0 Ξ Σ B 0 s Ξ Ξ B 0 s + p and B 0 s 0 n decays do not have any tree (T ) contribution [see Eqs. (A16) (A17) and (A18)] and consequently their rates idb in Table VI have vanishing second uncertainties. In particular the B 0 Ξ Σ and B 0 s Ξ Ξ decays are pure penguin modes which only have P DB and P iew DB contributions while B 0 s + p and B 0 s 0 n decays are pure exchange (E DB ) modes. The rates of some modes are related. Using the formulas in Appendices A and B we have [3] B(B Σ 0 p) = B(B Ξ 0 Σ + ) 3B(B Σ n) = 6B(B Ξ Σ 0 ) = B(B Ω Ξ 0 ) = B(B Ξ Λ) 3τ Bs B( B 0 s Σ Σ ) = 3τ Bs B( B 0 s Ξ Ξ ) = 3τ Bd B( B 0 Ξ Σ ) = τ Bd B( B 0 Ω Ξ ) B( B 0 s + p) = B( B 0 s 0 n) (0) where these relations are subjected to corrections from SU(3) breaking in p cm 3. Note that 17

TABLE VII: Same as Table I but with S = 0 B q DD modes. Mode B(10 8 ) Mode B(10 8 ) B + ++ 17.15 +8.31 B 0 + 6.47 +3.07 B 0 0.9 +0.41 B Σ 0 Σ + 3.0 +1.43 6.17 +1.61 +7.97.33 +1.0 +3.39 6.45 ± 0. B0 s + Σ + 5.18 +.51.67 ± 0.14 B0 s 0 Σ 0.93 +1.39.33 ± 0+0.39.3 +0.006 1.08 +0.47 +1.58.004 B 0 s Σ 0.83 +0.37 1.4 ± 0.07 B0 s Σ 0 Ξ 0.73 +1.9 1.86 +0.49 +.41 1.95 ± 0 1.05 +0.46 +1.53.9 ± 0+0.35 1.1 ± 0.9 ± 0.98 +0.43 +1.43 1.1 ± 0 B Σ Σ 0 0.57 +0.5.0 ± 0+0.4.0 ± 0.003 B0 s Σ Ξ 1.03 +0.46.37 ± 0+0.44.36 ± 0 B Ξ Ξ 0 0.6 +0.1.09 ± 0+0.11.09 +0.00.001 B 0 s Ξ Ω 0.71 +0.31.5 ± 0+0.30.5 ± 0 B 0 ++ ++ 0 ± 0 ± 0 ± 0 +0.004 B 0 Σ + Σ + 0 ± 0 ± 0 ± 0 +0.00 < 1.1 10 4 [38] B 0 + + 5.9 +.56 B 0 0 0 5.99 +.83 1.90 +0.50 +.46.15 +0.94 +3.14 < 1.5 10 5 [38].00 ± 0.16 B0 Σ 0 Σ 0 1.40 +0.66.47 ± 0.13 B0 Σ Σ 1.05 +0.47.50 +0. +0.73.37 ± 0+0.45.58 ± 0.06.37 ± 0.07 B 0.55 +1.13.91 ± 0+1.08.89 ± 0.11 B0 Ξ 0 Ξ 0 0 ± 0 ± 0 ± 0 +0.001 B 0 Ω Ω 0 ± 0 ± 0 ± 0 +0.001 B 0 Ξ Ξ 0.4 +0.11.09 ± 0+0.10.08 ± 0.03 the relations on B decay rates are new compared to those in [3]. These relations do not relay on the asymptotic relations. The rates in Table IV roughly satisfy the above relations and the agreement will be improved when the SU(3) breaking effects are taken into account. 4. Rates of B q DD decays Predictions on S = 0 B q DD decay rates are shown in Table VII. There are six modes that can cascadely decay to all charged final states. They are B 0 ++ ++ 0 0 Ω Ω Σ + Σ + Σ Σ and Ξ 0 Ξ 0 decays. However most of them have highly suppressed rates except B 0 0 0 and B 0 Σ Σ decays which have rates of order 10 8 and should be searchable. In particular the bound on B 0 0 0 rate has not been updated for almost three decades [38]. Note that with π 0 one can search for plenty of unsuppressed modes. These modes include B + ++ B 0 + B 0 s + Σ + B Σ 0 Σ + B 0 s 0 Σ 0 B 0 s Σ 0 Ξ 0 B 0 s Σ Ξ and B 0 s Ξ Ω decays. In particular the mode with the highest rate (the only one at order 10 7 ) can be searched through the 18

B + ++ pπ 0 pπ decay. With π 0 π 0 one can also search for B 0 + + and B 0 Σ 0 Σ 0 decays. Note that B 0 B Σ Σ 0 B Ξ Ξ 0 B 0 B 0 Σ Σ B 0 Ξ Ξ B 0 s Σ B 0 s Σ Ξ B 0 s Ξ Ω B 0 ++ ++ B 0 Σ + Σ + B 0 Ξ 0 Ξ 0 and B 0 Ω Ω decays do not have any tree (T DD ) contribution [see Eqs. (A19) (A0) and (A1)] and consequently their rates in Table VII have vanishing second uncertainties. In particular the B 0 B 0 Σ Σ B 0 Ξ Ξ B 0 s Σ B 0 s Σ Ξ and B 0 s Ξ Ω decays are pure penguin modes which only have P DD P EW DD and P A DD contributions the B 0 ++ ++ B 0 Σ + Σ + and B 0 Ξ 0 Ξ 0 decays are subleading modes which only have E DD and P DD contributions while the B 0 Ω Ω decay is a pure penguin annihilation (P A DD ) mode. The rates of some modes are related. Using formulas in Appendices A and B we have [3] B(B 0 ) = 3B(B Σ Σ 0 ) = 6B(B Ξ Ξ 0 ) B(B 0 + ) = B(B Σ 0 Σ + ) B(Bs 0 0 Σ 0 ) = B(Bs 0 Σ 0 Ξ 0 ) 4B(Bs 0 Σ ) = 4B(Bs 0 Ξ Ω ) = 3B(Bs 0 Σ Ξ ) (1) where these relations are subjected to SU(3) breaking from the phase space factors. These relations do not relay on the asymptotic relations. As shown in Table VII the rates of these modes roughly satisfy the above relations and the agreement will be improved when the SU(3) breaking effects are taken into account. Predictions on S = 1 B q DD decay rates are shown in Table VIII. There are eight modes that can cascadely decay to all charged final states with unsuppressed rates. They are B 0 s Ω Ω B Ξ 0 Σ + B 0 s Ξ 0 Ξ 0 B Σ + ++ B Ω Ξ 0 B 0 s Σ + Σ + B Σ 0 and B 0 s Σ Σ decays. It is interesting that many (the first five) of them have rates of order 10 7. In particular the B 0 s Ω Ω decay has the highest rate and good detectability. Note that with one π 0 one can search for B 0 Ξ Σ B 0 Ω Ξ B 0 s Ξ Ξ B Σ 0 + B 0 Σ 0 0 and B 0 Ξ 0 Σ 0 decays in the future. All of them have rates of order 10 7. With π 0 π 0 one can search for B Ξ Σ 0 B 0 s Σ 0 Σ 0 and B 0 Σ + + decays where the first one has rate of order 10 7. Note that B Σ 0 B Ξ Σ 0 B Ω Ξ 0 B 0 Σ B 0 Ξ Σ B 0 Ω Ξ B 0 s Σ Σ B 0 s Ω Ω B 0 s Ξ Ξ B 0 s ++ ++ B 0 s + + 19

TABLE VIII: Same as Table I but with S = 1 B q DD modes. Mode B(10 8 ) Mode B(10 8 ) B Σ + ++ 1.48 +8.9 B Σ 0 + 13.08 +5.63 7.65 +5.53 +8.86 4.66 +1.73 +5.3 7.15 ± 0.007 B0 Σ + + 6.63 +.75 4.39 ± 0.008 B0 Σ 0 0 1.11 +5.1.36 +1.71 +.73.1 ± 0 4.31 +1.60 +4.93 4.06 ± 0 B Σ 0 6.06 +.68.16 ± 0+.57.1 ± 0.005 B0 Σ 16.83 +7.46 5.99 ± 0+7.15 5.88 ± 0 B Ξ 0 Σ + 4.6 +10.44 8.64 +3.1 +9.87 8.13 ± 0.01 B0 Ξ 0 Σ 0 11.3 +4.83 4.00 +1.49 +4.57 3.77 ± 0 B Ξ Σ 0 11.4 +4.98 4.00 ± 0+4.77 3.93 ± 0.01 B0 Ξ Σ 0.79 +9.1 7.40 ± 0+8.83 7.7 ± 0 B Ω Ξ 0 15.49 +6.86 5.51 ± 0+6.57 5.41 ± 0.01 B0 Ω Ξ 14.3 +6.34 5.10 ± 0+6.08 5.01 ± 0 B 0 s ++ ++ 0 ± 0 ± 0 ± 0 +0.0 B 0 s Σ + Σ + 6.49 +.69 B 0 s + + 0 ± 0 ± 0 ± 0 +0.0 B 0 s Σ 0 Σ 0 5.93 +.55 B 0 s 0 0 0 ± 0 ± 0 ± 0 +0.0 B 0 s Σ Σ 5.48 +.43 B 0 s 0 ± 0 ± 0 ± 0 +0.0 B 0 s Ξ 0 Ξ 0 1.9 +9.44 B 0 s Ω Ω 41.95 +18.58 14.93 ± 0+17.81 14.67 +1.79 1.75 B 0 s Ξ Ξ 0.9 +8.99.31 +1.67 +.67.16 +0.77.7.11 +0.79 +.41 1.99 +0.74.70 1.95 ± 0+.33 1.9.67 +0.7 7.81 +.90 +8.9 7.35 +1.36 1.3 7. ± 0+8.6 7.10 1.6 +1.30 B 0 s 0 0 and B 0 s decays do not have any tree (T ) contribution [see Eqs. DD (A) (A3) and (A4)] and consequently their rates in Table VIII have vanishing second uncertainties. The B 0 Σ B 0 Ξ Σ B 0 Ω Ξ B 0 s Σ Σ B 0 s Ω Ω and B 0 s Ξ Ξ decays are pure penguin modes which only have P DD P EW DD and P A DD contributions the B 0 s ++ ++ B 0 s + + and B 0 s 0 0 decays are subleading modes which only have E DD and P DD contributions while the B0 s decay is a pure penguin annihilation (P A DD ) mode. The rates of some modes are related. Using formulas in Appendices A and B we have [3] 4 6B(B Σ 0 ) = B(B Ω Ξ 0 ) = 3B(B Ξ Σ 0 ) B(B Σ 0 + ) = B(B Ξ 0 Σ + ) B( B 0 Σ 0 0 ) = B( B 0 Ξ 0 Σ 0 ) 4B( B 0 Σ ) = 3B( B 0 Ξ Σ ) = 4B( B 0 Ω Ξ ). () where these relations are subjected to SU(3) breaking from the phase space factors. These relations do not relay on the asymptotic relations. The rates in Table VIII roughly satisfy 4 A typo in the last relation is corrected. 0

the above relations and the agreement will be improved when the SU(3) breaking effects are taken into account. We have shown that B 0 p p and B Λ p decays have better chance to be found experimentally. We have identified several modes that can cascadely decay to all charged final states with (relatively) unsuppressed rates. They are searchable in near future. In particular we note that the predicted B p ++ rate is close to the experimental bound which has not been updated in the last ten years [31]. Furthermore the bounds on B 0 p and B 0 Σ + p rates have not been updated in the last ten years [31 34] and the bound on B 0 0 0 rate was last given in 1989 [38] while their rates are predicted to be of the order of 10 8. Also note that the B 0 s Ω Ω rate is predicted to be the highest rate. It will also be interesting for Belle-II which will be turned on soon to search for baryonic modes as there are many modes having unsuppressed rates but require π 0 or γ for detection. We pointed out several modes without tree amplitudes. Some of them are pure penguin modes. As we shall see in the next subsection these will affect their CP asymmetries. We summarize our suggestions for experimental searches in Table IX. These modes have (relatively) unsuppressed rates and better detectability. Modes are arranged according to their quantum numbers detection and rates (in descending order) and they are assigned into Groups I II and III accordingly where Group I modes are modes that have unsuppressed rates and can cascadely decay to all charged final states Group II modes are modes that can be searched with π 0 or γ and Group III modes are modes that can be searched with π 0 π 0 π 0 γ or γγ. C. Numerical Results on Direct CP Asymmetries In this subsection we will give results of direct CP asymmetries of all modes and plot asymmetries of several interesting modes. Note that this study becomes possible as we now have the information of the tree-penguin ratio Eq. (15). 1. CP asymmetries of B q BB decays In Table X we give results of direct CP asymmetries of S = 0 B q BB modes. The central values are the asymmetries generated from tree-penguin interference where only the 1

TABLE IX: Modes with (relatively) unsuppressed rates and (relatively) good detectability are summarized. Group I Group II Group III All charged final states with single π 0 /γ with π 0 π 0 π 0 γ or γγ BB S = 0 B 0 pp; B 0 Σ 0 Λ B 0 s pσ + ; B 0 s Σ 0 Ξ 0 B Σ 0 Σ + B 0 Σ 0 Σ 0 BB S = 1 B Λp B 0 s Ξ Ξ ; B Ξ Σ 0 B 0 Ξ 0 Λ; B Ξ 0 Σ + B 0 s Ξ 0 Ξ 0 B 0 s ΛΛ B Ξ Λ; B 0 Σ + p B 0 Ξ 0 Σ 0 B 0 s Σ + Σ + B 0 s Σ 0 Σ 0 BD S = 0 B p ++ B 0 s pσ + ; B Σ 0 Σ + B 0 s Σ 0 Ξ 0 ; B 0 Σ 0 Σ 0 B 0 p + BD S = 1 B 0 Ξ Σ ; B Σ + ++ B 0 Σ 0 0 ; B Σ 0 + B 0 Σ + + B 0 s Σ + Σ + B Ξ 0 Σ + ; B 0 s Σ 0 Σ 0 B 0 Ξ 0 Σ 0 B 0 s Ξ Ξ B 0 s Ξ 0 Ξ 0 B Ξ Σ 0 DB S = 0 B 0 p B 0 s 0 Λ; B 0 + p B 0 Σ 0 Λ; B 0 s Σ 0 Ξ 0 B 0 s + Σ + B 0 s 0 Σ 0 B Σ 0 Σ + DB S = 1 B 0 Ω Ξ B 0 Ξ 0 Λ ; B Ω Ξ 0 B Ξ Λ; B 0 s Σ 0 Σ 0 B Ξ Σ 0 B 0 Σ + p B 0 s Ξ 0 Ξ 0 B 0 s Σ + Σ + B Ξ 0 Σ + B 0 s Ξ Ξ B Σ 0 p DD S = 0 B 0 0 0 B 0 Σ Σ ; B + ++ B 0 + ; B 0 + + B 0 Σ 0 Σ 0 B 0 s + Σ + B Σ 0 Σ + B 0 s 0 Σ 0 B 0 s Σ 0 Ξ 0 B 0 s Σ Ξ B 0 s Ξ Ω DD S = 1 B 0 s Ω Ω B Ξ 0 Σ + ; B 0 Ξ Σ B 0 Ω Ξ ; B 0 s Ξ Ξ B Σ 0 + B 0 s Ξ 0 Ξ 0 B Σ + ++ ; B 0 Σ 0 0 B 0 Ξ 0 Σ 0 ; B Ξ Σ 0 B 0 s Σ 0 Σ 0 B Ω Ξ 0 B 0 s Σ + Σ + ; B 0 Σ + + B Σ 0 B 0 s Σ Σ

TABLE X: Direct CP asymmetries (A in %) of S = 0 B q BB modes for φ = 0 ±π/4 and ±π/. The uncertainties are from varying the strong phases of r ( ) ti r( ) pi r( ) ewpi and η ijk (see Eqs. (11) and (13)). Mode φ = 0 φ = ±π/4 φ = ±π/ Mode φ = 0 φ = ±π/4 φ = ±π/ B np 0 ± 3.1 (74.0 5.6 +19.7 ) (97.9+.1 15.6 ) B0 s pσ + 0 ± 1. (36.0 +1. 17.7 ) (49.3+0. 15. ) B Σ 0 Σ + 0 ± 31. (59.3 7. +3.4 ) (79.5+16.1.4 ) B0 s nσ 0 0 ± 14.9 ±(6.7 +14.6 13.8 ) ±(38.8+13.9 1.8 ) B Σ Σ 0 0 ± 54.6 0 ± 54.6 0 ± 54.6 B 0 s nλ 0 ± 31.0 (71.6 +19.3 5.4 ) (94.9+5.1 16.3 ) B Σ Λ 0 ± 35.7 0 ± 35.7 0 ± 35.7 B 0 s Σ 0 Ξ 0 0 ± 15.4 (4.8 +14.1 13.9 ) (58.3+1.9 1.6 ) B Ξ Ξ 0 0 ± 35.7 0 ± 35.7 0 ± 35.7 B 0 s Σ Ξ 0 ± 30.0 0 ± 30.0 0 ± 30.0 B ΛΣ + 0 ± 70.7 0 ± 70.7 0 ± 70.7 B 0 s ΛΞ 0 0 ± 100 0 ± 100 0 ± 100 B 0 pp 0 ± 6.9 (36.0 +6.8.0 ) (49.3+5. 18.3 ) B0 Σ + Σ + 100 100 100 100 100 100 B 0 nn 0 ± 31.6 (59.3 7.9 +3.5 ) (79.5+16.0 3. ) B0 Σ 0 Σ 0 0 ± 33.6 (59.3 +4.8 9.6 ) (79.5+16.7 4.4 ) B 0 Ξ 0 Ξ 0 100 100 100 100 100 100 B 0 Σ Σ 0 ± 47.1 0 ± 47.1 0 ± 47.1 B 0 Ξ Ξ 0 ± 3.1 0 ± 3.1 0 ± 3.1 B 0 Σ 0 Λ 0 ± 13.3 (36.0 +1.8 1.1 ) (49.3+1.1 11. ) B 0 ΛΛ 100 100 100 100 100 100 B 0 ΛΣ 0 0 ± 70.7 0 ± 70.7 0 ± 70.7 asymptotic amplitudes and Eq. (16) are used. We show results for φ = 0 ±π/4 and ±π/. The uncertainties are from relaxing the asymptotic relation by using Eq. (11) and varying the strong phases of r ( ) ti r ( ) pi r ( ) ewpi and from sub-leading terms Eq. (13) with strong phases from η ijk. Note that to satisfy the experimental B 0 pp and B Λp rates the sizes of r ( ) are reduced by 60% in all B q BB modes. Now we discuss the CP asymmetries of Group I II III modes according to their rates and detectability as noted in Table IX. For the Group I mode the CP asymmetry of the B 0 pp decay can be as large as 49%. For the Group II modes A(B 0 Σ 0 Λ) and A(B 0 s pσ + ) are similar to A(B 0 pp). For the Group III modes the CP asymmetries of B Σ 0 Σ + and B 0 Σ 0 Σ 0 decays are similar and can reach 80% while the CP asymmetry of the B 0 s Σ 0 Ξ 0 decay is smaller then these two modes but can still reach 58%. The CP asymmetries of these modes basically all have the same sign when φ is large enough. As noted in the previous subsection B Σ Σ 0 B Σ Λ B Ξ Ξ 0 B 0 Ξ Ξ B 0 s Σ Ξ and B 0 Σ Σ decays do not have any tree (T ibb ) contribution. 3

TABLE XI: Same as Table X but with S = 1 B q BB modes. Mode φ = 0 φ = ±π/4 φ = ±π/ Mode φ = 0 φ = ±π/4 φ = ±π/ B Σ 0 p 0 ± 0.9 (8.9 18.6 +.1 ) (45.0+3.0 16.7 ) B0 Σ + p 0 ± 16.4 ±(7.1 +16.6 13.6 ) ±(35.3+15.5 11.9 ) B Σ n 0 ± 1.8 0 ± 1.8 0 ± 1.8 B 0 Σ 0 n 0 ± 9.0 ±(47.8 +3.8 3.1 ) ±(58.6+19.5 17.8 ) B Ξ 0 Σ + 0 ±.6 ±(5.8 +3.4.3 ) ±(8.1+3.9.8 ) B0 Ξ 0 Σ 0 0 ±.3 ±(5.8 +3.0.1 ) ±(8.1+3.6.5 ) B Ξ Σ 0 0 ± 1.8 0 ± 1.8 0 ± 1.8 B 0 Ξ 0 Λ 0 ± 4.5 ±(7.1 +8.3 16. ) ±(35.3+7. 13.0 ) B Ξ Λ 0 ± 5.5 0 ± 5.5 0 ± 5.5 B 0 Ξ Σ 0 ± 1.5 0 ± 1.5 0 ± 1.5 B Λp 0 ± 4.8 ±(9.6 4.0 +5.9 ) ±(13.+6.5 4.3 ) B0 Λn 0 ± 8.4 ±(18.7 +9.6 6.9 ) ±(4.9+9.8 6.6 ) B 0 s pp 100 100 100 100 100 100 B 0 s Σ + Σ + 0 ± 1.0 ±(7.1 +1.6 16.7 ) ±(35.3+0.1 14.1 ) B 0 s nn 100 100 100 100 100 100 B 0 s Σ 0 Σ 0 0 ± 11.1 ±(14. +1.5 8.9 ) ±(19.+1.7 8. ) B 0 s Ξ 0 Ξ 0 0 ± 8.4 ±(14. +10.4 6.5 ) ±(19.+11.1 6.3 ) B0 s Σ Σ 0 ± 1.6 0 ± 1.6 0 ± 1.6 B 0 s Ξ Ξ 0 ±.5 0 ±.5 0 ±.5 B 0 s Σ 0 Λ 0 ± 46.4 ±(74.3 +1.7 33.6 ) ±(84.7+13.8 19.7 ) B 0 s ΛΛ 0 ± 7.0 ±(14. 5.7 +8.3 ) ±(19.+8.8 5.6 ) B0 s ΛΣ 0 0 ± 4.8 ±(74.3 +19.7 33. ) ±(84.7+1.5.9 ) Although the B 0 Ξ Ξ B 0 Σ Σ and B 0 s Σ Ξ decays are pure penguins modes which only have P ibb P iew BB and P A BB contributions it is still possible for these modes to have sizable CP asymmetries. Since the sizes of u-penguin (P u ) and c-penguin (P c ) in S = 0 modes are not very different as their ratio can be estimated as P u P c V ub Vud V cb V 0.38 (3) and they can have different strong phases to produce CP asymmetries. 5 For subleading modes B 0 Ξ 0 Ξ 0 and B 0 Σ + Σ + decays which only have E ibb and P A BB contributions their CP asymmetries can be any value. Note that B ΛΣ + B 0 s ΛΞ 0 B 0 ΛΛ and B 0 ΛΣ 0 decay with tree amplitudes T ibb canceled in the asymptotic limits have have large uncertainties on CP asymmetries which mostly come from the corrections to the asymptotic relations. Measuring CP asymmetries of these modes can give information on the corrections to the asymptotic relations. In Table XI we give results of direct CP asymmetries of S = 1 B q BB modes. The CP asymmetries of the Group I modes are as following. The CP asymmetries of the cd 5 Similarly the sizable CP asymmetry of a pure penguin mode A(B 0 K 0 K 0 ) = 16.7 +4.7+4.5+1.5+4.6 3.7 5.1 1.7 3.6 % as predicted in a QCD factorization (QCDF) calculation [35] can be understood. 4

B Λp and B 0 s ΛΛ decays are similar reaching ±13% and ±19% respectively and their signs are opposite to the sign of A(B 0 pp). The CP asymmetries of B 0 s Ξ Ξ and B Ξ Λ decays are vanishingly small. We will return to these modes later. From the table we see that for the Group II modes A(B Ξ Σ 0 ) is vanishingly while A(B 0 Ξ 0 Λ) and A(B 0 Σ + p) are similar and can be as large as ±35%. For the Group III modes A(B 0 s Σ + Σ + ) is the largest one reaching ±35% A(B 0 s Ξ 0 Ξ 0 ) and A(B 0 s Σ 0 Σ 0 ) are similar and can be as large as ±19% while A(B Ξ 0 Σ + ) and A(B 0 Ξ 0 Σ 0 ) are similar and are not sizable but can reach ±8% Some of the above modes have rates of orders 10 7 (see Table II) and with unsuppressed CP asymmetries. From Table XI we see that the CP asymmetries of B Σ n B Ξ Σ 0 B Ξ Λ B 0 Ξ Σ B 0 s Σ Σ and B 0 s Ξ Ξ decays have vanishing central values as they do not have any tree (T ) contribution. In particular ibb B0 Ξ Σ B 0 s Σ Σ and B 0 s Ξ Ξ decays are pure penguin modes which only have P P and P ibb iew BB A BB terms. Their CP asymmetries are small. This can be understood as P u is much smaller than P c. We can estimate the CP asymmetry as following: P u A P c sin γ sin δ < V ub V us V cb Vcs sin γ 3.8% (4) where δ is the strong phase different of the penguins. 6 The smallness of the direct CP reflects the fact that in S = 1 decays the c-penguin is much larger than the u-penguin as their CKM factors are off by about a factor of 50. They can be tests of the Standard Model. In particular the B 0 s Ξ Ξ decay being a Group I mode can cascadely decays to all charged final states and with unsuppressed rate ( 10 7 see Table II). The CP asymmetry of this mode can be searched for. Nevertheless the search is quite demanding as it requires tagging of B s and suffers from low efficiency of Ξ recontruction. Subleading modes B 0 s pp and B 0 s nn decays with E ibb and P A BB value on their CP asymmetries. contributions can have any There are relations on direct CP asymmetries between S = 0 and S = 1 modes by using the so-called U-spin symmetry [36 37]. With CP (B q f) Γ(B q f) Γ(B q f) = τ(b q ) A(B q f)b(b q f) (5) 6 It is useful to compare to the CP asymmetry of the S = 1 pure penguin P P mode mode: A(B 0 s K 0 K 0 ) = 0.9...3 +0.+0.+0.1+0. % in a QCDF calculation [35]. 5