Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userReal-time lattice gauge theory actions: unitarity, convergence, and path integral contour deformations
74
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
Path integral contour deformations have been shown to mitigate sign and signal-to-noise problems associated with phase fluctuations in lattice field theories. We define a family of contour deformations applicable to $SU(N)$ lattice gauge theory that can reduce sign and signal-to-noise problems associated with complex actions and complex observables. For observables, these contours can be used to define deformed observables with identical expectation value but different variance. As a proof-of-principle, we apply machine learning techniques to optimize the deformed observables associated with Wilson loops in two dimensional $SU(2)$ and $SU(3)$ gauge theory. We study loops consisting of up to 64 plaquettes and achieve variance reduction of up to 4 orders of magnitude.
Monte Carlo studies of many quantum systems face exponentially severe signal-to-noise problems. We show that noise arising from complex phase fluctuations of observables can be reduced without introducing bias using path integral contour deformation techniques. A numerical study of contour deformations for correlation functions in Abelian gauge theory and complex scalar field theory demonstrates that variance can be reduced by orders of magnitude without modifying Monte Carlo sampling.
We investigate simulations for gauge theories on a Minkowskian space-time lattice. We employ stochastic quantization with optimized updating using stochastic reweighting or gauge fixing, respectively. These procedures do not affect the underlying theory but strongly improve the stability properties of the stochastic dynamics, such that simulations on larger real-time lattices can be performed.
Direct numerical evaluation of the real-time path integral has a well-known sign problem that makes convergence exponentially slow. One promising remedy is to use Picard-Lefschetz theory to flow the domain of the field variables into the complex plane, where the integral is better behaved. By Cauchys theorem, the final value of the path integral is unchanged. Previous analyses have considered the case of real scalar fields in thermal equilibrium, employing a closed Schwinger-Keldysh time contour, allowing the evaluation of the full quantum correlation functions. Here we extend the analysis by not requiring a closed time path, instead allowing for an initial density matrix for out-of-equilibrium initial value problems. We are able to explicitly implement Gaussian initial conditions, and by separating the initial time and the later times into a two-step Monte-Carlo sampling, we are able to avoid the phenomenon of multiple thimbles. In fact, there exists one and only one thimble for each sample member of the initial density matrix. We demonstrate the approach through explicitly computing the real-time propagator for an interacting scalar in 0+1 dimensions, and find very good convergence allowing for comparison with perturbation theory and the classical-statistical approximation to real-time dynamics.
Recently, a finite-temperature real-time static potential has been introduced via a Schrodinger-type equation satisfied by a certain heavy quarkonium Greens function. Furthermore, it has been pointed out that it possesses an imaginary part, which induces a finite width for the tip of the quarkonium peak in the thermal dilepton production rate. The imaginary part originates from Landau-damping of low-frequency gauge fields, which are essentially classical due to their high occupation number. Here we show how the imaginary part can be measured with classical lattice gauge theory simulations, accounting non-perturbatively for the infrared sector of finite-temperature field theory. We demonstrate that a non-vanishing imaginary part indeed exists non-perturbatively; and that its value agrees semi-quantitatively with that predicted by Hard Loop resummed perturbation theory.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
