We propose a new method to calculate electric dipole moments induced by the strong QCD $theta$-term. The method is based on the gradient flow for gauge fields and is free from renormalization ambiguities. We test our method by computing the nucleon e
lectric dipole moments in pure Yang-Mills theory at several lattice spacings, enabling a first-of-its-kind continuum extrapolation. The method is rather general and can be applied for any quantity computed in a $theta$ vacuum. This first application of the gradient flow has been successful and demonstrates proof-of-principle, thereby providing a novel method to obtain precise results for nucleon and light nuclear electric dipole moments.
Projection Monte Carlo calculations of lattice Chiral Effective Field Theory suffer from sign oscillations to a varying degree dependent on the number of protons and neutrons. Hence, such studies have hitherto been concentrated on nuclei with equal n
umbers of protons and neutrons, and especially on the alpha nuclei where the sign oscillations are smallest. Here, we introduce the symmetry-sign extrapolation method, which allows us to use the approximate Wigner SU(4) symmetry of the nuclear interaction to systematically extend the Projection Monte Carlo calculations to nuclear systems where the sign problem is severe. We benchmark this method by calculating the ground-state energies of the $^{12}$C, $^6$He and $^6$Be nuclei, and discuss its potential for studies of neutron-rich halo nuclei and asymmetric nuclear matter.
Luschers method is routinely used to determine meson-meson, meson-baryon and baryon-baryon s-wave scattering amplitudes below inelastic thresholds from Lattice QCD calculations - presently at unphysical light-quark masses. In this work we review the
formalism and develop the requisite expressions to extract phase-shifts describing meson-meson scattering in partial-waves with angular-momentum l<=6 and l=9. The implications of the underlying cubic symmetry, and strategies for extracting the phase-shifts from Lattice QCD calculations, are presented, along with a discussion of the signal-to-noise problem that afflicts the higher partial-waves.
The discrete energy-eigenvalues of two nucleons interacting with a finite-range nuclear force and confined to a harmonic potential are used to numerically reconstruct the free-space scattering phase shifts. The extracted phase shifts are compared to
those obtained from the exact continuum scattering solution and agree within the uncertainties of the calculations. Our results suggest that it might be possible to determine the amplitudes for the scattering of complex systems, such as n-d, n-t or n-alpha, from the energy-eigenvalues confined to finite volumes using ab-initio bound-state techniques.
This article presents differential equations and solution methods for the functions of the form $Q(x) = F^{-1}(G(x))$, where $F$ and $G$ are cumulative distribution functions. Such functions allow the direct recycling of Monte Carlo samples from one
distribution into samples from another. The method may be developed analytically for certain special cases, and illuminate the idea that it is a more precise form of the traditional Cornish-Fisher expansion. In this manner the model risk of distributional risk may be assessed free of the Monte Carlo noise associated with resampling. Examples are given of equations for converting normal samples to Student t, and converting exponential to hyperbolic, variance gamma and normal. In the case of the normal distribution, the change of variables employed allows the sampling to take place to good accuracy based on a single rational approximation over a very wide range of the sample space. The avoidance of any branching statement is of use in optimal GPU computations as it avoids the effect of {it warp divergence}, and we give examples of branch-free normal quantiles that offer performance improvements in a GPU environment, while retaining the best precision characteristics of well-known methods. We also offer models based on a low-probability of warp divergence. Comparisons of new and old forms are made on the Nvidia Quadro 4000, GTX 285 and 480, and Tesla C2050 GPUs. We argue that in single-precision mode, the change-of-variables approach offers performance competitive with the fastest existing scheme while substantially improving precision, and that in double-precision mode, this approach offers the most GPU-optimal Gaussian quantile yet, and without compromise on precision for Monte Carlo applications, working twice as fast as the CUDA 4 library function with increased precision.
I calculate finite-volume effects for three identical spin-1/2 fermions in a box assuming short-ranged repulsive interactions of `natural size. This analysis employs standard perturbation theory in powers of 1/L, where L^3 is the volume of the box. I
give results for the ground states in the A_1, T_1, and E cubic representations.
