Abstract: We present the first unquenched, continuum limit, lattice QCD results for the matrix elements of the operators describing neutral kaon oscillations in extensions of the Standard Model. Owing to the accuracy of our calculation on Delta S = 2 weak Hamiltonian matrix elements, we are able to provide a refined Unitarity Triangle analysis improving the bounds coming from model independent constraints on New Physics. In our non-perturbative computation we use a combination of N-f = 2 maximally twisted sea quarks and Osterwalder-Seiler valence quarks in order to achieve both O(a)-improvement and continuum-like renormalization properties for the relevant four-fermion operators. The calculation of the renormalization constants has been performed non-perturbatively in the RI-MOM scheme. Based on simulations at four values of the lattice spacing and a number of quark masses we have extrapolated/interpolated our results to the continuum limit and physical light/strange quark masses.

Abstract: The duality relation between one-loop integrals and phase-space integrals, developed in a previous work, is extended to higher-order loops. The duality relation is realized by a modification of the customary +i0 prescription of the Feynman propagators, which compensates for the absence of the multiple-cut contributions that appear in the Feynman tree theorem. We rederive the duality theorem at one-loop order in a form that is more suitable for its iterative extension to higher-loop orders. We explicitly show its application to two-and three-loop scalar master integrals, and we discuss the structure of the occurring cuts and the ensuing results in detail.

Abstract: In this paper we study the nonperturbative structure of the SU(3) four-gluon vertex in the Landau gauge, concentrating on contributions quadratic in the metric. We employ an approximation scheme where “one-loop” diagrams are computed using fully dressed gluon and ghost propagators, and tree-level vertices. When a suitable kinematical configuration depending on a single momentum scale p is chosen, only two structures emerge: the tree-level four-gluon vertex, and a tensor orthogonal to it. A detailed numerical analysis reveals that the form factor associated with this latter tensor displays a change of sign (zero-crossing) in the deep infrared, and finally diverges logarithmically. The origin of this characteristic behavior is proven to be entirely due to the masslessness of the ghost propagators forming the corresponding ghost-loop diagram, in close analogy to a similar effect established for the three-gluon vertex. However, in the case at hand, and under the approximations employed, this particular divergence does not affect the form factor proportional to the tree-level tensor, which remains finite in the entire range of momenta, and deviates moderately from its naive tree-level value. It turns out that the kinematic configuration chosen is ideal for carrying out lattice simulations, because it eliminates from the connected Green's function all one-particle reducible contributions, projecting out the genuine one-particle irreducible vertex. Motivated by this possibility, we discuss in detail how a hypothetical lattice measurement of this quantity would compare to the results presented here, and the potential interference from an additional tensorial structure, allowed by Bose symmetry, but not encountered within our scheme.

Abstract: We present an implementation of the relativistic three-particle quantization condition including both s- and d-wave two-particle channels. For this, we develop a systematic expansion of the three-particle K matrix, K-df,K-3, about threshold, which is the generalization of the effective range expansion of the two-particle K matrix, K-2. Relativistic invariance plays an important role in this expansion. We find that d-wave two-particle channels enter first at quadratic order. We explain how to implement the resulting multichannel quantization condition, and present several examples of its application. We derive the leading dependence of the threshold three-particle state on the two-particle d-wave scattering amplitude, and use this to test our implementation. We show how strong two-particle d-wave interactions can lead to significant effects on the finite-volume three-particle spectrum, including the possibility of a generalized three-particle Efimov-like bound state. We also explore the application to the 3 pi(+) system, which is accessible to lattice QCD simulations, where we study the sensitivity of the spectrum to the components of K-df,K-3. Finally, we investigate the circumstances under which the quantization condition has unphysical solutions.

Abstract: The next to leading order chiral corrections to the SU(2) x SU(2) Gell-Mann-Oakes- Renner (GMOR) relation are obtained using the pseudoscalar correlator to five-loop order in perturbative QCD, together with new finite energy sum rules (FESR) incorporating polynomial, Legendre type, integration kernels. The purpose of these kernels is to suppress hadronic contributions in the region where they are least known. This reduces considerably the systematic uncertainties arising from the lack of direct experimental information on the hadronic resonance spectral function. Three different methods are used to compute the FESR contour integral in the complex energy (squared) s-plane, i.e. Fixed Order Perturbation Theory, Contour Improved Perturbation Theory, and a fixed renormalization scale scheme. We obtain for the corrections to the GMOR relation, delta(pi), the value delta(pi) = (6.2 +/- 1.6)%. This result is substantially more accurate than previous determinations based on QCD sum rules; it is also more reliable as it is basically free of systematic uncertainties. It implies a light quark condensate < 0 vertical bar(u) over baru vertical bar 0 > similar or equal to < 0 vertical bar(d) over bard vertical bar 0 > < 0 vertical bar(q) over barq vertical bar 0 >vertical bar(2GeV) = (-267 +/- 5MeV)(3). As a byproduct, the chiral perturbation theory (unphysical) low energy constant H-2(r) is predicted to be H-2(r)(nu(X) = M-p) = -(5.1 +/- 1.8) x10(-3), or H-2(r) (nu(X) = M-eta) = -(5.7 +/- 2.0) x10(-3).