Center for Theoretical Physics Publications

We derive an effective field theory (EFT) for cosmological Lyman alpha forest fluctuations valid for the power spectrum at the one-loop order. The ``bottom-up'' EFT expansion at the level of the transmitted flux is identical to the line-of-sight dependent bias model first derived by Desjacques et al. We confirm this result by a ``top-down'' derivation based on the exponential map of the optical depth field. Specifically, we show that the combination of the exponential map and conditions of renormalizability generates the same EFT expansion as the ``bottom-up'' approach. In passing, we point out inconsistencies of the tree-level perturbative expansion of the exponential map without counterterms. To facilitate practical applications, we generalize the FFTLog method for efficient calculations of one-loop integrals from new line-of-sight dependent operators. Finally, we compare the one-loop EFT model against data from the Sherwood hydrodynamic simulation. The theory fits the data with sub percent accuracy up to \(k= \) 3 \(h\)Mpc\(^{-1}\) at \(z= 2.8\) for both 3D and 1D correlations. Our model can be readily used for cosmological full-shape analyses of Lyman alpha forest data.

We investigate the physics of one-dimensional symmetry protected topological (SPT) phases protected by symmetries whose symmetry generators exhibit spatial modulation. We focus in particular on phases protected by symmetries with linear (i.e., dipolar), quadratic and exponential modulations. We present a simple recipe for constructing modulated SPT models by generalizing the concept of decorated domain walls to spatially modulated symmetry defects, and develop several tools for characterizing and classifying modulated SPT phases. A salient feature of modulated symmetries is that they are generically only present for open chains, and are broken upon the imposition of periodic boundary conditions. Nevertheless, we show that SPT order is present even with periodic boundary conditions, a phenomenon we understand within the context of an object we dub a ``bundle symmetry''. In addition, we show that modulated SPT phases can avoid a certain no-go theorem, leading to an unusual algebraic structure in their matrix product state descriptions.

We study holographic entropy inequalities and their structural properties by making use of a judicious grouping of terms into certain multipartite information quantities. This allows us to recast cumbersome entropic expressions into much simpler ones which share interestingly rigid structures. By performing a systematic search over some of these structures, we are able to discover more than 300 novel entropy inequalities for six parties, thereby demonstrating that these recastings provide a fruitful generating technique for uncovering new holographic entropy inequalities. In attempting to interpret the corresponding sign-definite quantities as correlation measures, we also obtain a no-go result: the superbalance property of holographic entropy inequalities turns out to preclude them from being monotonic under partial tracing. In the process, we also comment on the geometrical significance of multipartite information quantities and present various structural relations amongst them.

We explore a self-interacting neutrino cosmology in which neutrinos experience a delayed onset of free-streaming. We use the effective field theory of large-scale structure (LSS) to model matter distribution on mildly non-linear scales within the self-interacting neutrino cosmology for the first time. We perform the first combined likelihood analysis of BOSS full-shape galaxy clustering, weak lensing, and Lyman-\(\alpha\) forest measurements, together with the cosmic microwave background (CMB) data from Planck. We find that the full data set strongly favors presence of a flavor-universal neutrino self-interaction, with a characteristic energy scale of order \(10\) MeV. The preference is at the \(>5\sigma\) level and is primarily driven by the Lyman-\(\alpha\) forest measurements and, to a lesser extent, the weak lensing data from DES. The self-interacting neutrino model reduces both the Hubble tension and the \(S_8\) tension between different cosmological data sets, but it does not fully resolve either. Finally, we note a preference for a non-zero sum of neutrino masses at the level of \(\sim 0.3\) eV under this model, consistent with previous bounds. These results call for further investigation in several directions, and may have significant implications for neutrino physics and for future new-physics searches with galaxy surveys.

We derive model-independent quantization conditions on the axion couplings (sometimes known as the anomaly coefficients) to the Standard Model gauge group \([SU(3)\times SU(2)\times U(1)_Y]/\mathbb{Z}_q\) with \(q=1,2,3,6\). Using these quantization conditions, we prove that any QCD axion model to the right of the \(E/N=8/3\) line on the \(g_{a\gamma\gamma}\)-\(m_a\) plot must necessarily face the axion domain wall problem in a post-inflationary scenario. We further demonstrate the higher-group and non-invertible global symmetries in the Standard Model coupled to a single axion. These generalized global symmetries lead to universal bounds on the axion string tension and the monopole mass. If axion were discovered in the future, our quantization conditions could be used to constrain the global form of the Standard Model gauge group.

The 3.5 keV line is a purported emission line observed in galaxies, galaxy clusters, and the Milky Way whose origin is inconsistent with known atomic transitions and has previously been suggested to arise from dark matter decay. We systematically re-examine the bulk of the evidence for the 3.5 keV line, attempting to reproduce six previous analyses that found evidence for the line. Surprisingly, we only reproduce one of the analyses; in the other five we find no significant evidence for a 3.5 keV line when following the described analysis procedures on the original data sets. For example, previous results claimed 4\(\sigma\) evidence for a 3.5 keV line from the Perseus cluster; we dispute this claim, finding no evidence for a 3.5 keV line. We find evidence for background mismodeling in multiple analyses. We show that analyzing these data in narrower energy windows diminishes the effects of mismodeling but returns no evidence for a 3.5 keV line. We conclude that there is little robust evidence for the existence of the 3.5 keV line. Some of the discrepancy of our results from those of the original works may be due to the earlier reliance on local optimizers, which we demonstrate can lead to incorrect results. For ease of reproducibility, all code and data are publicly available.

Consider the \(n!\) different unitary matrices that permute \(n\)\(d\)-dimensional quantum systems. If \(d\geq n\) then they are linearly independent. This paper discusses a sense in which they are approximately orthogonal (with respect to the Hilbert-Schmidt inner product) if \(d\gg n^2\), or, in a different sense, if \(d\gg n\). Previous work had shown pairwise approximate orthogonality of these matrices, but here we show a more collective statement, quantified in terms of the operator norm distance of the Gram matrix to the identity matrix. This simple point has several applications in quantum information and random matrix theory: (1) showing that random maximally entangled states resemble fully random states, (2) showing that Boson sampling output probabilities resemble those from Gaussian matrices, (3) improving the Eggeling-Werner scheme for multipartite data hiding, (4) proving that the product test of Harrow-Montanaro cannot be performed using LOCC without a large number of copies of the state to be tested, (5) proving that the purity of a quantum state also cannot be efficiently tested using LOCC, and (6, published separately) helping prove that poly-size random quantum circuits are poly-designs.

The two statistical methods, namely the frequentist and the Bayesian methods, are both commonly used for probabilistic inference in many scientific situations. However, it is not straightforward to interpret the result of one approach in terms of the concepts of the other. In this paper we explore the possibility of finding a Bayesian significance for the frequentist's main object of interest, the \(p\)-value, which is the probability assigned to the proposition -- which we call the {\it extremity proposition}-- that a measurement will result in a value that is at least as extreme as the value that was actually obtained. To make contact with the frequentist language, the Bayesian can choose to update probabilities based on the {\it extremity proposition}, which is weaker than the standard Bayesian update proposition, which uses the actual observed value. We then show that the posterior probability (or probability density) of a theory is equal to the prior probability (or probability density) multiplied by the ratio of the \(p\)-value for the data obtained, given that theory, to the mean \(p\)-value -- averaged over all theories weighted by their prior probabilities. Thus, we provide frequentist answers to Bayesian questions. Our result is generic -- it does not rely on restrictive assumptions about the situation under consideration or specific properties of the likelihoods or the priors.

The first stars are expected to form through molecular-hydrogen (H\(_2\)) cooling, a channel that is especially sensitive to the thermal and ionization state of gas, and can thus act as a probe of exotic energy injection from decaying or annihilating dark matter (DM). Here, we use a toy halo model to study the impact of DM-sourced energy injection on the H\(_2\) content of the first galaxies, and thus estimate the threshold mass required for a halo to form stars at high redshifts. We find that currently allowed DM models can significantly change this threshold, producing both positive and negative feedback. In some scenarios, the extra heating of the gas raises the halo mass required for collapse, whereas in others, energy injection lowers the threshold by increasing the free-electron fraction and catalyzing H\(_2\) formation. The direction of the effect can be redshift-dependent. We also bracket the uncertainties from self-shielding of halos from Lyman-Werner radiation. Hence, exotic energy injection can both delay and accelerate the onset of star formation; we show how this can impact the timing of 21cm signals at cosmic dawn. We encourage detailed simulation follow-ups in the most promising regions of parameter space identified in this work.

We present a comprehensive approach to quantum simulations at both zero and finite temperatures, employing a quantum information theoretic perspective and utilizing the Clifford + \(k\)Rz transformations. We introduce the "quantum magic ladder", a natural hierarchy formed by systematically augmenting Clifford transformations with the addition of Rz gates. These classically simulable similarity transformations allow us to reduce the quantumness of our system, conserving vital quantum resources. This reduction in quantumness is essential, as it simplifies the Hamiltonian and shortens physical circuit-depth, overcoming constraints imposed by limited error correction. We improve the performance of both digital and analog quantum computers on ground state and finite temperature molecular simulations, not only outperforming the Hartree-Fock solution, but also achieving consistent improvements as we ascend the quantum magic ladder. By facilitating more efficient quantum simulations, our approach enables near-term and early fault-tolerant quantum computers to address novel challenges in quantum chemistry.

Starting with the two-derivative limit of $D=2$ string theory, we explore the space of $T$-duality invariant ${\alpha }^{\prime }$ corrections, a space that contains a point representing the fully ${\alpha }^{\prime }$-corrected classical string theory. Using a parametrization introduced by Gasperini and Veneziano we obtain black hole solutions in this theory space. We prove that the dual of a solution with a regular horizon must have a curvature singularity. We find regions in the theory space where the black hole is deformed while preserving the horizon and the singularity, and regions where no black hole appears to exist. Furthermore, we find subregions in this theory space, probably not containing string theory, in which the black hole geometry exhibits a horizon leading to an interior that, having no singularity in the metric, curvature or dilaton, is a regular cosmology.

The influx of massive amounts of data from current and upcoming cosmological surveys necessitates compression schemes that can efficiently summarize the data with minimal loss of information. We introduce a method that leverages the paradigm of self-supervised machine learning in a novel manner to construct representative summaries of massive datasets using simulation-based augmentations. Deploying the method on hydrodynamical cosmological simulations, we show that it can deliver highly informative summaries, which can be used for a variety of downstream tasks, including precise and accurate parameter inference. We demonstrate how this paradigm can be used to construct summary representations that are insensitive to prescribed systematic effects, such as the influence of baryonic physics. Our results indicate that self-supervised machine learning techniques offer a promising new approach for compression of cosmological data as well its analysis.

Witten's open string field theory with a generalized version of stubs is reformulated as a cubic string field theory using an auxiliary string field. The gauge symmetries and equations of motion as well as the associative algebra of the resulting theory are investigated. Integrating out either the original or auxiliary field is shown to recover the conventional cubic theory. Our analysis demonstrates that deformations due to the stubs can be described as a homotopy transfer purely in the context of strong deformation retract. We also discuss to what extent the vertex regions resulting from stubs provide a model for the elementary interactions of closed string field theory.

Understanding the dynamics of bound state formation is one of the fundamental questions in confining quantum field theories such as Quantum Chromodynamics (QCD). One hadronization mechanism that has garnered significant attention is the breaking of a string initially connecting a fermion and an anti-fermion. Deepening our understanding of real-time string-breaking dynamics with simpler, lower dimensional models like the Schwinger model can improve our understanding of the hadronization process in QCD and other confining systems found in condensed matter and statistical systems. In this paper, we consider the string-breaking dynamics within the Schwinger model and investigate its modification inside a thermal medium, treating the Schwinger model as an open quantum system coupled to a thermal environment. Within the regime of weak coupling between the system and environment, the real-time evolution of the system can be described by a Lindblad evolution equation. We analyze the Liouvillian gaps of this Lindblad equation and the time dependence of the system's von Neumann entropy. We observe that the late-time relaxation rate decreases as the environment correlation length increases. Moreover, when the environment correlation length is infinite, the system exhibits two steady states, one in each of the sectors with definite charge-conjugation-parity (CP) quantum numbers. For parameter regimes where an initial string breaks in vacuum, we observe a delay of the string breaking in the medium, due to kinetic dissipation effects. Conversely, in regimes where an initial string remains intact in vacuum time evolution, we observe string breaking (melting) in the thermal medium. We further discuss how the Liouvillian dynamics of the open Schwinger model can be simulated on quantum computers and provide an estimate of the associated Trotter errors.

The theoretical description of observables at collider experiments relies on factorization theorems separating perturbative dynamics from universal non-perturbative matrix elements. Despite significant recent progress in extending these factorization theorems to increasingly differential jet substructure observables, the focus has been primarily on infrared safe observables sensitive only to correlations in the energy of final state hadrons. However, significant information about the dynamics of the underlying collision is encoded in how energy is correlated between hadrons of different quantum numbers. In this paper we extend the class of calculable correlations by deriving factorization theorems for a broad class of correlations, \(\langle \mathcal{E}_{R_1}(n_1) \cdots \mathcal{E}_{R_k}(n_k) \rangle\), between the energy flux carried by hadrons specified by quantum numbers, \(R_1, \cdots, R_k\). We show that these factorization theorems involve moments of a new class of universal non-perturbative functions, the "joint track functions", which extend the track function formalism to describe the fraction of energy carried by hadrons of multiple quantum numbers arising from the fragmentation of quarks or gluons. We study the general properties of these functions, and then apply this to the specific case of joint track functions for positive and negative electromagnetic charges. We extract these from parton shower Monte Carlo programs and use them to calculate correlations in electromagnetically charged energy flux. We additionally propose and study a C-odd \(\mathcal{E}_{\mathcal{Q}}\) detector, which results in a qualitatively distinct scaling behavior compared to the standard energy correlators. Our formalism significantly extends the class of observables that can be computed at hadron colliders, with a wide range of applications from particle to nuclear physics.

The confining transition from asymptotically free partons to hadrons remains one of the most mysterious aspects of Quantum Chromodynamics. With the wealth of high quality jet substructure data we can hope to gain new experimental insights into the details of its dynamics. Jet substructure has traditionally focused on correlations, \(\langle \mathcal{E}(n_1) \mathcal{E}(n_2) \cdots \mathcal{E}(n_k) \rangle\), in the energy flux of hadrons. However, significantly more information about the confinement transition is encoded in how energy is correlated between hadrons with different quantum numbers, for example electric charge. In this Letter we develop the field theoretic formalism to compute general correlations, \(\langle \mathcal{E}_{R_1}(n_1) \mathcal{E}_{R_2}(n_2) \cdots\mathcal{E}_{R_k}(n_k) \rangle\), between the energy flux carried by hadrons with quantum numbers \(R_i\), by introducing new universal non-perturbative functions, which we term joint track functions. Using this formalism we show that the strong interactions introduce enhanced small angle correlations between opposite-sign hadrons, relative to like-sign hadrons, identifiable as an enhanced scaling of \(\langle \mathcal{E}_+(n_1) \mathcal{E}_-(n_2) \rangle\) relative to \(\langle \mathcal{E}_+(n_1) \mathcal{E}_+(n_2) \rangle\). We are also able to compute the scaling of a \(C\)-odd three-point function, \(\langle \mathcal{E}_\mathcal{Q}(n_1) \mathcal{E}_\mathcal{Q}(n_2) \mathcal{E}_\mathcal{Q}(n_3) \rangle\). Our results greatly extend the class of systematically computable jet substructure observables, pushing perturbation theory deeper into the parton to hadron transition, and providing new observables to understand the dynamics of confinement.

The continued success of the jet substructure program will require widespread use of tracking information to enable increasingly precise measurements of a broader class of observables. The recent reformulation of jet substructure in terms of energy correlators has simplified the incorporation of universal non-perturbative matrix elements, so called "track functions", in jet substructure calculations. These advances make it timely to understand how these universal non-perturbative functions can be extracted from hadron collider data, which is complicated by the use jet algorithms. In this paper we introduce a new class of jet functions, which we call (semi-inclusive) track jet functions, which describe measurements of the track energy fraction in identified jets. These track jet functions can be matched onto the universal track functions, with perturbatively calculable matching coefficients that incorporate the jet algorithm dependence. We perform this matching, and present phenomenological results for the charged energy fraction in jets at the LHC and EIC/HERA at collinear next-to-leading logarithmic accuracy. We show that higher moments of the charged energy fraction directly exhibit non-linear Lorentzian renormalization group flows, allowing the study of these flows with collider data. Our factorization theorem enables the extraction of universal track functions from jet measurements, opening the door to their use for a precision jet substructure program.

We present an algorithm to compute correlation functions for systems with the quantum numbers of many identical mesons from lattice quantum chromodynamics (QCD). The algorithm is numerically stable and allows for the computation of \(n\)-pion correlation functions for \(n \in \{ 1, \dots, N\}\) using a single \(N \times N\) matrix decomposition, improving on previous algorithms. We apply the algorithm to calculations of correlation functions with up to 6144 \(\pi^+\)s using two ensembles of gauge field configurations generated with quark masses corresponding to a pion mass \(m_\pi = 170\) MeV and spacetime volumes of \((4.4^3\times 8.8)\ {\rm fm}^4\) and \((5.8^3\times 11.6)\ {\rm fm}^4\). We also discuss statistical techniques for the analysis of such systems, in which the correlation functions vary over many orders of magnitude. In particular, we observe that the many-pion correlation functions are well approximated by log-normal distributions, allowing the extraction of the energies of these systems. Using these energies, the large-isospin-density, zero-baryon-density region of the QCD phase diagram is explored. A peak is observed in the energy density at an isospin chemical potential \(\mu_I\sim 1.5 m_\pi\), signalling the transition into a Bose-Einstein condensed phase. The isentropic speed of sound in the medium is seen to exceed the ideal-gas (conformal) limit (\(c_s^2\leq 1/3\)) over a wide range of chemical potential before falling towards the asymptotic expectation at \(\mu_I\sim 15 m_\pi\). These, and other thermodynamic observables, indicate that the isospin chemical potential must be large for the system to be well described by an ideal gas or perturbative QCD.

A lattice QCD computation of the coupled channel \(\pi\Sigma\unicode{x2013}\bar{K}N\) scattering amplitudes in the \(\Lambda(1405)\) region is detailed. Results are obtained using a single ensemble of gauge field configurations with \(N_{\rm f} = 2+1\) dynamical quark flavors and \(m_{\pi} \approx 200\) MeV and \(m_K\approx487\) MeV. Hermitian correlation matrices using both single baryon and meson-baryon interpolating operators for a variety of different total momenta and irreducible representations are used. Several parametrizations of the two-channel scattering \(K\)-matrix are utilized to obtain the scattering amplitudes from the finite-volume spectrum. The amplitudes, continued to the complex energy plane, exhibit a virtual bound state below the \(\pi\Sigma\) threshold and a resonance pole just below the \(\bar{K}N\) threshold.

This work presents a determination of the quark Collins-Soper kernel, which relates transverse-momentum-dependent parton distributions (TMDs) at different rapidity scales, using lattice Quantum Chromodynamics (QCD). This is the first lattice QCD calculation of the kernel at quark masses corresponding to a close-to-physical value of the pion mass, with next-to-next-leading logarithmic matching to TMDs from the corresponding lattice-calculable distributions, and includes a complete analysis of systematic uncertainties arising from operator mixing. The kernel is extracted at transverse momentum scales \(240\,\text{MeV}\lesssim q_{T}\lesssim 1.6\,\text{GeV}\) with a precision sufficient to begin to discriminate between different phenomenological models in the non-perturbative region.