MIT Center for Theoretical Physics — a Leinweber Institute Publications

We introduce a new mechanism for the simultaneous generation of baryon and dark matter asymmetries through ultraviolet-dominated freeze-in scatterings. The mechanism relies on heavy Majorana neutrinos that connect the visible Standard Model sector to a dark sector through the neutrino portal. Following reheating of the visible sector to a temperature well below the heavy neutrino masses, we show that 2-to-2 scattering processes can populate the dark sector and generate both baryon and dark matter asymmetries. In some parameter regions, the dominant source of baryon asymmetry can be charge transfer from the dark sector, a process we call dark wash-in. We also demonstrate that annihilation of the dark matter to massless states within the dark sector can deplete the symmetric population without destroying the net baryon charge to leave only an asymmetric dark matter abundance today. Depending on the specific model parameters, the observed baryon and dark matter abundances can be attained with heavy neutrino masses $M_N \gtrsim 10^{10}$ GeV, and dark matter masses in the range 0.1 GeV $\lesssim m_\chi \lesssim 10^3$ GeV if the dark matter relic abundance is mainly asymmetric and even lower masses if it is symmetric.

We revisit perturbative unitarity in scalar field inflation with a nonminimal coupling, with Higgs inflation serving as the most prominent example. Although such models are phenomenologically successful, it is critical to examine whether or not unitarity violations spoil their theoretical self-consistency. The analysis of these issues has so far typically relied on order-of-magnitude estimates of scattering amplitudes, which are appropriate for generic parameters. It is not evident that these methods apply to scenarios relying on a near-critical inflationary potential, for which an interplay of both small scalar self-couplings and nonminimal couplings could partially alleviate the unitarity issues. To allow for an exploration of this possibility, we consider the full $S$-matrix for the relevant scattering processes, taking into account important phase space volume factors, leading to a precise evaluation of the cut-off scale. In the single-field case, we demonstrate that near-criticality raises the cut-off scale considerably, compared to previous estimates. In the multifield case, momentum-dependent self-interactions in the kinetic sector lower the cut-off compared to the single-field case to a value comparable to but slightly larger than previous estimates. We carefully study both the single-field and multifield cases in metric and metric-affine (Palatini) formulations of gravity, as well as introduce a new phenomenologically viable model with a canonical kinetic term and a significantly raised cut-off, and discuss the importance of background field effects.

Simulating generic quantum states and dynamics is practically intractable using classical computers. However, certain special classes -- namely Clifford and matchgate circuits -- permit efficient computation. They provide invaluable tools for studying many-body physics, quantum chemistry, and quantum computation. While both play foundational roles across multiple disciplines, the origins of their tractability seem disparate, and their relationship remain unclear. A deeper understanding of such tractable classes could expand their scope and enable a wide range of new applications. In this work, we make progress toward the unified understanding of the Clifford and matchgate -- these two classes are, in fact, distinct special cases of a single underlying structure. Specifically, we introduce the 2D Quon language, which combines Majorana worldlines with their underlying spacetime topology to diagrammatically represent quantum processes and tensor networks. In full generality, the 2D Quon language is universal -- capable of representing arbitrary quantum states, dynamics, or tensor networks -- yet they become especially powerful in describing Clifford and matchgate classes. Each class can be efficiently characterized in a visually recognizable manner using the Quon framework. This capability naturally gives rise to several families of efficiently computable tensor networks introduced in this work: punctured matchgates, hybrid Clifford-matchgate-MPS, and ansatze generated from factories of tractable networks. All of these exhibit high non-Cliffordness, high non-matchgateness, and large bipartite entanglement entropy. We discuss a range of applications of our approach, from recovering well-known results such as the Kramers-Wannier duality and the star-triangle relation of the Ising model, to enabling variational optimization with novel ansatz states.

We show that the honeycomb tight-binding model hosts an exact microscopic avatar of its low-energy SU(2) valley symmetry and parity anomaly. Specifically, the SU(2) valley symmetry arises from a collection of conserved, integer-quantized charge operators that obey the Onsager algebra. Along with lattice reflection and time-reversal symmetries, this Onsager symmetry has a Lieb-Schultz-Mattis (LSM) anomaly that matches the parity anomaly in the IR. Indeed, we show that any local Hamiltonian commuting with these symmetries cannot have a trivial unique gapped ground state. We study the phase diagram of the simplest symmetric model and survey various deformations, including Haldane's mass term, which preserves only the Onsager symmetry. Our results place the parity anomaly in 2+1D alongside Schwinger's anomaly in 1+1D and Witten's SU(2) anomaly in 3+1D as 't Hooft anomalies that can arise from the Onsager symmetry on the lattice.

Effective Field Theory (EFT) modeling is expected to be a useful tool in the era of future higher-redshift galaxy surveys such as DESI-II and Spec-S5 due to its robust description of various large-scale structure tracers. However, large values of EFT bias parameters of higher-redshift galaxies could jeopardize the convergence of the perturbative expansion. In this paper we measure the bias parameters and other EFT coefficients from samples of two types of star-forming galaxies in the state-of-the-art MilleniumTNG and Astrid hydrodynamical simulations. Our measurements are based on the field-level EFT forward model that allows for precision EFT parameter measurements by virtue of cosmic variance cancellation. Specifically, we consider approximately representative samples of Lyman-break galaxies (LBGs) and Lyman-alpha emitters (LAEs) that are consistent with the observed (angular) clustering and number density of these galaxies at $z=3$. Reproducing the linear biases and number densities observed from existing LAE and LBG data, we find quadratic bias parameters that are roughly consistent with those predicted from the halo model coupled with a simple halo occupation distribution model. We also find non-perturbative velocity contributions (Fingers of God) of a similar size for LBGs to the familiar case of Luminous Red Galaxies. However, these contributions are quite small for LAEs despite their large satellite fraction values of up to $\sim 30\%$. Our results indicate that the effective momentum reach $k_{\rm{Max}}$ at $z=3$ for LAEs (LBGs) will be in the range $0.3-0.6 ~h\rm{Mpc}^{-1}$ ($0.2-0.8~h\rm{Mpc}^{-1}$), suggesting that EFT will perform well for high redshift galaxy clustering. This work provides the first step toward obtaining realistic simulation-based priors on EFT parameters for LAEs and LBGs.

To date, holographic calculations in strongly coupled plasma have provided separate descriptions for the rates of energy loss either for ultrarelativistic massless quarks and gluons or for infinitely massive quarks, with the latter calculation valid for $\sqrt{\gamma} < M/(\sqrt{\lambda}T)$, where $\gamma$ is the Lorentz boost factor for a heavy quark with velocity $v$ and mass $M$ moving through plasma with 't Hooft coupling $\lambda$ and temperature $T$. These two calculations should apply sequentially in the description of the energy loss of a heavy quark that starts out ultrarelativistic, loses energy, slows down, becomes non-relativistic at later times, and ultimately comes to rest and diffuses in the strongly coupled plasma. We provide an ansatz for uniquely incorporating both regimes to give an approximate but unified description of how a heavy quark that is initially ultrarelativistic loses energy all the way until it comes to rest. We implement this ansatz in the Hybrid Strong/Weak Coupling Model. With this new, consistent, treatment of heavy quark energy loss at strong coupling, we confront our predictions for the suppression and azimuthal anisotropies of D- and B-mesons, as well as B-tagged jets, with available experimental data.

Scalar fields are ubiquitous in theories of high-energy physics. In the context of cosmic inflation, this suggests the existence of spectator fields, which provide a subdominant source of energy density. We show that spectator fields boost the inflationary production of primordial black holes, with single-field ultra-slow roll evolution supplanted by a phase of evolution along the spectator direction, and primordial perturbations amplified by the resulting multifield dynamics. This generic mechanism is largely free from the severe fine-tuning that afflicts single-field inflationary PBH models.

The far-from equilibrium dynamics of the pre-hydrodynamic quark-gluon plasma (QGP) formed in heavy ion collisions can be characterized by distinct stages, during each of which the system loses some memory of its initial condition, until only the hydrodynamic modes remain. However, even though it has been repeatedly observed, finding intuitive physical explanations of how and why attractor behavior occurs has remained a challenge. The Adiabatic Hydrodynamization (AH) framework provides exactly such an explanation, showing that the attractor solution can be thought of as the ground state of an analog to quantum mechanical adiabatic evolution, provided we identify appropriate coordinate rescalings. Using the example of a simplified QCD kinetic theory in the small-angle scattering limit, we show how AH can explain both the early pre-hydrodynamic attractor and the later hydrodynamizing attractor in a longitudinally expanding gluon gas in a unified framework. By doing this, we provide a unified description of, and intuition for, all the stages of what in QCD would be bottom-up thermalization, starting from a pre-hydrodynamic attractor and ending with hydrodynamization.

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.