Parallel software for lattice N=4 supersymmetric Yang--Mills theory

We present new parallel software, SUSY LATTICE, for lattice studies of four-dimensional $\mathcal N = 4$ supersymmetric Yang--Mills theory with gauge group SU(N). The lattice action is constructed to exactly preserve a single supersymmetry charge at non-zero lattice spacing, up to additional potential terms included to stabilize numerical simulations. The software evolved from the MILC code for lattice QCD, and retains a similar large-scale framework despite the different target theory. Many routines are adapted from an existing serial code, which SUSY LATTICE supersedes. This paper provides an overview of the new parallel software, summarizing the lattice system, describing the applications that are currently provided and explaining their basic workflow for non-experts in lattice gauge theory. We discuss the parallel performance of the code, and highlight some notable aspects of the documentation for those interested in contributing to its future development.Comment: Code available at https://github.com/daschaich/sus

Enhanced goal-oriented error assessment and computational strategies in adaptive reduced basis solver for stochastic problems

This work focuses on providing accurate low-cost approximations of stochastic ¿nite elements simulations in the framework of linear elasticity. In a previous work, an adaptive strategy was introduced as an improved Monte-Carlo method for multi-dimensional large stochastic problems. We provide here a complete analysis of the method including a new enhanced goal-oriented error estimator and estimates of CPU (computational processing unit) cost gain. Technical insights of these two topics are presented in details, and numerical examples show the interest of these new developments.Postprint (author's final draft

Multigrid Methods in Lattice Field Computations

The multigrid methodology is reviewed. By integrating numerical processes at all scales of a problem, it seeks to perform various computational tasks at a cost that rises as slowly as possible as a function of $n$, the number of degrees of freedom in the problem. Current and potential benefits for lattice field computations are outlined. They include: $O(n)$ solution of Dirac equations; just $O(1)$ operations in updating the solution (upon any local change of data, including the gauge field); similar efficiency in gauge fixing and updating; $O(1)$ operations in updating the inverse matrix and in calculating the change in the logarithm of its determinant; $O(n)$ operations per producing each independent configuration in statistical simulations (eliminating CSD), and, more important, effectively just $O(1)$ operations per each independent measurement (eliminating the volume factor as well). These potential capabilities have been demonstrated on simple model problems. Extensions to real life are explored.Comment: 4

Two-Qubit Separabilities as Piecewise Continuous Functions of Maximal Concurrence

The generic real (b=1) and complex (b=2) two-qubit states are 9-dimensional and 15-dimensional in nature, respectively. The total volumes of the spaces they occupy with respect to the Hilbert-Schmidt and Bures metrics are obtainable as special cases of formulas of Zyczkowski and Sommers. We claim that if one could determine certain metric-independent 3-dimensional "eigenvalue-parameterized separability functions" (EPSFs), then these formulas could be readily modified so as to yield the Hilbert-Schmidt and Bures volumes occupied by only the separable two-qubit states (and hence associated separability probabilities). Motivated by analogous earlier analyses of "diagonal-entry-parameterized separability functions", we further explore the possibility that such 3-dimensional EPSFs might, in turn, be expressible as univariate functions of some special relevant variable--which we hypothesize to be the maximal concurrence (0 < C <1) over spectral orbits. Extensive numerical results we obtain are rather closely supportive of this hypothesis. Both the real and complex estimated EPSFs exhibit clearly pronounced jumps of magnitude roughly 50% at C=1/2, as well as a number of additional matching discontinuities.Comment: 12 pages, 7 figures, new abstract, revised for J. Phys.

Status and Future Perspectives for Lattice Gauge Theory Calculations to the Exascale and Beyond

In this and a set of companion whitepapers, the USQCD Collaboration lays out a program of science and computing for lattice gauge theory. These whitepapers describe how calculation using lattice QCD (and other gauge theories) can aid the interpretation of ongoing and upcoming experiments in particle and nuclear physics, as well as inspire new ones.Comment: 44 pages. 1 of USQCD whitepapers

Practical solution to the Monte Carlo sign problem: Realistic calculations of 54Fe

We present a practical solution to the "sign problem" in the auxiliary field Monte Carlo approach to the nuclear shell model. The method is based on extrapolation from a continuous family of problem-free Hamiltonians. To demonstrate the resultant ability to treat large shell-model problems, we present results for 54Fe in the full fp-shell basis using the Brown-Richter interaction. We find the Gamow-Teller beta^+ strength to be quenched by 58% relative to the single-particle estimate, in better agreement with experiment than previous estimates based on truncated bases.Comment: 11 pages + 2 figures (not included

DNA electrophoresis studied with the cage model

The cage model for polymer reptation, proposed by Evans and Edwards, and its recent extension to model DNA electrophoresis, are studied by numerically exact computation of the drift velocities for polymers with a length L of up to 15 monomers. The computations show the Nernst-Einstein regime (v ~ E) followed by a regime where the velocity decreases exponentially with the applied electric field strength. In agreement with de Gennes' reptation arguments, we find that asymptotically for large polymers the diffusion coefficient D decreases quadratically with polymer length; for the cage model, the proportionality coefficient is DL^2=0.175(2). Additionally we find that the leading correction term for finite polymer lengths scales as N^{-1/2}, where N=L-1 is the number of bonds.Comment: LaTeX (cjour.cls), 15 pages, 6 figures, added correctness proof of kink representation approac