Neutral B mixing The Standard Model and Beyond E. Freeland, C. - - PowerPoint PPT Presentation

neutral b mixing the standard model and beyond
SMART_READER_LITE
LIVE PREVIEW

Neutral B mixing The Standard Model and Beyond E. Freeland, C. - - PowerPoint PPT Presentation

Neutral B mixing The Standard Model and Beyond E. Freeland, C. Bouchard, C. Bernard, A.X. El-Khadra, E. Gamiz, A.S. Kronfeld, J. Laiho, and R.S. Van de Water for the Fermilab Lattice and MILC Collaborations Neutral


slide-1
SLIDE 1

Neutral B mixing The Standard Model and Beyond

  • E. Freeland, C. Bouchard, C. Bernard, A.X. El-Khadra, E. Gamiz,

A.S. Kronfeld, J. Laiho, and R.S. Van de Water for the Fermilab Lattice and MILC Collaborations

slide-2
SLIDE 2

Neutral B mixing

In the Standard Model, mixing is suppressed loop process, diagram with dominates, but is Cabibbo suppressed B0

q

⇒ “Relatively” easy for new physics to cause observable effects.

mt

ɤ ɤ ɤ ɤ ɤ ɤ ɤ ɤ ɤ ɤ

B

q

B0

q

q b W ± W ∓

Vtq

u, c, t

⌇ ⌇ ⌇ ⌇ ⌇ ⌇ ɤ ɤ ɤ ɤ ɤ ɤ

q b

V ∗

tq

u, c, t

slide-3
SLIDE 3

Matrix Elements

⌇ ⌇ ⌇ ⌇

W ± W ∓ t, c, u t, c, u

B0

q

B

q

B0

q

B

q

Operators are

O1 = (¯ bαγµLqα) (¯ bβγµLqβ) O2 = (¯ bαLqα) (¯ bβLqβ) O3 = (¯ bαLqβ) (¯ bβLqα) O4 = (¯ bαLqα) (¯ bβRqβ) O5 = (¯ bαLqβ) (¯ bβRqα)

SM BSM

Oi

Heff =

5

  • i=1

CiOi

slide-4
SLIDE 4

¯ B0

q|Oi(µ)|B0 q ∝ f 2 BqBi(µ)

Matrix Elements

⌇ ⌇ ⌇ ⌇

W ± W ∓ t, c, u t, c, u

B0

q

B

q

B0

q

B

q

Operators are

O1 = (¯ bαγµLqα) (¯ bβγµLqβ) O2 = (¯ bαLqα) (¯ bβLqβ) O3 = (¯ bαLqβ) (¯ bβLqα) O4 = (¯ bαLqα) (¯ bβRqβ) O5 = (¯ bαLqβ) (¯ bβRqα)

Common parametrization

Oi

f 2

Bq

Bi(µ) Heff =

5

  • i=1

CiOi

slide-5
SLIDE 5

Experiment and SM: ∆Mq

“Tension” in the CKM matrix.

Lenz et al., arXiv: 1203:0238; Laiho et al.

  • PhysRevD. 81, 034503, and end-of-2011 update

Our ability to constrain , is limited by . |VtbV ∗

tq|2

¯ B0

q|O1(µ)|B0 q

known want need from lattice < 1%

∆Mq = G2

F M 2 W S0

4π2

  • ηB(µ)

|VtbV ∗

tq|2

¯ B0

q|O1(µ)|B0 q

experiment

slide-6
SLIDE 6

Experiment and SM: ∆Mq

∆Ms ∆Md =

  • Vts

Vtd

  • 2 ¯

B0

s|O1(µ)|B0 s

¯ B0

d|O1(µ)|B0 d ≡

  • Vts

Vtd

  • 2 MBs

MBd ξ2

experiment lattice want

slide-7
SLIDE 7

Experiment and SM: ∆Mq

SU(3)-breaking ratio

∆Ms ∆Md =

  • Vts

Vtd

  • 2 ¯

B0

s|O1(µ)|B0 s

¯ B0

d|O1(µ)|B0 d ≡

  • Vts

Vtd

  • 2 MBs

MBd ξ2

  • Some (lattice) errors cancel in the ratio of matrix elements,
  • In CKM matrix fits, use of can aid in minimizing correlations

between lattice inputs. ξ want lattice experiment

slide-8
SLIDE 8

Recent experimental results are putting focus on .

∆Γ

(Lenz et al., arXiv:1203.0238; Haisch, Moriond 2012 )

Experiment and SM: ∆Γq

Lenz, Nierste JHEP 0706:072, 2007 hep-ph/0612167 Beneke, Buchalla, Dunietz, PRD 54:4419, 1996, Erratum-ibid.D 83 119902 (2011); hep-ph/9605259v1

dominates also needed

∆Γq =

  • G1 ¯

B0

q|O1(µ)|B0 q + G3 ¯

B0

q|O3(µ)|B0 q

  • cos φq + O(1/mb, αs)
slide-9
SLIDE 9

Recent experimental results are putting focus on .

∆Γ

(Lenz et al., arXiv:1203.0238; Haisch, Moriond 2012 )

Experiment and SM: ∆Γq

∆Γq =

  • G1 ¯

B0

q|O1(µ)|B0 q + G3 ¯

B0

q|O3(µ)|B0 q

  • cos φq + O(1/mb, αs)

OR ≡ O2 + O3 + (1/2)O1

useful for estimating errors.

1/mb

yields

∆Γ/∆M O3/O1

slide-10
SLIDE 10

experiment

Including BSM contributions, takes the generic form above. ∆Mq Lattice values of (matrix elements of) through are needed to check that a given BSM model is consistent with experiment. O1 O5

Experiment and BSM: ∆Mq

model dependent need from lattice

∆Mq =

5

  • i=1

Ci(µ) B0

q|Oi(µ)|B q

< 1%

slide-11
SLIDE 11

Status of Lattice

slide-12
SLIDE 12

Status of Lattice:

FNAL-MILC 5.0% ξ = 1.268(63)

arXiv: 1205.7013, submitted to PRD

O1

slide-13
SLIDE 13

Albertus et al., Phys.Rev.D82:014505, 2010, arXiv:1001.2023

RBC

ξ = 1.13(12)

11%

domain-wall test calculation; one, 0.11 fm, lattice spacing

HPQCD

Gamiz et al., Phys.Rev.D80:014503, 2009, arXiv:0902.1815

fBd

  • ˆ

BBd = 216(15) MeV

6.8%

fBs

  • ˆ

BBs = 266(18) MeV

6.9% 2.6%

ξ = 1.258(33)

Status of Lattice: O1

slide-14
SLIDE 14

Albertus et al., Phys.Rev.D82:014505, 2010, arXiv:1001.2023

RBC

ξ = 1.13(12)

11%

domain-wall test calculation; one, 0.11 fm, lattice spacing

HPQCD

Gamiz et al., Phys.Rev.D80:014503, 2009, arXiv:0902.1815

fBd

  • ˆ

BBd = 216(15) MeV

6.8%

fBs

  • ˆ

BBs = 266(18) MeV

6.9% 2.6%

ξ = 1.258(33)

Status of Lattice: O1

HPQCD 2009 RBC 2010 FNAL-MILC 2012

1 1.2 1.4 1.6 1.8

slide-15
SLIDE 15

FNAL-MILC: Lattice 11 Proceedings (Dec 2011), arXiv:1112:5642

fBq

  • ˆ

BBq 9 to 6% on O1

Two ensembles:

  • E. Dalgic et al., PRD76:011501, 2007

Quenched: Becirevic et al., JEHP 0204 (2002) 0250

O1...5

Status of Lattice:

ETMC also working on this.

O2,3

Preliminary

Estimated

slide-16
SLIDE 16

Details of Our Calculation

slide-17
SLIDE 17

Previous vs current analysis

Ensembles more ensembles higher statistics smaller lattice spacing smaller light-quark mass Results full set of matrix elements bag parameters in conjunction with analysis able to do all ratios and combinations Use complete ChiPT expression This alone improves the error on from 5.0% to 3.8%. ξ fB

(Ethan Neil’s talk)

slide-18
SLIDE 18

Analysis Overview

  • Generate two- and three-points correlator data.
  • Fit 2pt+3pt correlators simultaneously for each meson.
  • Renormalize & match the matrix elements.
  • Do a chiral and continuum extrapolation for each matrix element.

B

q

B0

q

B0

q

B0

q

slide-19
SLIDE 19

Actions and Ensembles

gauge configurations

physical point

heavy valence quark

  • improved Wilson action
  • Fermilab interpretation

light valence quark

  • staggered action
  • mass from > to 0.05 .

ms ms

MILC (asqtad) gauge configurations

  • 2+1 asqtad sea quarks,
  • tadpole improved gluons
  • from 0.4 to 0.05

ml/ms

~2000 configurations ~500 configurations analyzed ~1000 configurations partially analyzed

0.02 0.04 0.06 0.08 0.1 0.12 0.14

lattice spacing in fm

0.1 0.2 0.3 0.4 0.5

ml / ms

arXiv: 1205.7013 ξ

slide-20
SLIDE 20

Actions and Ensembles

gauge configurations

physical point ~2000 configurations ~500 configurations analyzed ~1000 configurations partially analyzed

0.02 0.04 0.06 0.08 0.1 0.12 0.14

lattice spacing in fm

0.1 0.2 0.3 0.4 0.5

ml / ms

Lat11

heavy valence quark

  • improved Wilson action
  • Fermilab interpretation

MILC (asqtad) gauge configurations

  • 2+1 asqtad sea quarks,
  • tadpole improved gluons
  • from 0.4 to 0.05

ml/ms

light valence quark

  • staggered action
  • mass from > to 0.05 .

ms ms

slide-21
SLIDE 21

Actions and Ensembles

gauge configurations

physical point ~2000 configurations ~500 configurations analyzed ~1000 configurations partially analyzed

heavy valence quark

  • improved Wilson action
  • Fermilab interpretation

MILC (asqtad) gauge configurations

  • 2+1 asqtad sea quarks,
  • tadpole improved gluons
  • from 0.4 to 0.05

ml/ms

0.02 0.04 0.06 0.08 0.1 0.12 0.14

lattice spacing in fm

0.1 0.2 0.3 0.4 0.5

ml / ms

Lat12

light valence quark

  • staggered action
  • mass from > to 0.05 .

ms ms

slide-22
SLIDE 22

Three-points

✦ simultaneous fit of two-point + three-point; ✦ constrains energies and two-point amplitudes ✦ use constrained curve fitting

+Zp

mZp nOpp mn(−1)t1+t2e−Ep

mt1e−Ep nt2

  • +Zp

mZnOp mn(−1)t1e−Ep

mt1e−Ent2

+ZmZp

nOp mn(−1)t1e−Emt1e−Ep

nt2

C3pt(t1, t2) =

  • m,n
  • ZmZnOmne−Emt1e−Ent2

C2pt(t) =

  • m
  • Z2

me−Emt + (−1)(t+1)(Zp m)2e−Ep

mt

slide-23
SLIDE 23

Operators mix under renormalization (even in the continuum). E.g.

Renormalization and Matching

ζij are calculated using 1-loop perturbation theory.

αs = αv(2/a)

We use the “V” scheme as implemented by Q. Mason et al. with 4-loop running.

  • Q. Mason et al. [HPQCD Collaboration],
  • Phys. Rev. Lett. 95, 052002 (2005) hep-lat/0503005
  • T. van Ritbergen et al., Phys. Lett. B 400, 379 (1997) hep-ph/9701390

O1R = (1 + αsζ11)O1 + αsζ12O2

slide-24
SLIDE 24

Data O2

0.5 1 1.5 2

(r1 mπ )

2

  • 1.5
  • 1.25
  • 1
  • 0.75
  • 0.5

r1

3<O2> / MB

slide-25
SLIDE 25

Data O3

0.5 1 1.5 2

(r1 mπ )

2

0.05 0.1 0.15 0.2

r1

3<O3> / MB

slide-26
SLIDE 26

Status: ChiPT

continuum, PQ: Detmold and Lin, aXiv:0612028, hep-lat, 2006

We use SU(3), partially-quenched, heavy-meson, staggered ChiPT With staggered light quarks, matrix elements of wrong-spin

  • perators appear in the ChiPT.

Because the five matrix elements form a complete basis, wrong-spin contributions can be written in terms of them. O1...5

✦ Mixing occurs: ✦ No new LEC’s are introduced.

O1 O2 O3 O4 O5 and

Claude Bernard’s talk

slide-27
SLIDE 27

Status: ChiPT

We will do a simultaneous fit for each set of mixed operators.

✦ Mixing occurs. ✦ No new LEC’s are introduced.

E.g. B

q|Oq 1|B0 q = β1

  • 1 +

Wqb + Wbq 2 + Tq + Qq + ˜ T (a)

q

+ ˜ Q(a)

q

  • +(2β2 + 2β3) ˜

T (b)

q

+ (2β′

2 + 2β′ 3) ˜

Q(b)

q

+analytic terms O1

wrong-spin contributions

slide-28
SLIDE 28

Status of Lattice:

FNAL-MILC 5.0% ξ = 1.268(63)

arXiv: 1205.7013, submitted to PRD

O1

slide-29
SLIDE 29

✦ nearly done with three-point analysis ✦ fourteen ensembles across four lattice spacings ✦ corrected chiral form exists; extrapolation remains to be done ✦ results will include: ✦ complete set (5) of matrix elements and bag parameters ✦ ratios and and the combination ✦ “half-way point” (Lattice 11, arXiv:1112:5642): ✦ error on 9 to 6% (perturbation theory, continuum-ChiPT extrap.) ✦ Errors will be smaller for full-data set analysis.

ξ O3/O1 OR

fBq

  • ˆ

BBq

Conclusion

✦ arXiv: 1205.7013

ξ = 1.268(63)

F i n i s h e d I n p r

  • g

r e s s

slide-30
SLIDE 30

Backup Slides

slide-31
SLIDE 31

FNAL-MILC: Lattice 11 Proceedings (Dec 2011), arXiv:1112:5642

fBq

  • ˆ

BBq 9 to 6% on O1

O1...5

Status of Lattice:

slide-32
SLIDE 32

Data O1

0.5 1 1.5 2

(r1 mπ )

2

0.4 0.5 0.6 0.7 0.8 0.9 1 1.1

r1

3<O1>

slide-33
SLIDE 33

Data O4

0.05 0.1 0.15

r1 mq

1 1.5 2 2.5

r1

3<O4> / MB

slide-34
SLIDE 34

0.05 0.1 0.15

r1 mq

0.4 0.6 0.8 1

r1

3<O5> / MB

Data O5

slide-35
SLIDE 35

ChiPT systematic

0.5 1 1.5 2

(r1 mπ )

2

0.2 0.4 0.6 0.8 1 1.2

r1

3<O1>

coarse (0.005, 0.050) coarse (0.007, 0.050) coarse (0.010, 0.050) coarse (0.020, 0.050) fine (0.0031, 0.031) fine (0.00465, 0.031) coarse(0.005, 0.050) coarse (0.007, 0.050) coarse (0.010, 0.050) coarse (0.020 0.050) fine (0.0031, 0.031) fine (0.00465, 0.031) chiral-continuum extrapolation extrap from NO wrong-spin fit

<O1> versus (r1mπ)

2 Wrong Spin Included

4/13/12 NNLO7 fit to 4coarse+2fine Q = 1.0 (chisq/dof=0.16)

slide-36
SLIDE 36

PDG, J.Phys G37, 1 (2010) CDF, PRL 97, 242003 (2006)

∆Md = 0.507 ± 0.003(stat) ± 0.003(sys) ps−1 ∆Ms = 17.77 ± 0.10(stat) ± 0.07(sys) ps−1

Experiment and Lattice: ∆Mq

(LHCb-CONF-2011-050: 17.73(5))

“Tension” in the CKM matrix.

Lenz et al., arXiv: 1203:0238; Laiho et al.

  • PhysRevD. 81, 034503, and end-of-2011 update

Our ability to constrain , is limited by . |VtbV ∗

tq|2

¯ B0

q|O1(µ)|B0 q

known want need from lattice < 1%

∆Mq = G2

F M 2 W S0

4π2

  • ηB(µ)

|VtbV ∗

tq|2

¯ B0

q|O1(µ)|B0 q

slide-37
SLIDE 37

∆Γd Γd = 0.010 ± 0.037

HFAG/PDG 2011; BaBar, DELPHI

∆Γs = (0.163 ± 0.065) ps−1

D0 8 fb-1, arXiv 1109.3166

LHCb 1 fb-1, Moriond 2012 (0.37 fb-1 arXiv:1112.3183)

∆Γs = (0.116 ± 0.019) ps−1

‡Specifically, NP in versus .

M12 Γ12

∆Γs ∆Ms A scenario with new physics in yields a better fit to current data than a scenario with new physics in .‡ (Haisch, Moriond 2012 )

Lenz et al., arXiv:1203.0238

“We introduce a fourth scenario with NP in both and , which can accommodate all data.” Γd,s

12

M d,s

12

Experiment and Lattice: ∆Γq

Recent LHCb results are putting focus on ! ∆Γ

slide-38
SLIDE 38

green band = theory constraint on new physics Improved matrix elements may improve this band.

Experiment and Lattice: ∆Γq

  • L. Sabato, Lake Louise 2012

∆Γq = f 2

Bq [G1B1,q + G3B3,q] cos φq + O(1/mb, αs)

slide-39
SLIDE 39

Cross-checks

✦ We compare energies computed from 2pts to energies from the 2pt+3pt

fits.

✦ Two people fit . ✦ We use N=2,4,6 results to verify the plateau.

0.005 0.01 0.015 0.02

amq 1.235 1.24 1.245 1.25 1.255 1.26 1.265 1.27 aE0

from 2pt t = 5 - 30 O1 O2 O3

E0 for (0.0018, 0.0180)

6/16/12

O3

0.01 0.02 0.03 amq 0.12 0.14 0.16 0.18 0.2 0.22 β

(0.0124, 0.031) t=15 EDF (0.0124, 0.031) CMB

β versus mq for O3 fine ensembles

5/22/12