We point out that Monte Carlo simulations of theories with severe sign problems can be profitably performed over manifolds in complex space different from the one with fixed imaginary part of the action ("Lefschetz thimble"). We describe a family of such manifolds that interpolate between the tangent space at one critical point (where the sign problem is milder compared to the real plane but in some cases still severe) and the union of relevant thimbles (where the sign problem is mild but a multimodal distribution function complicates the Monte Carlo sampling). We exemplify this approach using a simple 0+1 dimensional fermion model previously used on sign problem studies and show that it can solve the model for some parameter values where a solution using Lefschetz thimbles was elusive.

Monte Carlo studies involving real time dynamics are severely restricted by the sign problem that emerges from highly oscillatory phase of the path integral. In this letter, we present a new method to compute real time quantities on the lattice using the Schwinger-Keldysh formalism via Monte Carlo simulations. The key idea is to deform the path integration domain to a complex manifold where the phase oscillations are mild and the sign problem is manageable. We use the previously introduced "contraction algorithm" to create a Markov chain on this alternative manifold. We substantiate our approach by analyzing the quantum mechanical anharmonic oscillator. Our results are in agreement with the exact ones obtained by diagonalization of the Hamiltonian. The method we introduce is generic and in principle applicable to quantum field theory albeit very slow. We discuss some possible improvements that should speed up the algorithm.

We present results of the numerical simulation of the two-dimensional Thirring model at finite density and temperature. The severe sign problem is dealt with by deforming the domain of integration into complex field space. This is the first example where a fermionic sign problem is solved in a quantum field theory by using the holomorphic gradient flow approach, a generalization of the Lefschetz thimble method.

Simulations of gauge theories on quantum computers require the digitization of continuous field variables. Digitization schemes that uses the minimum amount of qubits are desirable. We present a practical scheme for digitizing SU (3) gauge theories via its discrete subgroup S(1080). The S(1080) standard Wilson action cannot be used since a phase transition occurs as the coupling is decreased, well before the scaling regime. We proposed a modified action that allows simulations in the scaling window and carry out classical Monte Carlo calculations down to lattice spacings of order a ≈ 0.08 fm. We compute a set of observables with sub-percent precision at multiple lattice spacings and show that the continuum extrapolated value agrees with the full SU (3) results. This suggests that this digitization scheme provides sufficient precision for NISQ-era QCD simulations.

No abstract