Search

We propose a novel simulation strategy for Yang-Mills theories with a complex coupling, based on the Lefschetz thimble decomposition. We envisage, that the approach developed in the present work, can also be adapted to QCD at finite density, and real time simulations. Simulations with Lefschetz thimbles offer a potential solution to sign problems in Monte Carlo calculations within many different models with complex actions. We discuss the structure of Generalized Lefschetz thimbles for pure Yang-Mills theories with a complex gauge coupling $\beta$ and show how to incorporate the gauge orbits. We propose to simulate such theories on the union of the tangential manifolds to the relevant Lefschetz thimbles attached to the critical manifolds of the Yang-Mills action. We demonstrate our algorithm on a (1+1)-dimensional U(1) model and discuss how, starting from the main thimble result, successive subleading thimbles can be taken into account via a reweighting approach. While we face a residual sign problem, our novel approach performs exponentially better than the standard reweighting approach., Comment: 18 pages, 10 figures

Lefschetz thimbles have been proposed recently as a possible solution to the complex action problem (sign problem) in Monte Carlo simulations. Here we discuss pure abelian gauge theory with a complex coupling $\beta$ and apply the concept of Generalized Lefschetz thimbles. We propose to simulate the theory on the union of the tangential manifolds to the thimbles. We construct a local Metropolis-type algorithm, that is constrained to a specific tangential manifold. We also discuss how, starting from this result, successive subleading tangential manifolds can be taken into account via a reweighting approach. We demonstrate the algorithm on $U(1)$ gauge theory in 1+1 dimensions and investigate the residual sign problem., Comment: 7 pages, 3 figures, Proceedings of the 37th Annual International Symposium on Lattice Field Theory (Lattice 2019), 16-22 June 2019, Wuhan, China

We explore artificial neural networks as a tool for the reconstruction of spectral functions from imaginary time Green's functions, a classic ill-conditioned inverse problem. Our ansatz is based on a supervised learning framework in which prior knowledge is encoded in the training data and the inverse transformation manifold is explicitly parametrised through a neural network. We systematically investigate this novel reconstruction approach, providing a detailed analysis of its performance on physically motivated mock data, and compare it to established methods of Bayesian inference. The reconstruction accuracy is found to be at least comparable, and potentially superior in particular at larger noise levels. We argue that the use of labelled training data in a supervised setting and the freedom in defining an optimisation objective are inherent advantages of the present approach and may lead to significant improvements over state-of-the-art methods in the future. Potential directions for further research are discussed in detail., Comment: 20 pages, 16 figures

We present a novel reweighting technique to calculate the relative weights in the Lefschetz thimble decomposition of a path integral. Our method is put to work using a $U(1)$ one-link model providing for a suitable testing ground and sharing many features with realistic gauge theories with fermions at finite density. We discuss prospects and future challenges to our method., Comment: 7 pages, 2 figures, 1 table, Proceedings from the 36th International Symposium on Lattice Field Theory (Lattice 2018), July 22-28, 2018, East Lansing, USA

We employ the Complex Langevin method for simulation of complex-valued actions. First, we show how to test for convergence of the method by explicitely computing boundary terms and demonstrate this in a model. Then we investigate the deconfinement phase transition of QCD with $N_f=2$ Wilson-fermions using the Complex Langevin Method and. We give preliminary results for the transition temperatures up to $\mu/T_c(\mu=0)\approx 5$ and compute the curvature coefficient $\kappa_2$., Comment: Proceedings for The 36th Annual International Symposium on Lattice Field Theory - LATTICE2018; update: added some acknowledgements

As is well known the Complex Langevin (CL) method sometimes fails to converge or converges to the wrong limit. We identified one reason for this long ago: insufficient decay of the probability density either near infinity or near poles of the drift, leading to boundary terms that spoil the formal argument for correctness. To gain a deeper understanding of this phenomenon, we analyze the emergence of such boundary terms thoroughly in a simple model, where analytic results can be compared with numerics. We also show how some simple modification stabilizes the CL process in such a way that it can produce results agreeing with direct integration. Besides explicitly demonstrating the connection between boundary terms and correct convergence our analysis also suggests a correctness criterion which could be applied in realistic lattice simulations., Comment: 15 pages, 13 figures v2: changed format to two columns for journal submission, added references and a few comments on an old criterion for correctness v3: Added some plots and discussion on the criterion from citation [13]

One of the main challenges in simulations on Lefschetz thimbles is the computation of the relative weights of contributing thimbles. In this paper we propose a solution to that problem by means of computing those weights using a reweighting procedure. Besides we present recipes for finding parametrizations of thimbles and anti-thimbles for a given theory. Moreover, we study some approaches to combine the Lefschetz thimble method with the Complex Langevin evolution. Our numerical investigations are carried out by using toy models among which we consider a one-site z^4 model as well as a U(1) one-link model., Comment: 11 pages, 9 figures, 2 tables; in revised version added a few references

In the landscape of approaches toward the simulation of Lattice Models with complex action the Complex Langevin (CL) appears as a straightforward method with a simple, well defined setup. Its applicability, however, is controlled by certain specific conditions which are not always satisfied. We here discuss the procedures to meet these conditions and the estimation of systematic errors and present some actual achievements., Comment: 35th International Symposium on Lattice Field Theory

We present a procedure for reconstructing the decision function of an artificial neural network as a simple function of the input, provided the decision function is sufficiently symmetric. In this case one can easily deduce the quantity by which the neural network classifies the input. The procedure is embedded into a pipeline of machine learning algorithms able to detect the existence of different phases of matter, to determine the position of phase transitions and to find explicit expressions of the physical quantities by which the algorithm distinguishes between phases. We assume no prior knowledge about the Hamiltonian or the order parameters except Monte Carlo-sampled configurations. The method is applied to the Ising Model and SU(2) lattice gauge theory. In both systems we deduce the explicit expressions of the known order parameters from the decision functions of the neural networks.