• We present a first-principles lattice QCD+QED calculation at physical pion mass of the leading-order hadronic vacuum polarization contribution to the muon anomalous magnetic moment. The total contribution of up, down, strange, and charm quarks including QED and strong isospin breaking effects is found to be $a_\mu^{\rm HVP~LO}=715.4(16.3)(9.2) \times 10^{-10}$, where the first error is statistical and the second is systematic. By supplementing lattice data for very short and long distances with experimental R-ratio data using the compilation of Ref. [1], we significantly improve the precision of our calculation and find $a_\mu^{\rm HVP~LO} = 692.5(1.4)(0.5)(0.7)(2.1) \times 10^{-10}$ with lattice statistical, lattice systematic, R-ratio statistical, and R-ratio systematic errors given separately. This is the currently most precise determination of the leading-order hadronic vacuum polarization contribution to the muon anomalous magnetic moment. In addition, we present the first lattice calculation of the light-quark QED correction at physical pion mass.
  • We analyze commonly used expressions for computing the nucleon electric dipole form factors (EDFF) $F_3$ and moments (EDM) on a lattice and find that they lead to spurious contributions from the Pauli form factor $F_2$ due to inadequate definition of these form factors when parity mixing of lattice nucleon fields is involved. Using chirally symmetric domain wall fermions, we calculate the proton and the neutron EDFF induced by the CP-violating quark chromo-EDM interaction using the corrected expression. In addition, we calculate the electric dipole moment of the neutron using background electric field that respects time translation invariance and boundary conditions, and find that it decidedly agrees with the new formula but not the old formula for $F_3$. Finally, we analyze some selected lattice results for the nucleon EDM and observe that after the correction is applied, they either agree with zero or are substantially reduced in magnitude, thus reconciling their difference from phenomenological estimates of the nucleon EDM.
  • We present results for several light hadronic quantities ($f_\pi$, $f_K$, $B_K$, $m_{ud}$, $m_s$, $t_0^{1/2}$, $w_0$) obtained from simulations of 2+1 flavor domain wall lattice QCD with large physical volumes and nearly-physical pion masses at two lattice spacings. We perform a short, O(3)%, extrapolation in pion mass to the physical values by combining our new data in a simultaneous chiral/continuum `global fit' with a number of other ensembles with heavier pion masses. We use the physical values of $m_\pi$, $m_K$ and $m_\Omega$ to determine the two quark masses and the scale - all other quantities are outputs from our simulations. We obtain results with sub-percent statistical errors and negligible chiral and finite-volume systematics for these light hadronic quantities, including: $f_\pi$ = 130.2(9) MeV; $f_K$ = 155.5(8) MeV; the average up/down quark mass and strange quark mass in the $\bar {\rm MS}$ scheme at 3 GeV, 2.997(49) and 81.64(1.17) MeV respectively; and the neutral kaon mixing parameter, $B_K$, in the RGI scheme, 0.750(15) and the $\bar{\rm MS}$ scheme at 3 GeV, 0.530(11).
  • We present results for the leading hadronic contribution to the muon anomalous magnetic moment due to strange quark-connected vacuum polarisation effects. Simulations were performed using RBC--UKQCD's $N_f=2+1$ domain wall fermion ensembles with physical light sea quark masses at two lattice spacings. We consider a large number of analysis scenarios in order to obtain solid estimates for residual systematic effects. Our final result in the continuum limit is $a_\mu^{(2)\,{\rm had},\,s}=53.1(9)\left(^{+1}_{-3}\right)\times10^{-10}$.
  • In this document we address an error discovered in the ensemble generation for our calculation of the $I=0$ $K\to\pi\pi$ amplitude (Phys. Rev. Lett. 115, 212001 (2015), arXiv:1505.07863) whereby the same random numbers were used for the two independent quark flavors, resulting in small but measurable correlations between gauge observables separated by 12 units in the y-direction. We conclude that the effects of this error are negligible compared to the overall errors on our calculation.
  • We report the first lattice QCD calculation of the complex kaon decay amplitude $A_0$ with physical kinematics, using a $32^3\times 64$ lattice volume and a single lattice spacing $a$, with $1/a= 1.3784(68)$ GeV. We find Re$(A_0) = 4.66(1.00)(1.26) \times 10^{-7}$ GeV and Im$(A_0) = -1.90(1.23)(1.08) \times 10^{-11}$ GeV, where the first error is statistical and the second systematic. The first value is in approximate agreement with the experimental result: Re$(A_0) = 3.3201(18) \times 10^{-7}$ GeV while the second can be used to compute the direct CP violating ratio Re$(\varepsilon'/\varepsilon)=1.38(5.15)(4.59)\times 10^{-4}$, which is $2.1\sigma$ below the experimental value $16.6(2.3)\times 10^{-4}$. The real part of $A_0$ is CP conserving and serves as a test of our method while the result for Re$(\varepsilon'/\varepsilon)$ provides a new test of the standard-model theory of CP violation, one which can be made more accurate with increasing computer capability.
  • We report the first lattice QCD calculation of the hadronic vacuum polarization disconnected contribution to the muon anomalous magnetic moment at physical pion mass. The calculation uses a refined noise-reduction technique which enabled the control of statistical uncertainties at the desired level with modest computational effort. Measurements were performed on the $48^3 \times 96$ physical-pion-mass lattice generated by the RBC and UKQCD collaborations. We find $a_\mu^{\rm HVP~(LO)~DISC} = -9.6(3.3)(2.3)\times 10^{-10}$, where the first error is statistical and the second systematic.
  • We calculate the $B \to\pi\ell\nu$ and $B_s \to K \ell\nu$ form factors in dynamical lattice QCD. We use the (2+1)-flavor RBC-UKQCD gauge-field ensembles generated with the domain-wall fermion and Iwasaki gauge actions. For the $b$ quarks we use the anisotropic clover action with a relativistic heavy-quark interpretation. We analyze two lattice spacings $a \approx 0.11, 0.086$ fm and unitary pion masses as light as $M_\pi \approx 290$ MeV. We simultaneously extrapolate our numerical results to the physical light-quark masses and to the continuum and interpolate in the pion/kaon energy using SU(2) "hard-pion" chiral perturbation theory. We provide complete error budgets for the form factors $f_+(q^2)$ and $f_0(q^2)$ at three momenta that span the $q^2$ range accessible in our numerical simulations. We extrapolate these results to $q^2 = 0$ using a model-independent $z$-parametrization and present our final form factors as the $z$-coefficients and the matrix of correlations between them. Our results agree with other lattice determinations using staggered light quarks and provide important independent cross-checks. Both $B \to\pi\ell\nu$ and $B_s \to K \ell\nu$ decays enable a determination of the CKM matrix element $|V_{ub}|$. To illustrate this, we perform a combined $z$-fit of our numerical $B\to\pi\ell\nu$ form-factor data with the experimental branching-fraction measurements leaving the relative normalization as a free parameter; we obtain $|V_{ub}| = 3.61(32) \times 10^{-3}$, where the error includes statistical and systematic uncertainties. This approach can be applied to $B_s\to K \ell\nu$ decay to determine $|V_{ub}|$ once the process has been measured experimentally. Finally, in anticipation of future measurements, we make predictions for $B \to \pi\ell\nu$ and $B_s\to K \ell\nu$ Standard-Model differential branching fractions and forward-backward asymmetries.
  • We report initial nucleon structure results computed on lattices with 2+1 dynamical M\"obius domain wall fermions at the physical point generated by the RBC and UKQCD collaborations. At this stage, we evaluate only connected quark contributions. In particular, we discuss the nucleon vector and axial-vector form factors, nucleon axial charge and the isovector quark momentum fraction. From currently available statistics, we estimate the stochastic accuracy of the determination of $g_A$ and $<x>_{u-d}$ to be around 10%, and we expect to reduce that to 5% within the next year. To reduce the computational cost of our calculations, we extensively use acceleration techniques such as low-eigenmode deflation and all-mode-averaging (AMA). We present a method for choosing optimal AMA parameters.
  • We report on the first complete calculation of the $K_L-K_S$ mass difference, $\Delta M_K$, using lattice QCD. The calculation is performed on a 2+1 flavor, domain wall fermion ensemble with a 330MeV pion mass and a 575 MeV kaon mass. We use a quenched charm quark with a 949 MeV mass to implement Glashow-Iliopoulos-Maiani cancellation. For these heavier-than-physical particle masses, we obtain $\Delta M_K =3.19(41)(96)\times 10^{-12}$ MeV, quite similar to the experimental value. Here the first error is statistical and the second is an estimate of the systematic discretization error. An interesting aspect of this calculation is the importance of the disconnected diagrams, a dramatic failure of the OZI rule.
  • We present physical results for a variety of light hadronic quantities obtained via a combined analysis of three 2+1 flavour domain wall fermion ensemble sets. For two of our ensemble sets we used the Iwasaki gauge action with beta=2.13 (a^-1=1.75(4) GeV) and beta=2.25 (a^-1=2.31(4) GeV) and lattice sizes of 24^3 x 64 and 32^3 x 64 respectively, with unitary pion masses in the range 293(5)-417(10) MeV. The extent L_s for the 5^th dimension of the domain wall fermion formulation is L_s=16 in these ensembles. In this analysis we include a third ensemble set that makes use of the novel Iwasaki+DSDR (Dislocation Suppressing Determinant Ratio) gauge action at beta = 1.75 (a^-1=1.37(1) GeV) with a lattice size of 32^3 x 64 and L_s=32 to reach down to partially-quenched pion masses as low as 143(1) MeV and a unitary pion mass of 171(1) MeV, while retaining good chiral symmetry and topological tunneling. We demonstrate a significant improvement in our control over the chiral extrapolation, resulting in much improved continuum predictions for the above quantities. The main results of this analysis include the pion and kaon decay constants, f_\pi=127(3)_{stat}(3)_{sys} MeV and f_K = 152(3)_{stat}(2)_{sys} MeV respectively (f_K/f_\pi = 1.199(12)_{stat}(14)_{sys}); the average up/down quark mass and the strange-quark mass in the MSbar-scheme at 3 GeV, m_{ud}(MSbar, 3 GeV) = 3.05(8)_{stat}(6)_{sys} MeV and m_s(MSbar, 3 GeV) = 83.5(1.7)_{stat}(1.1)_{sys}; the neutral kaon mixing parameter in the MSbar-scheme at 3 GeV, B_K(MSbar,3 GeV) = 0.535(8)_{stat}(13)_{sys}, and in the RGI scheme, \hat B_K = 0.758(11)_{stat}(19)_{sys}; and the Sommer scales r_1 = 0.323(8)_{stat}(4)_{sys} fm and r_0 = 0.480(10)_{stat}(4)_{sys} (r_1/r_0 = 0.673(11)_{stat}(3)_{sys}). We also obtain values for the SU(2) ChPT effective couplings, \bar{l_3} = 2.91(23)_{stat}(7)_{sys}$ and \bar{l_4} = 3.99(16)_{stat}(9)_{sys}.
  • We develop and demonstrate techniques needed to compute the long distance contribution to the $K_{L}$-$K_{S}$ mass difference, $\Delta M_K$, in lattice QCD and carry out a first, exploratory calculation of this fundamental quantity. The calculation is performed on 2+1 flavor, domain wall fermion, $16^3\times32$ configurations with a 421 MeV pion mass and an inverse lattice spacing $1/a=1.73$ GeV. We include only current-current operators and drop all disconnected and double penguin diagrams. The short distance part of the mass difference in a 2+1 flavor calculation contains a quadratic divergence cut off by the lattice spacing. Here, this quadratic divergence is eliminated through the GIM mechanism by introducing a valence charm quark. The inclusion of the charm quark makes the complete calculation accessible to lattice methods provided the discretization errors associated with the charm quark can be controlled. The long distance effects are discussed for each parity channel separately. While we can see a clear signal in the parity odd channel, the signal to noise ratio in the parity even channel is exponentially decreasing as the separation between the two weak operators increases. We obtain a mass difference $\Delta M_K$ which ranges from $6.58(30)\times 10^{-12}$ MeV to $11.89(81)\times 10^{-12}$ MeV for kaon masses varying from 563 MeV to 839 MeV. Extensions of these methods are proposed which promise accurate results for both $\Delta M_K$ and $\epsilon_K$, including long distance effects.
  • After a brief self-contained introduction to the muon anomalous magnetic moment, (g-2), we review the status of lattice calculations of the hadronic vacuum polarization contribution and present first results from lattice QCD for the hadronic light-by-light scattering contribution. The signal for the latter is consistent with model calculations. While encouraging, the statistical error is large and systematic errors are mostly uncontrolled. The method is applied first to pure QED as a check.
  • We describe the computation of the amplitude A_2 for a kaon to decay into two pions with isospin I=2. The results presented in the letter Phys.Rev.Lett. 108 (2012) 141601 from an analysis of 63 gluon configurations are updated to 146 configurations giving Re$A_2=1.381(46)_{\textrm{stat}}(258)_{\textrm{syst}} 10^{-8}$ GeV and Im$A_2=-6.54(46)_{\textrm{stat}}(120)_{\textrm{syst}}10^{-13}$ GeV. Re$A_2$ is in good agreement with the experimental result, whereas the value of Im$A_2$ was hitherto unknown. We are also working towards a direct computation of the $K\to(\pi\pi)_{I=0}$ amplitude $A_0$ but, within the standard model, our result for Im$A_2$ can be combined with the experimental results for Re$A_0$, Re$A_2$ and $\epsilon^\prime/\epsilon$ to give Im$A_0/$Re$A_0= -1.61(28)\times 10^{-4}$ . Our result for Im\,$A_2$ implies that the electroweak penguin (EWP) contribution to $\epsilon^\prime/\epsilon$ is Re$(\epsilon^\prime/\epsilon)_{\mathrm{EWP}} = -(6.25 \pm 0.44_{\textrm{stat}} \pm 1.19_{\textrm{syst}}) \times 10^{-4}$.
  • We report on the first realistic \emph{ab initio} calculation of a hadronic weak decay, that of the amplitude $A_2$ for a kaon to decay into two \pi-mesons with isospin 2. We find Re$A_2=(1.436\pm 0.063_{\textrm{stat}}\pm 0.258_{\textrm{syst}})\,10^{-8}\,\textrm{GeV}$ in good agreement with the experimental result and for the hitherto unknown imaginary part we find {Im}$\,A_2=-(6.83 \pm 0.51_{\textrm{stat}} \pm 1.30_{\textrm{syst}})\,10^{-13}\,{\rm GeV}$. Moreover combining our result for Im\,$A_2$ with experimental values of Re\,$A_2$, Re\,$A_0$ and $\epsilon^\prime/\epsilon$, we obtain the following value for the unknown ratio Im\,$A_0$/Re\,$A_0$ within the Standard Model: $\mathrm{Im}\,A_0/\mathrm{Re}\,A_0=-1.63(19)_{\mathrm{stat}}(20)_{\mathrm{syst}}\times10^{-4}$. One consequence of these results is that the contribution from Im\,$A_2$ to the direct CP violation parameter $\epsilon^{\prime}$ (the so-called Electroweak Penguin, EWP, contribution) is Re$(\epsilon^\prime/\epsilon)_{\mathrm{EWP}} = -(6.52 \pm 0.49_{\textrm{stat}} \pm 1.24_{\textrm{syst}}) \times 10^{-4}$. We explain why this calculation of $A_2$ represents a major milestone for lattice QCD and discuss the exciting prospects for a full quantitative understanding of CP-violation in kaon decays.
  • We report a direct lattice calculation of the $K$ to $\pi\pi$ decay matrix elements for both the $\Delta I=1/2$ and 3/2 amplitudes $A_0$ and $A_2$ on 2+1 flavor, domain wall fermion, $16^3\times32\times16$ lattices. This is a complete calculation in which all contractions for the required ten, four-quark operators are evaluated, including the disconnected graphs in which no quark line connects the initial kaon and final two-pion states. These lattice operators are non-perturbatively renormalized using the Rome-Southampton method and the quadratic divergences are studied and removed. This is an important but notoriously difficult calculation, requiring high statistics on a large volume. In this paper we take a major step towards the computation of the physical $K\to\pi\pi$ amplitudes by performing a complete calculation at unphysical kinematics with pions of mass 422\,MeV at rest in the kaon rest frame. With this simplification we are able to resolve Re$(A_0)$ from zero for the first time, with a 25% statistical error and can develop and evaluate methods for computing the complete, complex amplitude $A_0$, a calculation central to understanding the $\Delta =1/2$ rule and testing the standard model of CP violation in the kaon system.
  • We determine the neutral kaon mixing matrix element $B_K$ in the continuum limit with 2+1 flavors of domain wall fermions, using the Iwasaki gauge action at two different lattice spacings. These lattice fermions have near exact chiral symmetry and therefore avoid artificial lattice operator mixing. We introduce a significant improvement to the conventional NPR method in which the bare matrix elements are renormalized non-perturbatively in the RI-MOM scheme and are then converted into the MSbar scheme using continuum perturbation theory. In addition to RI-MOM, we introduce and implement four non-exceptional intermediate momentum schemes that suppress infrared non-perturbative uncertainties in the renormalization procedure. We compute the conversion factors relating the matrix elements in this family of RI-SMOM schemes and MSbar at one-loop order. Comparison of the results obtained using these different intermediate schemes allows for a more reliable estimate of the unknown higher-order contributions and hence for a correspondingly more robust estimate of the systematic error. We also apply a recently proposed approach in which twisted boundary conditions are used to control the Symanzik expansion for off-shell vertex functions leading to a better control of the renormalization in the continuum limit. We control chiral extrapolation errors by considering both the NLO SU(2) chiral effective theory, and an analytic mass expansion. We obtain $B_K^{\msbar}(3 GeV) = 0.529(5)_{stat}(15)_\chi(2)_{FV}(11)_{NPR}$. This corresponds to $\hat{B}_K = 0.749(7)_{stat}(21)_\chi(3)_{FV}(15)_{NPR}$. Adding all sources of error in quadrature we obtain $\hat{B}_K = 0.749(27)_{combined}$, with an overall combined error of 3.6%.
  • The large mass of the ninth pseudoscalar meson, the $\eta^\prime$, is believed to arise from the combined effects of the axial anomaly and the gauge field topology present in QCD. We report a realistic, 2+1 flavor, lattice QCD calculation of the $\eta$ and $\eta^\prime$ masses and mixing which confirms this picture. The physical eigenstates show small octet-singlet mixing with a mixing angle of $\theta = -14.1(2.8)^\circ$. Extrapolation to physical light quark mass gives, with statistical errors only, $m_\eta=573(6)$ MeV and $m_{\eta^\prime}=947(142)$ MeV, consistent with the experimental values of 548 MeV and 958 MeV.
  • Results are presented for the electromagnetic mass splittings of the low lying hadrons. These are used to determine the non-degenerate light quark masses. It is found that m_u=2.24(10)(34), m_d=4.65(15)(32), and $m_s=97.6(2.9)(5.5)$ MeV (MSbar scheme, 2 GeV scale). The first error is statistical and the second systematic. We find the lowest order electromagnetic splitting (m_pi+-m_pi0)_QED=3.38(23) MeV, the splittings including next-to-leading order, (m_pi+-m_pi0)_QED=4.50(23) MeV, (m_K+-m_K0)_QED=1.87(10) MeV, and the m_u != m_d contribution to the kaon mass difference, (m_K+-m_K0)_(m_u-m_d)=-5.840(96) MeV. All errors are statistical only, and the next-to-leading order pion splitting is only approximate; it does not contain all next-to-leading order contributions. We also computed the proton-neutron mass difference, including for the first time, QED interactions in a realistic 2+1 flavor calculation. We find $(m_p-m_n)_{\rm QED}=0.383(68)$ MeV, (m_p-m_n)_(m_u-m_d)=-2.51(14) MeV, and the total m_p-m_n=-2.13(16)(70) MeV, where the first error is statistical, and the second, part of the systematic error. We use domain wall fermions and the Iwasaki gauge action (gauge coupling beta=2.13). We use two lattice sizes, 16^3 and 24^3, to address finite volume effects. Non-compact QED is treated in the quenched approximation. We present new results for the electromagnetic low energy constants in SU(3) and SU(2) partially-quenched chiral perturbation theory to the next-to-leading order, obtained from fits to our data. Detailed analysis of systematic errors in our results and methods for improving them are discussed. Finally, new analytic results for SU(2)_L x SU(2)_R-plus-kaon chiral perturbation theory, including the one-loop logs proportional to alpha_em*m, are given.
  • We present an evaluation of the quark mass renormalization factor for Nf=2+1 QCD. The Schroedinger functional scheme is employed as the intermediate scheme to carry out non-perturbative running from the low energy region, where renormalization of bare mass is performed on the lattice, to deep in the high energy perturbative region, where the conversion to the renormalization group invariant mass or the MS-bar scheme is safely carried out. For numerical simulations we adopted the Iwasaki gauge action and non-perturbatively improved Wilson fermion action with the clover term. Seven renormalization scales are used to cover from low to high energy regions and three lattice spacings to take the continuum limit at each scale. The regularization independent step scaling function of the quark mass for the Nf=2+1 QCD is obtained in the continuum limit. Renormalization factors for the pseudo scalar density and the axial vector current are also evaluated for the same action and the bare couplings as two recent large scale Nf=2+1 simulations; previous work of the CP-PACS/JLQCD collaboration, which covered the up-down quark mass range heavier than $m_\pi\sim 500$ MeV and that of PACS-CS collaboration for much lighter quark masses down to $m_\pi=155$ MeV. The quark mass renormalization factor is used to renormalize bare PCAC masses in these simulations.
  • We extend the Rome-Southampton regularization independent momentum-subtraction renormalization scheme(RI/MOM) for bilinear operators to one with a nonexceptional, symmetric subtraction point. Two-point Green's functions with the insertion of quark bilinear operators are computed with scalar, pseudoscalar, vector, axial-vector and tensor operators at one-loop order in perturbative QCD. We call this new scheme RI/SMOM, where the S stands for "symmetric". Conversion factors are derived, which connect the RI/SMOM scheme and the MSbar scheme and can be used to convert results obtained in lattice calculations into the MSbar scheme. Such a symmetric subtraction point involves nonexceptional momenta implying a lattice calculation with substantially suppressed contamination from infrared effects. Further, we find that the size of the one-loop corrections for these infrared improved kinematics is substantially decreased in the case of the pseudoscalar and scalar operator, suggesting a much better behaved perturbative series. Therefore it should allow us to reduce the error in the determination of the quark mass appreciably.
  • We present the results of the physical point simulation in 2+1 flavor lattice QCD with the nonperturbatively $O(a)$-improved Wilson quark action and the Iwasaki gauge action at $\beta=1.9$ on a $32^3 \times 64$ lattice. The physical quark masses together with the lattice spacing is determined with $m_\pi$, $m_K$ and $m_\Omega$ as physical inputs. There are two key algorithmic ingredients to make possible the direct simulation at the physical point: One is the mass-preconditioned domain-decomposed HMC algorithm to reduce the computational cost. The other is the reweighting technique to adjust the hopping parameters exactly to the physical point. The physics results include the hadron spectrum, the quark masses and the pseudoscalar meson decay constants. The renormalization factors are nonperturbatively evaluated with the Schr{\"o}dinger functional method. The results are compared with the previous ones obtained by the chiral extrapolation method.
  • We present an evaluation of the running coupling constant for Nf=2+1 QCD. The Schroedinger functional scheme is used as the intermediate scheme to carry out non-perturbative running from the low energy region, where physical scale is introduced, to deep in the high energy perturbative region, where conversion to the MS-bar scheme is safely performed. Possible systematic errors due to the use of perturbation theory occur only in the conversion from three-flavor to four-flavor running coupling constant near the charm mass threshold, where higher order terms beyond 5th order in the $\beta$ function may not be negligible. For numerical simulations we adopted Iwasaki gauge action and non-perturbatively improved Wilson fermion action with the clover term. Seven renormalization scales are used to cover from low to high energy region and three lattice spacings to take the continuum limit at each scale. A physical scale is introduced from the previous Nf=2+1 simulation of the CP-PACS/JL-QCD collaboration, which covered the up-down quark mass range heavier than $m_\pi\sim 500$ MeV.
  • We investigate the quark mass dependence of baryon masses in 2+1 flavor lattice QCD using SU(3) heavy baryon chiral perturbation theory up to one-loop order. The baryon mass data used for the analyses are obtained for the degenerate up-down quark mass of 3 MeV to 24 MeV and two choices of the strange quark mass around the physical value. We find that the SU(3) chiral expansion fails to describe both the octet and the decuplet baryon data if phenomenological values are employed for the meson-baryon couplings. The SU(2) case is also examined for the nucleon. We observe that higher order terms are controlled only around the physical point. We also evaluate finite size effects using SU(3) heavy baryon chiralperturbation theory, finding small values of order 1% even at the physical point.
  • We study heavy-heavy and heavy-light quark systems for charm with a relativistic heavy quark action in 2+1 flavor lattice QCD. Configurations are generated by the PACS-CS Collaboration at the lattice spacing is $a=0.09$ fm with the lattice size of $32^3\times 64$ employing the $O(a)$-improved Wilson quark action and the Iwasaki gauge action. We present preliminary results for the charmonium spectrum and the $D$ and $D_s$ meson decay constants evaluated at 3.5 MeV$< m_{\rm ud}<$ 12 MeV with $m_{\rm s}$ around the physical value. We investigate the dynamical quark mass dependences of the hyperfine and the orbital splittings. The decay constants are compared with the recent experimental values.