We introduce a general and accurate method for determining lattice phase shifts and mixing angles, which is applicable to arbitrary, non-cubic lattices. Our method combines angular momentum projection, spherical wall boundaries and an adjustable auxiliary potential. This allows us to construct radial lattice wave functions and to determine phase shifts at arbitrary energies. For coupled partial waves, we use a complex-valued auxiliary potential that breaks time-reversal invariance. We benchmark our method using a system of two spin-1/2 particles interacting through a finite-range potential with a strong tensor component. We are able to extract phase shifts and mixing angles for all angular momenta and energies, with precision greater than that of extant methods. We discuss a wide range of applications from nuclear lattice simulations to optical lattice experiments.
We study the breaking of rotational symmetry on the lattice for irreducible tensor operators and practical methods for suppressing this breaking. We illustrate the features of the general problem using an $alpha$ cluster model for $^{8}$Be. We focus on the lowest states with non-zero angular momentum and examine the matrix elements of multipole moment operators. We show that the physical reduced matrix element is well reproduced by averaging over all possible orientations of the quantum state, and this is expressed as a sum of matrix elements weighted by the corresponding Clebsch-Gordan coefficients. For our $alpha$ cluster model we find that the effects of rotational symmetry breaking can be largely eliminated for lattice spacings of $aleq 1.7$ fm, and we expect similar improvement for actual lattice Monte Carlo calculations.
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 numbers 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.
We explore the breaking of rotational symmetry on the lattice for bound state energies and practical methods for suppressing this breaking. We demonstrate the general problems associated with lattice discretization errors and finite-volume errors using an $alpha$ cluster model for $^8$Be and $^{12}$C. We consider the two and three $alpha$-particle systems and focus on the lowest states with non-zero angular momentum which split into multiplets corresponding to different irreducible representations of the cubic group. We examine the dependence of such splittings on the lattice spacing and box size. We find that lattice spacing errors are closely related to the commensurability of the lattice with the intrinsic length scales of the system. We also show that rotational symmetry breaking effects can be significantly reduced by using improved lattice actions, and that the physical energy levels are accurately reproduced by the weighted average of a given spin multiplets.
We present evidence, from Lattice Monte Carlo simulations of the phase diagram of graphene as a function of the Coulomb coupling between quasiparticles, that graphene in vacuum is likely to be an insulator. We find a semimetal-insulator transition at $alpha_g^text{crit} = 1.11 pm 0.06$, where $alpha_g^{} simeq 2.16$ in vacuum, and $alpha_g^{} simeq 0.79$ on a SiO$_2^{}$ substrate. Our analysis uses the logarithmic derivative of the order parameter, supplemented by an equation of state. The insulating phase disappears above a critical number of four-component fermion flavors $4 < N_f^{text{crit}} < 6$. Our data are consistent with a second-order transition.
