Search

Abstract Hadronic vacuum polarization (HVP) and light-by-light scattering (HLBL) are crucial for evaluating the Standard Model predictions concerning the muon’s anomalous magnetic moment. However, direct first-principle lattice gauge theory-based calculations of these observables in the timelike region remain challenging. Discrepancies persist between lattice quantum chromodynamics (QCD) calculations in the spacelike region and dispersive approaches relying on experimental data parametrization from the timelike region. Here, we introduce a methodology employing 1+1-dimensional quantum electrodynamics (QED), i.e. the Schwinger Model, to investigate the HVP and HLBL. To that end, we use both tensor network techniques, specifically matrix product states, and classical emulators of digital quantum computers. Demonstrating feasibility in a simplified model, our approach sets the stage for future endeavors leveraging digital quantum computers.

Abstract Anomaly detection is a vital technique for exploring signatures of new physics Beyond the Standard Model (BSM) at the Large Hadron Collider (LHC). The vast number of collisions generated by the LHC demands sophisticated deep learning techniques. Similarity learning, a self-supervised machine learning, detects anomalous signals by estimating their similarity to background events. In this paper, we explore the potential of quantum computers for anomaly detection through similarity learning, leveraging the power of quantum computing to enhance the known similarity learning method. In the realm of noisy intermediate-scale quantum (NISQ) devices, we employ a hybrid classical-quantum network to search for heavy scalar resonances in the di-Higgs production channel. In the absence of quantum noise, the hybrid network demonstrates improvement over the known similarity learning method. Moreover, we employ a clustering algorithm to reduce measurement noise from limited shot counts, resulting in 9% improvement in the hybrid network performance. Our analysis highlights the applicability of quantum algorithms for LHC data analysis, where improvements are anticipated with the advent of fault-tolerant quantum computers.

Abstract Flow-based architectures have recently proved to be an efficient tool for numerical simulations of Effective String Theories regularized on the lattice that otherwise cannot be efficiently sampled by standard Monte Carlo methods. In this work we use Stochastic Normalizing Flows, a state-of-the-art deep learning architecture based on non-equilibrium Monte Carlo simulations, to study different effective string models. After testing the reliability of this approach through a comparison with exact results for the Nambu-Gotō model, we discuss results on observables that are challenging to study analytically, such as the width of the string and the shape of the flux density. Furthermore, we perform a novel numerical study of Effective String Theories with terms beyond the Nambu-Gotō action, including a broader discussion on their significance for lattice gauge theories. The combination of these findings enables a quantitative description of the fine details of the confinement mechanism in different lattice gauge theories. The results presented in this work establish the reliability and feasibility of flow-based samplers for Effective String Theories and pave the way for future applications on more complex models.

Abstract It was recently suggested — based on general self-consistency arguments as well as results from the bootstrap [1–3] — that the CFT describing the Q-state Potts model is logarithmic for generic values of Q, with rank-two Jordan blocks for L 0 and L ¯ 0 $$ {\overline{L}}_0 $$ in many sectors of the theory. This is despite the well-known fact that the lattice transfer matrix (or Hamiltonian) is diagonalizable in (arbitrary) finite size. While the emergence of Jordan blocks only in the limit L → ∞ is perfectly possible conceptually, diagonalizability in finite size makes the measurement of logarithmic couplings (whose values are analytically predicted in [2, 3]) very challenging. This problem is solved in the present paper (which can be considered a companion to [2]), and the conjectured logarithmic structure of the CFT confirmed in detail by the study of the lattice model and associated “emerging Jordan blocks.”

Abstract We present a novel strategy to strongly reduce the severity of the sign problem, using line integrals along paths of changing imaginary action. Highly oscillating regions along these paths cancel out, decreasing their contributions. As a result, sampling with standard Monte-Carlo techniques becomes possible in cases that otherwise require methods taking advantage of complex analysis, such as Lefschetz-thimbles or Complex Langevin. We lay out how to write down an ordinary differential equation for the line integrals. As an example of its usage, we apply the results to a 1d quantum mechanical anharmonic oscillator with a x 4 potential in real time, finite temperature.

Abstract In this work we provide a non-perturbative solution to the theoretical problem of extracting scattering amplitudes from Euclidean correlators in infinite volume. We work within the solid axiomatic framework of the Haag-Ruelle scattering theory and derive formulae which can be used to approximate scattering amplitudes arbitrarily well in terms of linear combinations of Euclidean correlators at discrete time separations. Our result generalizes and extends the range of applicability of a result previously obtained by Barata and Fredenhagen [1]. We provide a concrete procedure to construct such approximations, making our formulae ready to be used in numerical calculations of non-perturbative QCD scattering amplitudes. A detailed numerical investigation is needed to assess whether the proposed strategy can lead to the calculation of scattering amplitudes with phenomenologically satisfactory precision with presently available lattice QCD data. This will be the subject of future work. Nevertheless, the numerical accuracy and precision of lattice simulations is systematically improvable, and we have little doubts that our approach will become useful in the future.

Abstract We present efficient data-driven approaches to predict the value of subdivergence-free Feynman integrals (Feynman periods) in ϕ 4-theory from properties of the underlying Feynman graphs, based on a statistical examination of almost 2 million graphs. We find that the numbers of cuts and cycles determines the period to better than 2% relative accuracy. Hepp bound and Martin invariant allow for even more accurate predictions. In most cases, the period is a multi-linear function of the properties in question. Furthermore, we investigate the usefulness of machine-learning algorithms to predict the period. When sufficiently many properties of the graph are used, the period can be predicted with better than 0.05% relative accuracy. We use one of the constructed prediction models for weighted Monte-Carlo sampling of Feynman graphs, and compute the primitive contribution to the beta function of ϕ 4-theory at L ∈ {13, … , 17} loops. Our results confirm the previously known numerical estimates of the primitive beta function and improve their accuracy. Compared to uniform random sampling of graphs, our new algorithm is 1000-times faster to reach a desired accuracy, or reaches 32-fold higher accuracy in fixed runtime. The dataset of all periods computed for this work, combined with a previous dataset, is made publicly available. Besides the physical application, it could serve as a benchmark for graph-based machine learning algorithms.

Abstract We present a flow-based method for simulating and calculating nucleation rates of first-order phase transitions in scalar field theory on a lattice. Motivated by recent advancements in machine learning tools, particularly normalizing flows for lattice field theory, we propose the “partitioning flow-based Markov chain Monte Carlo (PFMCMC) sampling” method to address two challenges encountered in normalizing flow applications for lattice field theory: the “mode-collapse” and “rare-event sampling” problems. Using a (2+1)-dimensional real scalar model as an example, we demonstrate the effectiveness of our PFMCMC method in modeling highly hierarchical order parameter probability distributions and simulating critical bubble configurations. These simulations are then used to facilitate the calculation of nucleation rates. We anticipate the application of this method to (3+1)-dimensional theories for studying realistic cosmological phase transitions.

Abstract We conjecture a multi-parameter generalization of the toric inequalities of [1]. We then extend their proof methods for the generalized toric inequalities in two ways. The first extension constructs the graph corresponding to the toric inequalities and the generalized toric conjectures by tiling the Euclidean space. An entanglement wedge nesting relation then determines the geometric structure of the tiles. In the second extension, we exploit the cyclic nature of the inequalities and conjectures to construct cycle graphs. Then, the graph can be obtained using graph Cartesian products of cycle graphs. In addition, we define a set of knots on the graph by following [1]. These graphs with knots then imply the validity of their associated inequality. We study the case where the graph can be decomposed into disjoint unions of torii. Under the specific case, we explore and prove the conjectures for some ranges of parameters. We also discuss ways to explore the conjectured inequalities whose corresponding geometries are d-dimensional torii (d > 2).

We present efficient data-driven approaches to predict the value of subdivergence-free Feynman integrals (Feynman periods) in ϕ4-theory from properties of the underlying Feynman graphs, based on a statistical examination of almost 2 million graphs. We find that the numbers of cuts and cycles determines the period to better than 2% relative accuracy. Hepp bound and Martin invariant allow for even more accurate predictions. In most cases, the period is a multi-linear function of the properties in question. Furthermore, we investigate the usefulness of machine-learning algorithms to predict the period. When sufficiently many properties of the graph are used, the period can be predicted with better than 0.05% relative accuracy. We use one of the constructed prediction models for weighted Monte-Carlo sampling of Feynman graphs, and compute the primitive contribution to the beta function of ϕ4-theory at L ∈ {13, … , 17} loops. Our results confirm the previously known numerical estimates of the primitive beta function and improve their accuracy. Compared to uniform random sampling of graphs, our new algorithm is 1000-times faster to reach a desired accuracy, or reaches 32-fold higher accuracy in fixed runtime. The dataset of all periods computed for this work, combined with a previous dataset, is made publicly available. Besides the physical application, it could serve as a benchmark for graph-based machine learning algorithms. [ABSTRACT FROM AUTHOR]

We conjecture a multi-parameter generalization of the toric inequalities of [1]. We then extend their proof methods for the generalized toric inequalities in two ways. The first extension constructs the graph corresponding to the toric inequalities and the generalized toric conjectures by tiling the Euclidean space. An entanglement wedge nesting relation then determines the geometric structure of the tiles. In the second extension, we exploit the cyclic nature of the inequalities and conjectures to construct cycle graphs. Then, the graph can be obtained using graph Cartesian products of cycle graphs. In addition, we define a set of knots on the graph by following [1]. These graphs with knots then imply the validity of their associated inequality. We study the case where the graph can be decomposed into disjoint unions of torii. Under the specific case, we explore and prove the conjectures for some ranges of parameters. We also discuss ways to explore the conjectured inequalities whose corresponding geometries are d-dimensional torii (d > 2). [ABSTRACT FROM AUTHOR]

We present a flow-based method for simulating and calculating nucleation rates of first-order phase transitions in scalar field theory on a lattice. Motivated by recent advancements in machine learning tools, particularly normalizing flows for lattice field theory, we propose the “partitioning flow-based Markov chain Monte Carlo (PFMCMC) sampling” method to address two challenges encountered in normalizing flow applications for lattice field theory: the “mode-collapse” and “rare-event sampling” problems. Using a (2+1)-dimensional real scalar model as an example, we demonstrate the effectiveness of our PFMCMC method in modeling highly hierarchical order parameter probability distributions and simulating critical bubble configurations. These simulations are then used to facilitate the calculation of nucleation rates. We anticipate the application of this method to (3+1)-dimensional theories for studying realistic cosmological phase transitions. [ABSTRACT FROM AUTHOR]

Abstract We investigate the phase structure of the (1+1)-dimensional U(1) gauge-Higgs model with a θ term, where the U(1) gauge action is constructed with Lüscher’s admissibility condition. Using the tensor renormalization group, both the complex action problem and topological freezing problem in the standard Monte Carlo simulation are avoided. We find the first-order phase transition with sufficiently large Higgs mass at θ = π, where the ℤ2 charge conjugation symmetry is spontaneously broken. On the other hand, the symmetry is restored with a sufficiently small mass. We determine the critical endpoint as a function of the Higgs mass parameter and show the critical behavior is in the two-dimensional Ising universality class.

Abstract We study the θ-dependent mass spectrum of the massive 2-flavor Schwinger model in the Hamiltonian formalism using the density-matrix renormalization group (DMRG). The masses of the composite particles, the pion and sigma meson, are computed by two independent methods. One is the improved one-point-function scheme, where we measure the local meson operator coupled to the boundary state and extract the mass from its exponential decay. Since the θ term causes a nontrivial operator mixing, we unravel it by diagonalizing the correlation matrix to define the meson operator. The other is the dispersion-relation scheme, a heuristic approach specific to Hamiltonian formalism. We obtain the dispersion relation directly by measuring the energy and momentum of the excited states. The sign problem is circumvented in these methods, and their results agree with each other even for large θ. We reveal that the θ-dependence of the pion mass at m/g = 0.1 is consistent with the prediction by the bosonized model. We also find that the mass of the sigma meson satisfies the semi-classical formula, M σ /M π = 3 $$ \sqrt{3} $$ , for almost all region of θ. While the sigma meson is a stable particle thanks to this relation, the eta meson is no longer protected by the G-parity and becomes unstable for θ ≠ 0.

Abstract The two-particle finite-volume scattering formalism derived by Lüscher and generalized in many subsequent works does not hold for energies far enough below the two-particle threshold to reach the nearest left-hand cut. The breakdown of the formalism is signaled by the fact that a real scattering amplitude is predicted in a regime where it should be complex. In this work, we address this limitation by deriving an extended formalism that includes the nearest branch cut, arising from single particle exchange. We focus on two-nucleon (NN → NN) scattering, for which the cut arises from pion exchange, but give expressions for any system with a single channel of identical particles. The new result takes the form of a modified quantization condition that can be used to constrain an intermediate K-matrix in which the cut is removed. In a second step, integral equations, also derived in this work, must be used to convert the K-matrix to the physical scattering amplitude. We also show how the new formalism reduces to the standard approach when the N → Nπ coupling is set to zero.

Abstract Krylov complexity is a measure of operator growth in quantum systems, based on the number of orthogonal basis vectors needed to approximate the time evolution of an operator. In this paper, we study the Krylov complexity of a PT-symmetric system of oscillators, which exhibits two phase transitions that separate a dissipative state, a Rabi-oscillation state, and an ultra-strongly coupled regime. We use a generalization of the su(1) algebra associated to the Bateman oscillator to describe the Hamiltonian of the coupled system, and construct a set of coherent states associated with this algebra. We compute the Krylov (spread) complexity using these coherent states, and find that it can distinguish between the PT-symmetric and PT symmetry-broken phases. We also show that the Krylov complexity reveals the ill-defined nature of the vacuum of the Bateman oscillator, which is a special case of our system. Our results demonstrate the utility of Krylov complexity as a tool to probe the properties and transitions of PT-symmetric systems.

Abstract The investigation of hadron interactions within lattice QCD has been facilitated by the well-known quantisation condition, linking scattering phase shifts to finite-volume energies. Additionally, the ability to utilise systems at finite total boosts has been pivotal in smoothly charting the energy-dependent behaviour of these phase shifts. The existing implementations of the quantization condition at finite boosts rely on momentum transformations between rest and moving frames, defined directly in terms of the energy eigenvalues. This energy dependence is unsuitable in the formulation of a Hamiltonian. In this work, we introduce a novel approach to generalise the three-momentum boost prescription, enabling the incorporation of energy-independent finite-volume Hamiltonians within moving frames. We demonstrate the application of our method through numerical comparisons, employing a phenomenological ππ scattering example.

Abstract We simulate N f = 2 + 1 QCD at the physical point combining open and periodic boundary conditions in a parallel tempering framework, following the original proposal by M. Hasenbusch for 2d CP N−1 models, which has been recently implemented and widely employed in 4d SU(N) pure Yang-Mills theories too. We show that using this algorithm it is possible to achieve a sizable reduction of the auto-correlation time of the topological charge in dynamical fermions simulations both at zero and finite temperature, allowing to avoid topology freezing down to lattice spacings as fine as a ∼ 0.02 fm. Therefore, this implementation of the Parallel Tempering on Boundary Conditions algorithm has the potential to substantially push forward the investigation of the QCD vacuum properties by means of lattice simulations.

Abstract We apply the previously-developed sub-volume method to study the θ-dependence of the four-dimensional SU(2) Yang-Mills theory at finite temperature. We calculate the first two coefficients, the topological susceptibility χ and the fourth cumulant b 2, in the θ-expansion of the free energy density around the critical temperature (T c ) for the confinement-deconfinement transition. Lattice calculations are performed with three different spatial sizes 243 , 323 , 483 to monitor finite size effects, while the temporal size is fixed to be 8. The systematic uncertainty associated with the sub-volume extrapolation is studied with special care. The sub-volume method allows us to determine the values of b 2 much more accurately than the standard full-volume method, and we successfully identify the temperature dependence of b 2 around T c . Our numerical results suggest that the θ-dependence of the free energy density near θ = 0 changes from 4χ(1 − cos(θ/2)) to χ(1 − cos θ) as the temperature crosses T c .

Abstract We reanalyze the origin of the large transverse logarithms associated with the QCD one loop β function coefficient in the NLO JIMWLK Hamiltonian. We show that some of these terms are not associated with the running of the QCD coupling constant but rather with the DGLAP evolution. The DGLAP-like resummation of these logarithms is mandatory within the JIMWLK Hamiltonian, as long as the color correlation length in the projectile is larger than that in the target. This regime in fact covers the whole range of rapidities at which JIMWLK evolution is supposed to be applicable. We derive the RG equation that resums these logarithms to all orders in α s in the JIMWLK Hamiltonian. This is a nonlinear equation for the eikonal scattering matrix S(x). We solve this equation, and perform the DGLAP resummation in two simple cases: the dilute limit, where both the projectile and the target are far from saturation, and the saturated regime, where the target correlation length also determines its saturation momentum.

Abstract The generalized Lefschetz thimble method is a promising approach that attempts to solve the sign problem in Monte Carlo methods by deforming the integration contour using the flow equation. Here we point out a general problem that occurs due to the property of the flow equation, which extends a region on the original contour exponentially to a region on the deformed contour. Since the growth rate for each eigenmode is governed by the singular values of the Hessian of the action, a huge hierarchy in the singular value spectrum, which typically appears for large systems, leads to various technical problems in numerical simulations. We solve this hierarchical growth problem by preconditioning the flow so that the growth rate becomes identical for every eigenmode. As an example, we show that the preconditioned flow enables us to investigate the real-time quantum evolution of an anharmonic oscillator with the system size that can hardly be achieved by using the original flow.

Abstract We propose an adaptation of Entanglement Renormalization for quantum field theories that, through the use of discrete wavelet transforms, strongly parallels the tensor network architecture of the Multiscale Entanglement Renormalization Ansatz (a.k.a. MERA). Our approach, called wMERA, has several advantages of over previous attempts to adapt MERA to continuum systems. In particular, (i) wMERA is formulated directly in position space, hence preserving the quasi-locality and sparsity of entanglers; and (ii) it enables a built-in RG flow in the implementation of real-time evolution and in computations of correlation functions, which is key for efficient numerical implementations. As examples, we describe in detail two concrete implementations of our wMERA algorithm for free scalar and fermionic theories in (1+1) spacetime dimensions. Possible avenues for constructing wMERAs for interacting field theories are also discussed.

We study the θ-dependent mass spectrum of the massive 2-flavor Schwinger model in the Hamiltonian formalism using the density-matrix renormalization group (DMRG). The masses of the composite particles, the pion and sigma meson, are computed by two independent methods. One is the improved one-point-function scheme, where we measure the local meson operator coupled to the boundary state and extract the mass from its exponential decay. Since the θ term causes a nontrivial operator mixing, we unravel it by diagonalizing the correlation matrix to define the meson operator. The other is the dispersion-relation scheme, a heuristic approach specific to Hamiltonian formalism. We obtain the dispersion relation directly by measuring the energy and momentum of the excited states. The sign problem is circumvented in these methods, and their results agree with each other even for large θ. We reveal that the θ-dependence of the pion mass at m/g = 0.1 is consistent with the prediction by the bosonized model. We also find that the mass of the sigma meson satisfies the semi-classical formula, Mσ/Mπ = 3 , for almost all region of θ. While the sigma meson is a stable particle thanks to this relation, the eta meson is no longer protected by the G-parity and becomes unstable for θ ≠ 0. [ABSTRACT FROM AUTHOR]

We investigate the phase structure of the (1+1)-dimensional U(1) gauge-Higgs model with a θ term, where the U(1) gauge action is constructed with Lüscher’s admissibility condition. Using the tensor renormalization group, both the complex action problem and topological freezing problem in the standard Monte Carlo simulation are avoided. We find the first-order phase transition with sufficiently large Higgs mass at θ = π, where the ℤ2 charge conjugation symmetry is spontaneously broken. On the other hand, the symmetry is restored with a sufficiently small mass. We determine the critical endpoint as a function of the Higgs mass parameter and show the critical behavior is in the two-dimensional Ising universality class. [ABSTRACT FROM AUTHOR]

We simulate Nf = 2 + 1 QCD at the physical point combining open and periodic boundary conditions in a parallel tempering framework, following the original proposal by M. Hasenbusch for 2d CPN−1 models, which has been recently implemented and widely employed in 4d SU(N) pure Yang-Mills theories too. We show that using this algorithm it is possible to achieve a sizable reduction of the auto-correlation time of the topological charge in dynamical fermions simulations both at zero and finite temperature, allowing to avoid topology freezing down to lattice spacings as fine as a ∼ 0.02 fm. Therefore, this implementation of the Parallel Tempering on Boundary Conditions algorithm has the potential to substantially push forward the investigation of the QCD vacuum properties by means of lattice simulations. [ABSTRACT FROM AUTHOR]

The investigation of hadron interactions within lattice QCD has been facilitated by the well-known quantisation condition, linking scattering phase shifts to finite-volume energies. Additionally, the ability to utilise systems at finite total boosts has been pivotal in smoothly charting the energy-dependent behaviour of these phase shifts. The existing implementations of the quantization condition at finite boosts rely on momentum transformations between rest and moving frames, defined directly in terms of the energy eigenvalues. This energy dependence is unsuitable in the formulation of a Hamiltonian. In this work, we introduce a novel approach to generalise the three-momentum boost prescription, enabling the incorporation of energy-independent finite-volume Hamiltonians within moving frames. We demonstrate the application of our method through numerical comparisons, employing a phenomenological ππ scattering example. [ABSTRACT FROM AUTHOR]

Krylov complexity is a measure of operator growth in quantum systems, based on the number of orthogonal basis vectors needed to approximate the time evolution of an operator. In this paper, we study the Krylov complexity of a PT-symmetric system of oscillators, which exhibits two phase transitions that separate a dissipative state, a Rabi-oscillation state, and an ultra-strongly coupled regime. We use a generalization of the su(1) algebra associated to the Bateman oscillator to describe the Hamiltonian of the coupled system, and construct a set of coherent states associated with this algebra. We compute the Krylov (spread) complexity using these coherent states, and find that it can distinguish between the PT-symmetric and PT symmetry-broken phases. We also show that the Krylov complexity reveals the ill-defined nature of the vacuum of the Bateman oscillator, which is a special case of our system. Our results demonstrate the utility of Krylov complexity as a tool to probe the properties and transitions of PT-symmetric systems. [ABSTRACT FROM AUTHOR]

The two-particle finite-volume scattering formalism derived by Lüscher and generalized in many subsequent works does not hold for energies far enough below the two-particle threshold to reach the nearest left-hand cut. The breakdown of the formalism is signaled by the fact that a real scattering amplitude is predicted in a regime where it should be complex. In this work, we address this limitation by deriving an extended formalism that includes the nearest branch cut, arising from single particle exchange. We focus on two-nucleon (NN → NN) scattering, for which the cut arises from pion exchange, but give expressions for any system with a single channel of identical particles. The new result takes the form of a modified quantization condition that can be used to constrain an intermediate K-matrix in which the cut is removed. In a second step, integral equations, also derived in this work, must be used to convert the K-matrix to the physical scattering amplitude. We also show how the new formalism reduces to the standard approach when the N → Nπ coupling is set to zero. [ABSTRACT FROM AUTHOR]

We apply the previously-developed sub-volume method to study the θ-dependence of the four-dimensional SU(2) Yang-Mills theory at finite temperature. We calculate the first two coefficients, the topological susceptibility χ and the fourth cumulant b2, in the θ-expansion of the free energy density around the critical temperature (Tc) for the confinement-deconfinement transition. Lattice calculations are performed with three different spatial sizes 243, 323, 483 to monitor finite size effects, while the temporal size is fixed to be 8. The systematic uncertainty associated with the sub-volume extrapolation is studied with special care. The sub-volume method allows us to determine the values of b2 much more accurately than the standard full-volume method, and we successfully identify the temperature dependence of b2 around Tc. Our numerical results suggest that the θ-dependence of the free energy density near θ = 0 changes from 4χ(1 − cos(θ/2)) to χ(1 − cos θ) as the temperature crosses Tc. [ABSTRACT FROM AUTHOR]

The generalized Lefschetz thimble method is a promising approach that attempts to solve the sign problem in Monte Carlo methods by deforming the integration contour using the flow equation. Here we point out a general problem that occurs due to the property of the flow equation, which extends a region on the original contour exponentially to a region on the deformed contour. Since the growth rate for each eigenmode is governed by the singular values of the Hessian of the action, a huge hierarchy in the singular value spectrum, which typically appears for large systems, leads to various technical problems in numerical simulations. We solve this hierarchical growth problem by preconditioning the flow so that the growth rate becomes identical for every eigenmode. As an example, we show that the preconditioned flow enables us to investigate the real-time quantum evolution of an anharmonic oscillator with the system size that can hardly be achieved by using the original flow. [ABSTRACT FROM AUTHOR]

We reanalyze the origin of the large transverse logarithms associated with the QCD one loop β function coefficient in the NLO JIMWLK Hamiltonian. We show that some of these terms are not associated with the running of the QCD coupling constant but rather with the DGLAP evolution. The DGLAP-like resummation of these logarithms is mandatory within the JIMWLK Hamiltonian, as long as the color correlation length in the projectile is larger than that in the target. This regime in fact covers the whole range of rapidities at which JIMWLK evolution is supposed to be applicable. We derive the RG equation that resums these logarithms to all orders in αs in the JIMWLK Hamiltonian. This is a nonlinear equation for the eikonal scattering matrix S(x). We solve this equation, and perform the DGLAP resummation in two simple cases: the dilute limit, where both the projectile and the target are far from saturation, and the saturated regime, where the target correlation length also determines its saturation momentum. [ABSTRACT FROM AUTHOR]

We propose an adaptation of Entanglement Renormalization for quantum field theories that, through the use of discrete wavelet transforms, strongly parallels the tensor network architecture of the Multiscale Entanglement Renormalization Ansatz (a.k.a. MERA). Our approach, called wMERA, has several advantages of over previous attempts to adapt MERA to continuum systems. In particular, (i) wMERA is formulated directly in position space, hence preserving the quasi-locality and sparsity of entanglers; and (ii) it enables a built-in RG flow in the implementation of real-time evolution and in computations of correlation functions, which is key for efficient numerical implementations. As examples, we describe in detail two concrete implementations of our wMERA algorithm for free scalar and fermionic theories in (1+1) spacetime dimensions. Possible avenues for constructing wMERAs for interacting field theories are also discussed. [ABSTRACT FROM AUTHOR]

Abstract Ab-initio simulations of multiple heavy quarks propagating in a Quark-Gluon Plasma are computationally difficult to perform due to the large dimension of the space of density matrices. This work develops machine learning algorithms to overcome this difficulty by approximating exact quantum states with neural network parametrisations, specifically Neural Density Operators. As a proof of principle demonstration in a QCD-like theory, the approach is applied to solve the Lindblad master equation in the 1 + 1d lattice Schwinger Model as an open quantum system. Neural Density Operators enable the study of in-medium dynamics on large lattice volumes, where multiple-string interactions and their effects on string-breaking and recombination phenomena can be studied. Thermal properties of the system at equilibrium can also be probed with these methods by variationally constructing the steady state of the Lindblad master equation. Scaling of this approach with system size is studied, and numerical demonstrations on up to 32 spatial lattice sites and with up to 3 interacting strings are performed.

Abstract We present a deterministic way of finding contraction maps for candidate holographic entanglement entropy inequalities modulo choices due to actual degeneracy. We characterize its complexity and give an argument for the completeness of the contraction map proof method as a necessary and sufficient condition for the validity of an entropy inequality for holographic entanglement.

Abstract The study of entanglement in gauge theories is expected to provide insights into many fundamental phenomena, including confinement. However, calculations of quantities related to entanglement in gauge theories are limited by ambiguities that stem from the non-factorizability of the Hilbert space. In this work we study lattice gauge theories that admit a dual description in terms of spin models, for which the replica trick and Rényi entropies are well defined. In the first part of this work, we explicitly perform the duality transformation in a replica geometry, deriving the structure of a replica space for a gauge theory. Then, in the second part, we calculate, by means of Monte Carlo simulations, the entropic c-function of the ℤ 2 gauge theory in three spacetime dimensions, exploiting its dual description in terms of the three-dimensional Ising model.

Abstract We propose an orbifold lattice formulation of QCD suitable for quantum simulations. We show explicitly how to encode gauge degrees of freedom into qubits using noncompact variables, and how to write down a simple truncated Hamiltonian in the coordinate basis. We show that SU(3) gauge group variables and quarks in the fundamental representation can be implemented straightforwardly on qubits, for arbitrary truncation of the gauge manifold.

Abstract We show how to construct a tensor network representation of the path integral for reduced staggered fermions coupled to a non-abelian gauge field in two dimensions. The resulting formulation is both memory and computation efficient because reduced staggered fermions can be represented in terms of a minimal number of tensor indices while the gauge sector can be approximated using Gaussian quadrature with a truncation. Numerical results obtained using the Grassmann TRG algorithm are shown for the case of SU(2) lattice gauge theory and compared to Monte Carlo results.

Abstract In this work, we establish a direct connection between generative diffusion models (DMs) and stochastic quantization (SQ). The DM is realized by approximating the reversal of a stochastic process dictated by the Langevin equation, generating samples from a prior distribution to effectively mimic the target distribution. Using numerical simulations, we demonstrate that the DM can serve as a global sampler for generating quantum lattice field configurations in two-dimensional ϕ 4 theory. We demonstrate that DMs can notably reduce autocorrelation times in the Markov chain, especially in the critical region where standard Markov Chain Monte-Carlo (MCMC) algorithms experience critical slowing down. The findings can potentially inspire further advancements in lattice field theory simulations, in particular in cases where it is expensive to generate large ensembles.

Abstract Motivated by the recently-established connection between Jarzynski’s equality and the theoretical framework of Stochastic Normalizing Flows, we investigate a protocol relying on out-of-equilibrium lattice Monte Carlo simulations to mitigate the infamous computational problem of topological freezing. We test our proposal on 2d CP N−1 models and compare our results with those obtained adopting the Parallel Tempering on Boundary Conditions proposed by M. Hasenbusch, obtaining comparable performances. Our work thus sets the stage for future applications combining our Monte Carlo setup with machine learning techniques.

Ab-initio simulations of multiple heavy quarks propagating in a Quark-Gluon Plasma are computationally difficult to perform due to the large dimension of the space of density matrices. This work develops machine learning algorithms to overcome this difficulty by approximating exact quantum states with neural network parametrisations, specifically Neural Density Operators. As a proof of principle demonstration in a QCD-like theory, the approach is applied to solve the Lindblad master equation in the 1 + 1d lattice Schwinger Model as an open quantum system. Neural Density Operators enable the study of in-medium dynamics on large lattice volumes, where multiple-string interactions and their effects on string-breaking and recombination phenomena can be studied. Thermal properties of the system at equilibrium can also be probed with these methods by variationally constructing the steady state of the Lindblad master equation. Scaling of this approach with system size is studied, and numerical demonstrations on up to 32 spatial lattice sites and with up to 3 interacting strings are performed. [ABSTRACT FROM AUTHOR]

The study of entanglement in gauge theories is expected to provide insights into many fundamental phenomena, including confinement. However, calculations of quantities related to entanglement in gauge theories are limited by ambiguities that stem from the non-factorizability of the Hilbert space. In this work we study lattice gauge theories that admit a dual description in terms of spin models, for which the replica trick and Rényi entropies are well defined. In the first part of this work, we explicitly perform the duality transformation in a replica geometry, deriving the structure of a replica space for a gauge theory. Then, in the second part, we calculate, by means of Monte Carlo simulations, the entropic c-function of the ℤ2 gauge theory in three spacetime dimensions, exploiting its dual description in terms of the three-dimensional Ising model. [ABSTRACT FROM AUTHOR]

We present a deterministic way of finding contraction maps for candidate holographic entanglement entropy inequalities modulo choices due to actual degeneracy. We characterize its complexity and give an argument for the completeness of the contraction map proof method as a necessary and sufficient condition for the validity of an entropy inequality for holographic entanglement. [ABSTRACT FROM AUTHOR]

We propose an orbifold lattice formulation of QCD suitable for quantum simulations. We show explicitly how to encode gauge degrees of freedom into qubits using noncompact variables, and how to write down a simple truncated Hamiltonian in the coordinate basis. We show that SU(3) gauge group variables and quarks in the fundamental representation can be implemented straightforwardly on qubits, for arbitrary truncation of the gauge manifold. [ABSTRACT FROM AUTHOR]

We show how to construct a tensor network representation of the path integral for reduced staggered fermions coupled to a non-abelian gauge field in two dimensions. The resulting formulation is both memory and computation efficient because reduced staggered fermions can be represented in terms of a minimal number of tensor indices while the gauge sector can be approximated using Gaussian quadrature with a truncation. Numerical results obtained using the Grassmann TRG algorithm are shown for the case of SU(2) lattice gauge theory and compared to Monte Carlo results. [ABSTRACT FROM AUTHOR]

In this work, we establish a direct connection between generative diffusion models (DMs) and stochastic quantization (SQ). The DM is realized by approximating the reversal of a stochastic process dictated by the Langevin equation, generating samples from a prior distribution to effectively mimic the target distribution. Using numerical simulations, we demonstrate that the DM can serve as a global sampler for generating quantum lattice field configurations in two-dimensional ϕ4 theory. We demonstrate that DMs can notably reduce autocorrelation times in the Markov chain, especially in the critical region where standard Markov Chain Monte-Carlo (MCMC) algorithms experience critical slowing down. The findings can potentially inspire further advancements in lattice field theory simulations, in particular in cases where it is expensive to generate large ensembles. [ABSTRACT FROM AUTHOR]

Motivated by the recently-established connection between Jarzynski's equality and the theoretical framework of Stochastic Normalizing Flows, we investigate a protocol relying on out-of-equilibrium lattice Monte Carlo simulations to mitigate the infamous computational problem of topological freezing. We test our proposal on 2d CPN−1 models and compare our results with those obtained adopting the Parallel Tempering on Boundary Conditions proposed by M. Hasenbusch, obtaining comparable performances. Our work thus sets the stage for future applications combining our Monte Carlo setup with machine learning techniques. [ABSTRACT FROM AUTHOR]

Abstract We investigate the entanglement and Rényi entropies for the (1+1)-dimensional O(3) nonlinear sigma model using the tensor renormalization group method. The central charge is determined from the asymptotic scaling properties of both entropies. We also examine the consistency between the entanglement entropy and the nth-order Rényi entropy with n → 1.

Abstract Effective String Theory (EST) represents a powerful non-perturbative approach to describe confinement in Yang-Mills theory that models the confining flux tube as a thin vibrating string. EST calculations are usually performed using the zeta-function regularization: however there are situations (for instance the study of the shape of the flux tube or of the higher order corrections beyond the Nambu-Goto EST) which involve observables that are too complex to be addressed in this way. In this paper we propose a numerical approach based on recent advances in machine learning methods to circumvent this problem. Using as a laboratory the Nambu-Goto string, we show that by using a new class of deep generative models called Continuous Normalizing Flows it is possible to obtain reliable numerical estimates of EST predictions.

Abstract We develop a transfer operator approach for the calculation of Rényi entanglement entropies in arbitrary (i.e. Abelian and non-Abelian) pure lattice gauge theory projected entangled pair states in 2+1 dimensions. It is explicitly shown how the long-range behavior of these quantities gives rise to an entanglement area law in both the thermodynamic limit and in the continuum. We numerically demonstrate the applicability of our method to the ℤ 2 lattice gauge theory and relate some entanglement properties to the confinement-deconfinement transition therein. We provide evidence that Rényi entanglement entropies in certain cases do not provide a complete probe of (de)confinement properties compared to Wilson loop expectation values as other genuine (nonlocal) observables.

Abstract We study the large N dynamics of two massless Yang-Mills coupled matrix quantum mechanics, by minimization of a loop truncated Jevicki-Sakita effective collective field Hamiltonian. The loop space constraints are handled by the use of master variables. The method is successfully applied directly in the massless limit for a range of values of the Yang-Mills coupling constant, and the scaling behaviour of different physical quantities derived from their dimensions are obtained with a high level of precision. We consider both planar properties of the theory, such as the large N ground state energy and multi-matrix correlator expectation values, and also the spectrum of the theory. For the spectrum, we establish that the U(N) traced fundamental constituents remain massless and decoupled from other states, and that bound states develop well defined mass gaps, with the mass of the two degenerate lowest lying bound states being determined with a particularly high degree of accuracy. In order to confirm, numerically, the physical interpretation of the spectrum properties of the U(N) traced constituents, we add masses to the system and show that, indeed, the U(N) traced fundamental constituents retain their “bare masses”. For this system, we draw comparisons with planar results available in the literature.