Abstract: Cosmic rays with energies up to a few PeV are known to be accelerated within the Milky Way(1,2). Traditionally, it has been presumed that supernova remnants were the main source of these very-high-energy cosmic rays(3,4), but theoretically it is difficult to accelerate protons to PeV energies(5,6) and observationally there simply is no evidence of the remnants being sources of hadrons with energies above a few tens of TeV7,8. One possible source of protons with those energies is the Galactic Centre region(9). Here, we report observations of 1-100 TeV gamma rays coming from the 'Cygnus Cocoon'(10), which is a superbubble that surrounds a region of massive star formation. These gamma rays are likely produced by 10-1,000 TeV freshly accelerated cosmic rays that originate from the enclosed star-forming region Cyg OB2. Until now it was not known that such regions could accelerate particles to these energies. The measured flux likely originates from hadronic interactions. The spectral shape and the emission profile of the Cocoon changes from GeV to TeV energies, which reveals the transport of cosmic particles and historical activity in the superbubble.

Abstract: The perturbative approach to quantum field theories has made it possible to obtain incredibly accurate theoretical predictions in high-energy physics. Although various techniques have been developed to boost the efficiency of these calculations, some ingredients remain specially challenging. This is the case of multiloop scattering amplitudes that constitute a hard bottleneck to solve. In this paper, we delve into the application of a disruptive technique based on the loop-tree duality theorem, which is aimed at an efficient computation of such objects by opening the loops to nondisjoint trees. We study the multiloop topologies that first appear at four loops and assemble them in a clever and general expression, the (NMLT)-M-4 universal topology. This general expression enables to open any scattering amplitude of up to four loops, and also describes a subset of higher order configurations to all orders. These results confirm the conjecture of a factorized opening in terms of simpler known subtopologies, which also determines how the causal structure of the entire loop amplitude is characterized by the causal structure of its subtopologies. In addition, we confirm that the loop-tree duality representation of the (NMLT)-M-4 universal topology is manifestly free of noncausal thresholds, thus pointing towards a remarkably more stable numerical implementation of multiloop scattering amplitudes.

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.