High Energy Physics - Lattice (hep-lat)

We discuss matters related to the point that topological quantization in the strong interaction is a consequence of an infinite spacetime volume. Because of the ensuing order of limits, i.e. infinite volume prior to summing over topological sectors, CP is conserved. Here, we show that this reasoning is consistent with the construction of the path integral from steepest-descent contours. We reply to some objections that aim to support the case for CP violation in the strong interactions that are based on the role of the CP-odd theta-parameter in three-form effective theories, the correct sampling of all configurations in the dilute instanton gas approximation and the volume dependence of the partition function. We also show that the chiral effective field theory derived from taking the volume to infinity first is in no contradiction with analyses based on partially conserved axial currents.

The spatial $z$-correlators of meson operators in $N_f=2+1+1$ lattice QCD with optimal domain-wall quarks at the physical point are studied for seven temperatures in the range of 190-1540 MeV. The meson operators include a complete set of Dirac bilinears (scalar, pseudoscalar, vector, axial vector, tensor vector, and axial-tensor vector), and each for six flavor combinations ($\bar u d$, $\bar u s$, $\bar s s$, $\bar u c$, $\bar s c$, and $\bar c c$). In Ref. \citeChiu:2023hnm, we focused on the meson correlators of $u$ and $d$ quarks, and discussed their implications for the effective restoration of $SU(2)_L \times SU(2)_R$ and $U(1)_A$ chiral symmetries, as well as the emergence of approximate $SU(2)_{CS}$ chiral spin symmetry. In this work, we extend our study to meson correlators of six flavor contents, and first observe the hierarchical restoration of chiral symmetries in QCD, from $SU(2)_L \times SU(2)_R \times U(1)_A $ to $SU(3)_L \times SU(3)_R \times U(1)_A $, and to $SU(4)_L \times SU(4)_R \times U(1)_A $, as the temperature is increased from 190 MeV to 1540 MeV. Moreover, we compare the temperature windows for the emergence of the approximate $SU(2)_{CS}$ symmetry in light and heavy vector mesons, and find that the temperature windows are dominated by the $(\bar u c, \bar s c, \bar c c)$ sectors.

The 4D compact U(1) gauge theory has a well-established phase transition between a confining and a Coulomb phase. In this paper, we revisit this model using state-of-the-art Monte Carlo simulations on anisotropic lattices. We map out the coupling-temperature phase diagram, and determine the location of the tricritical point, $T/K_0 \simeq 0.19$, below which the first-order transition is observed. We find the critical exponents of the high-temperature second-order transition to be compatible with those of the 3-dimensional $O(2)$ model. Our results at higher temperatures can be compared with literature results and are consistent with them. Surprisingly, below $T/K_0 \simeq 0.05$ we find strong indications of a second tricritical point where the first-order transition becomes continuous. These results suggest an unexpected second-order phase transition extending down to zero temperature, contrary to the prevailing consensus. If confirmed, these findings reopen the question of the detailed characterization of the transition including a suitable field theory description.

We present a spectroscopy scheme for the lattice field theory by using the tensor renormalization group method combining with the transfer matrix formalism. By using the scheme, we cannot only compute the energy spectrum for the lattice theory but also determine quantum numbers of the energy eigenstates. Furthermore, the wave function of the corresponding eigenstate can also be computed. The first step of the scheme is to coarse grain the tensor network of a given lattice model by using the higher order tensor renormalization group, and then after making a matrix corresponding to a transfer matrix from the coarse-grained tensors, its eigenvalues are evaluated to extract the energy spectrum. Second, the quantum number of the eigenstates can be identified by a selection rule that requires to compute matrix elements of an associated insertion operator. The matrix elements can be represented by an impurity tensor network and computed by the coarse-graining scheme. Moreover, we can compute the wave function of the energy eigenstate by putting the impurity tensor at each point in space direction of the network. Additionally, the momentum of the eigenstate can also be identified by computing appropriate matrix elements represented by the tensor network. As a demonstration of the new scheme, we show the spectroscopy of the $(1+1)$d Ising model and compare it with exact results. We also present a scattering phase shift obtained from two-particle state energy using Lüscher's formula.

The use of machine learning algorithms to investigate phase transitions in physical systems is a valuable way to better understand the characteristics of these systems. Neural networks have been used to extract information of phases and phase transitions directly from many-body configurations. However, one limitation of neural networks is that they require the definition of the model architecture and parameters previous to their application, and such determination is itself a difficult problem. In this paper, we investigate for the first time the relationship between the accuracy of neural networks for information of phases and the network configuration (that comprises the architecture and hyperparameters). We formulate the phase analysis as a regression task, address the question of generating data that reflects the different states of the physical system, and evaluate the performance of neural architecture search for this task. After obtaining the optimized architectures, we further implement smart data processing and analytics by means of neuron coverage metrics, assessing the capability of these metrics to estimate phase transitions. Our results identify the neuron coverage metric as promising for detecting phase transitions in physical systems.

We present a detailed study of the Bardeen-Cooper-Schrieffer (BCS) gap equation ``replicae'' or excited vacuum states, orthogonal to the ground-state one, in the chiral-quark sector of the Hamiltonian Coulomb-gauge model of chromodynamics. Analyzing the number of negative eigenmodes of the energy density's Hessian we believe that we have identified all of the (negative energy-density) vacua of this nonlinear system, namely the ground BCS state and two (or one) replicae for slightly massive (or massless) quarks, given the interaction strength typical of the strong interactions. The meson spectrum over each of the replicae looks similar, so the differences are not significant enough given model uncertainties, but matrix elements are more sensitive and allow to distinguish them. We propose to look for such excited vacua in lattice gauge theory by trying to identify excitations with scalar quantum numbers which have energies proportional to the lattice volume (unlike conventional mesons for which the mass stabilizes to a constant upon taking the infinite volume limit).

This note discusses a method for computing the energy spectra of quantum field theory utilizing digital quantum simulation. A quantum algorithm, called coherent imaging spectroscopy, quenches the vacuum with a time-oscillating perturbation and then reads off the excited energy levels from the loss in the vacuum-to-vacuum probability following the quench. As a practical demonstration, we apply this algorithm to the (1+1)-dimensional quantum electrodynamics with a topological term known as the Schwinger model, where the conventional Monte Carlo approach is practically inaccessible. In particular, on a classical simulator, we prepare the vacuum of the Schwinger model on a lattice by adiabatic state preparation and then apply various types of quenches to the approximate vacuum through Suzuki-Trotter time evolution. We discuss the dependence of the simulation results on the specific types of quenches and introduce various consistency checks, including the exact diagonalization and the continuum limit extrapolation. The estimation of the computational complexity required to obtain physically reasonable results implies that the method is likely efficient in the coming era of early fault-tolerant quantum computers.

We present a state-of-the-art calculation of the unpolarized pion valence-quark distribution in the framework of large-momentum effective theory (LaMET) with improved handling of systematic errors as well as two-loop perturbative matching. We use lattice ensembles generated by the MILC collaboration at lattice spacing $a\approx 0.09$~fm, lattice volume $64^3\times 96$, $N_f=2+1+1$ flavors of highly-improved staggered quarks and a physical pion mass. The LaMET matrix elements are calculated with pions boosted to momentum $P_z\approx 1.72$~GeV with high-statistics of $O(10^6)$ measurements. We study the pion PDF in both hybrid-ratio and hybrid-regularization-independent momentum subtraction (hybrid-RI/MOM) schemes and also compare the systematic errors with and without the addition of leading-renormalon resummation (LRR) and renormalization-group resummation (RGR) in both the renormalization and lightcone matching. The final lightcone PDF results are presented in the modified minimal-subtraction scheme at renormalization scale $\mu=2.0$~GeV. We show that the $x$-dependent PDFs are compatible between the hybrid-ratio and hybrid-RI/MOM renormalization with the same improvements. We also show that systematics are greatly reduced by the simultaneous inclusion of RGR and LRR and that these methods are necessary if improved precision is to be reached with higher-order terms in renormalization and matching.