Mixing in the $Sigma^0$-$Lambda^0$ system is a direct consequence of broken isospin symmetry and is a measure of both isospin-symmetry breaking as well as general SU(3)-flavour symmetry breaking. In this work we present a new scheme for calculating the extent of $Sigma^0$-$Lambda^0$ mixing using simulations in lattice QCD+QED and perform several extrapolations that compare well with various past determinations. Our scheme allows us to easily contrast the QCD-only mixing case with the full QCD+QED mixing.
In this Reply, we respond to the above Comment. Our computation [Phys. Rev. D 91 (2015) 074512] only took into account pure QCD effects, arising from quark mass differences, so it is not surprising that there are discrepancies in isospin splittings and in the Sigma - Lambda mixing angle. We expect that these discrepancies will be smaller in a full calculation incorporating QED effects.
Isospin breaking effects in baryon octet (and decuplet) masses are due to a combination of up and down quark mass differences and electromagnetic effects and lead to small mass splittings. Between the Sigma and Lambda this mass splitting is much larger, this being mostly due to their different wavefunctions. However when isospin is broken, there is a mixing between between these states. We describe the formalism necessary to determine the QCD mixing matrix and hence find the mixing angle and mass splitting between the Sigma and Lambda particles due to QCD effects.
Systems with the quantum numbers of up to twelve charged and neutral pseudoscalar mesons, as well as one-, two-, and three-nucleon systems, are studied using dynamical lattice quantum chromodynamics and quantum electrodynamics (QCD+QED) calculations and effective field theory. QED effects on hadronic interactions are determined by comparing systems of charged and neutral hadrons after tuning the quark masses to remove strong isospin breaking effects. A non-relativistic effective field theory, which perturbatively includes finite-volume Coulomb effects, is analyzed for systems of multiple charged hadrons and found to accurately reproduce the lattice QCD+QED results. QED effects on charged multi-hadron systems beyond Coulomb photon exchange are determined by comparing the two- and three-body interaction parameters extracted from the lattice QCD+QED results for charged and neutral multi-hadron systems.
Chiral perturbation theory makes definitive predictions for the extrinsic behavior of hadrons in external electric and magnetic fields. Near the chiral limit, the electric and magnetic polarizabilities of pions, kaons, and nucleons are determined in terms of a few well-known parameters. In this limit, hadrons become quantum mechanically diffuse as polarizabilities scale with the inverse square-root of the quark mass. In some cases, however, such predictions from chiral perturbation theory have not compared well with experimental data. Ultimately we must turn to first principles numerical simulations of QCD to determine properties of hadrons, and confront the predictions of chiral perturbation theory. To address the electromagnetic polarizabilities, we utilize the background field technique. Restricting our attention to calculations in background electric fields, we demonstrate new techniques to determine electric polarizabilities and baryon magnetic moments for both charged and neutral states. As we can study the quark mass dependence of observables with lattice QCD, the lattice will provide a crucial test of our understanding of low-energy QCD, which will be timely in light of ongoing experiments, such as at COMPASS and HIgamma S.
We present lattice QCD calculations of nucleon electromagnetic form factors using pion masses $m_pi$ = 149, 202, and 254 MeV and an action with clover-improved Wilson quarks coupled to smeared gauge fields, as used by the Budapest-Marseille-Wuppertal collaboration. Particular attention is given to removal of the effects of excited state contamination by calculation at three source-sink separations and use of the summation and generalized pencil-of-function methods. The combination of calculation at the nearly physical mass $m_pi$ = 149 MeV in a large spatial volume ($m_pi L_s$ = 4.2) and removal of excited state effects yields agreement with experiment for the electric and magnetic form factors $G_E(Q^2)$ and $G_M(Q^2)$ up to $Q^2$ = 0.5 GeV$^2$.