Variational Monte Carlo Calculations of n+3H Scattering

Abstract: A paramount goal in the field of nuclear physics is to unify ab-initio treatments of bound and unbound states. The position-space quantum Monte Carlo (QMC) methods have a long history of successful bound state calculations in light systems but so far have seen very little application to unbound systems. Here we introduce a numerical method to improve the efficiency and accuracy of unbound-state calculations in QMC, and as an initial application we compute scattering observables for the neutron-triton system using variational Monte Carlo (VMC) wave functions. The method consists of inferring long-range amplitudes in the wave function from integrals over the short-ranged region where all the particles interact. This approach using integral relations is well established in the literature; here we develop it for the QMC framework. We validate our code with a consistency check between short-range spectroscopic overlap functions computed from direct evaluation and from the integral relations; scattering amplitudes are long-range asymptotics of those overlaps. Comparison against published benchmark calculations using the same potential demonstrates that the integral method, when applied to the current VMC wave functions, produces more accurate scattering observables than direct evaluation from the same variational wave function. However, it still differs noticeably from exact results. We then present phase shifts and mixing parameters for the neutron-triton system using various interactions. Application of the integral method here paves the way for its use in Green's function Monte Carlo (GFMC) calculations. In GFMC the wave functions are more precise, but high-precision convergence of their tails is slow, and there are additional difficulties in reading out amplitudes. The integral methods will address both of those remaining problems.