Final state interactions in hadronic B decays

1. Final state interactions in hadronic B decays • Hai-Yang Cheng • Academia Sinica • FSIs • BRs & CPV in B decays • Polarization anomaly in BK* QCD & Hadronic Physics, Beijing, June 16-20, 2005

2. Importance of FSI in charm decays has long been recognized • some nearby resonances exist at energies  mD • charm is not very heavy • General folklore for B decays: • FSI plays a minor role due to large energy release in B decays • and the absence of nearby resonances • There are growing hints at some possible soft final-state • rescattering effects in B physics

3. 1. Direct CP violation B f ei(+) : strong phase : weak phase One needs at least two different B  f paths with distinct weak & strong phases _ first confirmed DCPV (5.7) in B decays (2004) _ _ Recall in kaon decays

4. Comparison with theory: pQCD & QCDF input • pQCD (Keum, Li, Sanda): A sizable strong phase from penguin-induced annihilation by introducing parton’s transverse momentum • QCD factorization (Beneke, Buchalla, Neubert, Sachrajda): • Because of endpoint divergences, QCD/mb power corrections due to annihilation and twist-3 spectator interactions can only be modelled • QCDF (S4 scenario): large annihilation with phase chosen so that a correct sign of A(K-+) is produced (A=1, A= -55 for PP, A=-20 for PV and A=-70 for VP)

5. Need sizable strong phases to explain the observed direct CPV • SD perturbative strong phases: penguin (BSS) vertex corrections (BBNS) annihilation (pQCD) • Nonperturbative LD strong phases induced from power correctionsespecially fromfinal-state rescattering If intermediate states are CKM more favored than final states, e.g. BDDsK  large phases  large corrections to rate strong weak

6. 2. Some color-suppressed or factorization-forbidden or penguin-dominated modes cannot be accommodatedin the naïve factorization approach • Some decay modes do not receive factorizable contributions e.g. B  Kc0 with sizable BR though c0|c(1-5)c|0=0. • Color-suppressed modes e.g. B0  D0 h0 (h0=0,,0,,’), 00, 00 have the measured rates larger than theoretical expectations. • Penguin-dominated modes such as BK*, K, K, K* predicted by QCDF are consistently lower than experiment by a factor of 2  3  importance of power corrections (inverse powers of mb) e.g. FSI, annihilation, EW penguin, New Physics, …

7. Our goal is to study FSI effect on BRs and CPV (direct & indirect) in B decays (Polarization anomaly in BK* will be briefly mentioned) LD rescattering can be incorporated in any SD approach but it requires modelling of 1/mb power corrections We would provide a specific model for FSI to compute strong phases so that we can predict (rather than accommodate) the sign and magnitude of direct CP asymmetries

8. Approaches for FSIs in charmless B decays • Regge approach [Donoghue,Golowich,Petrov,Soares] FSI phase is dominated by inelastic scattering and doesn’t vanish even in mb limit • QCDF [Beneke,Buchalla,Neubert,Sachrajda] strong phase is O(s, /mb): systematic cancellation of FSIs in mb • Charming penguin [Ciuchini et al.] [Colangelo et al.] [Isola et al.] long distance in nature, sources of strong phases, supported by SCET • One-particle-exchange model for LD rescattering • has been applied to charm and B decays [Du et al.][Lu,Zou,..] • Quasi elastic scattering model[Chua,Hou,Yang] Consider MMMM (M: octet meson) rescattering in BPP decays • Diagrammatic approach [Chiang, Gronau, Rosner et al.] …

9. Diagrammatic Approach All two-body hadronic decays of heavy mesons can be expressed in terms of six distinct quark diagrams [Chau, HYC(86)] (tree) (color-suppressed) (exchange) (vertical W loop) (penguin) (annihilation) All quark graphs are topological and meant to have all strong interactions included and hence they are not Feynman graphs. And SU(3) flavor symmetry is assumed.

10. Global fit to B, K data (BRs & DCPV) based on topological diagrammatic approach yields [Chiang et al.] consistent with that determined from B D decays

11. quark annihilation possible FSIs meson annihilation quark exchange relevant for e.g. B0 B0D00  W exchange  Color suppressed C At hadron level, FSIs manifest as resonant s-channel & OPE t-channel graphs

12. FSI as rescattering of intermediate two-body states [HYC, Chua, Soni] • FSIs via resonances are assumed to be suppressed in B decays due to the lack of resonances at energies close to B mass. • FSI is assumed to be dominated by rescattering of two-body intermediate states with one particle exchange in t-channel. Its absorptive part is computed via optical theorem: • Strong coupling is fixed on shell. For intermediate heavy mesons, • apply HQET+ChPT • Form factor or cutoff must be introduced as exchanged particle is • off-shell and final states are necessarily hard • Alternative: Regge trajectory [Nardulli,Pham][Falk et al.] [Du et al.] …

13. Form factor is introduced to render perturbative calculation meaningful  LD amp. vanishes in HQ limit •  = mexc + rQCD (r: of order unity) •  or r is determined by a fit to the measured rates •  r is process dependent • n=1 (monopole behavior), consistent with QCD sum rules Once cutoff is fixed  direct CPV can be predicted Dispersive part is obtained from the absorptive amplitude via dispersion relation subject to large uncertainties and will be ignored in the present work

14. Inputs Form factors: covariant light-front approach: relativistic QM for s-wave to s-wave and p-wave transitions (HYC,Chua,Hwang 2004) SD approach: QCD factorization (default scenario) with A=H=0 in double counting problem is circumvented

15. Theoretical uncertainties (SD) • variation of CKM parameters • =(6315) • quark masses: ms(2 GeV)= 9020 MeV • renormalization scale: from =2mb to mb/2 • heavy-to-light form factors: e.g. FB(0)=0.250.03 • meson distribution amplitudes

16. Theoretical uncertainties (LD) • 1. Model assumption • multi-body contributions • form-factor cutoff: • i). n=1 • ii).  = mexc + rQCD (15% error assigned for QCD) • rD=2.1, 1.6, 0.73, 0.67, respectively, for D, , , K modes • varies for penguin-dominated PV modes • dispersive contribution • 2. Input parameters • strong couplings of heavy mesons and their SU(3) breaking • g(D*D)=17.90.31.9 (CLEO) • heavy-to-heavy form factors

17. Br (10-6) • first error: SD, second error: LD • LD uncertainties are comparable to SD ones & SD errors are affected only slightly by FSIs. • No reliable estimate of LD rescattering effects for f0KS

18. B  K All rescattering diagrams contribute to penguin topology, dominated by charm intermediate states fit to rates  rD = rD*  0.67  predict direct CPV

19. _ _ _ _ • For simplicity only LD uncertainties are shown here • FSI yields correct sign and magnitude for A(+K-) ! • K anomaly: A(0K-)  A(+ K-), while experimentally they differ by 3.4 [See Fleischer’s talk]

20. B   B  ﹣ _ _ _ • Sign and magnitude for A(+-) are nicely predicted ! • DCPVs are sensitive to FSIs, but BRs are not (rD=1.6) • For 00, 1.40.7 BaBar Br(10-6)= 5.11.8 Belle 1.6+2.2-1.6 CLEO Discrepancy between BaBar and Belle should be clarified.

21. B  K* B  _ _ 6.92.4 _ _ BaBar K.F. Chen (CKM2005): BaBar, hep-ex/0504009  Br(B-0K*-)=(6.92.4)10-6 For 0K*0, Br(10-6)= 3.01.0 BaBar 0.4+1.9-1.7 Belle

22. Comparison with other approaches All known existing models fit the data of BRs and DCPV and then make predictions for mixing (indirect) CPV. e.g. 1. charming penguin (Ciuchini et al. and many others) Consider charmless B decay BKwith BDDsK  charming penguin is CKM doubly enhanced & gives dominant LD corrections a fit result CKM2005  S(0K0)=0.770.04

23. 2. Fit QCDF to data  fix unknown power correction • parameters A, H, A, H • Aleksan et al. (hep-ph/0301165) • S4 scenario of Beneke & Neubert (hep-ph/0308092) • Leitner, Guo, Thomas (hep-ph/0411392) • Cottingham et al. (hep-ph/0501040)

24. Mixing-induced CP violation [HYC,Chua,Soni] It is expected in SM that -fSf  sin2 = 0.726  0.037 with deviation at most O(0.1) in B0KS, KS, 0KS, ’KS, 0KS, f0KS, K+K-KS, KSKSKS [London,Soni; Grossman, Gronau, Ligeti, Nir, Rosner, Quinn,…

25. A current hot topic G. Kane (and others): The 2.7-3.7 anomaly seen in b→s penguin modes is the strongest hint of New Physics that has been searched in past many many years… It is extremely important to examine how much of the deviation is allowed in the SM and estimate the theoretical uncertainties as best as we can.

26. Time-dependent CP asymmetries: In general, Sf sin2eff sin(2+W). For bsqq modes, Since au is larger than ac, it is possible that S will be subject to significant “tree pollution”. However, au here is color-suppressed. • Penguin contributions to KS and 0KS are suppressed due to cancellation between two penguin terms (a4 & a6) • relative importance of tree contribution • large deviation of S from sin2

27. S(KS)>sin2, S(0KS)<sin2 • FSI can bring in additional weak phase via K*,K intermediate states (even when tree is absent at SD) (see also Beneke)

28. FSI effect is tiny due to small source (K*,K) amplitudes (Br~10-6) compared to Ds*D (Br~10-2,-3). It tends to alleviate the deviation from sin2 • For 0KS, S=S-sin2<0 at SD but it becomes positive after including FSI. • Sf is positive and less than 0.1 in SM, while experimentally Sf is always negative

29. Effective sin2 in K+K-KS & KSKSKS For K+K-KS, S= -(2f+-1)sin2eff (f+: CP-even fraction) For KSKSKS, S= -sin2eff • K+K-KS is subject to large tree pollution from color-allowed tree diagrams • KSKSKS is very clean for testing SM

30. Polarization anomaly in B  K* [HYC, Chua, Soni] Short-distance induced transverse polarization in B V1V2 (V: light vector meson) is expected to be suppressed Confirmed for B with fL 0.97 but for BK* fL 0.50, f|| 0.25, f 0.25 • Get large transverse polarization from B Ds*D* and then convey it to K* via FSI fL(Ds*D*)  0.51 contributes to f only f|| 0.41, f 0.08

31. Large cancellation occurs in B{Ds*D,DsD*}K* processes. This can be understood as CP & SU(3) symmetry  0 ! ＋  very small perpendicular polarization, f 2%, in sharp contrast to f 15% obtained by Colangelo, De FArzio, Pham While fT 0.50 is achieved, why is f not so small ? Cancellation in B{VP,PV}K* can be circumvented in B{SA,AS}K*. For S,A=D**,Ds**  f 0.22 It is very easy to explain why fL 0.50 by FSI, but it takes some efforts to understand why f f||

32. Conclusions • DCPV in charmless B decays is significantly affected by LD rescattering. Correct sign and right magnitude of DCPV in K-+ and +- are obtained after including FSI. • For penguin-dominated MKS modes, FSI tends to alleviate the deviation from sin2. • Large transverse polarization fT 0.50 can be obtained from final-state rescattering of B Ds*D* K*

33. Collinear vs. kT factorization • The subleading amplitudes in QCDF develop • end-point singularities in twist-3 nonspectator and in annihilation • An end-point singularity means breakdown • of simple collinear factorization • Use more conservative kT factorization • Include parton kT to smear the singularity Perturbative QCD approach [Keum, Li, Sanda; Lu, Yang, Ukai]

34. kT factorization Sudakov factors S describe the parton distribution in kT Parton-level diagrams Bound-state distribution amplitude KT accumulates after infinitely many gluon exchanges Similar to the DGLAP evolution up to kT~Q

35. Scales and penguin enhancement In PQCD this gluon is off-shell by Slow parton Fast parton For penguin-dominated modes, the branching ratios Wilson coefficients PQCD QCDF

36. Recent progress on PQCD • Nonfactorizable contributions are important for color-suppressed modes---explained B! D00, (J/,c0,c1,c )K (*) branching ratios, helicity amplitudes (Keum, Kurimoto et al.; Chen, Li). • Annihilation lowers longitudinal polarization in B! VV.Also predicted pure-annihilation modes, which cannot be done in FA (Lu et al.). • Predicted CP asymmetry, isospin breaking of B! K* (Matsumori et al.). • NLO PQCD enhances C, and resolves B! K puzzle (Li, Mishima, Sanda): LO NLO Data Acp(K+-) -0.13 -0.11 -0.1090.019 Acp(K+0) -0.09 +0.03 0.040.04 Annihilation generates large strong phase, which explains direct CP asymmetries.

37. Baryonic B decays 3-body baryonic B deacys were found to have larger BRs than 2-body decays There are extensive studies of baryonic B decays in Taiwan both experimentally and theoretically Theory Expt. Belle group at NTU Chua, Geng, Hou, Hsiao, Tsai, Yang, HYC,… B-→ppK-: first evidence of charmless baryonic B decay B→pp(K,K*,) →p(,K) →K B→pp, , p (stringent limits) Publication after 2001: (hep-ph) 0008079, 0107110, 0108068, 0110263, 0112245, 0112294, 0204185, 0204186, 0208185, 0210275, 0211240, 0302110,0303079, 0306092, 0307307, 0311035, 0503264 C.Q. Geng, this afternoon B→p: first evidence of b→s penguin in baryonic B decays 0201015, 0405283 first paper on radiative baryonic B decays

38. Back-up slides

39. Regge approach • In evaluating absorptive part via • replace Feynman-diagram strong scattering amplitude T by Regge • amplitude R(s,t): (t): residual function, linear Regge trajectory (t)=0+’t, intercept 0=0.45 for , K*,…, -1.8 for D, D*  suppression of FSI with increasing energy s  suppression of charming penguin relative to the light Regge exchanges • uncertainties: t-dep. of (t), BR sensitive to unknown ’ R(s,t) is valid at large s and t  0

40. Long-distance contributions to B  D 

41. Even if short-distance W-exchange vanishes (i.e. ESD=0), final-state rescattering does contribute to weak annihilation B0Ds+K- proceeds via W-exchange

42. B Cutoff scale is fixed by B K via SU(3) symmetry  too large +- ( 910-6) and too small 00 (0.410-6) An additional rescattering contribution unique to  but not available to K is needed to suppress +- and enhance 00 _ + D+(+) B0DD() has the same topology as vertical W-loop diagram V D-(-) -

43. Need to fit to rates and CPV of +- simultaneously • Charming penguin alone doesn’t suffice to explain 00 rate (rD=0.67) • Sign of A(00) can be used to discriminate between different models • W-exchange can receive LD contributions from FSI • Define Teff=T+E+V, Ceff=C-E-V  Ceff/Teff=(0.900.02) exp[-i(882)]