Search

We present lattice-QCD calculations of the hadronic form factors for the semileptonic decays $D\to\pi\ell\nu$, $D\to K\ell\nu$, and $D_s\to K\ell\nu$. Our calculation uses the highly improved staggered quark (HISQ) action for all valence and sea quarks and includes $N_f=2+1+1$ MILC ensembles with lattice spacings ranging from $a\approx0.12$ fm down to $0.042$ fm. At most lattice spacings, an ensemble with physical-mass light quarks is included. The HISQ action allows all the quarks to be treated with the same relativistic light-quark action, allowing for nonperturbative renormalization using partial conservation of the vector current. We combine our results with experimental measurements of the differential decay rates to determine $|V_{cd}|^{D\to\pi}=0.2238(11)^{\rm Expt}(15)^{\rm QCD}(04)^{\rm EW}(02)^{\rm SIB}[22]^{\rm QED}$ and $|V_{cs}|^{D\to K}=0.9589(23)^{\rm Expt}(40)^{\rm QCD}(15)^{\rm EW}(05)^{\rm SIB}[95]^{\rm QED}$ This result for $|V_{cd}|$ is the most precise to date, with a lattice-QCD error that is, for the first time for the semileptonic extraction, at the same level as the experimental error. Using recent measurements from BES III, we also give the first-ever determination of $|V_{cd}|^{D_s\to K}=0.258(15)^{\rm Expt}(01)^{\rm QCD}[03]^{\rm QED}$ from $D_s\to K \ell\nu$. Our results also furnish new Standard Model calculations of the lepton flavor universality ratios $R^{D\to\pi}=0.98671(17)^{\rm QCD}[500]^{\rm QED}$, $R^{D\to K}=0.97606(16)^{\rm QCD}[500]^{\rm QED}$, and $R^{D_s\to K}=0.98099(10)^{\rm QCD}[500]^{\rm QED}$, which are consistent within $2\sigma$ with experimental measurements. Our extractions of $|V_{cd}|$ and $|V_{cs}|$, when combined with a value for $|V_{cb}|$, provide the most precise test of second-row CKM unitarity, finding agreement with unitarity at the level of one standard deviation., Comment: 92 pages, V2 matches version accepted for publication in PRD. Expanded supplementary material for reconstructing our final results. An implementation of nonlinear shrinkage is also included

We study isovector unpolarized and helicity parton distribution functions (PDF) of the proton within the framework of Large Momentum Effective Theory. We use a gauge ensemble, generated by the MILC Collaboration, with a superfine lattice spacing of $0.042$ fm and a pion mass of $310$ MeV, enabling us to simultaneously reach sub-fermi spatial separations and larger nucleon momenta. We compare the spatial dependence of quasi-PDF matrix elements in different renormalization schemes with the corresponding results of the global fits, obtained using 1-loop perturbative matching. We present determinations of the first four moments of the unpolarized and helicity PDFs of proton from the Ioffe-time dependence of the isovector matrix elements, obtained by employing a ratio-based renormalization scheme., Comment: 28 pages, 34 figures. Accepted version for PRD

There have been rapid developments in the direct calculation in lattice QCD (LQCD) of the Bjorken-$x$ dependence of hadron structure through large-momentum effective theory (LaMET). LaMET overcomes the previous limitation of LQCD to moments (that is, integrals over Bjorken-$x$) of hadron structure, allowing LQCD to directly provide the kinematic regions where the experimental values are least known. LaMET requires large-momentum hadron states to minimize its systematics and allow us to reach small-$x$ reliably. This means that very fine lattice spacing to minimize lattice artifacts at order $(P_z a)^n$ will become crucial for next-generation LaMET-like structure calculations. Furthermore, such calculations require operators with long Wilson-link displacements (in finer lattice units), increasing the communication costs relative to that of the propagator inversion. In this work, we explore whether machine-learning (ML) algorithms can make correlator predictions to reduce the computational cost of these LQCD calculations. We consider two algorithms, gradient-boosting decision tree and linear models, applied to LaMET data, the matrix elements needed to determine the kaon and $\eta_s$ unpolarized parton distribution functions (PDFs), meson distribution amplitude (DA), and the nucleon gluon PDF. We find that both algorithms can reliably predict the target observables with different fit quality and systematic errors. The predictions from smaller displacement $z$ to larger ones work better than those for momentum $p$ due to the higher correlation among the data., Comment: 20 pages, 22 figures

We discuss preliminary results for the vector form factors $f_+^{\{\pi,K\}}$ at zero-momentum transfer for the decays $D\to\pi\ell\nu$ and $D\to K \ell\nu$ using MILC's $N_f = 2+1+1$ HISQ ensembles at four lattice spacings, $a \approx 0.042, 0.06, 0.09$, and 0.12 fm, and various HISQ quark masses down to the (degenerate) physical light quark mass. We use the kinematic constraint $f_+(q^2)= f_0(q^2)$ at $q^2 = 0$ to determine the vector form factor from our study of the scalar current, which yields $f_0(0)$. Results are extrapolated to the continuum physical point in the framework of hard pion/kaon SU(3) heavy-meson-staggered $\chi$PT and Symanzik effective theory. Our calculation improves upon the precision achieved in existing lattice-QCD calculations of the vector form factors at $q^2=0$. We show the values of the CKM matrix elements $|V_{cs}|$ and $|V_{cd}|$ that we would obtain using our preliminary results for the form factors together with recent experimental results, and discuss the implications of these values for the second row CKM unitarity., Comment: 10 pages, 3 figures, proceeding of The 36th Annual International Symposium on Lattice Field Theory

We report a state-of-the-art lattice calculation of the isovector quark transversity distribution of the proton at the physical pion mass. Within the framework of large-momentum effective theory (LaMET), we compute the transversity quasi-distributions using clover valence fermions on 2+1+1-flavor (up/down, strange, charm) HISQ-lattice configurations with boosted proton momenta as large as 3.0~GeV. The relevant lattice matrix elements are nonperturbatively renormalized in regularization-independent momentum-subtraction (RI/MOM) scheme and systematically matched to the physical transversity distribution. With high statistics, large proton momenta and meticulous control of excited-state contamination, we provide the best theoretical prediction for the large-$x$ isovector quark transversity distribution, with better precision than the most recent global analyses of experimental data. Our result also shows that the sea quark asymmetry in the proton transversity distribution is consistent with zero, which has been assumed in all current global analyses., Comment: 10 pages, 4 figures

We present a state-of-the-art calculation of the isovector quark helicity Bjorken-$x$ distribution in the proton using lattice-QCD ensembles at the physical pion mass. We compute quasi-distributions at proton momenta $P_z \in \{2.2, 2.6, 3.0\}$~GeV on the lattice, and match them systematically to the physical parton distribution using large-momentum effective theory (LaMET). We reach an unprecedented precision through high statistics in simulations, large-momentum proton matrix elements, and control of excited-state contamination. The resulting distribution with combined statistical and systematic errors is in agreement with the latest phenomenological analysis of the spin-dependent experimental data; in particular, $\Delta \bar{u}(x)>\Delta \bar{d}(x)$., Comment: 6 pages, 4 figures

With recent developments in parallel supercomputing architecture, many core, multi-core, and GPU processors are now commonplace, resulting in more levels of parallelism, memory hierarchy, and programming complexity. It has been necessary to adapt the MILC code to these new processors starting with NVIDIA GPUs, and more recently, the Intel Xeon Phi processors. We report on our efforts to port and optimize our code for the Intel Knights Landing architecture. We consider performance of the MILC code with MPI and OpenMP, and optimizations with QOPQDP and QPhiX. For the latter approach, we concentrate on the staggered conjugate gradient and gauge force. We also consider performance on recent NVIDIA GPUs using the QUDA library.

We review our work done to optimize the staggered conjugate gradient (CG) algorithm in the MILC code for use with the Intel Knights Landing (KNL) architecture. KNL is the second gener- ation Intel Xeon Phi processor. It is capable of massive thread parallelism, data parallelism, and high on-board memory bandwidth and is being adopted in supercomputing centers for scientific research. The CG solver consumes the majority of time in production running, so we have spent most of our effort on it. We compare performance of an MPI+OpenMP baseline version of the MILC code with a version incorporating the QPhiX staggered CG solver, for both one-node and multi-node runs., Comment: 8 pages, 4 figures

The nucleon axial form factor is a dominant contribution to errors in neutrino oscillation studies. Lattice QCD calculations can help control theory errors by providing first-principles information on nucleon form factors. In these proceedings, we present preliminary results on a blinded calculation of $g_A$ and the axial form factor using HISQ staggered baryons with 2+1+1 flavors of sea quarks. Calculations are done using physical light quark masses and are absolutely normalized. We discuss fitting form factor data with the model-independent $z$ expansion parametrization., Comment: 7 pages, 3 figures

The conjugate gradient (CG) algorithm is among the most essential and time consuming parts of lattice calculations with staggered quarks. We test the performance of CG and dslash, the key step in the CG algorithm, on the Intel Xeon Phi, also known as the Many Integrated Core (MIC) architecture. We try different parallelization strategies using MPI, OpenMP, and the vector processing units (VPUs)., Comment: 7 Pages, 1 figure, contribution to the 32nd International Symposium on Lattice Field Theory (Lattice 2014), 23-28 June 2014, Columbia University, New York, NY, USA