Center for Theoretical Physics Publications

A longstanding issue in the study of quantum chromodynamics (QCD) is its behavior at nonzero baryon density, which has implications for many areas of physics. The path integral has a complex integrand when the quark chemical potential is nonzero and therefore has a sign problem, but it also has a generalized \(\mathcal PT\) symmetry. We review some new approaches to \(\mathcal PT\)-symmetric field theories, including both analytical techniques and methods for lattice simulation. We show that \(\mathcal PT\)-symmetric field theories with more than one field generally have a much richer phase structure than their Hermitian counterparts, including stable phases with patterning behavior. The case of a \(\mathcal PT\)-symmetric extension of a \(\phi^4\) model is explained in detail. The relevance of these results to finite density QCD is explained, and we show that a simple model of finite density QCD exhibits a patterned phase in its critical region.

Algorithms based on normalizing flows are emerging as promising machine learning approaches to sampling complicated probability distributions in a way that can be made asymptotically exact. In the context of lattice field theory, proof-of-principle studies have demonstrated the effectiveness of this approach for scalar theories, gauge theories, and statistical systems. This work develops approaches that enable flow-based sampling of theories with dynamical fermions, which is necessary for the technique to be applied to lattice field theory studies of the Standard Model of particle physics and many condensed matter systems. As a practical demonstration, these methods are applied to the sampling of field configurations for a two-dimensional theory of massless staggered fermions coupled to a scalar field via a Yukawa interaction.

Direct searches for new particles at colliders have traditionally been factorized into model proposals by theorists and model testing by experimentalists. With the recent advent of machine learning methods that allow for the simultaneous unfolding of all observables in a given phase space region, there is a new opportunity to blur these traditional boundaries by performing searches on unfolded data. This could facilitate a research program where data are explored in their natural high dimensionality with as little model bias as possible. We study how the information about physics beyond the Standard Model is preserved by full phase space unfolding using an important physics target at the Large Hadron Collider (LHC): exotic Higgs boson decays involving hadronic final states. We find that if the signal cross section is high enough, information about the new physics is visible in the unfolded data. We will show that in some cases, quantifiably all of the information about the new physics is encoded in the unfolded data. Finally, we show that there are still many cases when the unfolding does not work fully or precisely, such as when the signal cross section is small. This study will serve as an important benchmark for enhancing unfolding methods for the LHC and beyond.

We propose an exact model of anyon ground states including higher Landau levels, and use it to obtain fractionally quantized Hall states at filling fractions \(\nu=p/(p(m-1)+1)\) with \(m\) odd, from integer Hall states at \(\nu=p\) through adiabatic localization of magnetic flux. For appropriately chosen two-body potential interactions, the energy gap remains intact during the process. The construction hence establishes the existence of incompressible states at these fillings.

We examine the vector-portal inelastic dark matter (DM) model with DM mass \(m_\chi\) and dark photon mass \(m_{A'}\), in the `forbidden dark matter' regime where \(1 \lesssim m_{A'}/m_\chi \lesssim 2\), carefully tracking the dark sector temperature throughout freezeout. The inelastic nature of the dark sector relaxes the stringent cosmic microwave background (CMB) and self-interaction constraints compared to symmetric DM models. We determine the CMB limits on both annihilations involving excited states and annihilation into \(e^+e^-\) through initial-state-radiation of an \(A'\), as well as limits on the DM self-scattering, which proceeds at the one-loop level. The unconstrained parameter space serves as an ideal target for accelerator \(A'\) searches, and provides a DM self-interaction cross section that is large enough to observably impact small-scale structure.

Massive scalar fields provide excellent dark matter candidates, whose dynamics are often explored analytically and numerically using nonrelativistic Schrödinger-Poisson (SP) equations in a cosmological context. In this paper, starting from the nonlinear and fully relativistic Klein-Gordon-Einstein (KGE) equations in an expanding universe, we provide a systematic framework for deriving the SP equations, as well as relativistic corrections to them, by integrating out `fast modes' and including nonlinear metric and matter contributions. We provide explicit equations for the leading-order relativistic corrections, which provide insight into deviations from the SP equations as the system approaches the relativistic regime. Upon including the leading-order corrections, our equations are applicable beyond the domain of validity of the SP system, and are simpler to use than the full KGE case in some contexts. As a concrete application, we calculate the mass-radius relationship of solitons in scalar dark matter and accurately capture the deviations of this relationship from the SP system towards the KGE one.

A high energy muon collider can provide new and complementary discovery potential to the LHC or future hadron colliders. Leptoquarks are a motivated class of exotic new physics models, with distinct production channels at hadron and lepton machines. We study a vector leptoquark model at a muon collider with \(\sqrt{s} = 3, 14\) TeV within a set of both UV and phenomenologically motivated flavor scenarios. We compute which production mechanism has the greatest reach for various values of the leptoquark mass and the coupling between leptoquark and Standard Model fermions. We find that we can probe leptoquark masses up to an order of magnitude beyond \(\sqrt{s}\) with perturbative couplings. Additionally, we can also probe regions of parameter space unavailable to flavor experiments. In particular, all of the parameter space of interest to explain recent low-energy anomalies in B meson decays would be covered even by a \(\sqrt{s} = 3\) TeV collider.

We study the contributions of supersymmetric models with a \(U(1)\) horizontal symmetry and only spontaneous CP breaking to various lepton flavor observables, such as \(\mu \to e\gamma\) and the electron electric dipole moment. We show that both a horizontal symmetry and a lack of explicit CP violation can alleviate the existing bounds from such observables. The undetermined \(\mathcal{O}(1)\) coefficients in such mass matrix models muddle the interpretation of the bounds from various flavor observables. To overcome this, we define a new fine-tuning measure for different observables in such setups. This allows us to study how naturally the observed IR flavor observables can emerge from a given mass matrix model. We use our flavor-naturalness measure in study of our supersymmetric models and quantify the degree of fine tuning required by the bounds from various lepton flavor observables at each mass scale of sleptons, neutralinos, and charginos. We also comment on implications of the observed discrepancy in the anomalous magnetic moment of the muon for our supersymmetric models.

We discuss why AI is hard and why physics is simple. We discuss how physical intuition and the approach of theoretical physics can be brought to bear on the field of artificial intelligence and specifically machine learning. We suggest that the underlying project of machine learning and the underlying project of physics are strongly coupled through the principle of sparsity, and we call upon theoretical physicists to work on AI as physicists. As a first step in that direction, we discuss an upcoming book on the principles of deep learning theory that attempts to realize this approach.

We study the effect of a first-order phase transition in a confining \(SU(N)\) dark sector with heavy dark quarks. The baryons of this sector are the dark matter candidate. During the confinement phase transition the heavy quarks are trapped inside isolated, contracting pockets of the deconfined phase, giving rise to a second stage of annihilation that dramatically suppresses the dark quark abundance. The surviving abundance is determined by the local accidental asymmetry in each pocket. The correct dark matter abundance is obtained for \(\mathcal{O}(1-100)\) PeV dark quarks, above the usual unitarity bound.

Parton distribution functions (PDFs) and light-cone distribution amplitudes (LCDAs) are central non-perturbative objects of interest in high-energy inelastic and elastic scattering, respectively. As a result, an ab-initio determination of these objects is highly desirable. In this paper we present theoretical details required for the calculation of the moments of the PDF and LCDA using a heavy-quark operator product expansion method. This article also provides results of the relevant Wilson coefficients for the unpolarized PDF, helicity PDF and meson LCDA which are an essential ingredient in this approach. In addition, we discuss the use of our method to determine the dependence on the parton momentum fraction \(\xi\) in the PDFs and LCDAs through a matching condition between these quantities and the relevant hadronic amplitudes.

We present an approach to measure the Milky Way (MW) potential using the angular accelerations of stars in aggregate as measured by astrometric surveys like Gaia. Accelerations directly probe the gradient of the MW potential, as opposed to indirect methods using e.g. stellar velocities. We show that end-of-mission Gaia stellar acceleration data may be used to measure the potential of the MW disk at approximately 3\(\sigma\) significance and, if recent measurements of the solar acceleration are included, the local dark matter density at ~2\(\sigma\) significance. Since the significance of detection scales steeply as \(t^{5/2}\) for observing time \(t\), future surveys that include angular accelerations in the astrometric solutions may be combined with Gaia to precisely measure the local dark matter density and shape of the density profile.

The Wilson action for Euclidean lattice gauge theory defines a positive-definite transfer matrix that corresponds to a unitary lattice gauge theory time-evolution operator if analytically continued to real time. Hoshina, Fujii, and Kikukawa (HFK) recently pointed out that applying the Wilson action discretization to continuum real-time gauge theory does not lead to this, or any other, unitary theory and proposed an alternate real-time lattice gauge theory action that does result in a unitary real-time transfer matrix. The character expansion defining the HFK action is divergent, and in this work we apply a path integral contour deformation to obtain a convergent representation for U(1) HFK path integrals suitable for numerical Monte Carlo calculations. We also introduce a class of real-time lattice gauge theory actions based on analytic continuation of the Euclidean heat-kernel action. Similar divergent sums are involved in defining these actions, but for one action in this class this divergence takes a particularly simple form, allowing construction of a path integral contour deformation that provides absolutely convergent representations for U(1) and SU(N) real-time lattice gauge theory path integrals. We perform proof-of-principle Monte Carlo calculations of real-time U(1) and SU(3) lattice gauge theory and verify that exact results for unitary time evolution of static quark-antiquark pairs in (1 + 1)D are reproduced.

A remarkable feature of chaos in many-body quantum systems is the existence of a bound on the quantum Lyapunov exponent. An important question is to understand what is special about maximally chaotic systems which saturate this bound. Here we provide further evidence for the `hydrodynamic' origin of chaos in such systems, and discuss hallmarks of maximally chaotic systems. We first provide evidence that a hydrodynamic effective field theory of chaos we previously proposed should be understood as a theory of maximally chaotic systems. We then emphasize and make explicit a signature of maximal chaos which was only implicit in prior literature, namely the suppression of exponential growth in commutator squares of generic few-body operators. We provide a general argument for this suppression within our chaos effective field theory, and illustrate it using SYK models and holographic systems. We speculate that this suppression indicates that the nature of operator scrambling in maximally chaotic systems is fundamentally different to scrambling in non-maximally chaotic systems. We also discuss a simplest scenario for the existence of a maximally chaotic regime at sufficiently large distances even for non-maximally chaotic systems.

We present predictions for the gluon-fusion Higgs \(p_T\) spectrum at N\(^3\)LL\('\)+N\(^3\)LO order including fiducial cuts as required by experimental measurements at the Large Hadron Collider. Integrating the spectrum, we predict for the first time the total fiducial cross section to third order (N\(^3\)LO) and improved by resummation. The N\(^3\)LO correction is enhanced by cut-induced logarithmic effects and is not reproduced by the inclusive N\(^3\)LO correction times a lower-order acceptance. These are the highest-order predictions of their kind achieved so far at a hadron collider.

Autoencoders have been proposed as a powerful tool for model-independent anomaly detection in high-energy physics. The operating principle is that events which do not belong to the space of training data will be reconstructed poorly, thus flagging them as anomalies. We point out that in a variety of examples of interest, the connection between large reconstruction error and anomalies is not so clear. In particular, for data sets with nontrivial topology, there will always be points that erroneously seem anomalous due to global issues. Conversely, neural networks typically have an inductive bias or prior to locally interpolate such that undersampled or rare events may be reconstructed with small error, despite actually being the desired anomalies. Taken together, these facts are in tension with the simple picture of the autoencoder as an anomaly detector. Using a series of illustrative low-dimensional examples, we show explicitly how the intrinsic and extrinsic topology of the dataset affects the behavior of an autoencoder and how this topology is manifested in the latent space representation during training. We ground this analysis in the discussion of a mock "bump hunt" in which the autoencoder fails to identify an anomalous "signal" for reasons tied to the intrinsic topology of \(n\)-particle phase space.

Pionless effective field theory in a finite volume (FVEFT\(_{\pi\!/}\)) is investigated as a framework for the analysis of multi-nucleon spectra and matrix elements calculated in lattice QCD (LQCD). By combining FVEFT\(_{\pi\!/}\) with the stochastic variational method, the spectra of nuclei with atomic number \(A\in\{2,3\}\) are matched to existing finite-volume LQCD calculations at heavier-than-physical quark masses corresponding to a pion mass \(m_\pi=806\) MeV, thereby enabling infinite-volume binding energies to be determined using infinite-volume variational calculations. Based on the variational wavefunctions that are constructed in this approach, the finite-volume matrix elements of various local operators are computed in FVEFT\(_{\pi\!/}\) and matched to LQCD calculations of the corresponding QCD operators in the same volume, thereby determining the relevant one and two-body EFT counterterms and enabling an extrapolation of the LQCD matrix elements to infinite volume. As examples, the scalar, tensor, and axial matrix elements are considered, as well as the magnetic moments and the isovector longitudinal momentum fraction.

We begin developing tools to compute off-shell string amplitudes with the recently proposed hyperbolic string vertices of Costello and Zwiebach. Exploiting the relation between a boundary value problem for Liouville's equation and a monodromy problem for a Fuchsian equation, we construct the local coordinates around the punctures for the generalized hyperbolic three-string vertex and investigate their various limits. This vertex corresponds to the general pants diagram with three boundary geodesics of unequal lengths. We derive the conservation laws associated with such vertex and perform sample computations. We note the relevance of our construction to the calculations of the higher-order string vertices using the pants decomposition of hyperbolic Riemann surfaces.

The axial charge of the triton is investigated using lattice quantum chromodynamics (QCD). Extending previous work at heavier quark masses, calculations are performed using three ensembles of gauge field configurations generated with quark masses corresponding to a pion mass of 450 MeV. Finite-volume energy levels for the triton, as well as for the deuteron and diproton systems, are extracted from analysis of correlation functions computed on these ensembles, and the corresponding energies are extrapolated to infinite volume using finite-volume pionless effective field theory (FVEFT). It is found with high likelihood that there is a compact bound state with the quantum numbers of the triton at these quark masses. The axial current matrix elements are computed using background field techniques on one of the ensembles and FVEFT is again used to determine the axial charge of the proton and triton. A simple quark mass extrapolation of these results and earlier calculations at heavier quark masses leads to a value of the ratio of the triton to proton axial charges at the physical quark masses of \(g_A^{^{3}{\rm H}}/g_A^p=0.91\substack{+0.07 \\ -0.09}\). This result is consistent with the ratio determined from experiment and prefers values less than unity (in which case the triton axial charge would be unmodified from that of the proton), thereby demonstrating that QCD can explain the modification of the axial charge of the triton.

I review recent applications of the open quantum system framework in the understanding of quarkonium suppression in heavy-ion collisions, which has been used as a probe of the quark-gluon plasma for decades. The derivation of the Lindblad equations for quarkonium in both the quantum Brownian motion and the quantum optical limits and their semiclassical counterparts is explained. The hierarchy of time scales assumed in the derivation is justified from the separation of energy scales in nonrelativistic effective field theories of QCD. Physical implications of the open quantum system approach are also discussed. Finally, I list some open questions for future studies.