Abstract: The three-particle K-matrix, K-df,K-3, is a scheme-dependent quantity that parametrizes short-range three-particle interactions in the relativistic-field-theory three particle finite-volume formalism. In this work, we compute its value for systems of three pions at maximal isospin through next-to-leading order (NLO) in Chiral Perturbation Theory (ChPT). We compare the values to existing lattice QCD results and find that the agreement between lattice QCD data and ChPT in the first two coefficients of the threshold expansion of K-df,K-3 is significantly improved with respect to leading order once NLO effects are incorporated.

Abstract: Recently, the formalism needed to relate the finite-volume spectrum of systems of nondegenerate spinless particles has been derived. In this work we discuss a range of issues that arise when implementing this formalism in practice, provide further theoretical results that can be used to check the implementation, and make available codes for implementing the three-particle quantization condition. Specifically, we discuss the need to modify the upper limit of the cutoff function due to the fact that the left-hand cut in the scattering amplitudes for two nondegenerate particles moves closer to threshold; we describe the decomposition of the three-particle amplitude K-df,K-3 into the matrix basis used in the quantization condition, including both s and p waves, with the latter arising in the amplitude for two nondegenerate particles; we derive the threshold expansion for the lightest three-particle state in the rest frame up to O(1/L-5); and we calculate the leading-order predictions in chiral perturbation theory for K-df,K-3 in the pi(+)pi K-+(+) and pi+K+K+ systems. We focus mainly on systems with two identical particles plus a third that is different (“2+1” systems). We describe the formalism in full detail, and present numerical explorations in toy models, in particular checking that the results agree with the threshold expansion, and making a prediction for the spectrum of pi(+)pi K-+(+) levels using the two- and three-particle interactions predicted by chiral perturbation theory.

Abstract: We present a general method for deriving the energy shift of an interacting system of N spinless particles in a finite volume. To this end, we use the nonrelativistic effective field theory (NREFT), and match the pertinent low-energy constants to the scattering amplitudes. Relativistic corrections are explicitly included up to a given order in the 1/L expansion. We apply this method to obtain the ground state of N particles, and the first excited state of two and three particles to order L-6 in terms of the threshold parameters of the two- and three-particle relativistic scattering amplitudes. We use these expressions to analyze the N-particle ground state energy shift in the complex phi (4) theory.

Abstract: We derive relations between finite-volume matrix elements and infinite-volume decay amplitudes, for processes with three spinless, degenerate and either identical or non-identical particles in the final state. This generalizes the Lellouch-Luscher relation for two-particle decays and provides a strategy for extracting three-hadron decay amplitudes using lattice QCD. Unlike for two particles, even in the simplest approximation, one must solve integral equations to obtain the physical decay amplitude, a consequence of the nontrivial finite-state interactions. We first derive the result in a simplified theory with three identical particles, and then present the generalizations needed to study phenomenologically relevant three-pion decays. The specific processes we discuss are the CP-violating K -> 3 pi weak decay, the isospin-breaking eta -> 3 pi QCD transition, and the electromagnetic gamma (*) -> 3 pi amplitudes that enter the calculation of the hadronic vacuum polarization contribution to muonic g – 2.

Abstract: We propose a modification of the Hybrid Monte Carlo (HMC) algorithm that overcomes the topological freezing of a two-dimensional U(1) gauge theory with and without fermion content. This algorithm includes reversible jumps between topological sectors – winding steps – combined with standard HMC steps. The full algorithm is referred to as winding HMC (wHMC), and it shows an improved behaviour of the autocorrelation time towards the continuum limit. We find excellent agreement between the wHMC estimates of the plaquette and topological susceptibility and the analytical predictions in the U(1) pure gauge theory, which are known even at finite beta. We also study the expectation values in fixed topological sectors using both HMC and wHMC, with and without fermions. Even when topology is frozen in HMC – leading to significant deviations in topological as well as non-topological quantities – the two algorithms agree on the fixed-topology averages. Finally, we briefly compare the wHMC algorithm results to those obtained with master-field simulations of size L similar to 8 x 10(3).