MIT Center for Theoretical Physics – a Leinweber Institute

We study the stochastic dynamics of heavy quarks propagating through the strongly coupled plasma of $\mathcal{N}=4$ supersymmetric Yang-Mills (SYM) theory at nonzero temperature in terms of the corresponding Kolmogorov equation, which correctly describes their kinetic equilibration and the non-Gaussian fluctuations in their momenta without having to restrict their velocity to the non-relativistic regime. Leveraging the heavy quark limit, we show that the evolution of the momentum space distribution function can be reformulated as a Hamilton-Jacobi problem, and therefore can be solved in terms of first-order ordinary differential equations. We solve these evolution equations in an infinite thermal plasma with a constant temperature for initial conditions specified by spherically symmetric heavy quark momentum distributions with a phenomenologically motivated shape that is steeply falling at large momentum. To highlight the distinctive features of the kinetic equilibration process, we compare their solutions with Fokker-Planck dynamics with the same drag coefficient and fluctuations chosen by hand to guarantee equilibration. We find qualitatively similar dynamics at small momentum, and very different dynamics at large momentum, where, much like in jet quenching phenomena, the steepness of the momentum distribution gives a larger relevance to unlikely events in which a heavy quark loses little momentum in a given time step. Such events are much less unlikely in the Kolmogorov evolution than in Fokker-Planck evolution with the same mean energy loss, meaning that equilibration at large momentum is significantly delayed. Our results provide a systematic description of heavy quarks propagating through strongly coupled plasma from the ultra-relativistic to the non-relativistic regime and point the way towards implementation in phenomenological studies.

We present the next-to-next-to-leading logarithmic (NNLL) collinear resummation of projected $N$-point energy correlators (ENCs) up to $N=6$, matched to fixed-order predictions at NLO, in both electron-positron annihilation and Higgs decay to gluons. The key new ingredient is the two-loop jet function for $N=4,5,6$, which we compute semi-analytically using Integration-by-Parts and differential equations. We further include the leading non-perturbative corrections for ENCs, described by two universal soft matrix elements $\overlineΩ_{1q},\overlineΩ_{1g}$ of order $Λ_{\rm QCD}$, whose evolution is governed by anomalous dimensions for $(N-1)$-point correlators. The matched distributions are compared with parton-shower simulations from Pythia8 and Herwig7, and we study the sensitivity of both the absolute spectra and their ratios to the two-point energy correlator under variations of $α_s$ and $\overlineΩ_{1q,1g}$. Our results show that higher-point projected energy correlators are now under quantitative control at NNLL accuracy, opening the door to future $α_s$ extractions with complementary systematics.

The reliable treatment of highly boosted hadrons is crucial for many lattice QCD applications. For all these cases the kinematically enhanced interpolators promise very significant improvements and are, therefore, ever more often used in recent calculations, especially for highly boosted mesons like the pion. Motivated by, e.g., the physics program of the future Electron-Ion Collider (EIC) in the US and Electron-Ion collider in China (EIcC), we systematically benchmark our code for the unpolarized isovector nucleon quark matrix elements extracted at large source-sink separations, where excited state artifacts are significantly suppressed. We find that the precision of the renormalized nucleon matrix elements is typically improved by an order of magnitude at momentum $P_z\sim2.5$ GeV. By comparing the results from three CLS ensembles with different lattice spacings $a$ but the same pion mass, we observe no statistically significant dependence on $a$ in the renormalized matrix elements at nearly identical values of $P_z$. These encouraging results suggest that the use of kinematically improved operators is highly advantageous for parton physics calculations and can be extended to a broader class of baryon observables, making them a promising candidate for a standard component of modern lattice QCD.

Primordial black holes (PBHs) in the asteroid-mass range, with typical masses $10^{17}\,{\rm g}\lesssim M \lesssim 10^{23}\,{\rm g}$, have drawn significant recent attention as a viable dark matter candidate. The peak frequencies of photons emitted via Hawking radiation from asteroid-mass PBHs range from infrared to $γ$-ray bands. We calculate expected local transit rates for extended PBH mass distributions which could comprise all the dark matter. We evaluate prospects for detecting Hawking-radiated photons from local PBH transits through the inner Solar System and from PBH explosions in the far outer edges of the Solar System. We consider several existing and proposed ground-based and space-based instruments sensitive to photons from the radio band to ultrahigh energy $γ$-rays. We find that proposed instruments, such as the AMEGO-X satellite, could reliably detect PBH transits within ${\it O} (0.1 \, {\rm AU})$ of the Earth, while the HAWC and LHAASO observatories are both sensitive to PBH explosions out to ${\it O}(0.1 \, {\rm pc})$ and ${\it O}(0.5 \, {\rm pc})$ respectively. We conclude by specifically considering potential companion electromagnetic signatures in the case of a PBH explosion about $10^3\,{\rm AU}$ from Earth, which has been suggested as a potential source for the $\sim 220 \, {\rm PeV}$ ultrahigh-energy KM3-230213A neutrino event observed by the KM3NeT collaboration in 2023. Whereas we find that the recent KM3NeT event would not have yielded detectable electromagnetic signals — due to its location on the sky, proposed distance from Earth, and the offline status of the HAWC observatory at that time — we demonstrate that future PBH explosions at comparable distances could yield measurable electromagnetic signals at Earth, depending on alignment of the PBH burst with detector fields of view.