MIT Center for Theoretical Physics — a Leinweber Institute Publications

Primordial black holes (PBHs) with masses between ${10}^{14}$ and ${10}^{20}$ kg are candidates to contribute a substantial fraction of the total dark matter abundance. When in orbit around the center of a star, which can possibly be a completely interior orbit, such objects would emit gravitational waves, as predicted by general relativity. In this work, we examine the gravitational wave signals emitted by such objects when they orbit typical stars, such as the Sun. We show that the magnitude of the waves that could eventually be detected on Earth from a possible PBH orbiting the Sun or a neighboring Sun-like star within our galaxy can be significantly stronger than those originating from a PBH orbiting a denser but more distant neutron star (NS). Such signals may be detectable by the LISA gravitational-wave detector. In addition, we estimate the contribution that a large collection of such PBH-star systems would make to the stochastic gravitational-wave background (SGWB) within a range of frequencies to which pulsar timing arrays are sensitive.

We present predictions for the DIS 1-jettiness event shape ${\tau }_1^b$, or DIS thrust, using the framework of Soft Collinear Effective Theory (SCET) for factorization, resummation of large logarithms, and rigorous treatment of nonperturbative power corrections, matched to fixed-order QCD away from the resummation region. Our predictions reach next-to-next-to-next-to-leading-logarithmic (N${}^{}$LL) accuracy in resummed perturbation theory, matched to $O\left({\alpha }_s^2\right)$ fixed-order QCD calculations obtained using the program NLOJet++. We include a rigorous treatment of hadronization corrections, which are universal across different event shapes and kinematic variables $x$ and $Q$ at leading power, and supplement them with a systematic scheme to remove $O\left({\Lambda }_{\text{QCD}}\right)$ renormalon ambiguities in their definition. The framework of SCET allows us to connect smoothly the nonperturbative, resummation, and fixed-order regions, whose relative importance varies with $x$ and $Q$, and to rigorously estimate theoretical uncertainties, across a broad range of $x$ and $Q$ covering existing experimental results from HERA as well as expected new measurements from the upcoming Electron-Ion-Collider (EIC). Our predictions will serve as an important benchmark for the EIC program, enabling the precise determination of the QCD strong coupling ${\alpha }_s$ and the universal nonperturbative first moment parameter ${\Omega }_1$.

We show that perturbative techniques inspired by effective field theory (EFT) can be used to "paint on" the 21~cm field during reionization using only the underlying linear density field. This procedure is accurate to within O(10%) on large scales and thus can be used to enlarge or "supersize" hydrodynamical simulations. In particular, the EFT provides a mapping between the linear density field and a nonlinear tracer field, both in real and redshift space. We show that this mapping can be reliably extracted from relatively small simulation volumes using the THESAN suite of simulations, which have a comoving volume of (95.5 Mpc)^3. Specifically, we show that if we fit the EFT coefficients in a small ~5% sub-volume of the simulation, we can accurately predict the 21cm field in the rest of the simulation given only the linear density field. We show that our technique is robust to different models of dark matter and differences in the sub-grid reionization modeling.

We explore the detection prospects for a minimal secluded dark matter model, where a fermionic dark matter particle interacts with the Standard Model (SM) via a kinetically mixed dark photon. We focus on scenarios where the dark photon decays visibly, making it a prime target for beam-dump experiments. In this model, the dark matter relic abundance can be achieved by a variety of mechanisms: freeze-in, out-of-equilibrium freeze-out, and secluded freeze-out. We demonstrate that the secluded freeze-out regime in the considered mass range is now entirely excluded by a combination of direct and indirect detection constraints. Moreover, we show that future direct detection and intensity frontier experiments offer complementary sensitivity to this minimal model in the parameter space where the hidden sector never enters equilibrium with the SM. In out-of-equilibrium freeze-out scenarios, nuclear-recoil direct detection experiments can still access signals above the neutrino fog that are mediated by dark photons that are too weakly coupled to be detected in future beam dump experiments. Meanwhile, future beam dump experiments provide a powerful probe of the freeze-in parameter space in this model, which is largely inaccessible to direct detection experiments. Notably, even in the absence of a future observation in direct detection experiments, a dark photon discovery remains possible at SHiP, DUNE, LHCb, and DarkQuest within this minimal dark matter model.

We pioneer the application of quantum information theory to experimentally distinguish between classes of hadronization models. We adapt the CHSH inequality to the fragmentation of a single parton to hadron pairs, a violation of which would rule out classical dynamics of hadronization altogether. Furthermore, we apply and extend the theory of quantum contextuality and local quantum systems to the neutral polarization of a single spin-1 hadronic system, namely the light constituents of excited Sigma baryons ${\Sigma }_{c,b}^{\ast }$ formed in the fragmentation of heavy quarks.

The ATLAS Collaboration has developed a method to analyze large-radius jets composed of skinny $R=0.2$ subjets in heavy-ion collisions. We first demonstrate that the measurements pioneered by ATLAS constrain the value of $L_{res}$, the resolution length of QGP -- and rule out any picture in which an entire parton shower loses energy coherently as a single entity. We then analyze the response of the medium to the passage of large-radius $R=2$ jets containing two skinny subjets in gamma-jet events. We introduce novel jet-shape observables that allow us to visualize how the internal structure of large-radius jets shapes the wakes they excite in the QGP. We find that even when the subjets are $\sim 0.8$ radians apart, the angular shape of the soft hadrons originating from their wake forms a single broad structure. Only when the two subjets are even farther apart are two sub-wakes revealed. We show that the way in which jet structure shapes the structure of jet-induced wakes can be visualized with similar clarity in experiments by using only low-$p_T$ hadrons. The observables we introduce offer a new and distinctive way of seeing jet-induced wakes -- and wake substructure -- in heavy-ion collision data.

The study of hadronic jets and their substructure at hadronic colliders is crucial for improving our understanding of QCD, and searching for new physics. As such, there has been a significant effort to improve their theoretical description. In the small radius limit, inclusive jet production exhibits a universal factorization, enabling the resummation of logarithms which greatly stabilizes theoretical predictions. In this paper, we show how to combine a recently introduced framework for small-$R$ resummation with the STRIPPER subtraction formalism for fragmentation, enabling next-to-next-to-leading order calculations of small-$R$ inclusive jet production for a wide variety of processes at the LHC. We extract the two-loop constants for the jet functions, enabling for the first time next-to-next-to-leading logarithmic resummation matched to next-to-next-to-leading order perturbative calculation. We compare with CMS data for small-$R$ jet production, and find that our results greatly improve the accuracy of the predictions at small-$R$, and stabilize the perturbative convergence and error estimates at larger $R$. Our approach is applicable to a wide class of jet substructure observables exhibiting similar factorization theorems, opening the door to an NNLO jet substructure program at the LHC.

The algebraic approach to quantum field theory focuses on the properties of local algebras, whereas the study of (possibly non-invertible) global symmetries emphasizes global aspects of the theory and spacetime. We study connections between these two perspectives by examining how either of two core algebraic properties -- "additivity" or "Haag duality" -- is violated in a 1+1D CFT or lattice model restricted to the symmetric sector of a general global symmetry. For the Verlinde symmetry of a bosonic diagonal RCFT, we find that additivity is violated whenever the symmetry algebra contains an invertible element, while Haag duality is violated whenever it contains a non-invertible element. We find similar phenomena for the Kramers-Wannier and Rep(D${}_8$) non-invertible symmetries on spin chains.

Recent work has shown that the (block) Lanczos algorithm can be used to extract approximate energy spectra and matrix elements from (matrices of) correlation functions in quantum field theory, and identified exact coincidences between Lanczos analysis methods and others. In this work, we note another coincidence: the Lanczos algorithm is equivalent to the well-known Rayleigh-Ritz method applied to Krylov subspaces. Rayleigh-Ritz provides optimal eigenvalue approximations within subspaces; we find that spurious-state filtering allows these optimality guarantees to be retained in the presence of statistical noise. We explore the relation between Lanczos and Prony's method, their block generalizations, generalized pencil of functions (GPOF), and methods based on the generalized eigenvalue problem (GEVP), and find they all fall into a larger "Prony-Ritz equivalence class", identified as all methods which solve a finite-dimensional spectrum exactly given sufficient correlation function (matrix) data. This equivalence allows simpler and more numerically stable implementations of (block) Lanczos analyses.

We present the first search for evidence of neutrino self-interaction with two new, state-of-the-art likelihoods for eBOSS Lyman-$\alpha$ data. These are an effective field theory (EFT) likelihood with priors from the Sherwood simulation suite, and a compressed likelihood derived from an emulator built using the PRIYA simulation suite. Previous analyses that combined Planck measurements with eBOSS Lyman-$\alpha$ likelihoods based on earlier simulations found a preference for neutrino self-interactions. In contrast, using either of the new eBOSS Lyman-$\alpha$ likelihoods, we find that a joint analysis with the cosmic microwave background (CMB) data from Planck prefers a negligible level of neutrino self-interaction, and derive new constraints on the neutrino self-coupling: ${log}_{10}\left(G_{eff}{MeV}^2\right)=-{5.57}_{-0.58}^{+0.98}$ for Planck + EFT Lyman-$\alpha$, and ${log}_{10}\left(G_{eff}{MeV}^2\right)=-{5.26}_{-1.49}^{+0.87}$ for Planck + PRIYA Lyman-$\alpha$, at 68% confidence. We also consider Planck in combination with DESI BAO data, and find that the latter does not provide significant constraining power for neutrino self-interactions.

We study the scaling of meson-meson scattering amplitudes with the number of colors, $N_{\text{c}}$. We use lattice calculations in a theory with $N_{\text{f}}=4$ degenerate flavors, with $N_{\text{c}}=3-6$ and pion mass $M_{\pi }\approx 560$ MeV. We focus on three different scattering channels, two of which have the same quantum numbers as some tetraquark candidates recently found at LHCb: the $T_{cs0}^0\left(2900\right)$, $T_{c\overline{s}0}^{++}\left(2900\right)$, $T_{c\overline{s}0}^0\left(2900\right)$ and $T_{c\overline{s}1}^0\left(2900\right)$ states. Finite-volume energies are extracted using a large set of operators, containing two-particle operators with the form of two pions or two vector mesons, and local tetraquark operators. The resulting energy spectra is used to constrain the infinite-volume scattering amplitude by means of Lüscher's quantization condition. We consider polynomial parametrizations of the phase shift, as well as one-loop chiral perturbation theory (ChPT) predictions. We find that our lattice results follow the expected $N_{\text{c}}$ scaling and are sensitive to subleading $N_{\text{c}}$ corrections. In addition, we constrain the scaling of different combinations of low-energy constants from matching to large $N_{\text{c}}$ ChPT. The results for the channel corresponding to a $({\pi }^{+}{D}_s^{+}-K^{+}{D}^{+})$ state show evidence of a virtual bound state with energy $E_{\text{virtual}}=1.63\left(10\right)M_{\pi }$ for $N_{\text{c}}=3$, while this pole disappears at $N_{\text{c}}>3$. This may be connected to the exotic states found in experiment.

The currently observing Dark Energy Spectroscopic Instrument (DESI) places sub-percent constraints on the Baryon Acoustic Oscillations (BAO) scaling parameters from the Lyman-$\alpha$ (Ly-$\alpha$) forest. However, no systematic error budget stemming from non-linearities in the 3D clustering of the Ly-$\alpha$ forest is included in the DESI-Ly-$\alpha$ analysis. In this work, we measure the size of the shift of the BAO peak using large Ly-$\alpha$ forest mocks produced on the $N$-body simulation suite \textsc{AbacusSummit}. Specifically, we measure the Ly-$\alpha$ auto-correlation and the Ly-$\alpha$-quasar cross-correlation functions. We use the DESI Ly-$\alpha$ forest fitting pipeline, \textsc{Vega}, with the publicly available covariance matrix from eBOSS DR16. To mitigate the noise, we adopt a linear control variates (LCV) technique, reducing the error bars by a factor of up to $\sim \sqrt{50}$ on large scales. From the auto-correlation, we detect a small positive shift in radial direction of $\Delta {\alpha }_{\parallel }=0.35\%$ at the 3$\sigma$ level and virtually no shift in the transverse direction, ${\alpha }_{\perp }$. From the cross-correlation, we see a similar shift to $\Delta {\alpha }_{\parallel }$, albeit with larger error bars, and a small negative shift, $\Delta {\alpha }_{\perp }=\sim$0.25%, at the 2$\sigma$ level. We also make a connection with the Ly-$\alpha$ forest effective field theory (EFT) framework and find that the one-loop EFT power spectrum yields unbiased measurements of the BAO shift parameters in radial and transverse direction for Ly-$\alpha$ auto and the Ly-$\alpha$-quasar cross-correlation measurements. When using the one-loop EFT framework, we find that we can recover the BAO parameters without a shift, which has important implications for future Ly-$\alpha$ forest analyses based on EFT. This work paves the way for the full-shape analysis of DESI and future surveys.

The pseudoscalar meson light-cone distribution amplitudes (LCDAs) are essential non-perturbative inputs for a range of high-energy exclusive processes in quantum chromodynamics. In this proceedings, progress towards a determination of the low Mellin moments of the pion and kaon LCDAs by the HOPE Collaboration is reported.

Effective field theory (EFT)-based full-shape analysis with simulation-based priors (SBPs) is a novel approach to galaxy clustering data analysis, which significantly boosts the constraining power by efficiently incorporating field-level simulation information from small scales. So far, SBPs have been mostly extracted from a large set of mock catalogs generated with the Halo Occupation Distribution (HOD) approach. We show that given a halo mass function model and assuming that the EFT parameters of halos depend only on the peak height, HOD-based priors can be computed analytically for the standard 7-parameter HOD model. We derive the relevant expressions for deterministic EFT parameters from the halo model. The halo model, however, fails to accurately describe stochastic EFT parameters, for which we use a physically motivated phenomenological prescription. We compare our analytic priors with simulations and find an excellent agreement. Our approach provides analytic insights into the physics behind the EFT parameters of galaxies, and allows one to reduce the computation time of SBPs to virtually nothing. As an application, we produce a set of 100,000 EFT parameters by sampling both HOD and cosmological models, and explicitly demonstrate that any cosmology dependence can be completely absorbed by relatively minor shifts of the HOD parameters. Our approach can be generalized to other variants of the HOD and cosmological models beyond $\Lambda$CDM.

The so-called metastability bound on the Higgs mass suggests that the smallness of the Higgs mass may be a byproduct of the metastability of the electroweak vacuum. A significantly strong bound requires new physics capable of lowering the scale where the Higgs quartic coupling turns negative through renormalization group effects, without destabilizing the electroweak vacuum entirely. We analyze in this context the low-scale Majoron model of neutrino masses, which automatically contains two key elements for a viable scenario: heavy fermions to lower the instability scale and a extended scalar sector to stabilize the potential and achieve realistic lifetimes for the electroweak vacuum. We show how the metastability bound can be generalized to theories with multiple scalars and present an efficient way of calculating the tunneling rate in such potentials. We also demonstrate that FCC will probe regions of the parameter space relevant for metastability: large regions of the fermionic sector at FCC-ee and some reach to the scalar sector at FCC-hh.

We provide a general prescription for gauging finite non-invertible symmetries in 1+1d lattice Hamiltonian systems. Our primary example is the Rep(D${}_8$) fusion category generated by the Kennedy-Tasaki transformation, which is the simplest anomaly-free non-invertible symmetry on a spin chain of qubits. We explicitly compute its lattice F-symbols and illustrate our prescription for a particular (non-maximal) gauging of this symmetry. In our gauging procedure, we introduce two qubits around each link, playing the role of "gauge fields" for the non-invertible symmetry, and impose novel Gauss's laws. Similar to the Kramers-Wannier transformation for gauging an ordinary $Z_2$, our gauging can be summarized by a gauging map, which is part of a larger, continuous non-invertible cosine symmetry.

Axion-like particles can couple to Standard Model gluons, electroweak gauge bosons, and massive fermions. A future multi-TeV muon collider provides a favorable environment to probe axion-like particles through multiple production channels, including vector boson fusion via electroweak gauge boson couplings and the top-associated production mediated by direct fermionic couplings. Motivated by the quality issue of the QCD axion, we focus on axion-like particles with masses and decay constants around the TeV scale. We explore how different axion-like particle couplings shape its production and decay modes, revealing a rich and intricate phenomenological landscape.

We study systems of two and three mesons composed of pions and kaons at maximal isospin using four CLS ensembles with $a\approx 0.063$fm, including one with approximately physical quark masses. Using the stochastic Laplacian-Heaviside method, we determine the energy spectrum of these systems including many levels in different momentum frames and irreducible representations. Using the relativistic two- and three-body finite-volume formalism, we constrain the two and three-meson K matrices, including not only the leading $s$ wave, but also $p$ and $d$ waves. By solving the three-body integral equations, we determine, for the first time, the physical-point scattering amplitudes for $3{\pi }^{+}$, $3K^{+}$, ${\pi }^{+}{\pi }^{+}{K}^{+}$ and $K^{+}{K}^{+}{\pi }^{+}$ systems. These are determined for total angular momentum $J^P=0^{-}$, $1^{+}$, and $2^{-}$. We also obtain accurate results for $2{\pi }^{+}$, ${\pi }^{+}{K}^{+}$, and $2K^{+}$ phase shifts. We compare our results to Chiral Perturbation Theory, and to phenomenological fits.

We use lattice QCD calculations of the finite-volume spectra of systems of two and three mesons to determine, for the first time, three-particle scattering amplitudes with physical quark masses. Our results are for combinations of ${\pi }^{+}$ and $K^{+}$, at a lattice spacing $a=0.063$fm, and in the isospin-symmetric limit. We also obtain accurate results for maximal-isospin two-meson amplitudes, with those for ${\pi }^{+}{K}^{+}$ and $2K^{+}$ being the first determinations at the physical point. Dense lattice spectra are obtained using the stochastic Laplacian-Heaviside method, and the analysis leading to scattering amplitudes is done using the relativistic finite-volume formalism. Results are compared to chiral perturbation theory and to phenomenological fits to experimental data, finding good agreement.