Center for Theoretical Physics Publications

Recently, the formalism needed to relate the finite-volume spectrum of systems of nondegenerate spinless particles has been derived. In this work we discuss a range of issues that arise when implementing this formalism in practice, provide further theoretical results that can be used to check the implementation, and make available codes for implementing the three-particle quantization condition. Specifically, we discuss the need to modify the upper limit of the cutoff function due to the fact that the left-hand cut in the scattering amplitudes for two nondegenerate particles moves closer to threshold; we describe the decomposition of the three-particle amplitude \(\mathcal{K}_\text{df,3}\) into the matrix basis used in the quantization condition, including both \(s\) and \(p\) waves, with the latter arising in the amplitude for two nondegenerate particles; we derive the threshold expansion for the lightest three-particle state in the rest frame up to \(\mathcal{O}(1/L^5)\); and we calculate the leading-order predictions in chiral perturbation theory for \(\mathcal{K}_\text{df,3}\) in the \(\pi^+\pi^+K^+\) and \(\pi^+K^+K^+\) systems. We focus mainly on systems with two identical particles plus a third that is different (``2+1'' systems). We describe the formalism in full detail, and present numerical explorations in toy models, in particular checking that the results agree with the threshold expansion, and making a prediction for the spectrum of \(\pi^+\pi^+K^+\) levels using the two- and three-particle interactions predicted by chiral perturbation theory.

The design of parametric quantum circuits (PQCs) for efficient use in variational quantum simulations (VQS) is subject to two competing factors. On one hand, the set of states that can be generated by the PQC has to be large enough to contain the solution state. Otherwise, one may at best find the best approximation of the solution restricted to the states generated by the chosen PQC. On the other hand, the PQC should contain as few parametric quantum gates as possible to minimize noise from the quantum device. Thus, when designing a PQC one needs to ensure that there are no redundant parameters. The dimensional expressivity analysis discussed in these proceedings is a means of addressing these counteracting effects. Its main objective is to identify independent and redundant parameters in the PQC. Using this information, superfluous parameters can be removed and the dimension of the space of states that are generated by the PQC can be computed. Knowing the dimension of the physical state space then allows us to deduce whether or not the PQC can reach all physical states. Furthermore, the dimensional expressivity analysis can be implemented efficiently using a hybrid quantum-classical algorithm. This implementation has relatively small overhead costs both for the classical and quantum part of the algorithm and could therefore be used in the future for on-the-fly circuit construction. This would allow for optimized circuits to be used in every loop of a VQS rather than the same PQC for the entire VQS. These proceedings review and extend work in [1, 2].

We study general freeze-out scenarios where an arbitrary number of initial and final dark matter particles participate in the number-changing freeze-out interaction. We consider a simple sector with two particle species undergoing such a thermal freeze-out; one of the relics is stable and gives rise to the dark matter today, while the other one decays to the Standard Model, injecting significant entropy into the thermal bath that dilutes the dark matter abundance. We show that this setup can lead to a stable relic population that reproduces the observed dark matter abundance without requiring weak scale masses or couplings. The final dark matter abundance is estimated analytically. We carry out this calculation for arbitrary temperature dependence in the freeze-out process and identify the viable dark matter mass and cross section ranges that explain the observed dark matter abundance. This setup can be used to open parameter space for both heavy (above the unitarity bound) or light (sub-GeV) dark matter candidates. We point out that the best strategy for probing most parts of our parameter space is to look for signatures of an early matter-dominant epoch.

Estimating the free energy, as well as other thermodynamic observables, is a key task in lattice field theories. Recently, it has been pointed out that deep generative models can be used in this context. Crucially, these models allow for the direct estimation of the free energy at a given point in parameter space. This is in contrast to existing methods based on Markov chains which generically require integration through parameter space. In this contribution, we will review this novel machine-learning-based estimation method. We will in detail discuss the issue of mode collapse and outline mitigation techniques which are particularly suited for applications at finite temperature.

We present an error mitigation scheme which corrects readout errors on Noisy Intermediate-Scale Quantum (NISQ) computers [1,2]. After a short review of applying the method to one qubit, we proceed to discuss the case when correlations between different qubits occur. We demonstrate how the readout error can be mitigated in this case. By performing experiments on IBMQ hardware, we show that such correlations do not have a strong effect on the results, justifying to neglect them.

Standard sampling algorithms for lattice QCD suffer from topology freezing (or critical slowing down) when approaching the continuum limit, thus leading to poor sampling of the distinct topological sectors. I will present a modified Hamiltonian Monte Carlo (HMC) algorithm that triggers topological sector jumps during the assembly of Markov chain of lattice configurations. We study its performance in the 2D Schwinger model and compare it to alternative methods, such as fixing topology or master field. We then briefly discuss the difficulties of the algorithm in a SU(2) gauge model in 4D.

We present results for \(B_{(s)}\)- and \(D_{(s)}\)-meson semileptonic decays from ongoing calculations by the Fermilab Lattice and MILC Collaborations. Our calculation employs the highly improved staggered quark (HISQ) action for both sea and valence quarks and includes several ensembles with physical-mass up, down, strange, and charm quarks and lattice spacings ranging from \(a\approx0.15\) fm down to 0.06 fm. At most lattice spacings, an ensemble with physical-mass light quarks is included. The use of the highly improved action, combined with the MILC Collaboration's gauge ensembles with lattice spacings down to \(a\approx0.042\) fm, allows heavy valence quarks to be treated with the same discretization as the light and strange quarks. This unified treatment of the valence quarks allows (in some cases) for absolutely normalized currents, bypassing the need for perturbative matching, which has been a leading source of uncertainty in previous calculations of \(B\)-meson decay form factors by our collaboration. All preliminary form-factor results are blinded.

Readout errors are among the most dominant errors on current noisy intermediate-scale quantum devices. Recently, an efficient and scaleable method for mitigating such errors has been developed, based on classical bit-flip correction. In this talk, we compare the performance of this method for IBM's and Rigetti's quantum devices, demonstrating how the method improves the noisy measurements of observables obtained on the quantum hardware. Moreover, we examine the variance amplification to the data after applying of our mitigation procedure, which is common to all mitigation strategies. We derive a new expression for the variance of the mitigated Pauli operators in terms of the corrected expectation values and the noisy variances.Our hardware results show good agreement with the theoretical prediction, and we demonstrate that the increase of the variance due to the mitigation procedure is only moderate.

We consider a nonlocal scalar field theory inspired by the tachyon action in open string field theory. The Lorentz-covariant action is characterized by a parameter \(\xi^2\) that quantifies the amount of nonlocality. Restricting to purely time-dependent configurations, we show that a field redefinition perturbative in \(\xi^2\) reduces the action to a local two-derivative theory with a \(\xi^2\)-dependent potential. This picture is supported by evidence that the redefinition maps the wildly oscillating rolling tachyon solutions of the nonlocal theory to conventional rolling in the new scalar potential. For general field configurations we exhibit an obstruction to a local Lorentz-covariant formulation, but we can still achieve a formulation local in time, as well as a light-cone formulation. These constructions provide an initial value formulation and a Hamiltonian. Their causality is consistent with a lack of superluminal behavior in the nonlocal theory.

Quantum and tensor network simulations have emerged as prominent sign-problem free approaches to lattice gauge theories. Unlike conventional Markov chain Monte Carlo methods, they are based on the Hamiltonian formulation. In this talk, we fill a gap in the literature and present the first derivation of the Hamiltonian 3+1D \(\theta\)-term -- which is an important sign-problem afflicted term -- for Abelian and non-Abelian lattice gauge theories. Furthermore, we perform exact diagonalization for a 3+1D U(1) lattice gauge theory including the \(\theta\)-term on a unit periodic cube. Our numerical results reveal a novel phase transition at fixed values of \(\theta\) in the strong-coupling regime. The transition is evidenced by an avoided level crossing in the ground state energy, as well as sudden changes in the plaquette expectation value, the electric energy density, and the topological charge density. Extensions of our work to larger lattices can be readily performed using state-of-the-art tensor network simulations. Moreover, our work provides a concrete starting point for an eventual quantum simulation of the \(\theta\)-dependent phase structure and dynamics of lattice gauge theories in 3+1D. This talk is mainly based on~\cite{kan2021investigating}. We expand beyond~\cite{kan2021investigating} by including a derivation of the (non-)Abelian fixed-length Higgs term in the Hamiltonian formulation for future studies of (non-)Abelian-Higgs models with a \(\theta\)-term.

Particles may be emitted efficiently from the solar interior if they are sufficiently light and weakly coupled to the solar plasma. In a narrow region of phase space, they are emitted with velocities smaller than the escape velocity of the solar system, thereby populating a gravitationally bound density that can accumulate over the solar lifetime, referred to as a "solar basin." Detection strategies that can succeed in spite of (or even be enhanced by) the low particle velocities are therefore poised to explore new regions of parameter space when taking this solar population into account. Here we identify "direct deflection" as a powerful method to detect such a population of millicharged particles. This approach involves distorting the local flow of gravitationally bound millicharges with an oscillating electromagnetic field and measuring these distortions with a resonant LC circuit. Since it is easier to distort the flow of slowly moving particles, the signal is parametrically enhanced by the small solar escape velocity near Earth. The proposed setup can probe couplings an order of magnitude smaller than other methods for millicharge masses ranging from 100 meV to 100 eV and can operate concurrently as a search for sub-GeV millicharged dark matter. The signal power scales as the millicharge coupling to the eighth power, meaning that even with conservative assumptions, direct deflection could begin to explore new regions of parameter space. We also highlight novel features of millicharge solar basins, including those associated with the phase space distribution and the possibility for the occupation number to vastly exceed that of a thermal distribution.

For any quantum system invariant under gauging a higher-form global symmetry, we construct a non-invertible topological defect by gauging in only half of spacetime. This generalizes the Kramers-Wannier duality line in 1+1 dimensions to higher spacetime dimensions. We focus on the case of a one-form symmetry in 3+1 dimensions, and determine the fusion rule. From modular invariance and a direct analysis of one-form symmetry protected topological phases, we show that the existence of certain kinds of duality defects is intrinsically incompatible with a trivially gapped phase. By further assuming time-reversal symmetry, we find that the presence of certain duality defects implies that the low-energy phase has to be gapless unless the one-form symmetry is spontaneously broken. We give an explicit realization of this duality defect in the free Maxwell theory, both in the continuum and in a modified Villain lattice model. The duality defect is realized by a Chern-Simons coupling between the gauge fields from the two sides. We further construct the duality defect in non-abelian gauge theories and the \(\mathbb{Z}_N\) lattice gauge theory.

We study the relation between approximate horizon symmetries of AdS black branes and approximately conserved currents in their dual hydrodynamic description. We argue that the existence of an approximately conserved enstrophy current unique to \(2+1\) dimensional fluid flow implies that AdS\({}_4\) black branes possess a special class of approximate supertranslations (which we identify).

A vanishing Yukawa coupling of the up quark could in principle solve the strong CP problem. To render this solution consistent with current algebra results, the up quark must receive an alternative mass contribution that conserves CP symmetry. Such a contribution could be provided by QCD through non-perturbative topological effects, including instantons. In this talk, we present the first direct lattice computation of this topological mass contribution, using gauge configurations generated by the Extended Twisted Mass collaboration. We use the Iwasaki gauge action, Wilson twisted mass fermions at maximal twist, and dynamical up, down, strange and charm quarks. Our result for the topological mass contribution is an order of magnitude too small to account for the phenomenologically required up-quark mass. This rules out the "massless" up-quark solution to the strong CP problem, in accordance with previous results relying on \(\chi\)PT fits to lattice data. The talk is based on [Alexandrou et al., PRL 125, 232001 (2020)], where more details can be found.

We numerically study the Hamiltonian lattice formulation of the two-flavor Schwinger model using matrix product states. Keeping the mass of the first flavor at a fixed positive value, we tune the mass of the second flavor through a range of negative values, thus exploring a regime where conventional Monte Carlo methods suffer from the sign problem and may run into instabilities due to zero modes. Our results indicate a phase transition at the point where the absolute value of the second flavor mass approaches the first flavor mass. The phase transition is accompanied by the formation of a fermion condensate, a steep drop of the average electric field, and a peak in the bipartite entanglement entropy. Our data hints at a second order transition, which is the 1+1D analog of the CP-violating Dashen phase transition in QCD.

We study the Large \(N_\text{c}\) scaling of pion-pion scattering amplitudes for \(N_\text{f}=4\) degenerate quark flavors. We focus on the standard isospin-2 channel and the adjoint-antisymmetric (AA) representation. The latter only exists for \(N_\text{f}\geq 4\) and a representative state is \(\frac{1}{\sqrt{2}}(|D_s^+\pi^+\rangle-|D^+ K^+\rangle)\). We compare the results obtained for two regularizations (Wilson and twisted-mass fermions) and three values of the lattice spacing, and observe significant discretization effects in the AA channel. Finally, we match our results to NLO SU(4) and NNLO U(4) Chiral Perturbation Theory and constrain the \(N_\text{c}\) scaling of the relevant low-energy couplings.

Finite-density QCD and many other field theories with sign problems have a \(\mathcal{PT}\)-type symmetry. After a brief introduction to \(\mathcal{PT}\)-symmetric field theories, a real dual representation for \(\mathcal{PT}\)-symmetric scalar field theories with complex actions is derived. We show that \(\mathcal{PT}\)-symmetric field theories can exhibit exotic behavior, including sinusoidally modulated propagators, disorder lines, and spatially inhomogeneous pattern-forming phases. We discuss the interplay of duality, \(\mathcal{PT}\)-symmetry and pattern formation using a \(\phi^4\) model and \(Z(N)\) spin model with sign problems as examples. These behaviors may occur in finite-density QCD and related models.

Axions are a theoretically promising dark matter (DM) candidate. In the presence of radiation from bright astrophysical sources at radio frequencies, nonrelativistic DM axions can undergo stimulated decay to two nearly back-to-back photons, meaning that bright sources of radio waves will have a counterimage (''gegenschein'') in nearly the exact opposite spatial direction. The counterimage will be spectrally distinct from backgrounds, taking the form of a narrow radio line centered at \(\nu = m_a/4\pi\) with a width determined by Doppler broadening in the DM halo, \(\Delta \nu/\nu \sim 10^{-3}\). In this work, we show that the axion decay-induced echoes of supernova remnants may be bright enough to be detectable. Their non-detection may be able to set the strongest limits to date on axion DM in the \(\sim 1-10 \, \mu\)eV mass range where there are gaps in coverage from existing experiments.

We introduce a global analysis of relativistic heavy ion collisions using Trajectum of a significantly higher precision and including a new option to vary the normalization of the centrality estimator. We use the posterior distribution of our parameters to generate a set of high statistics samples that allows us to make precise predictions including statistical and systematic uncertainties estimated from the model parameter distribution. The results are systematically compared with experiment whereby we also include many observables not included in the global analysis. This includes in particular (extremely) ultracentral anisotropic flow and mean transverse momentum, whereby we find satisfactory agreement with experiment where data is available. Lastly, we compute spectra and anisotropic flow for oxygen-oxygen collisions performed at RHIC and to be performed at the LHC and comment on how these collisions may inform us on properties of the Quark-Gluon-Plasma.

We study the impact of selection biases on jet structure and substructure observables and separate these effects from effects caused by jet quenching. We use the angular separation \(\Delta R\) of the hardest splitting in a jet as the primary example observable. We first conduct a simplified Monte Carlo study in which it is possible to identify the same jet after quenching in a heavy ion collision and as it would have been if it had formed in vacuum. We select a sample of jets by placing a cut on their quenched \(p_T\) and, as is possible only in a Monte Carlo study, compare to the same jets unquenched, and see that the \(\Delta R\) distribution seems to be unmodified. However, if we select a sample of jets formed in vacuum by placing a cut on their unquenched \(p_T\) and compare to those same jets after quenching, we see a significant enhancement in the number of jets with large \(\Delta R\), primarily due to the soft particles in the jet that originate from the wake in the droplet of quark-gluon plasma excited by the parton shower. We confirm that the jets contributing to this enhancement are those jets which lost the most energy, which were not included in the sample selected after quenching; jets selected after quenching are those which lose a small fraction of their energy. Next, we employ a method that is available to experimentalists: in a sample of jets with a recoiling \(Z\) boson, we show that selecting jets based on the jet \(p_T\) after quenching yields a \(\Delta R\) distribution that appears unmodified while selecting a sample of jets produced in association with a \(Z\) boson whose (unmodified) \(p_T\) is above some cut yields a significant enhancement in the number of jets with large \(\Delta R\). We again confirm that this is due to particles from the wake, and that the jets contributing to this enhancement are those which have lost a significant fraction of their energy.