HPQCD Collaboration 16 July 2015
Our approaches to b quarks In Glasgow, we take two complementary approaches to b quarks: Nonrelativistic QCD and heavy HISQ. Here I will focus exclusively on NRQCD (for b quarks). So why NRQCD? b quark can be simulated at its physical mass. NRQCD proceeds with a relatively straightforward evolution equation computationally inexpensive.
Gluon Field Configurations [1212.4768] Gluon field configurations provided by MILC Collaboration with 2 + 1 + 1 flavours of HISQ quarks in the sea. Those marked (*) are ensembles with physical light quark masses. Set β am l am s am c L/a T/a n cfg 1 5.80 0.013 0.065 0.838 16 48 1020 2 5.80 0.0064 0.064 0.828 24 48 1000 3 5.80 0.00235 0.0647 0.831 32 48 1000 4 6.00 0.0102 0.0509 0.635 24 64 1052 5 6.00 0.00507 0.0507 0.628 32 64 1000 6 6.00 0.00184 0.0507 0.628 48 64 1000 7 6.30 0.0074 0.037 0.440 32 96 1008 8 6.30 0.0012 0.0363 0.432 64 96 621 9 6.72 0.0048 0.024 0.286 48 144 1000
NRQCD We use the following NRQCD Hamiltonian: ( e ah = 1 aδh ) ( 1 ah ) n 0 U t 2 2n where we include terms up to O(v 4 ) ah 0 = (2) 2am b ( 1 ah 0 2n ( (2) ) 2 aδh = c 1 8(am b ) 3 + c i ( 2 8(am b ) 2 g ( ) c 3 8 (am b ) 2 σ Ẽ Ẽ ) ( 1 aδh ) 2 ) Ẽ Ẽ ( g a c 4 σ B 2 (4) a ) (2) 2 + c 5 c 6 2am b 24am b 16n (am b ) 2 with most c i coefficients O(α s ) improved. [1110.6887],[1105.5309]
Υ Decay Constant & Leptonic Width [1408.5768] Leptonic Width: Γ(Υ (n) e + e ) = 4π 3 α2 eme 2 b f 2 Υ (n) M Υ (n) Decay constant: 0 J V,i Υ (n) j = f Υ (n)m Υ (n)δ ij The Υ leptonic width is experimentally well measured, so a determination of this quantity on the lattice is a good test of QCD, and also of our approach to b quark physics.
Time moments [1408.5768] We need to renormalise our NRQCD currents, ( ) J V = Z V J (0) V,NRQCD + k 1J (1) V,NRQCD by determining both k 1 and overall normalisation Z V. We do this by comparing NRQCD to continuum QCD perturbation theory for time moments of the correlators: G V n = gv n (α s, µ/m b ) [am b (µ)] n 2 G V n = Z 2 V G V,NRQCD n
Υ Leptonic Width [1408.5768] 3 fυ M Υ (GeV 3/2 ) 2 1 0 m l /m s = 0.2 m l physical Experiment 0 0.005 0.01 0.015 0.02 0.025 a 2 (fm 2 ) Decay Constant f Υ = 0.649(31) GeV Leptonic Width Γ(Υ e + e ) = 1.19(11) kev
B Meson Decay Constants [1503.05762] We have updated our picture of B meson decay constants: B, B s and B c as well as B, B s and B c MB s R s = f Bs f Bs MBs R s 1.02 1.01 1 0.99 0.98 0.97 0.96 0.95 0.94 0.93 0.92 0 0.005 0.01 0.015 0.02 0.025 a 2 [fm 2 ]
f B q /f Bq [1503.05762] u, d, s, c sea B /B B s /B s HPQCD NRQCD this paper u, d sea twisted mass 1407.1091 Lucha et al 1410.6684 sum rules 0.8 0.9 1 1.1 1.2 f B q /f Bq Gelhausen et al 1305.5432
Decay Constants: summary DECAY CONSTANT [GeV] 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 Experiment : weak decays :emdecays Lattice QCD : predictions :postdictions π K B B D Bs B s φ D s Ds ψ η c ψ B c B c Υ η b Υ
Determination of m b [1408.5768] m b is an important parameter in the Standard Model, for example in accurately determining Higgs branching fraction to bb. Recall: G V n = gv n (α s, µ/m b ) [am b (µ)] n 2 so we can make a determination of m b using NRQCD time moments.
Determination of m b [1408.5768] mb(µ = 4.18, nf = 4) GeV 4.4 4.35 4.3 4.25 4.2 4.15 m l /m s = 0.2 m l physical 4.1 0 0.005 0.01 0.015 0.02 0.025 a 2 (fm 2 ) m b (µ = m b, n f = 5) = 4.196(23) GeV
m b : lattice summary HPQCD NRQCD JJ [1408.5768] HPQCD HISQ JJ n f = 3 [1004.4285] HPQCD HISQ ratio n f = 4 [1408.4169] HPQCD NRQCD E 0 [1302.3739] ETMC ratio [1311.2837] 4.0 4.1 4.2 4.3 4.4 4.5 m b (m b,n f = 5)(GeV)
Semileptonic decays We use NRQCD b quarks and HISQ light quarks in our calculation of B πlν. t J π B T We pick multiple values of T and fit B and π 2-point correlators simultaneously with 3-point correlators.
B π at zero recoil For B and π at rest, q 2 max 26.5 GeV 2, π V 0 B = f 0 ( q 2 max ) (mb + m π ) Soft pion theorem relates this to decay constants in m π 0 limit: f 0 (q 2 max) = f B /f π This relation was previously studied in lattice QCD, but it did not seem to hold. We now have the opportunity to study with light quarks at physical mass.
Soft Pion Theorem ) 1 + mπ mb 1.1 1 0.9 a 0.15 fm a 0.12 fm a 0.09 fm ( f0(q 2 max)/(fb/fπ) 0.8 0.7 0.6 0.5 0.4 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 m π (GeV)
A brief look at the future... We are optimistic about being able to explore the full q 2 range in B c η c lν decays on superfine lattices (a 0.06 fm). p ηc = 2.4 GeV in the B c rest frame at q 2 = 0
Summary Bottomonium: Using time moments from the lattice and continuum perturbation theory we have: calculated the Υ and Υ leptonic widths using NRQCD b quarks. made an accurate determination of m b (m b ). B Physics Decay constants have been calculated for B q in addition to B q We have calculated B π at q 2 max; we find results consistent with soft pion theorem. This program is now being extended and includes a full range of q 2 for B c η c on superfine ensembles.
Thank you!
Backup Slides
Υ Leptonic Width [1408.5768] fυ M Υ /f Υ M Υ Decay Constant f Υ 0.9 0.8 0.7 0.6 0.5 = 0.481(39) GeV 1 A = 0 J V Υ 0 J V Υ = f Υ f Υ MΥ M Υ Experiment 0 0.005 0.01 0.015 0.02 0.025 a 2 (fm 2 ) Leptonic Width Γ(Υ e + e ) = 0.69(9) kev
Z V [1408.5768] ZV 1.4 1.2 1.0 0.8 0.6 0.4 0.2 fine, J NRQCD, k 1 =-0.155, Z =1.019(8) 0.0 0 5 10 15 20 moment n
Time moments [1408.5768] Time moments of our NRQCD correlators defined as, G V,NRQCD n = 2 t (t/a) n C V,NRQCD (t)e (E0 M kin)t. for n = 4, 6, 8,...
Fit forms Bottomonium B π 3 h(a, m sea) =h phys 1 + b l δm sea/(10m s) + c j (aλ) 2j j=1 2 + (aλ) 2j [ ] c jb δx m + c jbb (δxm) 2 j=1 3 Γ(a, m π) =Γ phys 1 + c j (aλ) 2j + b j m j π j=1 ( ) Λ 2 [ + (aλ) 2 c b δx m + (aλ) 2 (δx m) 2] m b ( m + da 2 m 2 2 ) ] π l π log m 2 π 1.6
Determination of m b [1408.5768] [ ] n 2 [ ] 1/2 Rn V = rn V mb Rn r (α MS, µ/m b ) n 2 M kin m b (µ) R n 2 r n 2m b = M Υ,η b 2m b (µ) m b (µ) = M Υ,η b 2 [ ] 1/2 Rn 2 r n 2m b R n r n 2 M kin
Determination of m b [1408.5768] mb(µ = 4.18, nf = 4) (GeV) 4.6 4.5 4.4 4.3 4.2 4.1 fine very coarse 4 0 5 10 15 20 25 moment n