Source author record

Neill C. Warrington

Neill C. Warrington appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

8works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2020arXiv

Complex Paths Around The Sign Problem

The Monte Carlo evaluation of path integrals is one of a few general purpose methods to approach strongly coupled systems. It is used in all branches of Physics, from QCD/nuclear physics to the correlated electron systems. However, many systems of great importance (dense matter inside neutron stars, the repulsive Hubbard model away from half-filling, dynamical and non-equilibrium observables) are not amenable to the Monte Carlo method as it currently stands due to the so-called "sign-problem". We review a new set of ideas recently developed to tackle the sign problem based on the complexification of field space and the Picard-Lefshetz theory accompanying it. The mathematical ideas underpinning this approach, as well as the algorithms so far developed, are described together with non-trivial examples where the method has already been proved successful. Directions of future work, including the burgeoning use of machine learning techniques, are delineated.

preprint2020arXiv

Finite-Temperature Quarkyonic Matter with an Excluded Volume Model

We introduce a theoretical framework to analyze the thermal properties of a recently proposed "excluded volume" model of quarkyonic matter. This entails proposing finite temperature distribution functions and entropy functionals constrained by internal theoretical consistency. We use the formalism developed here to analyze the effect of temperature on the quark onset as a function of baryon density. We find that the baryon density at which the quarkyonic phase emerges decreases only slightly with temperature, an effect produced by large $N_c$. We furthermore find that a finite temperature increases the density of quarks produced. The framework presented here can be used to compute the equation of state of hot and dense quarkyonic matter which is of relevance to neutron star mergers.

preprint2020arXiv

Path integral contour deformations for noisy observables

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.

preprint2019arXiv

Structure Factors of The Unitary Gas Under Supernova Conditions

We compute with lattice field theory the vector and axial static structure factors of the unitary gas for arbitrary temperature above the superfluid transition and for fugacities 0.1 < z < 1.0. Using the lattice formulation, we calculate beyond the validity of the virial expansion, a commonly used technique in many-body physics. We find qualitative differences in the behavior of the structure factors at high fugacity compared to the predictions of the virial expansion. Due to the large scattering length of neutrons, we expect the unitary gas structure factors to approximate the structure factors of hot neutron gases, and we therefore expect our calculations to be useful in supernova simulations, where neutron gas structure factors are needed to compute in-medium neutrino-neutron scattering rates.

preprint2016arXiv

A study of symmetry breaking in a relativistic Bose gas using the contraction algorithm

A relativistic Bose gas at finite density suffers from a sign problem that makes direct numerical simulations not feasible. One possible solution to the sign problem is to re-express the path integral in terms of Lefschetz thimbles. Using this approach we study the relativistic Bose gas both in the symmetric phase (low-density) and the spontaneously broken phase (high-density). In the high-density phase we break explicitly the symmetry and determine the dependence of the order parameter on the breaking. We study the relative contributions of the dominant and sub-dominant thimbles in this phase. We find that the sub-dominant thimble only contributes substantially when the explicit symmetry breaking is small, a regime that is dominated by finite volume effects. In the regime relevant for the thermodynamic limit, this contribution is negligible.

preprint2016arXiv

Fast Estimator of jacobians in Monte Carlo Integration on Lefschetz Thimbles

A solution to the sign problem is the so-called "Lefschetz thimble approach" where the domain of integration for field variables in the path integral is deformed from the real axis to a sub-manifold in the complex space. For properly chosen sub-manifolds ("thimbles") the sign problem disappears or is drastically alleviated. The parametrization of the thimble by real coordinates require the calculation of a jacobian with a computational cost of order O(V^3), where V is proportional to the spacetime volume. In this note we propose two estimators for this jacobian with a computational cost of order O(V). We discuss analytically the regimes where we expect the estimator to work and show numerical examples in two different models.

preprint2016arXiv

Monte Carlo study of real time dynamics

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.

preprint2016arXiv

Sign problem and Monte Carlo calculations beyond Lefschetz thimbles

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. 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 Lefshetz thimbles was elusive.