# 1+1+1 flavor QCD+QED simulation at the physical point

###### Abstract

We present the results of 1+1+1 flavor QCD+QED simulation at the physical point, in which the dynamical quark effects in QED and the up-down quark mass difference are incorporated by the reweighting technique. The physical quark masses together with the lattice spacing are determined with , , and as physical inputs. Calculations are carried out using a set of 2+1 flavor QCD configurations near the physical point generated by the non-perturbatively -improved Wilson quark action and the Iwasaki gauge action at on a lattice. We evaluate the values of the up, down and strange quark masses individually with non-perturbative QCD renormalization.

###### pacs:

11.15.Ha, 12.38.-t, 12.38.GcPACS-CS Collaboration

## I Introduction

The ultimate goal of lattice QCD is to understand the non-perturbative dynamics of strong interactions at the physical values of the up, down and strange quark masses. Recently we have achieved this goal assuming artificial isospin symmetry with degenerate up and down quark masses PACS-CS_1 ; PACS-CS_2 . In nature, however, the isospin symmetry is broken due to the up-down quark mass difference and their electric charge difference. Their effects are observed in mass splittings among isospin multiplets of light hadrons, e.g. , . The magnitude of splittings is tiny yet important since, e.g., it is this difference which guarantees the stability of proton. Thus, our next step should be 1+1+1 QCD+QED lattice simulation at the physical point incorporating the isospin breaking effects.

A pioneering lattice QCD+QED simulation was carried out in Refs. Duncan_1 ; Duncan_2 employing quenched approximation both for QCD and QED, where the photon fields and the gluon fields are independently generated by conventional Monte Carlo methods and they are superimposed. This work was followed by several studies in quenched QED on quenched and dynamical SU(3) gauge configurationsNamekawa -Glaessle . These studies have shown that the mass splittings among the isospin multiplets of hadrons came out qualitatively consistent with the experimental ones. Nevertheless, dynamical quark effects in QED should be properly evaluated to obtain convincing results. An inclusion of dynamical quark effects in QED was first attempted in Ref. Duncan_3 by calculating the ratio of the quark determinant of QCD+QED to that of QCD on rather small lattices using the reweighting technique. Recently results for the low-energy constants including dynamical quark effects both in QCD and QED are reportedIshikawa .

In this article, we present results for 1+1+1 flavor QCD+QED simulation at the physical point. The photon fields are superimposed on the 2+1 flavor QCD configurations near the physical point generated by the PACS-CS CollaborationPACS-CS_2 employing the non-perturbatively -improved Wilson quark action and the Iwasaki gauge action at , corresponding to the lattice spacing of fm, on a lattice. We employ the reweighting technique to incorporate the dynamical quark effects in QEDDuncan_3 and adjust the up and down quark masses to their physical values independentlyPACS-CS_2 ; Hasenfratz ; RBC_UKQCD . We make several improvements on the generation of photon fields and evaluation of the reweighting factors.

This paper is organized as follows. Section II is devoted to explain the generation of the photon fields and the calculational method of the reweighting factors to incorporate the dynamical quark effects in QED. In Sec. III we present the results for 1+1+1 flavor QCD+QED simulation at the physical point. Our conclusions are summarized in Sec. IV.

## Ii Reweighting factors for full QED

Given a quenched U(1) gauge configuration superimposed on a 2+1 flavor QCD configuration, the reweighting factor is evaluated with a set of independent Gaussian random noises :

(1) | |||||

(2) |

where is the Wilson-Dirac matrix with a hopping parameter corresponding to the physical point and the physical electric charge multiplied by a fraction specific to each quark flavor. We employ the determinant breakup technique in terms of both the electric charge and the hopping parameters to reduce the noise fluctuations PACS-CS_2 ; Hasenfratz ; RBC_UKQCD . It is clear that the stochastic evaluation of det works efficiently when the ratio is close to the unit matrix, for which local fluctuations of the photon field should be suppressed. We employ the following procedure to prepare the photon field.

We first generate the photon field on a lattice, which is twice finer than the QCD lattice of , employing a non-compact pure gauge action in the Coulomb gauge and an appropriate treatment of zero modesKogut ; Duncan_1 ; Blum_1 ; Hayakawa :

(3) | |||||

with

(4) |

where is the unit vector in direction on the QCD lattice. The coefficient is chosen such that the averaged action density on the QCD lattice defined by

(5) |

with

(6) | |||||

(7) |

is numerically equal to that of the photon field generated with on the QCD lattice. To reduce local fluctuations we average the photon fields over independent paths inside the hypercube on the QED lattice embedded in the unit hypercube on the QCD lattice:

(8) | |||||

Finally, we construct the compact link variable with the physical electric charge of each quark.

We have checked that the averaging procedure does not cause variation of the electric charge by calculating the electromagnetic mass splittings among charged and neutral pseudoscalar mesons in the flavor non-singlet sector: Results obtained with the averaged photon field are consistent with those with the conventional ones within error bars. We have also found that this procedure sizably reduces the additive quark mass shift due to the QED effects for the Wilson-clover quark action we employ.

Further reduction of the fluctuations in the stochastic evaluation of the reweighting factors is achieved by calculating the combined reweighting factor for all three light quarks with a single set of noise as expressed in Eqs. (1) and (2)Ishikawa . It should be noted that the odd terms in for are not forbidden by the symmetries in the QCD+QED mixed system. The leading contribution in the fluctuations is expected to be of order of . Since the total charge for the up, down and strange quarks vanishes, , we can expect a cancellation of the leading contribution, , if the quark masses are identical. Even though the degeneracy is broken, the partial cancellation of the terms is expected to reduce the fluctuations in the stochastic evaluation. In order to illustrate the cancellation we plot the noise dependence of the odd part in for () with the first breakup interval to in Fig. 1, where the stochastic evaluation is applied to individual quarks. The odd part in is extracted from the calculations with . We observe that the fluctuation for the up quark is opposite to those for the down and strange quarks as a function of the noise, where we expect that the next-to-leading contribution of should be negligibly small.

## Iii 1+1+1 flavor QCD+QED simulation at the physical point

### iii.1 Parameters

Our QCD+QED calculation is based on a set of 2+1 flavor QCD configurations near the physical point generated with the non-perturbatively -improved Wilson quark action and the Iwasaki gauge action at on a lattice by the PACS-CS Collaboration PACS-CS_2 . The hopping parameters are . The number of configurations is 80 over 2000 MD time units. Using the U(1) link variables on a lattice, which are constructed as described in Sec. II, 2+1 flavor QCD partition function is reweighted into 1+1+1 flavor QCD+QED partition function with the hopping parameters () and the physical electric charge . The values of the hopping parameters are determined by , and masses as the physical inputs. We employ 400 determinant breakups for the QED part and 26 for the fine tuning of the hopping parameters to the physical point. Each piece of the divided determinant is evaluated stochastically with 12 sets of noises. The quark matrix inversions are obtained by a block solver algorithm developed in Ref. Nakamura , which reduce the computational cost by a factor of 3 to 4 compared with non-block solvers.

### iii.2 Results

In Fig. 2 we plot the configuration dependence of the reweighting factor from 2+1 flavor QCD to 1+1+1 flavor QCD+QED at the physical point, which are normalized by the configuration average. The largest value is about 10 which is the same order as found in Ref. Ishikawa .

Our results [MeV] | Experiment [MeV] | |
---|---|---|

137.7(8.0) | 139.57018(35) | |

492.3(4.7) | 493.677(16) | |

497.4(3.7) | 497.614(24) | |

input | 1672.45(29) |

This work | Refs. PACS-CS_3 ; PACS-CS_2 | Refs. PACS-CS_3 ; PACS-CS_1 | |
---|---|---|---|

[MeV] | 2.57(26)(07) | ||

[MeV] | 3.68(29)(10) | ||

[MeV] | 83.60(58)(2.23) | 86.7(2.3) | 87.7(3.1) |

[MeV] | 3.12(24)(08) | 2.78(27) | 3.05(12) |

0.698(51) | |||

26.8(2.0) | 31.2(2.7) | 28.78(40) |

We calculate the hadron correlators with point and smeared sources and a local sink employing the operators given in Refs. PACS-CS_1 ; PACS-CS_2 . The smeared source is constructed with an exponential smearing function with and . To reduce the statistical errors we employ 16 source points located at every 4 time slices with different spatial sites. Figure 3 shows effective mass plots for and propagators with the smeared source. The three horizontal lines are extended over the fit ranges, which are for the pseudoscalar mesons and for baryon, and show the resulting hadron masses with one standard deviation error band. The statistical errors are estimated with the jackknife method with a binsize of 100 MD time units. The fit results are summarized in Table 1. The scale is determined as GeV using the baryon mass as input.

In Fig. 4 we show the hadron mass ratios, , and , as a function of the number of noises employed for the stochastic evaluation of each determinant breakup. The ratios quickly reach a plateau around the experimental values (red broken lines) with small number of noises. This shows not only that twelve noises are sufficient for evaluation of each determinant breakup but also that we have achieved a successful tuning to the physical point.

Figure 5 shows the ratio of to propagators whose time dependence is expected as

(9) |

where we assume . The negative slope indicates . The fit result with gives 4.54(1.09) MeV, which is consistent with the experimental value 3.937(28) MeV Exp within the error bar.

In Table 2 we summarize the up, down and strange quark masses renormalized at GeV in the continuum scheme together with the values from our previous workPACS-CS_3 ; PACS-CS_2 ; PACS-CS_1 . Note that the present results employ the non-perturbative renormalization factor obtained with the Schrödinger functional schemePACS-CS_3 , and for a proper comparison, the previous results for the light quark masses in Refs. PACS-CS_1 ; PACS-CS_2 are updated with the use of the non-perturbative renormalization factor in Table 14 of Ref. PACS-CS_3 . The first errors are statistical and the second ones are associated with computation of the non-perturbative renormalization factors. We neglect the QED corrections to the renormalization factor, whose contribution would be at a level of 1% or less. A potentially large finite size effects in QED are discussed in Ref. Hayakawa . The authors in Ref. Blum_2 investigate the finite volume effects on the up, down and strange quark masses in (1.8 fm) and (2.7 fm) boxes, the latter of which is nearly equal to the spatial volume in this work, employing chiral perturbation theory including electromagnetic interactions. Based on their analyses we estimate the systematic errors for our results due to the finite size corrections of QED to be 13.50%, 2.48%, 0.07% for the up, down, strange quark masses, respectively.

## Iv Conclusions

We have presented the results for 1+1+1 flavor QCD+QED simulation at the physical point. The dynamical quark effects in QED and the up-down quark mass difference are incorporated by the reweighting technique. We have made several improvements to reduce the fluctuations in the stochastic evaluation of the reweighting factors. The physical values for the up, down and strange quark masses are determined by using and as input, and quoted in the scheme at GeV. A direct investigation of the finite size effects of QED by changing the spatial volume is left as a future work.

###### Acknowledgements.

Numerical calculations for the present work have been carried out under the “Interdisciplinary Computational Science Program” in Center for Computational Sciences, University of Tsukuba. This work is supported in part by Grants-in-Aid of the Ministry of Education, Culture, Sports, Science and Technology-Japan (Nos. 18104005, 20340047, 20540248, 21340049, 22244018, and 24540250), and Grant-in-Aid for Scientific Research on Innovative Areas (No. 2004: 20105001, 20105002, 20105003, and 20105005).## References

## References

- (1) S. Aoki et al. [PACS-CS Collaboration], Phys. Rev. D79, 034503 (2009).
- (2) S. Aoki et al. [PACS-CS Collaboration], Phys. Rev. D81, 074503 (2010).
- (3) A. Duncan, E. Eichten and H. Thacker, Phys. Rev. Lett. 76, 3894 (1996).
- (4) A. Duncan, E. Eichten and H. Thacker, Phys. Lett. B409, 387 (1997).
- (5) Y. Namekawa and Y. Kikukawa, PoS LATTICE2005, 090 (2006).
- (6) N. Yamada, T. Blum, M. Hayakawa and T. Izubuchi [RBC Collaboration], PoS LATTICE2005, 092 (2006).
- (7) T. Doi, T. Blum, M. Hayakawa, T. Izubuchi and N. Yamada, PoS LATTICE2006, 174 (2006).
- (8) T. Blum, T. Doi, M. Hayakawa, T. Izubuchi and N. Yamada, Phys. Rev. D76, 114508 (2007).
- (9) S. Basak et al. [MILC Collaboration], PoS LATTICE2008, 127 (2008).
- (10) R. Zhou, T. Blum, T. Doi, M. Hayakawa, T. Izubuchi and N. Yamada, PoS LATTICE2008, 131 (2008).
- (11) R. Zhou and S. Uno, PoS LATTICE2009, 182 (2009).
- (12) T. Blum, R. Zhou, T. Doi, M. Hayakawa, T. Izubuchi, S. Uno and N. Yamada, Phys. Rev. D82, 094508 (2010).
- (13) A. Portelli et al. [Budapest-Marseille-Wuppertal Collaboration], PoS LATTICE2010, 121 (2010).
- (14) A. Torok et al. [MILC Collaboration], PoS LATTICE2010, 127 (2010).
- (15) A. Portelli, PoS LATTICE2011, 136 (2011).
- (16) B. Gläßle, PoS LATTICE2011, 282 (2011).
- (17) A. Duncan, E. Eichten and R. Sedgewick, Phys. Rev. D71, 094509 (2005).
- (18) T. Ishikawa, T. Blum, M. Hayakawa, T. Izubuchi, C. Jung and R. Zhou, arXiv:1202.6018 [hep-lat].
- (19) A. Hasenfratz, R. Hoffmann and S Schaefer, Phys. Rev. D78, 014515 (2008).
- (20) Y. Aoki et al. [RBC Collaboration and UKQCD Collaboration], Phys. Rev. D83, 074508 (2011).
- (21) E. Dagotto, A. Kocic and J.B. Kogut, Nucl. Phys. B331, 500 (1990).
- (22) S. Uno and M. Hayakawa, Prog. Theor. Phys. 120, 413 (2008).
- (23) Y. Nakamura, K.-I. Ishikawa, Y. Kuramashi, T. Sakurai and H. Tadano, Comput. Phys. Commun. 183, 34 (2012).
- (24) K. Nakamura et al. [Particle Data Group], J. Phys. G 37, 075021 (2010).
- (25) S. Aoki et al. [PACS-CS Collaboration], JHEP 1008, 101 (2010).