Charmed and Light Pseudoscalar Meson Decay Constants from Four-Flavor Lattice QCD with Physical Light Quarks

INSPIRE: 1306303
arXiv: 1407.3772
DOI: 10.1103/PhysRevD.90.074509

Authors: Bazavov, A., Bernard, C., Bouchard, C.M., DeTar, C., Du, D., El-Khadra, A.X., Foley, J., Freeland, E.D., Gamiz, E., Gottlieb, Steven, Heller, U.M., Kim, J., Komijani, J., Kronfeld, A.S., Laiho, J., Levkova, L., Mackenzie, P.B., Neil, E.T., Simone, J.N., Sugar, R., Toussaint, D., Van de Water, R.S., Zhou, R.

Submitted: 14 July 2014

Subjects:

Journal reference: Phys.Rev.D 90 (7) 074509 (2014)

Citations: 207

Abstract

<p>We compute the leptonic decay constants <inline-formula><mml:math display="inline"><mml:msub><mml:mi>f</mml:mi><mml:msup><mml:mi>D</mml:mi><mml:mo>+</mml:mo></mml:msup></mml:msub></mml:math></inline-formula>, <inline-formula><mml:math display="inline"><mml:mrow><mml:msub><mml:mrow><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>D</mml:mi></mml:mrow><mml:mrow><mml:mi>s</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>, and <inline-formula><mml:math display="inline"><mml:msub><mml:mi>f</mml:mi><mml:msup><mml:mi>K</mml:mi><mml:mo>+</mml:mo></mml:msup></mml:msub></mml:math></inline-formula> and the quark-mass ratios <inline-formula><mml:math display="inline"><mml:msub><mml:mi>m</mml:mi><mml:mi>c</mml:mi></mml:msub><mml:mo stretchy="false">/</mml:mo><mml:msub><mml:mi>m</mml:mi><mml:mi>s</mml:mi></mml:msub></mml:math></inline-formula> and <inline-formula><mml:math display="inline"><mml:msub><mml:mi>m</mml:mi><mml:mi>s</mml:mi></mml:msub><mml:mo stretchy="false">/</mml:mo><mml:msub><mml:mi>m</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:math></inline-formula> in unquenched lattice QCD using the experimentally determined value of <inline-formula><mml:math display="inline"><mml:msub><mml:mi>f</mml:mi><mml:msup><mml:mi>π</mml:mi><mml:mo>+</mml:mo></mml:msup></mml:msub></mml:math></inline-formula> for normalization. We use the MILC highly improved staggered quark ensembles with four dynamical quark flavors—up, down, strange, and charm—and with both physical and unphysical values of the light sea-quark masses. The use of physical pions removes the need for a chiral extrapolation, thereby eliminating a significant source of uncertainty in previous calculations. Four different lattice spacings ranging from <inline-formula><mml:math display="inline"><mml:mi>a</mml:mi><mml:mo>≈</mml:mo><mml:mn>0.06</mml:mn></mml:math></inline-formula> to 0.15 fm are included in the analysis to control the extrapolation to the continuum limit. Our primary results are <inline-formula><mml:math display="inline"><mml:mrow><mml:msub><mml:mrow><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mi>D</mml:mi></mml:mrow><mml:mrow><mml:mo>+</mml:mo></mml:mrow></mml:msup></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>212.6</mml:mn><mml:mo stretchy="false">(</mml:mo><mml:mn>0.4</mml:mn><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mfrac linethickness="0"><mml:mrow><mml:mo>+</mml:mo><mml:mn>1.0</mml:mn></mml:mrow><mml:mrow><mml:mo>−</mml:mo><mml:mn>1.2</mml:mn></mml:mrow></mml:mfrac></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mtext> </mml:mtext><mml:mtext> </mml:mtext><mml:mi>MeV</mml:mi></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math display="inline"><mml:mrow><mml:msub><mml:mrow><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>D</mml:mi></mml:mrow><mml:mrow><mml:mi>s</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>249.0</mml:mn><mml:mo stretchy="false">(</mml:mo><mml:mn>0.3</mml:mn><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mfrac linethickness="0"><mml:mrow><mml:mo>+</mml:mo><mml:mn>1.1</mml:mn></mml:mrow><mml:mrow><mml:mo>−</mml:mo><mml:mn>1.5</mml:mn></mml:mrow></mml:mfrac></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mtext> </mml:mtext><mml:mtext> </mml:mtext><mml:mi>MeV</mml:mi></mml:mrow></mml:math></inline-formula>, and <inline-formula><mml:math display="inline"><mml:mrow><mml:msub><mml:mrow><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>D</mml:mi></mml:mrow><mml:mrow><mml:mi>s</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub><mml:mo stretchy="false">/</mml:mo><mml:msub><mml:mrow><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mi>D</mml:mi></mml:mrow><mml:mrow><mml:mo>+</mml:mo></mml:mrow></mml:msup></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>1.1712</mml:mn><mml:mo stretchy="false">(</mml:mo><mml:mn>10</mml:mn><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mfrac linethickness="0"><mml:mrow><mml:mo>+</mml:mo><mml:mn>29</mml:mn></mml:mrow><mml:mrow><mml:mo>−</mml:mo><mml:mn>32</mml:mn></mml:mrow></mml:mfrac></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:math></inline-formula>, where the errors are statistical and total systematic, respectively. The errors on our results for the charm decay constants and their ratio are approximately 2–4 times smaller than those of the most precise previous lattice calculations. We also obtain <inline-formula><mml:math display="inline"><mml:mrow><mml:msub><mml:mrow><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mi>K</mml:mi></mml:mrow><mml:mrow><mml:mo>+</mml:mo></mml:mrow></mml:msup></mml:mrow></mml:msub><mml:mo stretchy="false">/</mml:mo><mml:msub><mml:mrow><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mi>π</mml:mi></mml:mrow><mml:mrow><mml:mo>+</mml:mo></mml:mrow></mml:msup></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>1.1956</mml:mn><mml:mo stretchy="false">(</mml:mo><mml:mn>10</mml:mn><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mfrac linethickness="0"><mml:mrow><mml:mo>+</mml:mo><mml:mn>26</mml:mn></mml:mrow><mml:mrow><mml:mo>−</mml:mo><mml:mn>18</mml:mn></mml:mrow></mml:mfrac></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:math></inline-formula>, updating our previous result, and determine the quark-mass ratios <inline-formula><mml:math display="inline"><mml:mrow><mml:msub><mml:mrow><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">/</mml:mo><mml:msub><mml:mrow><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:mi>l</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>27.35</mml:mn><mml:mo stretchy="false">(</mml:mo><mml:mn>5</mml:mn><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mfrac linethickness="0"><mml:mrow><mml:mo>+</mml:mo><mml:mn>10</mml:mn></mml:mrow><mml:mrow><mml:mo>−</mml:mo><mml:mn>7</mml:mn></mml:mrow></mml:mfrac></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math display="inline"><mml:mrow><mml:msub><mml:mrow><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:mi>c</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">/</mml:mo><mml:msub><mml:mrow><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>11.747</mml:mn><mml:mo stretchy="false">(</mml:mo><mml:mn>19</mml:mn><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mfrac linethickness="0"><mml:mrow><mml:mo>+</mml:mo><mml:mn>59</mml:mn></mml:mrow><mml:mrow><mml:mo>−</mml:mo><mml:mn>43</mml:mn></mml:mrow></mml:mfrac></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:math></inline-formula>. When combined with experimental measurements of the decay rates, our results lead to precise determinations of the Cabibbo-Kobayashi-Maskawa matrix elements <inline-formula><mml:math display="inline"><mml:mo stretchy="false">|</mml:mo><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">|</mml:mo><mml:mo>=</mml:mo><mml:mn>0.22487</mml:mn><mml:mo stretchy="false">(</mml:mo><mml:mn>51</mml:mn><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mn>29</mml:mn><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mn>20</mml:mn><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mn>5</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula>, <inline-formula><mml:math display="inline"><mml:mo stretchy="false">|</mml:mo><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mi>c</mml:mi><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">|</mml:mo><mml:mo>=</mml:mo><mml:mn>0.217</mml:mn><mml:mo stretchy="false">(</mml:mo><mml:mn>1</mml:mn><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mn>5</mml:mn><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mn>1</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> and <inline-formula><mml:math display="inline"><mml:mo stretchy="false">|</mml:mo><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mi>c</mml:mi><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">|</mml:mo><mml:mo>=</mml:mo><mml:mn>1.010</mml:mn><mml:mo stretchy="false">(</mml:mo><mml:mn>5</mml:mn><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mn>18</mml:mn><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mn>6</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula>, where the errors are from this calculation of the decay constants, the uncertainty in the experimental decay rates, structure-dependent electromagnetic corrections, and, in the case of <inline-formula><mml:math display="inline"><mml:mo stretchy="false">|</mml:mo><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">|</mml:mo></mml:math></inline-formula>, the uncertainty in <inline-formula><mml:math display="inline"><mml:mo stretchy="false">|</mml:mo><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">|</mml:mo></mml:math></inline-formula>, respectively.</p>