We present a study of the IR behaviour of a three-dimensional super-renormalisable quantum field theory (QFT) consisting of a scalar field in the adjoint of $SU(N)$ with a $varphi^4$ interaction. A bare mass is required for the theory to be massless
at the quantum level. In perturbation theory the critical mass is ambiguous due to infrared (IR) divergences and we indeed find that at two-loops in lattice perturbation theory the critical mass diverges logarithmically. It was conjectured long ago in [Jackiw 1980, Appelquist 1981] that super-renormalisable theories are nonperturbatively IR finite, with the coupling constant playing the role of an IR regulator. Using a combination of Markov-Chain-Monte-Carlo simulations of the lattice-regularised theory, both frequentist and Bayesian data analysis, and considerations of a corresponding effective theory we gather evidence that this is indeed the case.
We compute the topological susceptibility $chi_t$ of 2+1-flavor lattice QCD with dynamical Mobius domain-wall fermions, whose residual mass is kept at 1 MeV or smaller. In our analysis, we focus on the fluctuation of the topological charge density in
a slab sub-volume of the simulated lattice, as proposed by Bietenholz et al. The quark mass dependence of our results agrees well with the prediction of the chiral perturbation theory, from which the chiral condensate is extracted. Combining the results for the pion mass $M_pi$ and decay constant $F_pi$, we obtain $chi_t$ = 0.227(02)(11)$M_pi^2 F_pi^2$ at the physical point, where the first error is statistical and the second is systematic.
We discuss practical methods to ensure near wirespeed performance from clusters with either one or two Intel(R) Omni-Path host fabric interfaces (HFI) per node, and Intel(R) Xeon Phi(TM) 72xx (Knights Landing) processors, and using the Linux operatin
g system. The study evaluates the performance improvements achievable and the required programming approaches in two distinct example problems: firstly in Cartesian communicator halo exchange problems, appropriate for structured grid PDE solvers that arise in quantum chromodynamics simulations of particle physics, and secondly in gradient reduction appropriate to synchronous stochastic gradient descent for machine learning. As an example, we accelerate a published Baidu Research reduction code and obtain a factor of ten speedup over the original code using the techniques discussed in this paper. This displays how a factor of ten speedup in strongly scaled distributed machine learning could be achieved when synchronous stochastic gradient descent is massively parallelised with a fixed mini-batch size. We find a significant improvement in performance robustness when memory is obtained using carefully allocated 2MB huge virtual memory pages, implying that either non-standard allocation routines should be used for communication buffers. These can be accessed via a LD_PRELOAD override in the manner suggested by libhugetlbfs. We make use of a the Intel(R) MPI 2019 library Technology Preview and underlying software to enable thread concurrency throughout the communication software stake via multiple PSM2 endpoints per process and use of multiple independent MPI communicators. When using a single MPI process per node, we find that this greatly accelerates delivered bandwidth in many core Intel(R) Xeon Phi processors.
We study spatial isovector meson correlators in $N_f=2$ QCD with dynamical domain-wall fermions on $32^3times 8$ lattices at temperatures up to 380 MeV with various quark masses. We measure the correlators of spin-one isovector operators including ve
ctor, axial-vector, tensor and axial-tensor. At temperatures above $T_c$ we observe an approximate degeneracy of the correlators in these channels, which is unexpected because some of them are not related under $SU(2)_L times SU(2)_R$ nor $U(1)_A$ symmetries. The observed approximate degeneracy suggests emergent $SU(2)_{CS}$ (chiral-spin) and $SU(4)$ symmetries at high $T$.
We compute the chiral condensate in 2+1-flavor QCD through the spectrum of low-lying eigenmodes of Dirac operator. The number of eigenvalues of the Dirac operator is evaluated using a stochastic method with an eigenvalue filtering technique on the ba
ckground gauge configurations generated by lattice QCD simulations including the effects of dynamical up, down and strange quarks described by the Mobius domain-wall fermion formulation. The low-lying spectrum is related to the chiral condensate, which is one of the leading order low-energy constants in chiral effective theory, as dictated by the Banks-Casher relation. The spectrum shape and its dependence on the sea quark masses calculated in numerical simulations are consistent with the expectation from one-loop chiral perturbation theory. After taking the chiral limit as well as the continuum limit using the data at three lattice spacings ranging 0.080-0.045 fm, we obtain $Sigma^{1/3}$(2 GeV) = 270.0(4.9) MeV, with the error combining those from statistical and from various sources of systematic errors. Finite volume effect is confirmed to be under control by a direct comparison of the results from two different volumes at the lightest available sea quarks corresponding to 230 MeV pions.
In this proceedings we discuss the motivation, implementation details, and performance of a new physics code base called Grid. It is intended to be more performant, more general, but similar in spirit to QDP++cite{QDP}. Our approach is to engineer th
e basic type system to be consistently fast, rather than bolt on a few optimised routines, and we are attempt to write all our optimised routines directly in the Grid framework. It is hoped this will deliver best known practice performance across the next generation of supercomputers, which will provide programming challenges to traditional scalar codes. We illustrate the programming patterns used to implement our goals, and advances in productivity that have been enabled by using new features in C++11.
We explore the phase structure and symmetry breaking in four-dimensional SU(3) gauge theory with one spatial compact dimension on the lattice in the presence of fermions in the adjoint and fundamental representations with general boundary conditions.
The eigenvalue phases of Polyakov loops and the associated susceptibility are measured on 16^3 x 4 lattice. We establish a correspondence between the phases found on the lattice and the gauge symmetry breaking by the Hosotani mechanism.
We examine the axial U(1) symmetry near and above the finite temperature phase transition in two-flavor QCD using lattice QCD simulations. Although the axial U(1) symmetry is always violated by quantization, (i.e.) the chiral anomaly, the correlation
functions may manifest effective restoration of the symmetry in the high temperature phase. We explicitly study this possibility by calculating the meson correlators as well as the Dirac operator spectral density near the critical point. Our numerical simulations are performed on a $16^3times 8$ lattice with two flavors of dynamical quarks represented by the overlap fermion formalism. Chiral symmetry and its violation due to the axial anomaly is manifestly realized with this formulation, which is a prerequisite for the study of the effective restoration of the axial U(1) symmetry. In order to avoid discontinuity in the gauge configuration space, which occurs for the exactly chiral lattice fermions, the simulation is confined in a fixed topological sector. It induces finite volume effect, which is well described by a formula based on the Fourier transform from the $theta$-vacua. We confirm this formula at finite temperature by calculating the topological susceptibility in the quenched theory. Our two flavor simulations show degeneracy of the meson correlators and a gap in the Dirac operator spectral density, which implies that the axial U(1) symmetry is effectively restored in the chirally symmetric phase.
