
We present a firstprinciples lattice QCD+QED calculation at physical pion
mass of the leadingorder 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 Rratio 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, Rratio statistical,
and Rratio systematic errors given separately. This is the currently most
precise determination of the leadingorder hadronic vacuum polarization
contribution to the muon anomalous magnetic moment. In addition, we present the
first lattice calculation of the lightquark 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 CPviolating quark
chromoEDM 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 nearlyphysical
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 subpercent statistical errors and
negligible chiral and finitevolume 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 quarkconnected vacuum polarisation
effects. Simulations were performed using RBCUKQCD'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
ydirection. 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
standardmodel 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 noisereduction 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$ physicalpionmass 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 RBCUKQCD gaugefield ensembles
generated with the domainwall fermion and Iwasaki gauge actions. For the $b$
quarks we use the anisotropic clover action with a relativistic heavyquark
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 lightquark masses and to the
continuum and interpolate in the pion/kaon energy using SU(2) "hardpion"
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 modelindependent $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 crosschecks. 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$ formfactor data with the experimental branchingfraction
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$ StandardModel
differential branching fractions and forwardbackward 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 axialvector
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>_{ud}$ 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 loweigenmode
deflation and allmodeaveraging (AMA). We present a method for choosing
optimal AMA parameters.

We report on the first complete calculation of the $K_LK_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
GlashowIliopoulosMaiani cancellation. For these heavierthanphysical
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
partiallyquenched 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 strangequark mass in the MSbarscheme 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
MSbarscheme 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 currentcurrent 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 selfcontained introduction to the muon anomalous magnetic
moment, (g2), we review the status of lattice calculations of the hadronic
vacuum polarization contribution and present first results from lattice QCD for
the hadronic lightbylight 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 \pimesons
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 socalled 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
CPviolation 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, fourquark
operators are evaluated, including the disconnected graphs in which no quark
line connects the initial kaon and final twopion states. These lattice
operators are nonperturbatively renormalized using the RomeSouthampton 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 nonperturbatively in the
RIMOM scheme and are then converted into the MSbar scheme using continuum
perturbation theory. In addition to RIMOM, we introduce and implement four
nonexceptional intermediate momentum schemes that suppress infrared
nonperturbative uncertainties in the renormalization procedure. We compute the
conversion factors relating the matrix elements in this family of RISMOM
schemes and MSbar at oneloop order. Comparison of the results obtained using
these different intermediate schemes allows for a more reliable estimate of the
unknown higherorder 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 offshell 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 octetsinglet 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 nondegenerate 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 nexttoleading 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_um_d)=5.840(96)
MeV.
All errors are statistical only, and the nexttoleading order pion splitting
is only approximate; it does not contain all nexttoleading order
contributions. We also computed the protonneutron mass difference, including
for the first time, QED interactions in a realistic 2+1 flavor calculation. We
find $(m_pm_n)_{\rm QED}=0.383(68)$ MeV, (m_pm_n)_(m_um_d)=2.51(14) MeV,
and the total m_pm_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. Noncompact QED is treated in
the quenched approximation.
We present new results for the electromagnetic low energy constants in SU(3)
and SU(2) partiallyquenched chiral perturbation theory to the nexttoleading
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)_Rpluskaon chiral perturbation theory, including the oneloop 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 nonperturbative 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 MSbar scheme is safely carried out. For numerical
simulations we adopted the Iwasaki gauge action and nonperturbatively 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
CPPACS/JLQCD collaboration, which covered the updown quark mass range heavier
than $m_\pi\sim 500$ MeV and that of PACSCS 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 RomeSouthampton regularization independent
momentumsubtraction renormalization scheme(RI/MOM) for bilinear operators to
one with a nonexceptional, symmetric subtraction point. Twopoint Green's
functions with the insertion of quark bilinear operators are computed with
scalar, pseudoscalar, vector, axialvector and tensor operators at oneloop
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 oneloop 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 masspreconditioned
domaindecomposed 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
nonperturbative running from the low energy region, where physical scale is
introduced, to deep in the high energy perturbative region, where conversion to
the MSbar scheme is safely performed. Possible systematic errors due to the
use of perturbation theory occur only in the conversion from threeflavor to
fourflavor 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
nonperturbatively 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
CPPACS/JLQCD collaboration, which covered the updown 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 oneloop
order. The baryon mass data used for the analyses are obtained for the
degenerate updown 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 mesonbaryon 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 heavyheavy and heavylight quark systems for charm with a
relativistic heavy quark action in 2+1 flavor lattice QCD. Configurations are
generated by the PACSCS 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.