“…Φ → e iq Φ a/fa Φ (Φ = ψ, H u,d ), (2.5) and subsequently integrate out all massive fields heavier than µ = O(1) GeV, where q Φ is the PQ charge of Φ (in the normalization convention with q σ = 1) for the linearly realized U (1) PQ , and ψ stands for all chiral fermions in the model. Then the axion-gluon coupling c G , which arises as a consequence of the axion-dependent field redefinition of ψ, corresponds to the coefficient of the U (1) PQ -SU (3) c -SU (3) c anomaly, while the couplings C u,d to the light quarks are determined by (i) a contribution from the axion-dependent fieldredefinition of {Q 1 , u c 1 , d c 1 }, (ii) the tree-level threshold correction from the axion mixing with the Z boson which is induced by the field redefinition of H u,d , and finally (iii) the radiative corrections caused by the gauge and Yukawa couplings in the model [33]. Putting these together, one finds…”