Lattice calculation of heavy-light decay constants with two flavors of dynamical quarks

We present results for f(B), f(Bs), f(D), f(Ds) and their ratios in the presence of two flavors of light sea quarks (N-f=2). We use Wilson light valence quarks and Wilson and static heavy valence quarks; the sea quarks are simulated with staggered fermions. Additional quenched simulations with nonperturbatively improved clover fermions allow us to improve our control of the continuum extrapolation. For our central values the masses of the sea quarks are not extrapolated to the physical u, d masses; that is, the central values are "partially quenched." A calculation using "fat-link clover" valence fermions is also discussed but is not included in our final results. We find, for example, f(B)=190(7)((+24)(-17))((+11)(-2))((+8)(-0)) MeV, f(B s)/f(B)=1.16(1)(2)(2)((+4)(-0)), f(D s)=241(5)((+27)(-26))((+9)(-4))((+5)(-0)) MeV, and f(B)/f(D s)=0.79(2) ((+5)(-4))(3)((+5)(-0)), where in each case the first error is statistical and the remaining three are systematic: the error within the partially quenched N-f=2 approximation, the error due to the missing strange sea quark and to partial quenching, and an estimate of the effects of chiral logarithms at small quark mass. The last error, though quite significant in decay constant ratios, appears to be smaller than has been recently suggested by Kronfeld and Ryan, and Yamada. We emphasize, however, that as in other lattice computations to date, the lattice u,d quark masses are not very light and chiral log effects may not be fully under control.