• We present a determination of the corrections to Dashen's theorem and of the individual up and down quark masses from a lattice calculation based on quenched QED and $N_f=2+1$ QCD simulations with 5 lattice spacings down to 0.054 fm. The simulations feature lattice sizes up to 6 fm and average up-down quark masses all the way down to their physical value. For the parameter which quantifies violations to Dashens's theorem we obtain $\epsilon=0.73(2)(5)(17)$, where the first error is statistical, the second is systematic, and the third is an estimate of the QED quenching error. For the light quark masses we obtain, $m_u=2.27(6)(5)(4) \, MeV$ and $m_d=4.67(6)(5)(4) \, MeV$ in the $\overline{MS}$ scheme at $2 \, GeV$ and the isospin breaking ratios $m_u/m_d=0.485(11)(8)(14)$, $R=38.2(1.1)(0.8)(1.4)$ and $Q=23.4(0.4)(0.3)(0.4)$. Our results exclude the $m_u=0$ solution to the strong CP problem by more than 24 standard deviations.
  • We compute the slope and curvature, at vanishing four-momentum transfer squared, of the leading order hadron vacuum polarization function, using lattice QCD. Calculations are performed with 2+1+1 flavors of staggered fermions directly at the physical values of the quark masses and in volumes of linear extent larger than 6fm. The continuum limit is carried out using six different lattice spacings. All connected and disconnected contributions are calculated, up to and including those of the charm.
  • We present a QCD calculation of the $u$, $d$ and $s$ scalar quark contents of nucleons based on $47$ lattice ensembles with $N_f = 2+1$ dynamical sea quarks, $5$ lattice spacings down to $0.054\,\text{fm}$, lattice sizes up to $6\,\text{fm}$ and pion masses down to $120\,\text{MeV}$. Using the Feynman-Hellmann theorem, we obtain $f^N_{ud} = 0.0405(40)(35)$ and $f^N_s = 0.113(45)(40)$, which translates into $\sigma_{\pi N}=38(3)(3)\,\text{MeV}$, $\sigma_{sN}=105(41)(37)\,\text{MeV}$ and $y_N=0.20(8)(8)$ for the sigma terms and the related ratio, where the first errors are statistical and the second are systematic. Using isospin relations, we also compute the individual up and down quark contents of the proton and neutron (results in the main text).
  • In a previous letter (arXiv:1306.2287) we determined the isospin mass splittings of the baryon octet from a lattice calculation based on quenched QED and $N_f{=}2{+}1$ QCD simulations with 5 lattice spacings down to $0.054~\mathrm{fm}$, lattice sizes up to $6~\mathrm{fm}$ and average up-down quark masses all the way down to their physical value. Using the same data we determine here the corrections to Dashen's theorem and the individual up and down quark masses. For the parameter which quantifies violations to Dashens's theorem, we obtain $\epsilon=0.73(2)(5)(17)$, where the first error is statistical, the second is systematic, and the third is an estimate of the QED quenching error. For the light quark masses we obtain, $m_u=2.27(6)(5)(4)~\mathrm{MeV}$ and $m_d=4.67(6)(5)(4)~\mathrm{MeV}$ in the $\bar{\mathrm{MS}}$ scheme at $2~\mathrm{GeV}$ and the isospin breaking ratios $m_u/m_d=0.485(11)(8)(14)$, $R=38.2(1.1)(0.8)(1.4)$ and $Q=23.4(0.4)(0.3)(0.4)$. Our results exclude the $m_u=0$ solution to the strong CP problem by more than $24$ standard deviations.
  • By using lattice QCD computations we determine the sigma terms and strangeness content of all octet baryons by means of an application of the Hellmann-Feynman theorem. In addition to polynomial and rational expressions for the quark mass dependence of octet members, we use SU(3) covariant baryon chiral perturbation theory to perform the extrapolation to the physical up and down quark masses. Our N_f=2+1 lattice ensembles include pion masses down to about 190 MeV in large volumes (M_\pi L > 4), and three values of the lattice spacing. Our main results are the nucleon sigma term \sigma_{\pi N} = 39(4)(^{+18}_{-7}) and the strangeness content y_{N} = 0.20(7)(^{+13}_{-17}). Under the assumption of validity of covariant baryon \chi PT in our range of masses one finds y_{N} = 0.276(77)(^{+90}_{-62}).
  • The existence and stability of atoms rely on the fact that neutrons are more massive than protons. The measured mass difference is only 0.14\% of the average of the two masses. A slightly smaller or larger value would have led to a dramatically different universe. Here, we show that this difference results from the competition between electromagnetic and mass isospin breaking effects. We performed lattice quantum-chromodynamics and quantum-electrodynamics computations with four nondegenerate Wilson fermion flavors and computed the neutron-proton mass-splitting with an accuracy of $300$ kilo-electron volts, which is greater than $0$ by $5$ standard deviations. We also determine the splittings in the $\Sigma$, $\Xi$, $D$ and $\Xi_{cc}$ isospin multiplets, exceeding in some cases the precision of experimental measurements.
  • 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 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$.
  • We report the first Lattice QCD calculation using the almost physical pion mass mpi=149 MeV that agrees with experiment for four fundamental isovector observables characterizing the gross structure of the nucleon: the Dirac and Pauli radii, the magnetic moment, and the quark momentum fraction. The key to this success is the combination of using a nearly physical pion mass and excluding the contributions of excited states. An analogous calculation of the nucleon axial charge governing beta decay has inconsistencies indicating a source of bias at low pion masses not present for the other observables and yields a result that disagrees with experiment.
  • Recent results for moments of multiplicity distributions of net-protons and net-electric charge from the STAR collaboration are compared to lattice QCD results for higher order fluctuations of baryon number and electric charge by the Wuppertal-Budapest collaboration, with the purpose of extracting the freeze-out temperature and chemical potential. All lattice simulations are performed for a system of 2+1 dynamical quark flavors, at the physical mass for light and strange quarks; all results are continuum extrapolated. We show that it is possible to extract an upper value for the freeze-out temperature, as well as precise baryo-chemical potential values corresponding to the four highest collision energies of the experimental beam energy scan. Consistency between the freeze-out parameters obtained from baryon number and electric charge fluctuations is found. The freeze-out chemical potentials are now in agreement with the statistical hadronization model.
  • While electromagnetic and up-down quark mass difference effects on octet baryon masses are very small, they have important consequences. The stability of the hydrogen atom against beta decay is a prominent example. Here we include these effects by adding them to valence quarks in a lattice QCD calculation based on $N_f=2+1$ simulations with 5 lattice spacings down to 0.054 fm, lattice sizes up to 6 fm and average up-down quark masses all the way down to their physical value. This allows us to gain control over all systematic errors, except for the one associated with neglecting electromagnetism in the sea. We compute the octet baryon isomultiplet mass splittings, as well as the individual contributions from electromagnetism and the up-down quark mass difference. Our results for the total splittings are in good agreement with experiment.
  • We present continuum extrapolated lattice results for the higher order fluctuations of conserved charges in high temperature Quantum Chromodynamics. Through the matching of the grand canonical ensemble on the lattice to the net charge and net baryon distribution realized in heavy ion experiments the temperature and the chemical potential may be estimated at the time of chemical freeze-out
  • In lattice QCD computations a substantial amount of work is spent in solving linear systems arising in Wilson's discretization of the Dirac equations. We show first numerical results of the extension of the two-level DD-\alpha AMG method to a true multilevel method based on our parallel MPI-C implementation. Using additional levels pays off, allowing to cut down the core minutes spent on one system solve by a factor of approximately 700 compared to standard Krylov subspace methods and yielding another speed-up of a factor of 1.7 over the two-level approach.
  • We present our results for ratios of higher order fluctuations of electric charge as functions of the temperature. These results are obtained in a system of 2+1 quark flavors at physical quark masses and continuum extrapolated. We compare them to preliminary data on higher order moments of the net electric charge distribution from the STAR collaboration. This allows us to determine the freeze-out temperature and chemical potential from first principles. We also show continuum-extrapolated results for ratios of higher order fluctuations of baryon number. These will allow to test the consistency of the approach, by comparing them to the corresponding experimental data (once they become available) and thus extracting the freeze-out parameters in an independent way.
  • We consider the effect of strong external electromagnetic fields on thermodynamic observables in QCD, through lattice simulations with 1+1+1 flavors of staggered quarks at physical quark masses. Continuum extrapolated results are presented for the light quark condensates and for their tensor polarizations, as functions of the temperature and the magnetic field. We find the light condensates to undergo inverse magnetic catalysis in the transition region, in a manner that the transition temperature decreases with growing magnetic field. We also compare the results to other approaches and lattice simulations. Furthermore, we relate the tensor polarization to the spin part of the magnetic susceptibility of the QCD vacuum, and show that this contribution is diamagnetic.
  • Scale setting is of central importance in lattice QCD. It is required to predict dimensional quantities in physical units. Moreover, it determines the relative lattice spacings of computations performed at different values of the bare coupling, and this is needed for extrapolating results into the continuum. Thus, we calculate a new quantity, $w_0$, for setting the scale in lattice QCD, which is based on the Wilson flow like the scale $t_0$ (M. Luscher, JHEP 1008 (2010) 071). It is cheap and straightforward to implement and compute. In particular, it does not involve the delicate fitting of correlation functions at asymptotic times. It typically can be determined on the few per-mil level. We compute its continuum extrapolated value in 2+1-flavor QCD for physical and non-physical pion and kaon masses, to allow for mass-independent scale setting even away from the physical mass point. We demonstrate its robustness by computing it with two very different actions (one of them with staggered, the other with Wilson fermions) and by showing that the results agree for physical quark masses in the continuum limit.
  • We present 2+1 flavor Lattice QCD calculations of the nucleon scalar and tensor charges. Using the BMW clover-improved Wilson action with pion masses between 150 and 350 MeV and three source-sink separations between 0.9 and 1.4 fm, we achieve good control over excited-state contamination and extrapolation to the physical pion mass. As a consistency check, we also present results from calculations using unitary domain wall fermions with pion masses between 300 and 400 MeV, and using domain wall valence quarks and staggered sea quarks with pion masses between 300 and 600 MeV.
  • QCD thermodynamics is considered using Wilson fermions in the fixed scale approach. The temperature dependence of the renormalized chiral condensate, quark number susceptibility and Polyakov loop is measured at four lattice spacings allowing for a controlled continuum limit. The light quark masses are fixed to heavier than physical values in this first study. Finite volume effects are ensured to be negligible by using approriately large box sizes. The final continuum results are compared with staggered fermion simulations performed in the fixed N_t approach. The same continuum renormalization conditions are used in both approaches and the final results agree perfectly.
  • We determine the equation of state of QCD for nonzero chemical potentials via a Taylor expansion of the pressure. The results are obtained for N_f=2+1 flavors of quarks with physical masses, on various lattice spacings. We present results for the pressure, interaction measure, energy density, entropy density, and the speed of sound for small chemical potentials. At low temperatures we compare our results with the Hadron Resonance Gas model. We also express our observables along trajectories of constant entropy over particle number. A simple parameterization is given (the Matlab/Octave script parameterization.m, submitted to the arXiv along with the paper), which can be used to reconstruct the observables as functions of T and mu, or as functions of T and S/N.
  • We study QCD thermodynamics using two flavors of dynamical overlap fermions with quark masses corresponding to a pion mass of 350 MeV. We determine several observables on N_t=6 and 8 lattices. All our runs are performed with fixed global topology. Our results are compared with staggered ones and a nice agreement is found.
  • We use the Wilson flow to define the gauge anisotropy at a given physical scale. We demonstrate the use of the anisotropic flow by performing the tuning of the bare gauge anisotropy in the tree-level Symanzik action for several lattice spacings and target anisotropies. We use this method to tune the anisotropy parameters in full QCD, where we also exploit the diminishing effect of a well chosen smearing on the renormalization of the fermion anisotropy.
  • The effect of an external (electro)magnetic field on the finite temperature transition of QCD is studied. We generate configurations at various values of the quantized magnetic flux with $N_f=2+1$ flavors of stout smeared staggered quarks, with physical masses. Thermodynamic observables including the chiral condensate and susceptibility, and the strange quark number susceptibility are measured as functions of the field strength. We perform the renormalization of the studied observables and extrapolate the results to the continuum limit using $N_t=6,8$ and 10 lattices. We also check for finite volume effects using various lattice volumes. We find from all of our observables that the transition temperature $T_c$ significantly decreases with increasing magnetic field. This is in conflict with various model calculations that predict an increasing $T_c(B)$. From a finite volume scaling analysis we find that the analytic crossover that is present at B=0 persists up to our largest magnetic fields $eB \approx 1 \textmd{GeV}^2$, and that the transition strength increases mildly up to this $eB\approx1 \textmd{GeV}^2$.
  • The effect of an external magnetic field on the finite temperature transition of QCD is studied. We measure thermodynamic observables including the quark condensates and susceptibilities and the strange quark number susceptibility. We generate configurations at various values of the quantized magnetic flux with N_f=2+1 flavors of stout smeared staggered quarks at physical quark masses. We perform the renormalization of our observables and approach the continuum limit with N_t=6,8 and 10 lattices. We also check finite volume effects using various lattice volumes. Our main result is that the transition temperature significantly decreases with growing magnetic field, and that the transition remains an analytic crossover up to our largest external field eB ~ 1 GeV^2.
  • Indirect CP violation in K \rightarrow {\pi}{\pi} decays plays a central role in constraining the flavor structure of the Standard Model (SM) and in the search for new physics. For many years the leading uncertainty in the SM prediction of this phenomenon was the one associated with the nonperturbative strong interaction dynamics in this process. Here we present a fully controlled lattice QCD calculation of these effects, which are described by the neutral kaon mixing parameter B_K . We use a two step HEX smeared clover-improved Wilson action, with four lattice spacings from a\approx0.054 fm to a\approx0.093 fm and pion masses at and even below the physical value. Nonperturbative renormalization is performed in the RI-MOM scheme, where we find that operator mixing induced by chiral symmetry breaking is very small. Using fully nonperturbative continuum running, we obtain our main result B_K^{RI}(3.5GeV)=0.531(6)_{stat}(2)_{sys}. A perturbative 2-loop conversion yields B_K^{MSbar-NDR}(2GeV)=0.564(6)_{stat}(3)_{sys}(6)_{PT}, which is in good agreement with current results from fits to experimental data.
  • We study the spectra of heavy-light and heavy-heavy mesons containing charm quarks, including higher spin states. We use two sets of $N_f = 2 + 1$ gauge configurations, one set from QCDSF using the SLiNC action, and the other configurations from the Budapest-Marseille-Wuppertal collaboration, using the HEX smeared clover action. To extract information about the excited states, we choose a suitable basis of operators to implement the variational method.