Source author record

Daping Du

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Catalog footprint

Research overview

22works
5topics
4close collaborators

Published work

22 publications

preprint2016arXiv

B→Kl+l−B\to Kl^+l^- decay form factors from three-flavor lattice QCD

We compute the form factors for the B→Kl+l−B \to Kl^+l^- semileptonic decay process in lattice QCD using gauge-field ensembles with 2+1 flavors of sea quark, generated by the MILC Collaboration. The ensembles span lattice spacings from 0.12 to 0.045 fm and have multiple sea-quark masses to help control the chiral extrapolation. The asqtad improved staggered action is used for the light valence and sea quarks, and the clover action with the Fermilab interpretation is used for the heavy bb quark. We present results for the form factors f+(q2)f_+(q^2), f0(q2)f_0(q^2), and fT(q2)f_T(q^2), where q2q^2 is the momentum transfer, together with a comprehensive examination of systematic errors. Lattice QCD determines the form factors for a limited range of q2q^2, and we use the model-independent zz expansion to cover the whole kinematically allowed range. We present our final form-factor results as coefficients of the zz expansion and the correlations between them, where the errors on the coefficients include statistical and all systematic uncertainties. We use this complete description of the form factors to test QCD predictions of the form factors at high and low q2q^2. We also compare a Standard-Model calculation of the branching ratio for B→Kl+l−B \to Kl^+l^- with experimental data.

preprint2016arXiv

B(s)0B^0_{(s)}-mixing matrix elements from lattice QCD for the Standard Model and beyond

We calculate---for the first time in three-flavor lattice QCD---the hadronic matrix elements of all five local operators that contribute to neutral B0B^0- and BsB_s-meson mixing in and beyond the Standard Model. We present a complete error budget for each matrix element and also provide the full set of correlations among the matrix elements. We also present the corresponding bag parameters and their correlations, as well as specific combinations of the mixing matrix elements that enter the expression for the neutral BB-meson width difference. We obtain the most precise determination to date of the SU(3)-breaking ratio ξ=1.206(18)(6)ξ= 1.206(18)(6), where the second error stems from the omission of charm sea quarks, while the first encompasses all other uncertainties. The threefold reduction in total uncertainty, relative to the 2013 Flavor Lattice Averaging Group results, tightens the constraint from BB mixing on the Cabibbo-Kobayashi-Maskawa (CKM) unitarity triangle. Our calculation employs gauge-field ensembles generated by the MILC Collaboration with four lattice spacings and pion masses close to the physical value. We use the asqtad-improved staggered action for the light valence quarks, and the Fermilab method for the bottom quark. We use heavy-light meson chiral perturbation theory modified to include lattice-spacing effects to extrapolate the five matrix elements to the physical point. We combine our results with experimental measurements of the neutral BB-meson oscillation frequencies to determine the CKM matrix elements ∣Vtd∣=8.00(34)(8)×10−3|V_{td}| = 8.00(34)(8) \times 10^{-3}, ∣Vts∣=39.0(1.2)(0.4)×10−3|V_{ts}| = 39.0(1.2)(0.4) \times 10^{-3}, and ∣Vtd/Vts∣=0.2052(31)(10)|V_{td}/V_{ts}| = 0.2052(31)(10), which differ from CKM-unitarity expectations by about 2σσ. These results and others from flavor-changing-neutral currents point towards an emerging tension between weak processes that are mediated at the loop and tree levels.

preprint2016arXiv

Phenomenology of semileptonic B-meson decays with form factors from lattice QCD

We study the exclusive semileptonic BB-meson decays B→K(π)ℓ+ℓ−B\to K(π)\ell^+\ell^-, B→K(π)ννˉB\to K(π)ν\barν, and B→πτνB\toπτν, computing observables in the Standard model using the recent lattice-QCD results for the underlying form factors from the Fermilab Lattice and MILC Collaborations. These processes provide theoretically clean windows into physics beyond the Standard Model because the hadronic uncertainties are now under good control for suitably binned observables. For example, the resulting partially integrated branching fractions for B→πμ+μ−B\toπμ^+μ^- and B→Kμ+μ−B\to Kμ^+μ^- outside the charmonium resonance region are 1-2σσ higher than the LHCb Collaboration's recent measurements, where the theoretical and experimental errors are commensurate. The combined tension is 1.7σσ. Combining the Standard-Model rates with LHCb's measurements yields values for the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements ∣Vtd∣=7.45(69)×10−3|V_{td}|=7.45{(69)}\times10^{-3}, ∣Vts∣=35.7(1.5)×10−3|V_{ts}|=35.7(1.5)\times10^{-3}, and ∣Vtd/Vts∣=0.201(20)|V_{td}/V_{ts}|=0.201{(20)}, which are compatible with the values obtained from neutral B(s)B_{(s)}-meson oscillations and have competitive uncertainties. Alternatively, taking the CKM matrix elements from unitarity, we constrain new-physics contributions at the electroweak scale. The constraints on the Wilson coefficients Re(C9){\rm Re}(C_9) and Re(C10){\rm Re}(C_{10}) from B→πμ+μ−B\toπμ^+μ^- and B→Kμ+μ−B\to Kμ^+μ^- are competitive with those from B→K∗μ+μ−B\to K^* μ^+μ^-, and display a 2.0σσ tension with the Standard Model. Our predictions for B→K(π)ννˉB\to K(π)ν\barν and B→πτνB\toπτν are close to the current experimental limits.

preprint2015arXiv

∣Vub∣|V_{ub}| from B→πℓνB\toπ\ellν decays and (2+1)-flavor lattice QCD

We present a lattice-QCD calculation of the B→πℓνB\toπ\ellν semileptonic form factors and a new determination of the CKM matrix element ∣Vub∣|V_{ub}|. We use the MILC asqtad 2+1-flavor lattice configurations at four lattice spacings and light-quark masses down to 1/20 of the physical strange-quark mass. We extrapolate the lattice form factors to the continuum using staggered chiral perturbation theory in the hard-pion and SU(2) limits. We employ a model-independent zz parameterization to extrapolate our lattice form factors from large-recoil momentum to the full kinematic range. We introduce a new functional method to propagate information from the chiral-continuum extrapolation to the zz expansion. We present our results together with a complete systematic error budget, including a covariance matrix to enable the combination of our form factors with other lattice-QCD and experimental results. To obtain ∣Vub∣|V_{ub}|, we simultaneously fit the experimental data for the B→πℓνB\toπ\ellν differential decay rate obtained by the BaBar and Belle collaborations together with our lattice form-factor results. We find ∣Vub∣=(3.72±0.16)×10−3|V_{ub}|=(3.72\pm 0.16)\times 10^{-3} where the error is from the combined fit to lattice plus experiments and includes all sources of uncertainty. Our form-factor results bring the QCD error on ∣Vub∣|V_{ub}| to the same level as the experimental error. We also provide results for the B→πℓνB\toπ\ellν vector and scalar form factors obtained from the combined lattice and experiment fit, which are more precisely-determined than from our lattice-QCD calculation alone. These results can be used in other phenomenological applications and to test other approaches to QCD.

preprint2015arXiv

B→πℓℓB\toπ\ell\ell form factors for new-physics searches from lattice QCD

The rare decay B→πℓ+ℓ−B\toπ\ell^+\ell^- arises from b→db\to d flavor-changing neutral currents and could be sensitive to physics beyond the Standard Model. Here, we present the first abab-initioinitio QCD calculation of the B→πB\toπ tensor form factor fTf_T. Together with the vector and scalar form factors f+f_+ and f0f_0 from our companion work [J. A. Bailey et al.et~al., Phys. Rev. D 92, 014024 (2015)], these parameterize the hadronic contribution to B→πB\toπ semileptonic decays in any extension of the Standard Model. We obtain the total branching ratio BR(B+→π+μ+μ−)=20.4(2.1)×10−9{\text{BR}}(B^+\toπ^+μ^+μ^-)=20.4(2.1)\times10^{-9} in the Standard Model, which is the most precise theoretical determination to date, and agrees with the recent measurement from the LHCb experiment [R. Aaij et al.et~al., JHEP 1212, 125 (2012)]. Note added: after this paper was submitted for publication, LHCb announced a new measurement of the differential decay rate for this process [T. Tekampe, talk at DPF 2015], which we now compare to the shape and normalization of the Standard-Model prediction.

preprint2015arXiv

The B→DℓνB \to D \ell ν form factors at nonzero recoil and ∣Vcb∣|V_{cb}| from 2+12+1-flavor lattice QCD

We present the first unquenched lattice-QCD calculation of the hadronic form factors for the exclusive decay B‾→Dℓν‾\overline{B} \rightarrow D \ell \overlineν at nonzero recoil. We carry out numerical simulations on fourteen ensembles of gauge-field configurations generated with 2+1 flavors of asqtad-improved staggered sea quarks. The ensembles encompass a wide range of lattice spacings (approximately 0.045 to 0.12 fm) and ratios of light (up and down) to strange sea-quark masses ranging from 0.05 to 0.4. For the bb and cc valence quarks we use improved Wilson fermions with the Fermilab interpretation, while for the light valence quarks we use asqtad-improved staggered fermions. We extrapolate our results to the physical point using rooted staggered heavy-light meson chiral perturbation theory. We then parameterize the form factors and extend them to the full kinematic range using model-independent functions based on analyticity and unitarity. We present our final results for f+(q2)f_+(q^2) and f0(q2)f_0(q^2), including statistical and systematic errors, as coefficients of a series in the variable zz and the covariance matrix between these coefficients. We then fit the lattice form-factor data jointly with the experimentally measured differential decay rate from BaBar to determine the CKM matrix element, ∣Vcb∣=(39.6±1.7QCD+exp±0.2QED)×10−3|V_{cb}|=(39.6 \pm 1.7_{\rm QCD+exp} \pm 0.2_{\rm QED})\times 10^{-3}. As a byproduct of the joint fit we obtain the form factors with improved precision at large recoil. Finally, we use them to update our calculation of the ratio R(D)R(D) in the Standard Model, which yields R(D)=0.299(11)R(D) = 0.299(11).

preprint2014arXiv

Update of ∣Vcb∣|V_{cb}| from the Bˉ→D∗ℓνˉ\bar{B}\to D^*\ell\barν form factor at zero recoil with three-flavor lattice QCD

We compute the zero-recoil form factor for the semileptonic decay Bˉ0→D∗+ℓ−νˉ\bar{B}^0\to D^{*+}\ell^-\barν (and modes related by isospin and charge conjugation) using lattice QCD with three flavors of sea quarks. We use an improved staggered action for the light valence and sea quarks (the MILC \asqtad\ configurations), and the Fermilab action for the heavy quarks. Our calculations incorporate higher statistics, finer lattice spacings, and lighter quark masses than our 2008 work. As a byproduct of tuning the new data set, we obtain the DsD_s and BsB_s hyperfine splittings with few-MeV accuracy. For the zero-recoil form factor, we obtain F(1)=0.906(4)(12)\mathcal{F}(1)=0.906(4)(12), where the first error is statistical and the second is the sum in quadrature of all systematic errors. With the latest HFAG average of experimental results and a cautious treatment of QED effects, we find ∣Vcb∣=(39.04±0.49expt±0.53QCD±0.19QED)×10−3|V_{cb}| = (39.04 \pm 0.49_\text{expt} \pm 0.53_\text{QCD} \pm 0.19_\text{QED})\times10^{-3}. The QCD error is now commensurate with the experimental error.

preprint2013arXiv

B→πℓνB\toπ\ellν and B→πℓ+ℓ−B\toπ\ell^+\ell^- semileptonic form factors from unquenched lattice QCD

We update the lattice calculation of the B→πB\toπ semileptonic form factors, which have important applications to the CKM matrix element ∣Vub∣|V_{ub}| and the B→πℓ+ℓ−B\toπ\ell^+\ell^- rare decay. We use MILC asqtad ensembles with Nf=2+1N_f=2+1 sea quarks and over a range of lattice spacings a≈0.045a \approx 0.045--0.120.12 fm. We perform a combined chiral and continuum extrapolation of our lattice data using SU(2) staggered chiral perturbation theory in the hard pion limit. To extend the results for the form factors to the full kinematic range, we take a functional approach to parameterize the form factors using the Bourrely-Caprini-Lellouch formalism in a model-independent way. Our analysis is still blinded with an unknown off-set factor which will be disclosed when we present the final results.

preprint2013arXiv

Heavy-meson semileptonic decays for the Standard Model and beyond

We calculate the form factors for the semileptonic decays Bs→KℓνB_s\to K\ellν and B→KℓℓB\to K\ell\ell with lattice QCD. We work at several lattice spacings and a range of light quark masses, using the MILC 2+1-flavor asqtad ensembles. We use the Fermilab method for the bb quark. We obtain chiral-continuum extrapolations for EKE_K up to ∼1.2\sim1.2 GeV and then extend to the entire kinematic range with the model-independent zz expansion.

preprint2013arXiv

Kaon semileptonic vector form factor and determination of |V_{us}| using staggered fermions

Using staggered fermions and twisted boundary conditions, we calculate the K meson semileptonic decay vector form factor at zero momentum transfer. The HISQ formulation is used for the valence quarks, while the sea quarks are simulated with the asqtad action (MILC N_f=2+1 configurations). For the chiral and continuum extrapolation we use two-loop continuum CHPT, supplemented by partially quenched staggered CHPT at one loop. Our result is f_+^{Kπ}(0) = 0.9667+-0.0023+-0.0033, where the first error is statistical and the second is the sum in quadrature of the systematic uncertainties. This result is the first N_f=2+1 calculation with two lattice spacings and a controlled continuum extrapolation. It is also the most precise result to date for the vector form factor and, although the central value is larger than previous unquenched lattice calculations, it is compatible with them within errors. Combining our value for f_+^{Kπ}(0) with the latest experimental measurements of K semileptonic decays, we obtain |V_{us}| = 0.2238+-0.0009+-0.0005, where the first error is from f_+^{Kπ}(0) and the second one is experimental. As a byproduct of our calculation, we obtain the combination of low-energy constants [C_{12}^r+C_{34}^r-(L_5^r)^2](M_ρ) = (3.62+-1.00)x10^{-6}.

preprint2012arXiv

B_s->D_s/B->D Semileptonic Form-Factor Ratios and Their Application to BR(B^0_s->μ^+μ^-)

We calculate form-factor ratios between the semileptonic decays \bar{B}->D^+\ell^-\barν and \bar{B}_s->D_s^+\ell^-\barν with lattice QCD. These ratios are a key theoretical input in a new strategy to determine the fragmentation fractions of the neutral B decays, which are needed for measurements of BR(B^0_s-> μ^+μ^-). We use the MILC ensembles of gauge configurations with 2+1 flavors of sea quarks at two lattice spacings of approximately 0.12 fm and 0.09 fm. We use the model-independent z-parametrization to extrapolate our simulation results at small recoil toward maximum recoil. Our results for the form-factor ratios are f0(s)(Mπ2)/f0(d)(MK2)=1.046(44)stat.(15)syst.f_0^{(s)}(M^2_π)/f_0^{(d)}(M^2_K) =1.046(44)_{stat.}(15)_{syst.} and f0(s)(Mπ2)/f0(d)(Mπ2)=1.054(47)stat.(17)syst.f_0^{(s)}(M^2_π)/f_0^{(d)}(M^2_π)=1.054(47)_{stat.}(17)_{syst.}. In contrast to a QCD sum-rule calculation, no significant departure from U-spin (d<->s) symmetry is observed.

preprint2012arXiv

Density of states and Fisher's zeros in compact U(1) pure gauge theory

We present high-accuracy calculations of the density of states using multicanonical methods for lattice gauge theory with a compact gauge group U(1) on 4^4, 6^4 and 8^4 lattices. We show that the results are consistent with weak and strong coupling expansions. We present methods based on Chebyshev interpolations and Cauchy theorem to find the (Fisher's) zeros of the partition function in the complex beta=1/g^2 plane. The results are consistent with reweighting methods whenever the latter are accurate. We discuss the volume dependence of the imaginary part of the Fisher's zeros, the width and depth of the plaquette distribution at the value of beta where the two peaks have equal height. We discuss strategies to discriminate between first and second order transitions and explore them with data at larger volume but lower statistics. Higher statistics and even larger lattices are necessary to draw strong conclusions regarding the order of the transition.

preprint2012arXiv

Fisher zeros and conformality in lattice models

Fisher zeros are the zeros of the partition function in the complex beta=2N_c/g^2 plane. When they pinch the real axis, finite size scaling allows one to distinguish between first and second order transition and to estimate exponents. On the other hand, a gap signals confinement and the method can be used to explore the boundary of the conformal window. We present recent numerical results for 2D O(N) sigma models, 4D U(1) and SU(2) pure gauge and SU(3) gauge theory with N_f=4 and 12 flavors. We discuss attempts to understand some of these results using analytical methods. We discuss the 2-lattice matching and qualitative aspects of the renormalization group (RG) flows in the Migdal-Kadanoff approximation, in particular how RG flows starting at large beta seem to move around regions where bulk transitions occur. We consider the effects of the boundary conditions on the nonperturbative part of the average energy and on the Fisher zeros for the 1D O(2) model.

preprint2012arXiv

Refining new-physics searches in B -> D tau nu decay with lattice QCD

The semileptonic decay channel B -> D tau nu is sensitive to the presence of a scalar current, such as that mediated by a charged-Higgs boson. Recently the BaBar experiment reported the first observation of the exclusive semileptonic decay B -> D tau nu, finding an approximately 2-sigma disagreement with the Standard-Model prediction for the ratio R(D)=BR(B->D tau nu)/BR(B->D l nu), where l=e,mu. We compute this ratio of branching fractions using hadronic form factors computed in unquenched lattice QCD and obtain R(D) = 0.316(12)(7), where the errors are statistical and total systematic, respectively. This result is the first Standard-Model calculation of R(D) from ab initio full QCD. Its error is smaller than that of previous estimates, primarily due to the reduced uncertainty in the scalar form factor f_0(q^2). Our determination of R(D) is approximately 1-sigma higher than previous estimates and, thus, reduces the tension with experiment. We also compute R(D) in models with electrically charged scalar exchange, such as the type II two-Higgs doublet model. Once again, our result is consistent with, but approximately 1-sigma higher than, previous estimates for phenomenologically relevant values of the scalar coupling in the type II model. As a byproduct of our calculation, we also present the Standard-Model prediction for the longitudinal polarization ratio P_L (D)= 0.325(4)(3).

preprint2012arXiv

Semileptonic B to D decays at nonzero recoil with 2+1 flavors of improved staggered quarks. An update

The Fermilab Lattice and MILC collaborations are completing a comprehensive program of heavy-light physics on MILC (2+1)-flavor asqtad ensembles with lattice spacings as small as 0.045 fm and light-to-strange-quark mass ratios as low as 1/20. We use the Fermilab interpretation of the clover action for heavy valence quarks and the asqtad action for the light valence quarks. The central goal of the program is to provide ever more exacting tests of the unitarity of the CKM matrix. We present preliminary results for one part of the program, namely the analysis of the semileptonic decay B -> D l nu at nonzero recoil.

preprint2012arXiv

Semileptonic form factor ratio B_s->D_s/B->D and its application to BR(B^0_s->μ^+μ^-)

We present a (2+1)-flavor lattice QCD calculation of the form factor ratio between the semileptonic decays Bˉs0→Ds+l−νˉ\bar{B}^0_s \to D^+_sl^-\barν and Bˉ0→D+l−νˉ\bar{B}^0 \to D^+l^-\barν . This ratio is an important theoretical input to the hadronic determination of the BB meson fragmentation fraction ratio fs/fdf_s/f_d which enters in the measurement of BR(Bs0→μ+μ−)\mathrm{BR}(B^0_s\to μ^+μ^-). Small lattice spacings and high statistics enable us to simulate the decays with a dynamic final DD meson of small momentum and reliably extract the hadronic matrix elements at nonzero recoil. We report our preliminary result for the form factor ratio at the corresponding momentum transfer of the two decays f0(s)(Mπ2)/f0(d)(MK2)f_0^{(s)}(M^2_π)/f_0^{(d)}(M^2_K).

preprint2011arXiv

Fisher's zeros, complex RG flows and confinement in LGT models

The zeros of the partition function in the complex beta plane (Fisher's zeros) play an important role in our understanding of phase transitions and RG flows. Recently, we argued that they act as gates or separatrices for complex RG flows. Using histogram reweighting to construct the density of states, we calculate the Fisher's zeros for pure gauge SU(2) and U(1) on L^4 lattices. For SU(2), these zeros appear to move almost horizontally when the volume increases. They stay away from the real axis which indicates a confining theory at zero temperature. We discuss the effect of an adjoint term on these results. In contrast, using recent multicanonical simulations for the U(1) model for L up to 8 we find that the zeros pinch the real axis near beta =1.0113. Preliminary results concerning U(1) at larger volumes, SU(3) with 3 light flavors and plans to delimit the boundary of the conformal window are briefly discussed.

preprint2010arXiv

Fisher's zeros as boundary of renormalization group flows in complex coupling spaces

We propose new methods to extend the renormalization group transformation to complex coupling spaces. We argue that the Fisher's zeros are located at the boundary of the complex basin of attraction of infra-red fixed points. We support this picture with numerical calculations at finite volume for two-dimensional O(N) models in the large-N limit and the hierarchical Ising model. We present numerical evidence that, as the volume increases, the Fisher's zeros of 4-dimensional pure gauge SU(2) lattice gauge theory with a Wilson action, stabilize at a distance larger than 0.15 from the real axis in the complex beta=4/g^2 plane. We discuss the implications for proofs of confinement and searches for nontrivial infra-red fixed points in models beyond the standard model.

preprint2010arXiv

Fisher's zeros as boundary of RG flows in complex coupling space

We discuss the possibility of extending the RG flows to complex coupling spaces. We argue that the Fisher's zeros are located at the boundary of the complex basin of attraction of IR fixed points. We support this picture with numerical calculations at finite volume for2D O(N) models in the large-N limit and the hierarchical Ising model using the two-lattice matching method. We present numerical evidence supporting the idea that, as the volume increases, the Fisher's zeros of 4-dimensional pure gauge SU(2) lattice gauge theory with a Wilson action, stabilize at a distance larger than 0.1 from the real axis in the complex beta=4/g^2 plane. We show that when a positive adjoint term is added, the zeros get closer to the real axis. We compare the situation with the U(1) case. We discuss the implications of this new framework for proofs of confinement and searches for nontrivial IR fixed points in models beyond the standard model.

preprint2009arXiv

Dyson's Instability in Lattice Gauge Theory

We discuss Dyson's argument that the vacuum is unstable under a change g^2 -> - g^2, in the context of lattice gauge theory. For compact gauge groups, the partition function is well defined at negative g^2, but the average plaquette P has a discontinuity when g^2 changes sign. This reflects a change of vacuum rather than a loss of vacuum. In addition, P has poles in the complex g^2 plane, located at the complex zeros of the partition function (Fisher's zeros). We discuss the relevance of these singularities for lattice perturbation theory. We present new methods to locate Fisher's zeros using numerical values for the density of state in SU(2) and U(1) pure gauge theory. We briefly discuss similar issues for O(N) nonlinear sigma models where the local integrals are also over compact spaces.

preprint2008arXiv

Volume dependence of Fisher's zeros

We study the location of the partition function zeros in the complex beta plane (Fisher's Zeros) for SU(2) lattice gauge theory on L^4 lattices. We discuss recent attempts to locate complex zeros for L=4 and 6. We compare results obtained using various polynomial approximations of the logarithm of the density of states and a straightforward MC reweighting. We conclude that the method based on a combination of discrete Chebyshev orthogonality and patching plaquette distributions at different beta provides the more reliable estimates.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.