Search

We present an exact solution for the linearized junction array modes of the superconducting qubit fluxonium in the absence of array disorder. This solution holds for arrays of any length and ground capacitance, and for both differential and grounded devices. Array mode energies are determined by roots of convex combinations of Chebyshev polynomials, and their spatial profiles are plane waves. We also provide a simple, approximate solution, which estimates array mode properties over a wide range of circuit parameters, and an accompanying Mathematica file that implements both the exact and approximate solutions., Comment: 11 pages, 3 figures, 1 ancillary Mathematica notebook

We construct a lattice field theory method for computing the real-time dynamics of spin systems in a thermal bath. This is done by building on previous work of Takano with Schwinger-Keldysh and functional differentiation techniques. We derive a Schwinger-Keldysh path integral for generic spin Hamiltonians, then demonstrate the method on a simple system. Our path integral has a sign problem, which generally requires exponential run time in the system size, but requires only linear storage. The latter may place this method at an advantage over exact diagonalization, which is exponential in both. Our path integral is amenable to contour deformations, a technique for reducing sign problems., Comment: 7 pages, 2 figures, 1 table

We examine the behavior of charged pions in neutron-rich matter using heavy-baryon chiral perturbation theory. This study is motivated by the prospect that pions, or pion-like excitations, may be relevant in neutron-rich matter encountered in core-collapse supernovae and neutron star mergers. We find, as previously expected, that the $\pi^-$ mass increases with density and precludes s-wave condensation at $n_B\lesssim n_{\text{sat}}$, where $n_{\text{sat}} \simeq 0.16 \text{fm}^{-3}$ is the nuclear saturation density, and the mass of the $\pi^+$ mode decreases with density. The uncertainty in these predictions increases rapidly for $n_B \gtrsim n_{\text{sat}} $ because low energy constants associated with the two-pion-two-nucleon operators in chiral perturbation theory are poorly constrained. We find that these uncertainties are especially large in symmetric nuclear matter and should be included in the analysis of pion-nucleus interactions at low energy and pionic atoms. In neutron-rich matter, accounting for the self-energy difference between neutrons and protons related to the nuclear symmetry energy has several effects. It alters the power counting of certain higher-order contributions to the pion self-energy. Previously unimportant but attractive diagrams are enhanced and result in a modest reduction of the pion masses. Furthermore, in the low-wavelength limit, a collective mode with the quantum numbers of the $\pi^+$ appears., Comment: 16 pages, 6 figures, 1 appendix

Monte Carlo simulations are useful tools for modeling quantum systems, but in some cases they suffer from a sign problem, leading to an exponential slow down in their convergence to a value. While solving the sign problem is generically NP-hard, many techniques exist for mitigating the sign problem in specific cases; in particular, the technique of deforming the Monte Carlo simulation's plane of integration onto Lefschetz thimbles (complex hypersurfaces of stationary phase) has seen significant success in the context of quantum field theories. We extend this methodology to spin systems by utilizing spin coherent state path integrals to re-express the spin system's partition function in terms of continuous variables. Using some toy systems, we demonstrate its effectiveness at lessening the sign problem in this setting, despite the fact that the initial mapping to spin coherent states introduces its own sign problem. The standard formulation of the spin coherent path integral is known to make use of uncontrolled approximations; despite this, for large spins they are typically considered to yield accurate results, so it is somewhat surprising that our results show significant systematic errors. Therefore, possibly of independent interest, our use of Lefschetz thimbles to overcome the intrinsic sign problem in spin coherent state path integral Monte Carlo enables a novel numerical demonstration of a breakdown in the spin coherent path integral., Comment: 11 pages, 5 figures

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., Comment: 25 pages, 12 figures

We compute continuum and infinite volume limit extrapolations of the structure factors of neutron matter at finite temperature and density. Using a lattice formulation of leading-order pionless effective field theory, we compute the momentum dependence of the structure factors at finite temperature and at densities beyond the reach of the virial expansion. The Tan contact parameter is computed and the result agrees with the high momentum tail of the vector structure factor. All errors, statistical and systematic, are controlled for. This calculation is a first step towards a model-independent understanding of the linear response of neutron matter at finite temperature, a realm until now little explored., Comment: 5 pages, 3 figures, data included in submission; further experimental data included

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., Comment: 35 pages, 21 figures

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., Comment: 7 pages, 2 figures

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., Comment: 14 pages, 6 figures

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., Comment: 10 pages, 6 figures, 1 table

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\approx 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., Comment: 5 pages, 5 figures

We present Monte Carlo calculations of the thermodynamics of the (2+1) dimensional Thirring model at finite density. We bypass the sign problem by deforming the domain of integration of the path integral into complex space in such a way as to maximize the average sign within a parameterized family of manifolds. We present results for lattice sizes up to $10^3$ and we find that at high densities and/or temperatures the chiral condensate is abruptly reduced., Comment: 5 pages, 4 figures

We study by Monte Carlo methods the thermodynamics of a spin polarized gas of non-relativistic fermions in 1+1 dimensions. The main result of this work is that our action suffers no significant sign problem for any spin polarization in the region relevant for dilute degenerate fermi gases. This lack of sign problem allows us to study attractive spin polarized fermions non-perturbatively at spin polarizations not previously explored. For some parameters values we verify results previously obtained by methods which include an uncontrolled step like complex Langevin and/or analytical continuation from imaginary chemical potential. For others, larger values of the polarization, we deviate from these previous results., Comment: 9 pages, 5 figures

We calculate the differential scattering rate for thermal neutrinos in a hot and dilute gas of interacting neutrons using linear response theory. The dynamical structure factors for density and spin fluctuations of the strongly interacting neutron matter, expected in the neutrino decoupling regions of supernovae and neutron star mergers, are calculated in the virial expansion for the first time. Correlations due to nucleon-nucleon interactions are taken into account using a pseudo-potential that reproduces measured nucleon-nucleon phase shifts, and we find that attractive s-wave interactions enhance the density response and suppress the spin response of neutron matter. The net effect of neutron correlations is to strongly suppress backscattering. Moreover, we find nearly exact scaling laws for the response functions, valid for the range $T = 5 - 10$ MeV and q < 30 MeV, allowing us to obtain analytic results for the dynamic structure factors at second-order in the fugacity of the neutron gas. We find that the modification of scattering rates depends on the energy and momentum exchanged, implying that dynamical structure factors are essential to describe neutrino decoupling in supernovae and neutron star mergers.

Quantum field theories with complex actions cannot be investigated using importance sampling due to the sign problem. One possible solution is to use the holomorphic gradient flow, a method we introduced related to the Lefschetz thimbles idea. In many cases the probability distribution generated by this method is multi-modal and standard Monte-Carlo sampling fails. We propose an algorithm that incorporates tempered proposals to solve this problem. We apply this algorithm to the 0+1 dimensional Thirring model at finite density for a parameter set where standard sampling fails and show that tempered proposals cure this problem., Comment: 10 pages, 6 figures

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., Comment: 6 pages, 5 figures

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., Comment: 12 pages, 6 figures, 1 table

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., Comment: 5 pages, 3 figures

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., Comment: 10 pages, 3 figures

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., Comment: 12 pages, 7 figures, figure added, corrected typos